Hans Fangohr

Skyrmion-skyrmion and skyrmion-edge repulsions in skyrmion-based racetrack memory

Magnetic skyrmions are promising for building next-generation magnetic memories and spintronic devices due to their stability, small size and the extremely low currents needed to move them. In particular, skyrmion-based racetrack memory is attractive for information technology, where skyrmions are used to store information as data bits instead of traditional domain walls. Here we numerically demonstrate the impacts of skyrmion-skyrmion and skyrmion-edge repulsions on the feasibility of skyrmion-based racetrack memory. The reliable and practicable spacing between consecutive skyrmionic bits on the racetrack as well as the ability to adjust it are investigated. Clogging of skyrmionic bits is found at the end of the racetrack, leading to the reduction of skyrmion size. Further, we demonstrate an effective and simple method to avoid the clogging of skyrmionic bits, which ensures the elimination of skyrmionic bits beyond the reading element. Our results give guidance for the design and development of future skyrmion-based racetrack memory.

Ground state search, hysteretic behaviour, and reversal mechanism of skyrmionic textures in confined helimagnetic nanostructures

Magnetic skyrmions have the potential to provide solutions for low-power, high-density data storage and processing. One of the major challenges in developing skyrmion-based devices is the skyrmions’ magnetic stability in confined helimagnetic nanostructures. Through a systematic study of equilibrium states, using a full three-dimensional micromagnetic model including demagnetisation effects, we demonstrate that skyrmionic textures are the lowest energy states in helimagnetic thin film nanostructures at zero external magnetic field and in absence of magnetocrystalline anisotropy. We also report the regions of metastability for non-ground state equilibrium configurations. We show that bistable skyrmionic textures undergo hysteretic behaviour between two energetically equivalent skyrmionic states with different core orientation, even in absence of both magnetocrystalline and demagnetisation-based shape anisotropies, suggesting the existence of Dzyaloshinskii-Moriya-based shape anisotropy. Finally, we show that the skyrmionic texture core reversal dynamics is facilitated by the Bloch point occurrence and propagation.

Joule heating in nanowires

We study the effect of Joule heating from electric currents flowing through ferromagnetic nanowires on the temperature of the nanowires and on the temperature of the substrate on which the nanowires are grown. The spatial current density distribution, the associated heat generation, and diffusion of heat is simulated within the nanowire and the substrate. We study several different nanowire and constriction geometries as well as different substrates: (thin) silicon nitride membranes, (thick) silicon wafers, and (thick) diamond wafers. The spatially resolved increase in temperature as a function of time is computed. For effectively three-dimensional substrates (where the substrate thickness greatly exceeds the nanowire length), we identify three different regimes of heat propagation through the substrate: regime (i), where the nanowire temperature increases approximately logarithmically as a function of time. In this regime, the nanowire temperature is well-described analytically by You et al. [APL89, 222513 (2006)]. We provide an analytical expression for the time tc that marks the upper applicability limit of the You model. After tc, the heat flow enters regime (ii), where the nanowire temperature stays constant while a hemispherical heat front carries the heat away from the wire and into the substrate. As the heat front reaches the boundary of the substrate, regime (iii) is entered where the nanowire and substrate temperature start to increase rapidly. For effectively two-dimensional substrates (where the nanowire length greatly exceeds the sub- strate thickness), there is only one regime in which the temperature increases logarithmically with time for large times. We provide an analytical expression, valid for all pulse durations, that allows one to accurately compute this temperature increase in the nanowire on thin substrates

BioSimGrid: grid-enabled biomolecular simulation data storage and analysis

In computational biomolecular research, large amounts of simulation data are generated to capture the motion of proteins. These massive simulation data can be analysed in a number of ways to reveal the biochemical properties of the proteins. However, the legacy way of storing these data (usually in the laboratory where the simulations have been run) often hinders a wider sharing and easier cross-comparison of simulation results. The data is commonly encoded in a way specific to the simulation package that produced the data and can only be analysed with tools developed specifically for that simulation package. The BioSimGrid platform seeks to provide a solution to these challenges by exploiting the potential of the Grid in facilitating data sharing. By using BioSimGrid either in a scripting or web environment, users can deposit their data and reuse it for analysis. BioSimGrid tools manage the multiple storage locations transparently to the users and provide a set of retrieval and analysis tools for processing the data in a convenient and efficient manner. This paper details the usage and implementation of BioSimGrid using a combination of commercial databases, the Storage Resource Broker and Python scripts, gluing the building blocks together. It introduces a case study of how BioSimGrid can be used for better storage, retrieval and analysis of biomolecular simulation data.

The stability of buoyant bubbles in the atmospheres of galaxy clusters

The buoyant rise of hot plasma bubbles inflated by active galactic nuclei outflows in galaxy clusters can heat the cluster gas and thereby compensate radiative energy losses of this material. Numerical simulations of this effect often show the complete disruption of the bubbles followed by the mixing of the bubble material with the surrounding cluster gas due to fluid instabilities on the bubble surface. This prediction is inconsistent with the observations of apparently coherent bubble structures in clusters. We derive a general description in the linear regime of the growth of instabilities on the surface between two fluids under the influence of a gravitational field, viscosity, surface tension provided by a magnetic field and relative motion of the two fluids with respect to each other. We demonstrate that Kelvin–Helmholtz instabilities are always suppressed, if the fluids are viscous. They are also suppressed in the inviscid case for fluids of very different mass densities. We show that the effects of shear viscosity as well as a magnetic field in the cluster gas can prevent the growth of Rayleigh–Taylor instabilities on relevant scalelengths. Rayleigh–Taylor instabilities on parsec scales are suppressed even if the kinematic viscosity of the cluster gas is reduced by two orders of magnitude compared to the value given by Spitzer for a fully ionized, unmagnetized gas. Similarly, magnetic fields exceeding a few ?G result in an effective surface tension preventing the disruption of bubbles. For more massive clusters, instabilities on the bubble surface grow faster. This may explain the absence of thermal gas in the north-west bubble observed in the Perseus cluster compared to the apparently more disrupted bubbles in the Virgo cluster.

A Comparison of C, MATLAB, and Python as Teaching Languages in Engineering

We describe and compare the programming languages C, MATLAB and Python as teaching languages for engineering students. We distinguish between two distinct phases in the process of converting a given problem into a computer program that can provide a solution: (i) finding an algorithmic solution and (ii) implementing this in a particular programming language. It is argued that it is most important for the understanding of the students to perform the first step whereas the actual implementation in a programming language is of secondary importance for the learning of problem-solving techniques. We therefore suggest to chose a well-structured teaching language that provides a clear and intuitive syntax and allows students to quickly express their algorithms. In our experience in engineering computing we find that MATLAB is much better suited than C for this task but the best choice in terms of clarity and functionality of the language is provided by Python.

Proposal for a Standard Micromagnetic Problem: Spin Wave Dispersion in a Magnonic Waveguide

In this paper, we propose a standard micromagnetic problem, of a nanostripe of permalloy. We study the magnetization dynamics and describe methods of extracting features from simulations. Spin wave dispersion curves, relating frequency and wave vector, are obtained for wave propagation in different directions relative to the axis of the waveguide and the external applied field. Simulation results using both finite element (Nmag) and finite difference (OOMMF) methods are compared against analytic results, for different ranges of the wave vector.

Proposal for a Standard Problem for Micromagnetic Simulations Including Spin-Transfer Torque

of micromagnetic simulation tools. The work is based on the micromagnetic model extended by the spin-transfer torque in continuously varying magnetizations as proposed by Zhang and Li. The standard problem geometry is a permalloy cuboid of 100 nm edge length and 10 nm thickness, which contains a Landau pattern with a vortex in the center of the structure. A spin-polarized dc current density of 10 12 A/m 2 ows laterally through the cuboid and moves the vortex core to a new steady-state position. We show that the new vortex-core position is a sensitive measure for the correctness of micromagnetic simulators that include the spin-transfer torque. The suitability of the proposed problem as a standard problem is tested by numerical results from four dierent nite-dierence and nite-element-based simulation tools.

The effects of thermal conduction on the intracluster medium of the Virgo cluster

Thermal conduction has been suggested as a possible mechanism by which sufficient energy is supplied to the central regions of galaxy clusters to balance the effect of radiative cooling. Here we present the results of a simulated, high-resolution, 3-d Virgo cluster for different values of thermal conductivity (1, 1/10, 1/100, 0 times the full Spitzer value). Starting from an initially isothermal cluster atmosphere we allow the cluster to evolve freely over timescales of roughly

Heating rate profiles in galaxy clusters

1.3-4.7 \times 10^{9}