Magnetic Direct-Write Skyrmion Nanolithography.

Magnetic skyrmions are stable spin textures with quasi-particle behavior and attract significant interest in fundamental and applied physics. The metastability of magnetic skyrmions at zero magnetic field is particularly important to enable, for instance, a skyrmion racetrack memory. Here, the results of the nucleation of stable skyrmions and formation of ordered skyrmion lattices by magnetic force microscopy in (Pt/CoFeSiB/W)n multilayers, exploiting the additive effect of the interfacial Dzyaloshinskii-Moriya interaction, are presented. The appropriate conditions under which skyrmion lattices are confined with a dense two-dimensional liquid phase are identified. A crucial parameter to control the skyrmion lattice characteristics and the number of scans resulting in the complete formation of a skyrmion lattice is the distance between two adjacent scanning lines of a magnetic force microscopy probe. The creation of skyrmion patterns with complex geometry is demonstrated, and the physical mechanism of direct magnetic writing of skyrmions is comprehended by micromagnetic simulations. This study shows a potential of a direct-write (maskless) skyrmion (topological) nanolithography with sub-100 nm resolution, where each skyrmion acts as a pixel in the final topological image.

Controlling the wetting properties of the Asakura-Oosawa model and applications to spherical confinement.

We demonstrate for the Asakura-Oosawa model and an extension of this model that uses continuous rather than hard potentials, how wetting properties at walls can be easily controlled. By increasing the interaction range of the repulsive wall potential acting on the colloids (while keeping the polymer-wall interactions constant) polymers begin to substitute colloids at walls and the system can be driven from complete wetting of colloids via partial wetting to complete wetting of polymers. As an application, we discuss the morphology and wetting behavior of colloid-polymer mixtures in spherical confinement. We apply the recently developed 'ensemble switch method' where the Hamiltonian is extended to a combination of a system with walls and of a system without walls to calculate the surface excess free energies of colloid-rich and polymer-rich phases. The contact angle then is inferred from Young's equation.

Effects of confinement and external fields on structure and transport in colloidal dispersions in reduced dimensionality.

In this work, we focus on low-dimensional colloidal model systems, via simulation studies and also some complementary experiments, in order to elucidate the interplay between phase behavior, geometric structures and transport properties. In particular, we try to investigate the (nonlinear!) response of these very soft colloidal systems to various perturbations: uniform and uniaxial pressure, laser fields, shear due to moving boundaries and randomly quenched disorder. We study ordering phenomena on surfaces or in monolayers by Monte Carlo computer simulations of binary hard-disk mixtures, the influence of a substrate being modeled by an external potential. Weak external fields allow a controlled tuning of the miscibility of the mixture. We discuss the laser induced de-mixing for the three different possible couplings to the external potential. The structural behavior of hard spheres interacting with repulsive screened Coulomb or dipolar interaction in 2D and 3D narrow constrictions is investigated using Brownian dynamics simulations. Due to misfits between multiples of the lattice parameter and the channel widths, a variety of ordered and disordered lattice structures have been observed. The resulting local lattice structures and defect probabilities are studied for various cross sections. The influence of a self-organized order within the system is reflected in the velocity of the particles and their diffusive behavior. Additionally, in an experimental system of dipolar colloidal particles confined by gravity on a solid substrate we investigate the effect of pinning on the dynamics of a two-dimensional colloidal liquid. This work contains sections reviewing previous work by the authors as well as new, unpublished results. Among the latter are detailed studies of the phase boundaries of the de-mixing regime in binary systems in external light fields, configurations for shear induced effects at structured walls, studies on the effect of confinement on the structures and defect densities in three-dimensional systems, the effect of confinement and barriers on two-dimensional flow and diffusion, and the effect of pinning sites on the diffusion.

Langevin dynamics simulations of a two-dimensional colloidal crystal under confinement and shear.

Langevin dynamics simulations are used to study the effect of shear on a two-dimensional colloidal crystal (with implicit solvent) confined by structured parallel walls. When walls are sheared very slowly, only two or three crystalline layers next to the walls move along with them, while the inner layers of the crystal are only slightly tilted. At higher shear velocities, this inner part of the crystal breaks into several pieces with different orientations. The velocity profile across the slit is reminiscent of shear banding in flowing soft materials, where liquid and solid regions coexist; the difference, however, is that in the latter case the solid regions are glassy while here they are crystalline. At even higher shear velocities, the effect of the shearing becomes smaller again. Also the effective temperature near the walls (deduced from the velocity distributions of the particles) decreases again when the wall velocity gets very large. When the walls are placed closer together, thereby introducing an incommensurability between the periodicity of the confined crystal and the walls, a structure containing a soliton staircase arises in simulations without shear. Introducing shear increases the disorder in these systems until no solitons are visible anymore. Instead, similar structures like in the case without mismatch result. At high shear rates, configurations where the incommensurability of the crystalline structure is compensated by the creation of holes become relevant.

Numerical approaches to determine the interface tension of curved interfaces from free energy calculations.

A recently proposed method to obtain the surface free energy σ(R) of spherical droplets and bubbles of fluids, using a thermodynamic analysis of two-phase coexistence in finite boxes at fixed total density, is reconsidered and extended. Building on a comprehensive review of the basic thermodynamic theory, it is shown that from this analysis one can extract both the equimolar radius R(e) as well as the radius R(s) of the surface of tension. Hence the free energy barrier that needs to be overcome in nucleation events where critical droplets and bubbles are formed can be reliably estimated for the range of radii that is of physical interest. It is found that the conventional theory of nucleation, where the interface tension of planar liquid-vapor interfaces is used to predict nucleation barriers, leads to a significant overestimation, and this failure is particularly large for bubbles. Furthermore, different routes to estimate the effective radius-dependent Tolman length δ(R(s)) from simulations in the canonical ensemble are discussed. Thus we obtain an instructive exemplification of the basic quantities and relations of the thermodynamic theory of metastable droplets/bubbles using simulations. However, the simulation results for δ(R(s)) employing a truncated Lennard-Jones system suffer to some extent from unexplained finite size effects, while no such finite size effects are found in corresponding density functional calculations. The numerical results are compatible with the expectation that δ(R(s) → ∞) is slightly negative and of the order of one tenth of a Lennard-Jones diameter, but much larger systems need to be simulated to allow more precise estimates of δ(R(s) → ∞).

Active nonlinear microrheology in a glass-forming Yukawa fluid.

A molecular dynamics computer simulation of a glass-forming Yukawa mixture is used to study the anisotropic dynamics of a single particle pulled by a constant force. Beyond linear response, a scaling regime is found where a force-temperature superposition principle of a Peclet number holds. In the latter regime, the diffusion dynamics perpendicular to the force can be mapped on the equilibrium dynamics in terms of an effective temperature, whereas parallel to the force a superdiffusive behavior is seen in the long-time limit. This behavior is associated with a hopping motion from cage to cage and can be qualitatively understood by a simple trap model.

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects.

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.

Slow process in confined polymer melts: layer exchange dynamics at a polymer solid interface.

Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.

Coarse-graining dipolar interactions in simple fluids and polymer solutions: Monte Carlo studies of the phase behavior.

In this paper we investigate the phase diagram of pure dipolar substances and their mixtures with short alkanes, using grand canonical Monte Carlo simulations of simplified coarse-grained models. Recently, an efficient coarse-grained model for simple quadrupolar molecules, based on a Lennard-Jones (LJ) interaction plus a spherically averaged quadrupolar potential, has been shown to be successful in predicting single-component and mixture phase diagrams. Motivated by these results, we investigate the phase diagrams of simple dipolar molecules (and their mixtures with alkanes) using a spherically averaged potential. First, we test the model on pure components. A generalized (state-dependent) mapping procedure allows us to recycle Monte Carlo results of the simple Lennard-Jones (LJ) potential. Considering ammonia, nitrous oxide, and hydrogen sulfide, we generally observe improvements in the single-component phase diagram compared to a pure LJ description, but also some discrepancies in the coexistence pressure near the critical point and in the liquid branch of the coexistence densities well below criticality. In addition, we present results for mixtures. We consider mixtures of ammonia (NH3) with methane (CH4), nonane (C9H20) and hexadecane (C16H34)--for which experimental results are available--and compare the predictions from this modeling ansatz with predictions from simple LJ models. We also present results for the hydrogen sulfide-pentane mixture (H2S and C5H12) for which big discrepancies between simulations and experiments are present. Possible explanations for these discrepancies and limitations of the modeling are discussed.

Coarse-grained models for fluids and their mixtures: Comparison of Monte Carlo studies of their phase behavior with perturbation theory and experiment.

The prediction of the equation of state and the phase behavior of simple fluids (noble gases, carbon dioxide, benzene, methane, and short alkane chains) and their mixtures by Monte Carlo computer simulation and analytic approximations based on thermodynamic perturbation theory is discussed. Molecules are described by coarse grained models, where either the whole molecule (carbon dioxide, benzene, and methane) or a group of a few successive CH(2) groups (in the case of alkanes) are lumped into an effective point particle. Interactions among these point particles are fitted by Lennard-Jones (LJ) potentials such that the vapor-liquid critical point of the fluid is reproduced in agreement with experiment; in the case of quadrupolar molecules a quadrupole-quadrupole interaction is included. These models are shown to provide a satisfactory description of the liquid-vapor phase diagram of these pure fluids. Investigations of mixtures, using the Lorentz-Berthelot (LB) combining rule, also produce satisfactory results if compared with experiment, while in some previous attempts (in which polar solvents were modeled without explicitly taking into account quadrupolar interaction), strong violations of the LB rules were required. For this reason, the present investigation is a step towards predictive modeling of polar mixtures at low computational cost. In many cases Monte Carlo simulations of such models (employing the grand-canonical ensemble together with reweighting techniques, successive umbrella sampling, and finite size scaling) yield accurate results in very good agreement with experimental data. Simulation results are quantitatively compared to an analytical approximation for the equation of state of the same model, which is computationally much more efficient, and some systematic discrepancies are discussed. These very simple coarse-grained models of small molecules developed here should be useful, e.g., for simulations of polymer solutions with such molecules as solvent.

Spherically averaged versus angle-dependent interactions in quadrupolar fluids.

Employing simplified models in computer simulation is on the one hand often enforced by computer time limitations but on the other hand it offers insights into the molecular properties determining a given physical phenomenon. We employ this strategy to the determination of the phase behavior of quadrupolar fluids, where we study the influence of omitting angular degrees of freedom of molecules via an effective spherically symmetric potential obtained from a perturbative expansion. Comparing the liquid-vapor coexistence curve, vapor pressure at coexistence, interfacial tension between the coexisting phases, etc., as obtained from both the models with the full quadrupolar interactions and the (approximate) isotropic interactions, we find discrepancies in the critical region to be typically (such as in the case of carbon dioxide) of the order of 4%. However, when the Lennard-Jones parameters are rescaled such that critical temperatures and critical densities of both models coincide with the experimental results, almost perfect agreement between the above-mentioned properties of both models is obtained. This result justifies the use of isotropic quadrupolar potentials. We also present a detailed comparison of our simulations with a combined integral equation-density functional approach and show that the latter provides an accurate description except for the vicinity of the critical point.

Efficient prediction of thermodynamic properties of quadrupolar fluids from simulation of a coarse-grained model: the case of carbon dioxide.

Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter q=Q(*4)T*, where Q* is the reduced quadrupole moment of the molecule and T* the reduced temperature. The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density, and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of perturbation theory in the mean spherical approximation. As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures.