Properties of a soft-sphere liquid from non-Newtonian molecular dynamics

A soft-sphere, inverse-12 liquid is simulated in both the isokinetic-isochoric and the isokinetic-isobaric ensemble using nonequilibrium molecular dynamics. The simulation for the isobaric ensemble is discussed in detail. The non-Newtonian characteristics of the liquid are clearly demonstrated; namely, the shear-rate-dependent pressure and density (shear dilatancy), the shear-rate-dependent shear viscosity (shear thinning), and evidence of normal pressure differences. For the first time, it is clearly shown that a significant component of isobaric shear thinning is due to shear dilatancy. The isochoric and isobaric results are checked for consistency. Simple empirical relations for the equation of state and transport properties of the fluid are presented.This is a publication in part of the National Institute of Standards and Technology (formerly NBS) and is not subject to copyright.     

A simple effective pair potential for the molecular simulation of the thermodynamic properties of ammonia

A new optimized effective pair potential model is proposed, which is appropriate for the prediction of thermodynamic properties of fluid ammonia including vapour—liquid coexistence data. The phase behaviour is determined using a recently developed version of the Gibbs ensemble Monte Carlo method. Furthermore, liquid structure characteristics, the dielectric constant and supercritical properties are determined by Monte Carlo simulations in the isothermal—isobaric ensemble. The second virial coefficient of the pair potential model is calculated over a broad range of temperature. All properties are compared with experimental data or results of a multi-parameter equation of state for ammonia. The new model is found to yield coexistence properties and second virial coefficients in good agreement with experimental data and the results of the equation of state, respectively.     

氢在Nd晶体中行为的分子动力学模拟

由三维Mobius反演变换所得的金属Nd原子和H原子间的相互作用势和组合规则的方法得到的Nd H原子间的相互作用势 ,利用正则系统分子动力学算法研究了在一定加载应力强度因子K =0 .6MPam下 ,氢在Nd晶体中的行为。模拟结果表明 ,氢在Nd晶体裂尖富集成许多氢原子团或氢气团。这可用来部分地解释NdFeB稀土永磁体吸氢后的氢爆行为。     

Dynamics of hydrogen bonds in water and consequences for the unusual behaviour of supercooled water

The dynamics of liquid water is evaluated by the coherent quasi-elastic scattering at two different momentum transfers, in order to discriminate hydrogen bond lifetime from molecular dynamics. The results indicate a possible issue for the puzzle of the behaviour of supercooled water.      

Chaotic behavior in nonlinear Hamiltonian systems and equilibrium statistical mechanics

The relation between chaotic dynamics of nonlinear Hamiltonian systems and equilibrium statistical mechanics in its canonical ensemble formulation has been investigated for two different nonlinear Hamiltonian systems. We have compared time averages obtained by means of numerical simulations of molecular dynamics type with analytically computed ensemble averages. The numerical simulation of the dynamic counterpart of the canonical ensemble is obtained by considering the behavior of a small part of a given system, described by a microcanonical ensemble, in order to have fluctuations of the energy of the subsystem. The results for the Fermi-Pasta-Ulam model (i.e., a one-dimensional anharmonic solid) show a substantial agreement between time and ensemble averages independent of the degree of stochasticity of the dynamics. On the other hand, a very different behavior is observed for a chain of weakly coupled rotators, where linear exchange effects are absent. In the high-temperature limit (weak coupling) we have a strong disagreement between time and ensemble averages for the specific heat even if the dynamics is chaotic. This behavior is related to the presence of spatially localized chaos, which prevents the complete filling of the accessible phase space of the system. Localized chaos is detected by the distribution of all the characteristic Liapunov exponents.     

Structure of a liquid-vapor interface

The structure of the interface of an argonlike fluid in equilibrium with its own vapor at low temperature is studied using molecular dynamics. The longitudinal pair correlations in the interface are found to be consistent with a simply defined ensemble of local thermodynamic states. However, the transverse correlations exhibit very long-range behavior not predicted by straightforward local thermodynamics. These results strongly suggest that the interface is made up of an ensemble of configurations in each of which the transition from liquid to vapor is locally sharp, but that the transition surface fluctuates strongly in space and time.Supported in part by ERDA, Contract No. EY-76-C-02-3077, and by NSF, Grant DMR-72-03029, through Materials Science Center, Cornell University.The classical theoretical density profile dates back to van der Waals. (1)     

Zeeman effect in the N.Q.R. of 81Br in single crystals of 2,4-dibromoanisole and 2,4,6-tribromoanisole

The influence of the intramolecular degrees of freedom on the vapour–liquid equilibrium properties of ammonia is studied for vapour pressure, saturated densities and enthalpy of vaporization. Molecular force fields with and without intramolecular degrees of freedom, keeping all other parameters unchanged, show significantly different phase envelopes. For ammonia, the angle potential is particularly important, because the hydrogen sites are more aligned in the liquid than in the vapour, leading to a significantly enhanced molecular dipole moment in the condensed phase. Based on a rigid force field for ammonia from prior work of our group [Eckl et al., Mol. Phys. 106, 1039 (2008)], a new accurate force field with intramolecular degrees of freedom is developed.     

Hydrogen bonding in liquid water: An ab initio QM/MM MD simulation study

Pattern and dynamics of hydrogen bonds in liquid water were investigated by a quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulation at Hartree–Fock (HF) level of theory. A large subregion of the whole system comprising two complete coordination shells was treated quantum mechanically in order to include all polarization and charge transfer effects and to obtain accurate data about structure and dynamics of the intermolecular bonds. The results of this investigation are in agreement with recent experimental findings and suggest that in liquid water every molecule forms in average 2.8, but almost as a rule less than four intermolecular hydrogen bonds.     

Thermodynamic and transport properties of nitrogen and butane mixtures

A force field has been developed to describe the phase behaviour, interfacial, and transport properties of nitrogen and hydrocarbon mixtures under conditions relevant to those found in the high pressure extraction of oil from underground reservoirs. A Gibbs ensemble Monte Carlo method is used to parametrize intermolecular potentials for the pure components by matching experimental and simulated liquid and vapour coexisting densities. Also the surface tension, diffusion coefficient and shear viscosity of nitrogen and its mixtures with butane have been determined. The latter properties were obtained by canonical molecular dynamics simulations. The diffusion coefficient and shear viscosity were calculated by a Green-Kubo method. Results for pure nitrogen are given for temperatures ranging from 70 K to 110K. For mixtures of nitrogen with butane, results are presented at 339.4 K and 380.2 K. Good agreement is found between the results of simulations and available experimental data.     

On the material frame indifference controversy: Some results from group theory and computer simulation.

Group theory and computer simulation are used to show that molecular and atomic scale time cross correlation functions have different symmetry properties in different frames of reference. Constitutive equations, written in any frame of reference, are approximations to classical molecular dynamics, and must therefore reproduce the results from group theory applied to ensemble averages such as time correlation functions, and corroborated by computer simulation. This may help to resolve some of the contemporary literature claims regarding material frame indifference and the implementation of “co- rotating” frames.     

固壁加热的分子动力学模拟研究

本文应用分子动力学模拟方法,讨论了NVE系综中不同初始晶格类型对模拟系统平衡温度和截断半径对模拟系统平衡总能的影响.建立了一种新的固壁加热模型,模拟研究了由铂组成固壁对流体氩快速加热过程.在本文的模拟条件下,由于固壁与液体之间的快速加热作用,固壁与液体之间产生的蒸汽推动液体向上移动,可以明显地观察到Leidenfrost现象的产生.     

从常温常压到超临界乙醇的分子动力学模拟

采用分子动力学方法系统地研究了从常温常压到超临界状态乙醇的热力学性质、结构性质和动力学性质.模拟发现随着温度的升高,体系焓值增大,乙醇分子间的氢键作用减弱,自扩散系数增大;随着压强的增大,乙醇分子间的氢键作用增强,自扩散系数减小;乙醇自扩散系数在液相区随温度变化不明显,在气相区随压强增大很快减小,超临界区乙醇的自扩散系数比液相区大十几倍.温度和压强对乙醇自扩散系数的影响可通过密度来体现.与常温常压相比,超临界条件下的乙醇体系因密度涨落存在分子聚集现象,且在低密度区域更显著;乙醇分子间的氢键作用明显减弱,结     

R32热力学性质的蒙特卡罗模拟

建立在统计热力学和分子力学理论基础上的分子模拟方法逐渐运用于计算制冷剂的热力学性质。文中首先在NVT系综条件下,采用吉布斯蒙特卡罗模拟方法(GEMC),模拟了R32的气液相平衡的密度、饱和蒸汽压及蒸发焓;其次在NPT系综条件下,采用蒙特卡罗模拟方法(MC),模拟了R32在4MPa条件下的过冷液态密度。模拟结果同美国国家标准研究所(NIST)的Refprop 8.0相比,有很好的一致性。结果表明,运用该方法预测单一制冷剂的热力学性质是可行的。     

Canonical ensemble and nonequilibrium states by molecular dynamics

We present a new technique to simulate the contact of a molecular dynamics system with a thermal wall. A canonical ensemble is obtained, and its statistical and thermodynamic fluctuations are studied. The values of the specific heat found by simulation agree with the experimental data. By means of thermal walls at different temperatures, thermal gradients are obtained. The values of the thermal conductivity are consistent with the experimental data.This work was supported in part by the Commissariat à l'Énergie Atomique (France).     

Molecular dynamics simulations of PVP kinetic inhibitor in liquid water and hydrate/liquid water systems

The possible effects of PVP (poly(N-vinylpyrrolidone)) on the properties of liquid and water in clathrate hydrate has been investigated using NVT molecular dynamics simulations. A model for a monomer of the PVP polymer is immersed in three systems, liquid water, a unit cell of a hydrate in liquid water with a hydrate former and a third system where some of the liquid water molecules of this last system are replaced by a PVP monomer. Both molecular dynamics simulation and integral equation theory predict hydrogen bonding between the double bonded oxygen in the PVP ring and hydrogens in water. For the composite system, the PVP monomer has a preference for hydrogen bonding to hydrogens from the water molecules at the surface of the hydrate lattice. The simulations indicate that the PVP monomer tends to orient perpendicular to the hydrate surface. For the model systems in this study PVP may form hydrogen bonds with liquid water through the double bonded oxygen in the ring. When a hydrate crystal is immersed in the liquid water phase this hydrogen bonding is shifted towards the hydrate due to a more favourable Coulomb interaction involving hydrogens from more than one water molecule at the hydrate surface. The PVP monomer has a preference for perpendicular orientation with respect to the hydrate surface. A scheme is suggested for the characterization of kinetic hydrate inhibitors based on molecular dynamics simulations and on three basic properties. In addition to the energy between the active groups of the inhibitor and hydrate water another point of focus is the free energy changes in the interactions between the inhibitor and water as the charges are changed from zero to the original model charges. In particular the difference between this integral for the (hydrate water)–(PVP monomer) interaction and the (liquid water)–(PVP inhibitor) interaction should reflect the driving forces in freezing the inhibitor out from the liquid water phase and onto the hydrate surface. The third property in the characterization scheme is the diffusivities of groups connecting to the hydrate crystal, relative to the diffusivities of the hydrate crystal. Results are presented from simulations where a small cavity with a methane model as a guest is immersed in liquid water with free methane molecules at a temperature of 150K. Changes in structure, diffusivities and energy indicate a tendency towards a more solid–like structurde around the cavity.     

Dynamic properties of hydrogen-bonded networks in supercritical water

Dynamic properties of supercritical water at temperatures between 573 and 773 K and densities between 0.49 and 0.83 g/cm(3) have been investigated by molecular dynamics simulation and compared to states located on the vapor-liquid coexistence curve. A flexible simple point charge potential has been assumed for interactions in the subcritical states, whereas a reparameterization of that model has been performed to model the supercritical states. The hydrogen bonding structure and the diffusion coefficients in an ensemble of simulated states were previously analyzed and a good agreement with available experimental data was found. Dynamic properties of hydrogen bonding like the bond lifetimes and the influence of hydrogen bonds in the vibrational spectra are discussed along a range of simulation conditions. A nonlinear behavior of the hydrogen-bond lifetime as a function of temperature is observed in subcritical water whereas a linear dependence is found in supercritical water. Special attention is paid to the intermolecular vibrational spectrum (10-400 cm(-1)). It has been observed that the mode centered at 200 cm(-1), attributed to the intermolecular O-O stretching vibration in the ambient state remains active in the supercritical states.     

Atomistic simulations of liquid water using Lekner electrostatics

Equilibrium molecular dynamics simulations have been performed for liquid water using three different potential models in the NVT and NPT ensembles. The flexible SPC model, the rigid TIP4P model and the rigid/polarizable TIP4P-FQ potential were studied. The Lekner method was used to handle long range electrostatic interactions, and an efficient trivariate cubic spline interpolation method was devised for this purpose. A partitioning of the electrostatic interactions into medium and long range parts was performed, and the concomitant use of multiple timestep techniques led to substantially enhanced computation speeds. The simulations were carried out using 256 molecules in the NVT ensemble at 25°C and 997 kg m^?3 and in the NPT ensemble at 25°C and 1 bar. Various dynamic, structural, dielectric, rotational and thermodynamic properties were calculated, and it was found that the simulation methodologies performed satisfactorily vis-à-vis previous simulation results and experimental observations.     

Thermodynamically guided nonequilibrium Monte Carlo method for generating realistic shear flows in polymeric systems

A thermodynamically guided atomistic Monte Carlo methodology is presented for simulating systems beyond equilibrium by expanding the statistical ensemble to include a tensorial variable accounting for the overall structure of the system subjected to flow. For a given shear rate, the corresponding tensorial conjugate field is determined iteratively through independent nonequilibrium molecular dynamics simulations. Test simulations for the effect of flow on the conformation of a C50H102 polyethylene liquid show that the two methods (expanded Monte Carlo and nonequilibrium molecular dynamics) provide identical results.     

Rotational dynamics of a molecular ensemble in the presence of a strong laser field

The rotational dynamics of a molecular ensemble in the presence of a strong laser field is investigated in the framework of the density-matrix approach. The results obtained are compared with the rotation of molecules in a classical ensemble. The thermal molecular motion in the ensemble is taken into account in a model of random collisions, and various values of the relaxation time are considered. The effect of the relaxation process on molecular alignment is analyzed for the laser-pulse action and in the afterpulse regime. The similarity between quantum and classical rotational dynamics of a molecular ensemble is examined.     

Structure and dynamics of liquid formamide at high pressure and high temperature: comparison of X-ray diffraction and molecular dynamics results

New molecular dynamics simulations with optimised potentials for liquid simulation are presented for liquid formamide at high pressures and high temperatures. The structural results are compared to those found by X-ray diffraction measurements, and the H-bonding structure is analysed in detail. While it is ambiguous from purely experimental data, the simulation results support the idea that pressure increase can enhance ring dimer formation at the expense of linear chain structure, but increasing temperature results in an opposite effect. Dynamical results (reorientational correlation times and lifetimes of hydrogen bonds) are in agreement with these findings.