Search

Sodium chloride (NaCl) is one of the key components of atmospheric aerosols. The surface tension of aqueous NaCl solution (σNaCl,sol) and its concentration dependence are essential to determine the equilibrium water vapor pressure of aqueous NaCl droplets. Supersaturated NaCl solution droplets are observed in laboratory experiments and under atmospheric conditions, but the experimental data for σNaCl,sol are mostly limited up to subsaturated solutions. In this study, the surface tension of aqueous NaCl is investigated by molecular dynamics (MD) simulations and the pressure tensor method from dilute to highly supersaturated solutions. We show that the linear approximation of concentration dependence of σNaCl,sol at molality scale can be extended to the supersaturated NaCl solution until a molality of ∼10.7 mol kg−1 (i.e., solute mass fraction (xNaCl) of ∼0.39). Energetic analyses show that this monotonic increase in surface tension is driven by the increase in excess surface enthalpy (ΔH) as the solution becomes concentrated. After that, the simulated σNaCl,sol remains almost unchanged until xNaCl of ∼0.47 (near the concentration upon efflorescence). The existence of the “inflection point” at xNaCl of ∼0.39 and the stable surface tension of xNaCl between ∼0.39 and ∼0.47 can be attributed to the nearly unchanged excess surface entropy term (T⋅ΔS) and the excess surface enthalpy term (ΔH). After a “second inflection point” at xNaCl of ∼0.47, the simulated σNaCl,sol gradually regains the growing momentum with a tendency to approach the surface tension of molten NaCl (∼175.58 mN m−1 at 298.15 K, MD simulation-based extrapolation). This fast increase in σNaCl,sol at xNaCl>0.47 is a process driven by excess surface enthalpy and excess surface entropy. Our results reveal different regimes of concentration dependence of the surface tension of aqueous NaCl at 298.15 K: a water-dominated regime (xNaCl from 0 to ∼0.39), a transition regime (xNaCl from ∼0.39 to ∼0.47) and a molten NaCl-dominated regime (xNaCl from ∼0.47 to 1).

When physicians are asked to determine the positive predictive value from the a priori probability of a disease and the sensitivity and false positive rate of a medical test (Bayesian reasoning), it often comes to misjudgments with serious consequences. In daily clinical practice, however, it is not only important that doctors receive a tool with which they can correctly judge—the speed of these judgments is also a crucial factor. In this study, we analyzed accuracy and efficiency in medical Bayesian inferences. In an empirical study we varied information format (probabilities vs. natural frequencies) and visualization (text only vs. tree only) for four contexts. 111 medical students participated in this study by working on four Bayesian tasks with common medical problems. The correctness of their answers was coded and the time spent on task was recorded. The median time for a correct Bayesian inference is fastest in the version with a frequency tree (2:55 min) compared to the version with a probability tree (5:47 min) or to the text only versions based on natural frequencies (4:13 min) or probabilities (9:59 min).The score diagnostic efficiency (calculated by: median time divided by percentage of correct inferences) is best in the version with a frequency tree (4:53 min). Frequency trees allow more accurate and faster judgments. Improving correctness and efficiency in Bayesian tasks might help to decrease overdiagnosis in daily clinical practice, which on the one hand cause cost and on the other hand might endanger patients’ safety.

Semiflexible polymers bound to planar substrates by a short-range surface potential are studied by Molecular Dynamics simulations to clarify the extent to which these chain molecules can be considered as strictly two-dimensional. Applying a coarse-grained bead-spring model, the chain length N and stiffness &kappa, as well as the strength of the adsorption potential ϵ w a l l are varied over a wide range. The excluded-volume (EV) interactions inherent in this model can also be &ldquo, switched off&rdquo, to provide a discretized version of the Kratky&ndash, Porod wormlike chain model. We study both local order parameters (fraction f of monomers within the range of the potential, bond-orientational order parameter &eta, ) and the mean square gyration radius parallel, &lang, R g 2 &rang, | , and perpendicular, &lang, &perp, to the wall. While for strongly adsorbed chains EV has negligible effect on f and &eta, we find that &lang, | is strongly affected when the chain contour length exceeds the persistence length. Monomer coordinates in perpendicular (&perp, ) direction are correlated over the scale of the deflection length which is estimated. It is found that f , &eta, and &lang, converge to their asymptotic values with 1 / N corrections. For both weakly and strongly adsorbed chains, the distribution functions of &ldquo, loops&rdquo, &ldquo, trains&rdquo, and &ldquo, tails&rdquo, are analyzed. Some consequences pertaining to the analysis of experiments on adsorbed semiflexible polymers are pointed out.

Conformations of partially or fully adsorbed semiflexible polymer chains are studied varying both contour length $L$, chain stiffness, $\ensuremath{\kappa}$, and the strength of the adsorption potential over a wide range. Molecular dynamics simulations show that partially adsorbed chains (with ``tails,'' surface attached ``trains,'' and ``loops'') are not described by the Kratky-Porod wormlike chain model. The crossover of the persistence length from its three-dimensional value (${\ensuremath{\ell}}_{p}$) to the enhanced value in two dimensions ($2{\ensuremath{\ell}}_{p}$) is analyzed, and excluded volume effects are identified for $L\ensuremath{\gg}{\ensuremath{\ell}}_{p}$. Consequences for the interpretation of experiments are suggested. We verify the prediction that the adsorption threshold scales as ${\ensuremath{\ell}}_{p}^{\ensuremath{-}1/3}$.

Changing the information format from probabilities into frequencies as well as employing appropriate visualizations such as tree diagrams or 2 × 2 tables are important tools that can facilitate people’s statistical reasoning. Previous studies have shown that despite their widespread use in statistical textbooks, both of those visualization types are only of restricted help when they are provided with probabilities, but that they can foster insight when presented with frequencies instead. In the present study, we attempt to replicate this effect and also examine, by the method of eye tracking, why probabilistic 2 × 2 tables and tree diagrams do not facilitate reasoning with regard to Bayesian inferences (i.e., determining what errors occur and whether they can be explained by scan paths), and why the same visualizations are of great help to an individual when they are combined with frequencies. All ten inferences of N = 24 participants were based solely on tree diagrams or 2 × 2 tables that presented either the famous “mammography context” or an “economics context” (without additional textual wording). We first asked participants for marginal, conjoint, and (non-inverted) conditional probabilities (or frequencies), followed by related Bayesian tasks. While solution rates were higher for natural frequency questions as compared to probability versions, eye-tracking analyses indeed yielded noticeable differences regarding eye movements between correct and incorrect solutions. For instance, heat maps (aggregated scan paths) of distinct results differed remarkably, thereby making correct and faulty strategies visible in the line of theoretical classifications. Moreover, the inherent structure of 2 × 2 tables seems to help participants avoid certain Bayesian mistakes (e.g., “Fisherian” error) while tree diagrams seem to help steer them away from others (e.g., “joint occurrence”). We will discuss resulting educational consequences at the end of the paper.

Semiflexible polymers in concentrated lyotropic solution are studied within a bead-spring model by molecular dynamics simulations, focusing on the emergence of a smectic A phase and its properties. We systematically vary the density of the monomeric units for several contour lengths that are taken smaller than the chain persistence length. The difficulties concerning the equilibration of such systems and the choice of appropriate ensemble (constant volume versus constant pressure, where all three linear dimensions of the simulation box can fluctuate independently) are carefully discussed. Using HOOMD-blue on graphics processing units, systems containing more than a million monomeric units are accessible, making it possible to distinguish the order of the phase transitions that occur. While in this model the nematic-smectic transition is continuous, the transition from the smectic phase to a related crystalline structure with true three-dimensional long-range order is clearly of first order. Further, both orientational and positional correlations of monomeric units are studied as well as the order parameters characterizing the nematic, smectic A, and crystalline phases. The analogy between smectic order and one-dimensional harmonic crystals with respect to the behavior of the structure factor is also explored. Finally, the results are put in perspective with pertinent theoretical predictions and possible experiments.

The ordering of semiflexible polymers with persistence length lp and contour length L confined in a sphere of radius R is studied by molecular dynamics simulations of a coarse-grained model. Monomer densities are chosen where the corresponding bulk lyotropic solution or melt is a well-ordered nematic, and purely repulsive walls of the rigid confining sphere are considered. It is found that polymers close to the walls are bent according to the curvature of the confining spheres with all their monomers in a few layers parallel to the sphere surface, whereas the remaining macromolecules closer to the sphere center have one chain end and their center of mass far from the surface. The latter chains are responsible for the average nematic order of the system for sufficiently stiff chains. If L exceeds R, a single transition from isotropic to nematic states is found when lp increases. However, if R is close to an even multiple of L, a second transition occurs where chain ends are enriched near the equatorial plane (and two further planes parallel to it, if R 2L); i.e., the system develops a distorted smectic structure. The order in the surface shell then is characterized by topological defects. Fil: Milchev, Andrey. Bulgarian Academy Of Sciences; Bulgaria Fil: Egorov, Sergei A.. University of Virginia; 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: Binder, Kurt. Johannes Gutenberg Universitat Mainz; Alemania Fil: Nikoubashman, Arash. Johannes Gutenberg Universitat Mainz; Alemania

We study active Brownian particles as a paradigm for a genuine nonequilibrium phase transition requiring steady driving. Access to the critical point in computer simulations is obstructed by the fact that the density is conserved. We propose a method based on arguments from finite-size scaling to determine critical points and successfully test it for the two-dimensional (2D) Ising model. Using this method allows us to accurately determine the critical point of two-dimensional active Brownian particles at ${\mathrm{Pe}}_{\text{cr}}=40(2), {\ensuremath{\phi}}_{\text{cr}}=0.597(3)$. Based on this estimate, we study the corresponding critical exponents $\ensuremath{\beta}, \ensuremath{\gamma}/\ensuremath{\nu}$, and $\ensuremath{\nu}$. Our results are incompatible with the 2D-Ising exponents, thus raising the question whether there exists a corresponding nonequilibrium universality class.

We performed extensive simulations of the random-bond Ising model confined between walls where competitive surface fields act. By properly taking the thermodynamic limit we unambiguously determined wetting transition points of the system, as extrapolation of localization-delocalization transitions of the interface between domains of different orientation driven by the respective fields. The finite-size scaling theory for wetting with short-range fields establishes that the average magnetization of the sample, with critical exponent β, is the proper order parameter for the study of wetting. While the hyperscaling relationship given by γ+2β=ν +ν requires β=1/2 (γ=4, ν =3, and ν =2), the thermodynamic scaling establishes that Δs=γ+β, which in contrast requires β=0 (Δs=4), where γ, ν, ν, and Δs are the critical exponents of the susceptibility, the correlation lengths parallel and perpendicular to the interface, and the gap exponent, respectively. So, we formulate a finite-size scaling theory for wetting without hyperscaling and perform numerical simulations that provide strong evidence of hyperscaling violation (i.e., β=0) and a direct measurement of the susceptibility critical exponent γ/ν =2.0±0.2, in agreement with theoretical results for the strong fluctuation regime of wetting transitions with quenched noise., Instituto de Física de Líquidos y Sistemas Biológicos

The two-dimensional Blume–Capel model with free surfaces where a surface field H1 acts and the “crystal field” (controlling the density of the vacancies) takes a value Ds different from the value D in the bulk, is studied by Monte Carlo methods. Using a recently developed finite size scaling method that studies thin films in a L × M geometry with antisymmetric surface fields (HL = −H1) and keeps a generalized aspect ratio c = L2/M constant, surface phase diagrams are computed for several typical choices of the parameters. It is shown that both second order and first order wetting transitions occur, separated by tricritical wetting behavior. The special role of vacancies near the surface is investigated in detail., Instituto de Física de Líquidos y Sistemas Biológicos

We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length $L$ is comparable to the persistence length $\ell_p$ and the case $L\gg \ell_p$. Our theory captures the early time monomer dynamics of a stiff chain characterized by $t^{3/4}$ dependence for the mean square displacement(MSD) of the monomers, but predicts a first crossover to the Rouse regime of $t^{2\nu/{1+2\nu}}$ for $\tau_1 \sim \ell_p^3$, and a second crossover to the purely diffusive dynamics for the entire chain at $\tau_2 \sim L^{5/2}$. We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16 - 2048 and persistence length $\ell_p = 1 - 500$ Lennard-Jones (LJ) units. These BD simulation results further confirm the absence of Gaussian regime for a 2d swollen chain from the slope of the plot of $\langle R_N^2 \rangle/2L \ell_p \sim L/\ell_p$ which around $L/\ell_p \sim 1$ changes suddenly from $\left(L/\ell_p \right) \rightarrow \left(L/\ell_p \right)^{0.5} $, also manifested in the power law decay for the bond autocorrelation function disproving the validity of the WLC in 2d. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness $\sqrt{\langle l_{\bot}^2\rangle}/L$ as a function of renormalized contour length $L/\ell_p$ collapse on the same master plot and exhibits power law scaling $\sqrt{\langle l_{\bot}^2\rangle}/L \sim (L/\ell_p)^\eta $ at extreme limits, where $\eta = 0.5$ for extremely stiff chains ($L/\ell_p \gg 1$), and $\eta = -0.25$ for fully flexible chains., Comment: 14 pages, 18 figures

We present a unified scaling theory for the dynamics of monomers for dilute solutions of semiflexible polymers under good solvent conditions in the free draining limit. Our theory encompasses the well-known regimes of mean square displacements (MSDs) of stiff chains growing like t^{3/4} with time due to bending motions, and the Rouse-like regime t^{2 \nu / (1+ 2\nu)} where \nu is the Flory exponent describing the radius R of a swollen flexible coil. We identify how the prefactors of these laws scale with the persistence length l_p, and show that a crossover from stiff to flexible behavior occurs at a MSD of order l^2_p (at a time proportional to l^3_p). A second crossover (to diffusive motion) occurs when the MSD is of order R^2. Large scale Molecular Dynamics simulations of a bead-spring model with a bond bending potential (allowing to vary l_p from 1 to 200 Lennard-Jones units) provide compelling evidence for the theory, in D=2 dimensions where \nu=3/4. Our results should be valuable for understanding the dynamics of DNA (and other semiflexible biopolymers) adsorbed on substrates., Comment: 4-page paper with 5 figures. 3-page supplemental information with 3 figures

A brief review focusing on low-dimensional colloidal model systems is given describing both simulation studies and complementary experiments, elucidating the interplay between phase behavior, geometric structures, and transport phenomena. These studies address the response of these very soft colloidal systems to perturbations such as uniform or uniaxial compression, laser fields, randomly quenched disorder, and shear deformation caused by moving boundaries.Binary hard-disk mixtures are studied by Monte Carlo simulation, to investigate ordering on surfaces or in monolayers, modeling the effect of a substrate by an external potential. By weak external laser fields the miscibility of the mixture can be controlled, and the underlying mechanism (laser-induced demixing) is clarified. The stability of various space-filling structures is discussed only for the case where no laser fields are present.Hard spheres interacting with repulsive screened Coulomb or dipolar interaction confined in 2D and 3D narrow constrictions are investigated by Brownian Dynamics simulation. With respect to the structural behavior, it is found that layers or planes throughout the microchannel are formed. The arrangement of the particles is disturbed by diffusion, and can also be modified by an external driving force causing a density gradient along the channel. Then the number of layers or planes gets reduced, adjusting to the density gradient, and this self-organized change of order also shows up in the particle velocities.The experimental work that is reviewed here addresses dipolar colloidal particles confined by gravity on a solid substrate on which a set of pinning sites has been randomly distributed. The dynamics of the system is studied by tracking the trajectories of individual particles, and it is found that the mean square displacements of particles that are nearest neighbors of pinned particles are strongly affected by these defects. The influence of the pinning sites on the order and microscopic mechanism of phase transitions in two dimensions is investigated.

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.

The Asakura-Oosawa model for colloid-polymer mixtures is studied by Monte Carlo simulations at densities inside the two-phase coexistence region of fluid and solid. Choosing a geometry where the system is confined between two flat walls, and a wall-colloid potential that leads to incomplete wetting of the crystal at the wall, conditions can be created where a single nanoscopic wall-attached crystalline cluster coexists with fluid in the remainder of the simulation box. Following related ideas that have been useful to study heterogeneous nucleation of liquid droplets at the vapor-liquid coexistence, we estimate the contact angles from observations of the crystalline clusters in thermal equilibrium. We find fair agreement with a prediction based on Young's equation, using estimates of interface and wall tension from the study of flat surfaces. It is shown that the pressure versus density curve of the finite system exhibits a loop, but the pressure maximum signifies the "droplet evaporation-condensation" transition and thus has nothing in common with a van der Waals-like loop. Preparing systems where the packing fraction is deep inside the two-phase coexistence region, the system spontaneously forms a "slab state", with two wall-attached crystalline domains separated by (flat) interfaces from liquid in full equilibrium with the crystal in between; analysis of such states allows a precise estimation of the bulk equilibrium properties at phase coexistence.

A model for two-dimensional colloids confined laterally by ``structured boundaries'' (i.e., ones that impose a periodicity along the slit) is studied by Monte Carlo simulations. When the distance $D$ between the confining walls is reduced at constant particle number from an initial value ${D}_{0}$, for which a crystalline structure commensurate with the imposed periodicity fits, to smaller values, a succession of phase transitions to imperfectly ordered structures occur. These structures have a reduced number of rows parallel to the boundaries (from $n$ to $n\ensuremath{-}1$ to $n\ensuremath{-}2$, etc.) and are accompanied by an almost periodic strain pattern, due to ``soliton staircases'' along the boundaries. Since standard simulation studies of such transitions are hampered by huge hysteresis effects, we apply the phase switch Monte Carlo method to estimate the free energy difference between the structures as a function of the misfit between $D$ and ${D}_{0}$, thereby locating where the transitions occur in equilibrium. For comparison, we also obtain this free energy difference from a thermodynamic integration method: The results agree, but the effort required to obtain the same accuracy as provided by phase switch Monte Carlo would be at least three orders of magnitude larger. We also show for a situation where several ``candidate structures'' exist for a phase, that phase switch Monte Carlo can clearly distinguish the metastable structures from the stable one. Finally, applying the method in the conjugate statistical ensemble (where the normal pressure conjugate to $D$ is taken as an independent control variable), we show that the standard equivalence between the conjugate ensembles of statistical mechanics is violated.

By large-scale Monte Carlo simulations of semiflexible polymers in $d=2$ dimensions the applicability of the Kratky-Porod model is tested. This model is widely used as "standard model" for describing conformations and force versus extension curves of stiff polymers. It is shown that semiflexible polymers in $d=2$ show a crossover from hard rods to self-avoiding walks, the intermediate Gaussian regime (implied by the Kratky-Porod model) is completely absent. Hence the latter can also describe force versus extension curves only if the contour length is only a few times larger than the persistence length. Consequences for experiments on biopolymers at interfaces are briefly discussed., 6 pages, 4 figures

A comparative simulation study of polymer brushes formed by grafting at a planar surface either flexible linear polymers (chain length $N_L$) or (non-catenated) ring polymers (chain length $N_R=2 N_L$) is presented. Two distinct off-lattice models are studied, one by Monte Carlo methods, the other by Molecular Dynamics, using a fast implementation on graphics processing units (GPUs). It is shown that the monomer density profiles $\rho(z)$ in the $z$-direction perpendicular to the surface for rings and linear chains are practically identical, $\rho_R(2 N_L, z)=\rho_L(N_L, z)$. The same applies to the pressure, exerted on a piston at hight z, as well. While the gyration radii components of rings and chains in $z$-direction coincide, too, and increase linearly with $N_L$, the transverse components differ, even with respect to their scaling properties: $R_{gxy}^{(L)} \propto N_L^{1/2}$, $R_{gxy}^{(R)} \propto N_L^{0.4}$. These properties are interpreted in terms of the statistical properties known for ring polymers in dense melts., Comment: 13 pages, 4 figures

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.

Phase transitions of systems confined in long cylindrical pores (capillary condensation, wetting, crystallization, etc.) are intrinsically not sharply defined but rounded. The finite size of the cross section causes destruction of long range order along the pore axis by spontaneous nucleation of domain walls. This rounding is analyzed for two models (Ising/lattice gas and Asakura-Oosawa model for colloid-polymer mixtures) by Monte Carlo simulations and interpreted by a phenomenological theory. We show that characteristic differences between the behavior of pores of finite length and infinitely long pores occur. In pores of finite length a rounded transition occurs first, from phase coexistence between two states towards a multi-domain configuration. A second transition to the axially homogeneous phase follows near pore criticality., PRL, accepted

9 pages; International audience; Molecular dynamics simulations of a dense melt of short bead-spring polymer chains containing N = 5, 10, or 25 effective monomers are presented and analyzed. Parts of our simulations include also a single dumbbell (N = 2) of the same type, which is interpreted to represent a coarse-grained model for a fluorescent probe molecule as used in corresponding experiments. We obtain the mean-square displacements of monomers and chains center of mass, and intermediate incoherent scattering functions of both monomers in the chains and particles in the dumbbells as function of time for a broad regime of temperatures above the critical temperature Tc of mode-coupling theory. For both the chains and the dumbbell, also orientational autocorrelation functions are calculated and for the dumbbell time series for the time evolution of linear dichroism and its autocorrelation function are studied. From both sets of data we find that both the mode-coupling critical temperature Tc (representing the “cage effect”) and the Vogel–Fulcher temperature T0 (representing the caloric glass transition temperature) systematically increase with chain length. Furthermore, the dumbbell dynamics yields detailed information on the differences in the matrix dynamics that are caused by the chain length variation. Deviations from the Stokes–Einstein relation are discussed, and an outlook to related experiments is given.

8 pages; International audience; We present molecular dynamics simulations of coarse-grained model systems of a glassforming polymer matrix containing fluorescent probe molecules. These probe molecules are either dispersed in the matrix or covalently attached to the center or the end of a dilute fraction of the polymer chains. We show that in all cases the translational and rotational relaxation of the probe molecules is a faithful sensor for the glass transition of the matrix as determined from a mode-coupling analysis or Vogel-Fulcher analysis of their R-relaxation behavior. Matrix and dumbbell related relaxation processes show a clear violation of the Stokes-Einstein-Debye laws. In accordance with recent experimental results, the long time behavior of single molecule spectroscopy observables like the linear dichroism is not susceptible to distinguish between center-attached and end-attached fluorophores. However we show that it is different from the behavior of dispersed fluorophores. We also show that the difference between the two attachment forms does show up in the caging regime of the relaxation functions and that this difference increases upon supercooling the melt toward its glass transition.

Monte Carlo Simulation in Statistical Physics deals with the computer simulation of many-body systems in condensed-matter physics and related fields of physics, chemistry and beyond, to traffic flows, stock market fluctuations, etc.). Using random numbers generated by a computer, probability distributions are calculated, allowing the estimation of the thermodynamic properties of various systems. This book describes the theoretical background to several variants of these Monte Carlo methods and gives a systematic presentation from which newcomers can learn to perform such simulations and to analyze their results. The fifth edition covers Classical as well as Quantum Monte Carlo methods. Furthermore a new chapter on the sampling of free-energy landscapes has been added. To help students in their work a special web server has been installed to host programs and discussion groups (http://wwwcp.tphys.uni-heidelberg.de). Prof. Binder was awarded the Berni J. Alder CECAM Award for Computational Physics 2001 as well as the Boltzmann Medal in 2007.

We report Monte Carlo studies of a two-dimensional soft colloidal crystal confined in a strip geometry by parallel walls. The wall-particle interaction has corrugations along the length of the strip. Compressing the crystal by decreasing the distance between the walls induces a structural transition characterized by the sudden appearance of a one-dimensional array of extended defects each of which span several lattice parameters, a "soliton staircase". We obtain the effective interaction between these defects. A Lindemann criterion shows that the reduction of dimensionality causes a finite periodic chain of these defects to readily melt as the temperature is raised. We discuss possible experimental realizations and speculate on potential applications., 4 pages 5 embedded figures

A simple fluid, described by point-like particles interacting via the Lennard-Jones potential, is considered under confinement in a slit geometry between two walls at distance Lz apart for densities inside the vapor-liquid coexistence curve. Equilibrium then requires the coexistence of a liquid "bridge" between the two walls, and vapor in the remaining pore volume. We study this equilibrium for several choices of the wall-fluid interaction (corresponding to the full range from complete wetting to complete drying, for a macroscopically thick film), and consider also the kinetics of state changes in such a system. In particular, we study how this equilibrium is established by diffusion processes, when a liquid is inserted into an initially empty capillary (partial or complete evaporation into vacuum), or when the volume available for the vapor phase increases. We compare the diffusion constants describing the rates of these processes in such inhomogeneous systems to the diffusion constants in the corresponding bulk liquid and vapor phases., accepted for publication in Physical Review E

Confining a colloidal crystal within a long narrow channel produced by two parallel walls can be used to impose a meso-scale superstructure of a predominantly mechanical elastic character [Chui et al., EPL 2008, 83, 58004]. When the crystal is compressed in the direction perpendicular to the walls, we obtain a structural transition when the number of rows of particles parallel to the walls decreases by one. All the particles of this vanishing row are distributed throughout the crystal. If the confining walls are structured (say with a corrugation along the length of the walls), then these extra particles are distributed neither uniformly nor randomly; rather, defect structures are created along the boundaries resembling "soliton staircases", inducing a non-uniform strain pattern within the crystal. Here we study the conditions of stability, formation and annihilation of these solitons using a coarse grained description of the dynamics. The processes are shown by comparing superimposed configurations as well as molecular animations obtained from our simulations. Also the corresponding normal and shear stresses during the transformation are calculated. A study of these dynamical processes should be useful for controlling strain wave superstructures in the self-assembly of various nano- and meso scaled particles., 26 pages, 6 figures

As a generic model system of an asymmetric binary fluid mixture, hexadecane dissolved in carbon dioxide is considered, using a coarse-grained bead-spring model for the short polymer, and a simple spherical particle with Lennard-Jones interactions for the carbon dioxide molecules. In previous work, it has been shown that this model reproduces the real phase diagram reasonable well, and also the initial stages of spinodal decomposition in the bulk following a sudden expansion of the system could be studied. Using the parallelized simulation package ESPResSo on a multiprocessor supercomputer, phase separation of thin fluid films confined between parallel walls that are repulsive for both types of molecules are simulated in a rather large system (1356 x 1356 x 67.8 A^3, corresponding to about 3.2 million atoms). Following the sudden system expansion, a complicated interplay between phase separation in the directions perpendicular and parallel to the walls is found: in the early stages the hexadecane molecules accumulate mostly in the center of the slit pore, but as the coarsening of the structure in the parallel direction proceeds, the inhomogeneity in the perpendicular direction gets much reduced. Studying then the structure factors and correlation functions at fixed distances from the wall, the densities are essentially not conserved at these distances, and hence the behavior differs strongly from spinodal decomposition in the bulk. Some of the characteristic lengths show a nonmonotonic variation with time, and simple coarsening described by power-law growth is only observed if the domain sizes are much larger than the film thickness., accepted for publication in PRE

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 behaviour 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 present also 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., Phys. Rev. E, accepted

Ordering phenomena in colloidal dispersions exposed to external one-dimensional, periodic fields or under confinement are studied systematically by Monte Carlo computer simulations. Such systems are useful models for the study of monolayers on a substrate. We find that the interaction with a substrate potential completely changes the miscibility of a binary, hard disc mixture at low external field amplitudes. The underlying ordering mechanisms leading to this laser-induced de-mixing differ, depending on which components interact with the substrate potential. Generic effects of confinement on crystalline order in two dimensions are studied in a model system of point particles interacting via a potential ∝ r −12 . The state of the system (a strip of width D) depends very sensitively on the precise boundary conditions at the two confining walls. Commensurate, corrugated boundary conditions enhance both orientational order and positional order. In contrast, smooth repulsive boundaries enhance only the orientational order and destroy positional (quasi-)long range order. As external fields have a strong impact on the elastic behaviour of colloidal crystals there is a need to analyse the elastic response in such systems for the field-free case first. To this aim we study the strain–strain correlation functions in a two-dimensional crystal formed by super-paramagnetic colloids, as monitored by standard video microscopy. (Some figures in this article are in colour only in the electronic version)

Scaling predictions and results from self-consistent field calculations for bottle-brush polymers with a rigid backbone and flexible side chains under good solvent conditions are summarized and their validity and applicability is assessed by a comparison with Monte Carlo simulations of a simple lattice model. It is shown that under typical conditions, as they are also present in experiments, only a rather weak stretching of the side chains is realized, and then the scaling predictions based on the extension of the Daoud-Cotton blob picture are not applicable. Also two-component bottle brush polymers are considered, where two types (A,B) of side chains are grafted, assuming that monomers of different kind repel each other. In this case, variable solvent quality is allowed for, such that for poor solvent conditions rather dense cylinder-like structures result. Theories predict ``Janus Cylinder''-type phase separation along the backbone in this case. The Monte Carlo simulations, using the pruned-enriched Rosenbluth method (PERM) then are restricted to rather short side chain length. Nevertheless, evidence is obtained that the phase separation between an A-rich part of the cylindrical molecule and a B-rich part can only occur locally. The correlation length of this microphase separation can be controlled by the solvent quality. This lack of a phase transition is interpreted by an analogy with models for ferromagnets in one space dimension., 27pages, 33 figures

In this work we compare and characterize the behavior of Langevin and Dissipative Particle Dynamics (DPD) thermostats in a broad range of non-equilibrium simulations of polymeric systems. Polymer brushes in relative sliding motion, polymeric liquids in Poiseuille and Couette flows, and brush-melt interfaces are used as model systems to analyze the efficiency and limitations of different Langevin and DPD thermostat implementations. Widely used coarse-grained bead-spring models under good and poor solvent conditions are employed to assess the effects of the thermostats. We considered equilibrium, transient, and steady state examples for testing the ability of the thermostats to maintain constant temperature and to reproduce the underlying physical phenomena in non-equilibrium situations. The common practice of switching-off the Langevin thermostat in the flow direction is also critically revisited. The efficiency of different weight functions for the DPD thermostat is quantitatively analyzed as a function of the solvent quality and the non-equilibrium situation., 12 pages, introduction improved, references added, to appear in Phys. Rev. E

A lattice model for a symmetrical copolymer "bottle brush" molecule, where two types (A,B) of flexible side chains are grafted with one chain end to a rigid backbone, is studied by a variant of the pruned-enriched Rosenbluth method (PERM), allowing for simultaneous growth of all side chains in the Monte Carlo sampling. Choosing repulsive binary interactions between unlike monomers and varying the solvent quality, it is found that phase separation into an $A$-rich part of the cylindrical molecule and a $B$-rich part can occur only locally. Long range order (in the direction of the backbone) does not occur, and hence the transition from the randomly mixed state of the bottle brush to the phase separated structure is strongly rounded, in contrast to corresponding mean field predictions. This lack of a phase transition can be understood from an analogy with spin models in one space dimension. We predict that the range of microphase separation along the bottle brush backbone can be controlled on the nanoscale by varying the solvent quality., 5 pages, 7 figures. Edited for typos

A combination of classical molecular dynamics (MD) and ab initio Car-Parrinello molecular dynamics (CPMD) simulations is used to investigate the adsorption of water on a free amorphous silica surface. From the classical MD SiO_2 configurations with a free surface are generated which are then used as starting configurations for the CPMD.We study the reaction of a water molecule with a two-membered ring at the temperature T=300K. We show that the result of this reaction is the formation of two silanol groups on the surface. The activation energy of the reaction is estimated and it is shown that the reaction is exothermic., 12 pages, 6 figures, to be published in J. Phys.: Condens. Matter

Conformational properties of polymer melts confined between two hard structureless walls are investigated by Monte Carlo simulation of the bond-fluctuation model. Parallel and perpendicular components of chain extension, bond-bond correlation function and structure factor are computed and compared with recent theoretical approaches attempting to go beyond Flory's and Silberberg's hypotheses. We demonstrate that for ultrathin films where the thickness, $H$, is smaller than the excluded volume screening length (blob size), $\xi$, the chain size parallel to the walls diverges logarithmically, $R^2/2N \approx b^2 + c \log(N)$ with $c \sim 1/H$. The corresponding bond-bond correlation function decreases like a power law, $C(s) = d/s^{\omega}$ with $s$ being the curvilinear distance between bonds and $\omega=1$. % Upon increasing the film thickness, $H$, we find -- in contrast to Flory's hypothesis -- the bulk exponent $\omega=3/2$ and, more importantly, an {\em decreasing} $d(H)$ that gives direct evidence for an {\em enhanced} self-interaction of chain segments reflected at the walls. Systematic deviations from the Kratky plateau as a function of $H$ are found for the single chain form factor parallel to the walls in agreement with the {\em non-monotonous} behaviour predicted by theory. This structure in the Kratky plateau might give rise to an erroneous estimation of the chain extension from scattering experiments. For large $H$ the deviations are linear with the wave vector, $q$, but are very weak. In contrast, for ultrathin films, $H, Comment: 16 pages, 7 figures. Dedicated to L. Sch\"afer on the occasion of his 60th birthday

Using a coarse-grained bead-spring model of flexible polymer chains, the structure of a polymeric nanodroplet adsorbed on a chemically decorated flat wall is investigated by means of Molecular Dynamics simulation. We consider sessile drops on a lyophilic (attractive for the monomers) region of circular shape with radius R_D while the remaining part of the substrate is lyophobic. The variation of the droplet shape, including its contact angle, with R_D is studied, and the density profiles across these droplets also are obtained. In addition, the interaction of droplets adsorbed on two walls forming a slit pore with two lyophilic circular regions just opposite of one another is investigated, paying attention to the formation of a liquid bridge between both walls. A central result of our study is the measurement of the force between the two substrate walls at varying wall separation as well as the kinetics of droplet merging. Our results are compared to various phenomenological theories developed for liquid droplets of mesoscopic rather than nanoscopic size., 8 pages, 9 figures, accepted in J. Chem. Phys. 2004

We present the results of molecular dynamics computer simulations of a binary Lennard-Jones liquid confined between two parallel rough walls. These walls are realized by frozen amorphous configurations of the same liquid and therefore the structural properties of the confined fluid are identical to the ones of the bulk system. Hence this setup allows us to study how the relaxation dynamics is affected by the pure effect of confinement, i.e. if structural changes are completely avoided. We find that the local relaxation dynamics is a strong function of z, the distance of the particles from the wall, and that close to the surface the typical relaxation times are orders of magnitude larger than the ones in the bulk. Because of the cooperative nature of the particle dynamics, the slow dynamics also affects the dynamics of the particles for large values of z. Using various empirical laws, we are able to parameterize accurately the z-dependence of the generalized incoherent intermediate scattering function F_s(q,z,t) and also the spatial dependence of structural relaxation times. These laws allow us to determine various dynamical length scales and we find that their temperature dependence is compatible with an Arrhenius law. Furthermore, we find that at low temperatures time and space dependent correlation function fulfill a generalized factorization property similar to the one predicted by mode-coupling theory for bulk systems. For thin films and/or at sufficiently low temperatures, we find that the relaxation dynamics is influenced by the two walls in a strongly non-linear way in that the slowing down is much stronger than the one expected from the presence of only one confining wall. ...., Comment: 22 pages of Latex

A simple modification of the ``Wang-Landau sampling'' algorithm removes the systematic error that occurs at the boundary of the range of energy over which the random walk takes place in the original algorithm., 3 pages, Latex, 2 Figures

Polymers are good glass formers and allow for the study of melts near the glass transition in (meta-)stable equilibrium. Theories of the glass transition imply such an equilibrium and can, hence, be tested by the study of polymer melts. After a brief summary of the basic experimental facts about the glass transition in polymers, the main theoretical concepts are reviewed: mode coupling theory (MCT), entropy theory, free-volume theory, the idea of a growing length describing the size of cooperative regions, etc. Then, two basic coarse-grained models of polymers are described, which have been developed aiming at a test of these concepts. The first model is the bond-fluctuation model on the simple cubic lattice; the second is a bead-spring model in the continuum. While the first model is studied by kinetic Monte Carlo methods, the second is studied by molecular dynamics simulation: the issue is addressed which aspects of the results are model-dependent and which aspects are generic and should universally apply, including real materials. Attempts to include chemical detail in order to describe real materials are discussed, too. It is shown that idealized MCT is a good description for the onset of slow relaxation (‘cage effect’) in the discussed model polymer melts, but the singularities at the critical temperature of MCT are rounded, and the correct description of the dynamics close to the calorimetric glass transition remains controversial. While the configurational entropy of the polymer melt does decrease strongly when the glass transition is approached, evidence is presented that the ‘entropy catastrophe’ of the theory of Gibbs and Di Marzio is an artifact of inaccurate approximations. Simulations of polymer melts confined in a thin film geometry are also reviewed, emphasizing the question to which extent they shed light on the issue of a growing glass correlation length. Finally, we discuss the implications of the glass transition on the motion of polymers at larger length scales. It is shown that for short (non-entangled) chains the Rouse model stays essentially valid, the friction coefficient reflecting the slowing down as described by the α-relaxation, while for small-scale motions (e.g. β-relaxation regime of MCT) the connectivity of the polymer is not very relevant. We conclude that the theories describe some aspects of vitrification phenomena correctly, but a unified description that puts all these phenomena into one coherent framework, does not yet exist.

Results of large-scale molecular-dynamics simulations of a supercooled polymer film are presented (F. Varnik, J. Baschnagel, K. Binder, J. Phys. IV 10, 239 (2000)). The dynamic and static properties of the system are studied for a wide range of film thicknesses (from 3 to about 55 times the bulk radius of gyration) and temperatures (from the normal liquid state to the supercooled region). The system is confined between two completely smooth and purely repulsive walls. Motivated by the previous results on the enhancement of the local relaxation dynamics due to the confinement (F. Varnik, J. Baschnagel, K. Binder, Eur. Phys. J. E 8, 175 (2002); Phys. Rev. E. 65, 021507 (2002)), we now study the effect of the walls on the dynamics of the Rouse modes. It has been reported from Monte Carlo studies of the Bond Fluctuation Model (BFM) that, contrary to the enhancement of the "cage dynamics" (exemplified by a faster relaxation of the incoherent scattering function in the film), Rouse modes exhibit a slower relaxation in the confined system (C. Mischler, J. Baschnagel, K. Binder, Adv. Colloid Interface Sci. 94, 197 (2001)). However, we do not observe such a discrepancy for our continuum model: At a given temperature, the relaxation of a given Rouse mode is faster in the film than in the bulk in accordance with the acceleration of the dynamics around the cage.