The nature of the calculation of the pressure in molecular simulations of continuous models from volume perturbations

Published versio

New apparatus for DTA at 2000 bar: thermodynamic studies on Au, Ag, Al and HTSC oxides

A new DTA (Differential Thermal Analysis) device was designed and installed in a Hot Isostatic Pressure (HIP) furnace in order to perform high-pressure thermodynamic investigations up to 2 kbar and 1200C. Thermal analysis can be carried out in inert or oxidising atmosphere up to p(O2) = 400 bar. The calibration of the DTA apparatus under pressure was successfully performed using the melting temperature (Tm) of pure metals (Au, Ag and Al) as standard calibration references. The thermal properties of these metals have been studied under pressure. The values of DV (volume variation between liquid and solid at Tm), ROsm (density of the solid at Tm) and ALPHAm (linear thermal expansion coefficient at Tm) have been extracted. A very good agreement was found with the existing literature and new data were added. This HP-DTA apparatus is very useful for studying the thermodynamics of those systems where one or more volatile elements are present, such as high TC superconducting oxides. DTA measurements have been performed on Bi,Pb(2223) tapes up to 2 kbar under reduced oxygen partial pressure (p(O2) = 0.07 bar). The reaction leading to the formation of the 2223 phase was found to occur at higher temperatures when applying pressure: the reaction DTA peak shifted by 49C at 2 kbar compared to the reaction at 1 bar. This temperature shift is due to the higher stability of the Pb-rich precursor phases under pressure, as the high isostatic pressure prevents Pb from evaporating.Comment: 6 figures, 3 tables, Thermodynamics, Thermal property, Bi-2223, fundamental valu

The ELBA Force Field for Coarse-Grain Modeling of Lipid Membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities

Computer simulation study of gas–liquid nucleation in a Lennard-Jones system

We report a computer-simulation study of homogeneous gas–liquid nucleation in a Lennard-Jones system. Using umbrella sampling, we compute the free energy of a cluster as a function of its size. A thermodynamic integration scheme is employed to determine the height of the nucleation barrier as a function of supersaturation. Our simulations illustrate that the mechanical and the thermodynamical surfaces of tension and surface tension differ significantly. In particular, we show that the mechanical definition of the surface tension cannot be used to compute this barrier height. We find that the relations recently proposed by McGraw and Laaksonen [J. Chem. Phys. 106, 5284 (1997)] for the height of the barrier and for the size of the critical nucleus are obeyed

Evaluation of pressure tensor in constant-volume simulations of hard and soft convex bodies

A method for calculating the pressure tensor in constant-volume Monte Carlo simulations of convex bodies is presented. In contrast to other approaches, the method requires only an isotropic scaling of the simulation box and the counting of simple geometric quantities characterizing overlapping pairs. Nonsphericity presents no special difficulties. The result is expressed as a sum of pairwise contributions and can therefore be used to compute pressure tensor profiles in a conventional way. (c) 2006 American Institute of Physics

Anisotropic surface stresses of a solid/fluid interface: Molecular dynamics calculations for the copper/methane interface

International audienc