Search

We study a passive probe immersed in a fluid of active particles. Despite the system's non-equilibrium nature, the trajectory of the probe does not exhibit non-equilibrium signatures: its velocity distribution remains Gaussian, the second fluctuation dissipation theorem is not fundamentally violated, and the motion does not indicate breaking of time reversal symmetry. To tell that the probe is out of equilibrium requires examination of its behavior in tandem with that of the active fluid: the kinetic temperature of the probe does not equilibrate to that of the surrounding active particles. As a strategy to diagnose non-equilibrium from probe trajectories alone, we propose to examine their response to a small perturbation which reveals a non-equilibrium signature through a violation of the first fluctuation dissipation theorem.

We study the lateral domain structure in a coarse-grained molecular model for multicomponent lipid bilayers by semi-grandcanonical Monte Carlo simulations. The membranes are filled with liquid ordered (lo) domains surrounded by a liquid disordered (ld) matrix. Depending on the membrane composition and temperature, we identify different morphological regimes: one regime (I) where the lo domains are small and relatively compact, and two regimes (II, II') where they are larger and often interconnected. In the latter two regimes, the ld matrix forms a network of disordered trenches separating the lo domains, with a relatively high content of interdigitated line defects. Since such defects are also a structural element of the modulated ripple phase in one component membranes, we argue that the regimes II, II' may be amorphous equivalents of the ripple phase in multicomponent membranes. We also analyze the local structure and provide evidence that the domains in regime I are stabilized by a monolayer curvature mechanism postulated in earlier work [S. Meinhardt et al., PNAS, 2013, 110, 4476].

We propose a dynamic coarse-graining (CG) scheme for mapping heterogeneous polymer fluids onto extremely CG models in a dynamically consistent manner. The idea is to use as target function for the mapping a wave-vector dependent mobility function derived from the single-chain dynamic structure factor, which is calculated in the microscopic reference system. In previous work, we have shown that dynamic density functional calculations based on this mobility function can accurately reproduce the order/disorder kinetics in polymer melts, thus it is a suitable starting point for dynamic mapping. To enable the mapping over a range of relevant wave vectors, we propose to modify the CG dynamics by introducing internal friction parameters that slow down the CG monomer dynamics on local scales, without affecting the static equilibrium structure of the system. We illustrate and discuss the method using the example of infinitely long linear Rouse polymers mapped onto ultrashort CG chains. We show that our method can be used to construct dynamically consistent CG models for homopolymers with CG chain length N=4, whereas for copolymers, longer CG chain lengths are necessary, Comment: to appear in the J. Phys.: Cond. Matter - update of references and small changes in Sec. 2.1

The generalized Langevin equation is a model for the motion of coarse-grained particles where dissipative forces are represented by a memory term. The numerical realization of such a model requires the implementation of a stochastic delay-differential equation and the estimation of a corresponding memory kernel. Here we develop a new approach for computing a data-driven Markov model for the motion of the particles, given equidistant samples of their velocity autocorrelation function. Our method bypasses the determination of the underlying memory kernel by representing it via up to about twenty auxiliary variables. The algorithm is based on a sophisticated variant of the Prony method for exponential interpolation and employs the positive real lemma from model reduction theory to extract the associated Markov model. We demonstrate the potential of this approach for the test case of anomalous diffusion, where data are given analytically, and then apply our method to velocity autocorrelation data of molecular dynamics simulations of a colloid in a Lennard-Jones fluid. In both cases, the velocity autocorrelation function and the memory kernel can be reproduced very accurately. Moreover, we show that the algorithm can also handle input data with large statistical noise. We anticipate that it will be a very useful tool in future studies that involve dynamic coarse-graining of complex soft matter systems.

Soft matter covers a wide range of materials based on linear or branched polymers, gels and rubbers, amphiphilic (macro)molecules, colloids, and self-assembled structures. These materials have applications in various industries, all highly important for our daily life, and they control all biological functions; therefore, controlling and tailoring their properties is crucial. One way to approach this target is defect engineering, which aims to control defects in the material's structure, and/or to purposely add defects into it to trigger specific functions. While this approach has been a striking success story in crystalline inorganic hard matter, both for mechanical and electronic properties, and has also been applied to organic hard materials, defect engineering is rarely used in soft matter design. In this review, we present a survey on investigations on defects and/or defect engineering in nine classes of soft matter composed of liquid crystals, colloids, linear polymers with moderate degree of branching, hyperbranched polymers and dendrimers, conjugated polymers, polymeric networks, self-assembled amphiphiles and proteins, block copolymers and supramolecular polymers. This overview proposes a promising role of this approach for tuning the properties of soft matter.

Using a mixed classical Molecular dynamics (MD)/ab initio simulation scheme combined with a quasi-harmonic approximation, we calculate the linear thermal expansion coefficient αL(T) in vitreous silica glasses. The systems are first cooled down by classical MD simulations. Then they are structurally relaxed by ab initio DFT calculations. The vibrational properties are calculated employing the frozen phonon method, and these results are finally used to calculate the Helmholtz free energy as a function of volume. In agreement with experiments, our simulations predict that αL(T) is negative at low temperatures up to T ≈ 150 K. In this low-temperature regime, the simulation results are in quantitative agreement with experiments. To elucidate the origin of the negative thermal expansion, we analyze in detail the microscopic mode Gruneisen parameters in the system and show that the anomalous behavior of αL(T) can be related to the fact that the Gruneisen parameters for the lowest modes become negative at low temperatures – i.e., the lowest eigenfrequencies become stiffer with increasing volume.

We present a multiscale hybrid particle-field scheme for the simulation of relaxation and diffusion behavior of soft condensed matter systems. It combines particle-based Brownian dynamics and field-based local dynamics in an adaptive sense such that particles can switch their level of resolution on the fly. The switching of resolution is controlled by a tuning function which can be chosen at will according to the geometry of the system. As an application, the hybrid scheme is used to study the kinetics of interfacial broadening of a polymer blend, and is validated by comparing the results to the predictions from pure Brownian dynamics and pure local dynamics calculations., Comment: 10 Pages, 5 Figures

Correction for ‘Generalized Langevin dynamics: construction and numerical integration of non-Markovian particle-based models’ by Gerhard Jung et al., Soft Matter, 2018, DOI: 10.1039/c8sm01817k.

We propose a generalized Langevin dynamics (GLD) technique to construct non-Markovian particle-based coarse-grained models from fine-grained reference simulations and to efficiently integrate them. The proposed GLD model has the form of a discretized generalized Langevin equation with distance-dependent two-particle contributions to the self- and pair-memory kernels. The memory kernels are iteratively reconstructed from the dynamical correlation functions of an underlying fine-grained system. We develop a simulation algorithm for this class of non-Markovian models that scales linearly with the number of coarse-grained particles. Our GLD method is suitable for coarse-grained studies of systems with incomplete time scale separation, as is often encountered, e.g., in soft matter systems. We apply the method to a suspension of nanocolloids with frequency-dependent hydrodynamic interactions. We show that the results from GLD simulations perfectly reproduce the dynamics of the underlying fine-grained system. The effective speedup of these simulations amounts to a factor of about $10^4$. Additionally, the transferability of the coarse-grained model with respect to changes of the nanocolloid density is investigated. The results indicate that the model is transferable to systems with nanocolloid densities that differ by up to one order of magnitude from the density of the reference system., Comment: 16 pages, 6 figures, accepted by Soft Matter (https://pubs.rsc.org/en/content/articlelanding/2018/sm/c8sm01817k#!divAbstract)

Experimental data on thin films of cylinder-forming block copolymers (BC)—free-standing BCmembranes as well as supported BC films—strongly suggest that the local orientation of the BC patternsis coupled to the geometry in which the patterns are embedded. We analyze this phenomenon using generalsymmetry considerations and numerical self-consistent field studies of curved BC films in cylindricalgeometry. The stability of the films against curvature-induced dewetting is also analyzed. In goodagreement with experiments, we find that the BC cylinders tend to align along the direction of curvature athigh curvatures. At low curvatures, we identify a transition from perpendicular to parallel alignment insupported films, which is absent in free-standing membranes. Hence both experiments and theory show thatcurvature can be used to manipulate and align BC patterns. Fil: Vu, Giang Thi. Johannes Gutenberg Universitat Mainz; Alemania Fil: Abate, Anabella Angela. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Física del Sur. Universidad Nacional del Sur. Departamento de Física. Instituto de Física del Sur; Argentina Fil: Gomez, Leopoldo Raimundo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Física del Sur. Universidad Nacional del Sur. Departamento de Física. Instituto de Física del Sur; Argentina Fil: Pezzutti, Aldo Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Física del Sur. Universidad Nacional del Sur. Departamento de Física. Instituto de Física del Sur; Argentina Fil: Register, Richard A.. University of Princeton; Estados Unidos Fil: Vega, Daniel Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Física del Sur. Universidad Nacional del Sur. Departamento de Física. Instituto de Física del Sur; Argentina Fil: Schmid, Friederike. Johannes Gutenberg Universitat Mainz; Alemania

The co-solvency method is a method for the size controlled preparation of nanoparticles like polymersomes, where a poor co-solvent is mixed into a homogeneous copolymer solution to trigger precipitation of the polymer. The size of the resulting particles is determined by the rate of co-solvent addition. We use the Cahn-Hilliard equation with a Flory-Huggins free energy model to describe the precipitation of a polymer under changing solvent quality by applying a time dependent Flory-Huggins interaction parameter. The analysis focuses on the characteristic size R of polymer aggregates that form during the initial spinodal decomposition stage, and especially on how R depends on the rate s of solvent quality change. Both numerical results and a perturbation analysis predict a power law dependence R ∼ s−⅙, which is in agreement with power laws for the final particle sizes that have been reported from experiments and molecular dynamics simulations. Hence, our model results suggest that the nanoparticle size in size-controlled precipitation is essentially determined during the spinodal decomposition stage.

We have extended an existing hybrid MD-SCF simulation technique that employs a coarsening step to enhance the computational efficiency of evaluating non-bonded particle interactions. This technique is conceptually equivalent to the single chain in mean-field (SCMF) method in polymer physics, in the sense that non-bonded interactions are derived from the non-ideal chemical potential in self-consistent field (SCF) theory, after a particle-to-field projection. In contrast to SCMF, however, MD-SCF evolves particle coordinates by the usual Newton's equation of motion. Since collisions are seriously affected by the softening of non-bonded interactions that originates from their evaluation at the coarser continuum level, we have devised a way to reinsert the effect of collisions on the structural evolution. Merging MD-SCF with multi-particle collision dynamics (MPCD), we mimic particle collisions at the level of computational cells and at the same time properly account for the momentum transfer that is important for a realistic system evolution. The resulting hybrid MD-SCF/MPCD method was validated for a particular coarse-grained model of phospholipids in aqueous solution, against reference full-particle simulations and the original MD-SCF model. We additionally implemented and tested an alternative and more isotropic finite difference gradient. Our results show that efficiency is improved by merging MD-SCF with MPCD, as properly accounting for hydrodynamic interactions considerably speeds up the phase separation dynamics, with negligible additional computational costs compared to efficient MD-SCF. This new method enables realistic simulations of large-scale systems that are needed to investigate the applications of self-assembled structures of lipids in nanotechnologies.

We review recent computer simulation studies of undulating lipid bilayers. Theoretical interpretations of such fluctuating membranes are most commonly based on generalized Helfrich-type elastic models, with additional contributions of local "protrusions" and/or density fluctuations. Such models provide an excellent basis for describing the fluctuations of tensionless bilayers in the fluid phase at a quantitative level. However, this description is found to fail for membranes in the gel phase and for membranes subject to high tensions. The fluctuations of tilted gel membranes show a signature of the modulated ripple structure, which is a nearby phase observed in the pretransition regime between the fluid and tilted gel state. This complicates a quantitative analysis on mesoscopic length scales. In the case of fluid membranes under tension, the large-wavelength fluctuation modes are found to be significantly softer than predicted by theory. In the latter context, we also address the general problem of the relation between frame tension and the fluctuation tension, which has been discussed somewhat controversially in recent years. Simulations of very simple model membranes with fixed area show that the fluctuations should be controlled by the frame tension, and not by the internal tension., Comment: Invited Review article for Biophysical Reviews and Letters (written in 2012)

We review and compare recent work on the properties of fluctuating interfaces between nematic and isotropic liquid-crystalline phases. Molecular dynamics and Monte Carlo simulations have been carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, and the fluctuation spectrum of interface positions (the capillary wave spectrum) has been analyzed. In addition, the capillary wave spectrum has been calculated analytically within the Landau-de Gennes theory. The theory predicts that the interfacial fluctuations can be described in terms of a wave vector dependent interfacial tension, which is anisotropic at small wavelengths (stiff director regime) and becomes isotropic at large wavelengths (flexible director regime). After determining the elastic constants in the nematic phase, theory and simulation can be compared quantitatively. We obtain good agreement for the stiff director regime. The crossover to the flexible director regime is expected at wavelengths of the order of several thousand particle diameters, which was not accessible to our simulations.

Fibrinogen, a blood glycoprotein of vertebrates, plays an essential role in blood clotting by polymerizing into fibrin when activated. Upon adsorption on material surfaces, it also contributes to determine their biocompatibility and has been implicated in the onset of thrombosis and inflammation at medical implants. Here we present the first fully atomistic simulations of the initial stages of the adsorption process of fibrinogen on mica and graphite surfaces. The simulations reveal a weak adsorption on mica that allows frequent desorption and reorientation events. This adsorption is driven by electrostatic interactions between the protein and the silicate surface as well as the counterion layer. Preferred adsorption orientations for the globular regions of the protein are identified. The adsorption on graphite is found to be stronger with fewer reorientation and desorption events and shows the onset of denaturation of the protein.

Applying electric fields is an attractive way to control and manipulate single particles or molecules, e.g., in lab-on-a-chip devices. However, the response of nanosize objects in electrolyte solution to external fields is far from trivial. It is the result of a variety of dynamical processes taking place in the ion cloud surrounding charged particles and in the bulk electrolyte, and it is governed by an intricate interplay of electrostatic and hydrodynamic interactions. Already systems composed of one single particle in electrolyte solution exhibit a complex dynamical behaviour. In this review, we discuss recent coarse-grained simulations that have been performed to obtain a molecular-level understanding of the dynamic and dielectric response of single particles and single macromolecules to external electric fields. We address both the response of charged particles to constant fields (DC fields), which can be characterized by an electrophoretic mobility, and the dielectric response of both uncharged and charged particles to alternating fields (AC fields), which is described by a complex polarizability. Furthermore, we give a brief survey of simulation algorithms and highlight some recent developments., Comment: short review, submitted to Soft Matter

The effect of rough walls on the phase behaviour of a confined liquid crystal (LC) fluid is studied using constant pressure Monte Carlo simulations. The LC is modelled as a fluid of soft ellipsoidal molecules and the rough walls are represented as a hard wall with a number of molecules randomly embedded in them. It is found that the isotropic-nematic (IN) transition is shifted to higher pressures for rougher walls., 4 pages, 4 figures Accepted in Chemical Physics Letters

Hydrodynamic interactions play an important role in many areas of soft matter science. In simulations with implicit solvent, various techniques such as Brownian or Stokesian dynamics explicitly include hydrodynamic interactions a posteriori by using hydrodynamic diffusion tensors derived from the Stokes equation. However, this equation assumes the interaction to be instantaneous which is an idealized approximation and only valid on long time scales. In the present paper, we go one step further and analyze the time-dependence of hydrodynamic interactions in a compressible fluid on the basis of the linearized Navier-Stokes equation. The theoretical results show that the compressibility of the fluid has a significant impact on frequency-dependent pair interactions. The predictions of the hydrodynamic theory are compared to molecular dynamics simulations of two solid spheres in a Lennard-Jones fluid. For this system we reconstruct memory functions by extending the inverse Volterra technique. The simulation data agree very well with the theory, therefore, the theory can be used to implement dynamically consistent hydrodynamic interactions in the increasingly popular field of non-Markovian modeling., Comment: 12 pages, 15 figures, published by Physics of Fluids

Using Brownian Dynamics, we study the dynamical behavior of a polymer grafted onto an adhesive surface close to the mechanically induced adsorption-stretching transition. Even though the transition is first order, (in the infinite chain length limit, the stretching degree of the chain jumps discontinuously), the characteristic relaxation time is found to grow according to a power law as the transition point is approached. We present a dynamic effective interface model which reproduces these observations and provides an excellent quantitaive description of the simulations data. The generic nature of the theoretical model suggests that the unconventional mixing of features that are characteristic for first-order transitions (a jump in an order parameter) and features that are characteristic of critical points (anomalous slowdown) may be a common phenomenon in force-driven phase transitions of macromolecules., 5 figures

We briefly review recent theoretical and simulation studies of charged colloidal dispersions in alternating electric fields (AC fields). The response of single colloid to an external field can be characterized by a complex polarizability, which describes the dielectric properties of the colloid and its surrounding electrical double layer. We present computer simulation studies of single spherical colloid, using a coarse-grained mesoscale approach that accounts in full for hydrodynamic and electrostatic interactions as well as for thermal fluctuations. We investigate systematically a number of controlling parameters, such as the amplitude and frequency of the AC-fields. The results are in good agreement with recent theoretical predictions., Comment: 12 pages, 4 figures

Self-assembled monolayers of organic molecules on solid substrates build well-defined surfaces, which can be structured efficiently on a nanometre scale. The molecular order of the molecules, in particular their tilt with respect to the surface, is believed to be important for the stability of those patterns. We discuss Monte Carlo simulations of a simple, idealized model for self-assembled monolayers, with particular emphasis on the different types of tilt order and tilting transitions.

We present Monte Carlo simulations of surface induced disordering in a model of a binary alloy on a bcc lattice which undergoes a first order bulk transition from the ordered DO3 phase to the disordered A2 phase. The data are analyzed in terms of an effective interface Hamiltonian for a system with several order parameters in the framework of the linear renormalization approach due to Brezin, Halperin and Leibler. We show that the model provides a good description of the system in the vicinity of the interface. In particular, we recover the logarithmic divergence of the thickness of the disordered layer as the bulk transition is approached, we calculate the critical behavior of the maxima of the layer susceptibilities, and demonstrate that it is in reasonable agreement with the simulation data. Directly at the (110) surface, the theory predicts that all order parameters vanish continuously at the surface with a nonuniversal, but common critical exponent. However, we find different exponents for the order parameter of the DO3 phase and the order parameter of the B2 phase. Using the effective interface model, we derive the finite size scaling function for the surface order parameter and show that the theory accounts well for the finite size behavior of the DO3 ordering but not for that of B2 ordering. The situation is even more complicated in the neighborhood of the (100) surface, due to the presence of an ordering field which couples to the B2 order., To appear in Physical Review B

Symmetric binary polymer mixtures are studied by Monte Carlo simulation of the bond fluctuation model, considering both interfaces between coexisting bulk phases and interfaces confined in thin films. It is found that the critical behavior of interfacial tension and width is compatible with that of the Ising model, as expected from the universality principle. In the strong segregation limit, only qualitative but not quantitative agreement with the self-consistent field (SCF) theory is found. It is argued that the SCF theory requires √ 6 X √D for short-range forces, in agreement with experiment.

We consider an interface between two demixed A and B phases, confined in a thin film between two antisymmetric walls, one of which prefers A and the other B. Above the wetting transition, the interface is stabilized in the center of the film. Based on a suitable extension of the Cahn-Hilliard gradient-square theory in combination with the capillary wave theory, we argue that the confinement influences the interfacial structure in two ways: It squeezes the intrinsic structure and cuts off the capillary wave spectrum. As a result, the interfacial width is proportional to the film thickness D in very thin films, and proportional to the square root of D in thicker films. These effects are then discussed in detail for the special case of an interface between demixed homopolymer phases. The width of the intrinsic profile is calculated analytically as a function of film thickness in the Cahn-Hilliard approximation (weak segregation limit) and in the Helfand theory (strong segregation limit), and numerically in the self-consistent field approximation. The results are confirmed by Monte Carlo simulations of a lattice model for a polymer mixture.

A generalized self-consistent field approach for polymer networks with fixed topology is developed. It is shown that the theory reproduces the localization of crosslinks which is characteristic for gels. The theory is then used to study the order-disorder transition in regular networks of endlinked diblock copolymers. Compared to diblock copolymer melts, the transition is shifted towards lower values of the incompatibility parameter $\chi$ (the Flory- Huggins parameter). Moreover, the transition becomes strongly first order already at the mean-field level. If stress is applied, the transition is further shifted and finally vanishes in a critical point., Comment: Includes supplemental information

Wetting transitions in balanced amphiphilic systems containing only short-range forces are studied. In the absence of fluctuations, a weak continuous wetting transition would take place at three-phase coexistence at the disorder line. We argue that capillary wave fluctuations cause the transition to become first order and to occur on the structured side of the disorder line. With reasonable model parameters, the transition occurs in the vicinity of the Lifshitz line where it is observed experimentally.

We study the response of single nanosized spherical colloids in electrolyte solution to an alternating electric field (AC field) by computer simulations. We use a coarse-grained mesoscopic simulation approach that accounts in full for hydrodynamic and electrostatic interactions as well as for thermal fluctuations. The solvent is modeled as a fluid of single Dissipative Particle Dynamics (DPD) beads, and the colloidal particle is modeled as a rigid body made of DPD beads. We compute the mobility and the polarizability of a single colloid and investigate systematically the effect of amplitude and frequency of the AC-fields. Even though the thickness of the Debye layer is not "thin" compared to the radius of the colloid, and the thermal fluctuations are significant, the results are in good agreement with the theoretical prediction of the Maxwell-Wagner-O'Konski theory, especially for uncharged colloids., Comment: 22 pages, 11 figure, for proceedings of CODEF III conference

Free (100)-surfaces of body-centered cubic binary alloys are studied in a parameter range where the bulk turns from the ordered B2-phase to the disordered A2-phase. A model is chosen that describes iron-aluminium alloys in a fairly realistic way. Mean field treatments and Monte Carlo investigations both show that under certain circumstances the surface remains ordered far above the bulk disordering temperatureT c, though the surface order parameter and the surface susceptibility exhibit a singularity atT c with critical exponents characteristic for the ordinary transition. One finds, that if the surface is nonstoechiometric and different layers are not equivalent with respect to perfect bulk ordering, this induces an effective ordering field at the surface. Finally surface potentials are introduced into the model. Realistic values are determined by mean field considerations from experimental data and their effect on the surface is discussed.

The structure of antiphase domain boundaries in the ordered phase of body-centered-cubic alloys is investigated by Monte Carlo simulations, studying a fairly realistic model of iron-aluminum alloys that includes a description of the magnetic properties of iron within a nearest-neighbor Heisenberg Hamiltonian. An interface is enforced by employing in one lattice direction an odd number of lattice planes and choosing periodic boundary conditions throughout. At the temperatures of interest, the antiphase domain boundary is not localized in the system and is rough, which implies a logarithmic dependence of the interfacial width on the linear dimension ${\mathit{L}}_{\mathrm{\ensuremath{\parallel}}}$ parallel to the interface. Methods are introduced to define an instantaneous interface position for each configuration that is analyzed. Choosing a coordinate frame whose origin is fixed at this (moving) interface position, profiles of order parameter, concentration, etc., can be estimated. It is found that capillary wave theory accounts for the interface widths reasonably well. Near the bulk critical temperature the lengths characterizing the intrinsic interfacial widths agree with the bulk correlation length, which is obtained from an independent investigation of correlation functions in the bulk. Finally, the interfacial enrichment of the majority species in the domain wall is studied and it is shown that the critical behavior of the interfacial excess concentration is consistent with scaling predictions.

We reconsider the mobility of hydrogen in bcc metals, in particular H and D in Nb, and its influence on the coherent and incoherent quasielastic neutron and x-ray scattering. An approach to describe time-dependent lattice distortions around highly mobile protons is presented, using the frequency-dependent Green's function of the host lattice so that effects of the so-called quasielastic effects and resonances with host lattice phonons are included. We discuss within this formalism how noncubic components of the long-ranged displacement field can fade away in the presence of a high defect mobility and we demonstrate by a detailed calculation that the experimentally observed disappearance of the coherent quasielastic diffuse (x-ray and neutron) scattering in ${\mathrm{NbH}}_{\mathit{x}}$ (${\mathrm{NbD}}_{\mathit{x}}$) from these Fourier components of the displacement field points to an unusual, hitherto neglected rapid local diffusion process between adjacent tetrahedral sites with a short-time constant of ${\mathrm{\ensuremath{\tau}}}_{\mathrm{loc}}$\ensuremath{\simeq}${10}^{\mathrm{\ensuremath{-}}13}$--${10}^{\mathrm{\ensuremath{-}}14}$ s. The consequences of such a diffusion mechanism on the incoherent quasielastic neutron scattering are discussed in detail. We show that the elastic incoherent structure factor of the interstitial proton exhibits, in addition to the usual Debye-Waller factor behavior, distinct oscillations as a function of the neutron momentum transfer.A neutron-scattering study of ${\mathrm{NbH}}_{0.02}$ at various temperatures between room temperature and 573 K is presented. The experimental results of the elastic incoherent structure factor displays, after a proper correction for the proton Debye-Waller factor, a pronounced oscillatory behavior in the [100] and [111] direction as predicted for a rapid local diffusion on neighboring tetrahedral sites. It is shown that this property of the proton form factor is temperature dependent. We present measurements of the quasielastic broadening of the coherent diffuse neutron scattering from the time-dependent long-ranged displacement field (``Huang linewidth'') and discuss further implications of our findings on the interpretation of structural and dynamical properties of H in Nb, as the local modes and the coherent diffuse neutron scattering from the lattice distortions around the local defect neighborhood.

Information from recent X-ray and neutron diffuse scattering investigations of Fe-Al alloys and the resulting effective interaction parameters are used to construct a model in which the lattice configurational degrees of freedom of these alloys are described by Ising spins, and the magnetic moments of iron atoms are described by classical Heisenberg spins. Starting from the exact ground state of the model, information on the phase diagram at non-zero temperature is obtained from both spin wave theory, mean-field approximations and extensive Monte Carlo simulations. It is shown that use of the unrenormalized interaction parameters from the various experiments yields a phase diagram which disagrees severely with experiment even at the qualitative level. However, with a suitable renormalization of the ratio between magnetic and non-magnetic interactions, a qualitatively reasonable phase diagram of the model is obtained. The consequences for the authors' understanding of effective interaction parameters in Fe-Al alloys are spelt out, and the implications for experimental work are discussed.

We analyze the interplay between wetting and anchoring of nematic liquid crystals on disordering, e.g., rough substrates in the framework of the Landau-de Gennes theory, in situations of competing homeotropic and planar easy axes on the substrate and the nematic-isotropic (NI) interface. The phase diagram for azimuthally symmetric substrates is calculated. We identify two regimes - a strongly coupled regime, where the wetting transition coincides with an anchoring transition, and a weakly coupled regime, where the two are separated. The anchoring transition is first order and switches between homeotropic and planar anchoring. The two competing orientations are metastable over a broad parameter range. Hence such surfaces can be used to generate bistable surfaces.

A thermostat of the Nose-Hoover type, based on relative velocities and a local definition of the temperature, is presented. The thermostat is momentum-conserving and Galilean-invariant, which should make it suitable for use in Dissipative Particle Dynamics simulations, as well as nonequilibrium molecular dynamics simulations., Comment: Symposium on Progress and Future Prospects in Molecular Dynamics Simulation - In Memory of Professor Shuichi Nose

Nematic liquid crystals at rough and fluctuating interfaces are analyzed within the Frank elastic theory and the Landau-de Gennes theory. We study specifically interfaces that locally favor planar anchoring. In the first part we reconsider the phenomenon of Berreman anchoring on fixed rough surfaces, and derive new simple expressions for the corresponding azimuthal anchoring energy. Surprisingly, we find that for strongly aligning surfaces, it depends only on the geometrical surface anisotropy and the bulk elastic constants, and not on the precise values of the chemical surface parameters. In the second part, we calculate the capillary waves at nematic-isotropic interfaces. If one neglects elastic interactions, the capillary wave spectrum is characterized by an anisotropic interfacial tension. With elastic interactions, the interfacial tension, i.e. the coefficient of the leading $q^2$ term of the capillary wave spectrum, becomes isotropic. However, the elastic interactions introduce a strongly anisotropic {\em cubic} $q^3$ term. The amplitudes of capillary waves are largest in the direction perpendicular to the director. These results are in agreement with previous molecular dynamics simul, to appear in European Physical Journal E

The interfacial width of polyolefins blends has been probed as a function of distance away from the critical point by using neutron reflectivity. For strongly immiscible polymer pairs, the width of the interface increases slowly when the degree of immiscibility is decreased and the interfacial width varies with the interaction parameter $\ensuremath{\chi}$ of the polymers. Closer to the critical point the dependence on the degree of miscibility becomes stronger and the way in which the interfacial width diverges, as criticality is approached, is related to both the chain length and $\ensuremath{\chi}$. The self-consistent field theory numerical calculations, with the additional contribution due to capillary waves, provides a good description of the width of the interface between two polymer bulk phases in particular at intermediate values of the degree of immiscibility.

We present a large-scale molecular dynamics study of nematic-paranematic interfaces under shear. We use a model of soft repulsive ellipsoidal particles with well-known equilibrium properties, and consider interfaces which are oriented normal to the direction of the shear gradient (common stress case). The director at the interface is oriented parallel to the interface (planar). A fixed average shear rate is imposed with Lees-Edwards boundary conditions, and the heat is dissipated with a profile-unbiased thermostat. First we study the properties of the interface at one particular shear rate in detail. The local interfacial profiles and the capillary wave fluctuations of the interfaces are calculated and compared with those of the corresponding equilibrium interface. Under shear, the interfacial width broadens and the capillary wave amplitudes at large wavelengths increase. The strain is distributed inhomogeneously in the system (shear banding), the local shear rate in the nematic region being distinctly higher than in the paranematic region. Surprisingly, we also observe (symmetry breaking) flow in the {\em vorticity} direction, with opposite direction in the nematic and the paranematic state. Finally, we investigate the stability of the interface for other shear rates and construct a nonequilibrium phase diagram., Comment: to appear in J. Chem. Phys

By computer simulations of systems of ellipsoids, we study the influence of the isotropic/nematic phase transition on the direct correlation functions (DCF) in anisotropic fluids. The DCF is determined from the pair distribution function by solving the full Ornstein-Zernike equation, without any approximations. Using a suitable molecular-fixed reference frame, we can distinguish between two qualitatively different contributions to the DCF: One which preserves rotational invariance, and one which breaks it and vanishes in the isotropic phase. We find that the symmetry preserving contribution is barely affected by the phase transition. However, symmetry breaking contributions emerge in the nematic phase and may become quite substantial. Thus the DCF in a nematic fluid is not rotationally invariant. In the isotropic fluid, the DCF is in good agreement with the prediction of the Percus-Yevick theory., Comment: to appear in J. Chem. Phys

Systems that exhibit a first-order phase transition in the bulk, such as binary alloys where the order parameter vanishes discontinuously at some critical value of a control parameter, may show a continuous vanishing of the local order parameter at the surface. This “surface-induced disordering” is described theoretically as a variant of critical wetting, where an interface between the locally disordered surface and the ordered bulk gradually moves towards the bulk. We test this description by Monte Carlo simulations for a body centered cubic model alloy, with interactions between nearest and next nearest neighbors, for which the phase diagram in the bulk has been calculated very accurately. A critical vanishing of the order parameter components is found both for the (110) and the (100) surface. In contrast to the theory, the simulations indicate different critical exponents v from the order parameter of the B2 phase and the order parameter of the DO3 phase. Observations from simulations for the face centered cubic lattice and from experiments are also briefly discussed.

We propose a density functional for anisotropic fluids of hard body particles. It interpolates between the well-established geometrically based Rosenfeld functional for hard spheres and the Onsager functional for elongated rods. We test the new approach by calculating the location of the the nematic-isotropic transition in systems of hard spherocylinders and hard ellipsoids. The results are compared with existing simulation data. Our functional predicts the location of the transition much more accurately than the Onsager functional, and almost as good as the theory by Parsons and Lee. We argue that it might be suited to study inhomogeneous systems., Comment: To appear in J. Physics: Condensed Matter

We present large-scale molecular dynamics simulations of a nematic-isotropic interface in a system of repulsive ellipsoidal molecules, focusing in particular on the capillary wave fluctuations of the interfacial position. The interface anchors the nematic phase in a planar way, i.e., the director aligns parallel to the interface. Capillary waves in the direction parallel and perpendicular to the director are considered separately. We find that the spectrum is anisotropic, the amplitudes of capillary waves being larger in the direction perpendicular to the director. In the long wavelength limit, however, the spectrum becomes isotropic and compares well with the predictions of a simple capillary wave theory., to appear in Phys. Rev. E

We present the first calculations of the pressure tensor profile in the vicinity of the planar interface between isotropic liquid and nematic liquid crystal, using Onsager's density functional theory and computer simulation. When the liquid crystal director is aligned parallel to the interface, the situation of lowest free energy, there is a large tension on the nematic side of the interface and a small compressive region on the isotropic side. By contrast, for perpendicular alignment, the tension is on the isotropic side. There is excellent agreement between theory and simulation both in the forms of the pressure tensor profiles, and the values of the surface tension., Comment: Minor changes; to appear in Phys. Rev. E

We present recent Monte Carlo results on surfaces of bcc-structured binary alloys which undergo an order-disorder phase transformation in the bulk. In particular, we discuss surface order and surface induced disorder at the bulk transition between the ordered (DO3) phase and the disordered (A2) phase. An intricate interplay between different ordering and segregation phenomena leads to a complex surface behavior, which depends on the orientation of the surface under consideration.

Recent developments in off-lattice self-consistent field theories for inhomogeneous complex fluids are reviewed. Particular emphasis is given to the treatment of intermolecular interactions and compressibility, to the role of fluctuations, and to the discussion of the coarse-graining length which is inherent to the theory. Valuable insight can be gained from the comparison of self-consistent field calculations with Monte Carlo simulations. Finally, some applications of the theory to orientational properties of polymers and copolymers at interfaces, and to the phase behavior of amphiphiles at surfaces, are presented., Comment: Review article, revised version with changes at the end of sctn. 2.4

The structure of polymer coils near interfaces between coexisting phases of symmetrical polymer mixtures (AB) is discussed, as well as the structure of symmetric diblock copolymers of the same chain length N adsorbed at the interface. The problem is studied by Monte Carlo simulations of the bond fluctuation model on the simple cubic lattice, using massively parallel computers (CRAY T3D). While homopolymer coils in the strong segregation limit are oriented parallel to the interface, the diblocks form ``dumbbells'' oriented perpendicular to the interface. However, in the dilute case (``mushroom regime'' rather than ``brush regime''), the diblocks are only weakly stretched. Distribution functions for monomers at the chain ends and in the center of the polymer are obtained, and a comparison to the self consistent field theory is made., to appear in Physica A

A three-component model representing a balanced, amphiphilic system at a three-phase coexistence of oil-, water-rich, and disordered phases is simulated, and interfacial tensions between all phases are obtained from histogram methods. We find that the disordered phase wets the oil-water interface, and that its tension is reduced due to the presence of the amphiphile by a factor of 93. The determination of the structure function shows that the system is between disorder and Lifshitz lines. These three factors are all consistent with the simulated amphiphile being weak. Fluctuations do not reduce the tension significantly from mean-field estimates, but do alter the location of the wetting transition. The latter is in accord with the expected behavior of a system in which the interfacial tension is small.

We illustrate the use of Monte Carlo simulations in the study of order-disorder phenomena in metallic alloys by presenting detailed work on a fairly realistic lattice model for iron aluminum. The model has been constructed based on recent measurements of effective interaction parameters and includes a description of the magnetism of iron within a Heisenberg Hamiltonian. We show that it reproduces the bulk phase diagram in a qualitatively correct way. Then internal antiphase boundaries and free surfaces in the (100)-direction are studied. An interfacial roughening transition is predicted as well as critical broadening of the profiles as the bulk approaches a second order transition. Nonstoechiometry with respect to perfect bulk ordering is shown to induce interfacial enrichment and under certain circumstances an ordering surface field.

The effect of confinement on a microemulsion, a disordered, fluid phase with two characteristic length scales, \ensuremath{\xi} the correlation length, and d the wavelength of oil and water density variations, is studied within a simple Ginzburg-Landau theory. We find a qualitative difference between those systems which contain a strong amphiphile as opposed to a weak one. As the distance between walls is increased in the former, a series of first-order transitions occurs which, in principle, can continue without limit, while in the latter, only a small number of such transitions are expected. Such transitions will be manifest in discontinuities in the forces between the walls. In our model calculation, the boundary between strong and weak amphiphiles is 2\ensuremath{\pi}\ensuremath{\xi}/d= \ensuremath{\surd}3 , a value easily accessible experimentally.