scholarly journals Effective Determination of Coexistence Curves using Reversible-Scaling Molecular Dynamics Simulations

2000 ◽  
Vol 653 ◽  
Author(s):  
Maurice de Koning ◽  
Alex Antonelli ◽  
Sidney Yip

AbstractWe present a simulation technique that allows the calculation of a phase coexistence curve from a single nonequilibrium molecular dynamics (MD) simulation. The approach is based on the simultaneous simulation of two coexisting phases, each in its own computational cell, and the integration of the relevant Clausius-Clapeyron equation starting from a known coexistence point. As an illustration of the effectiveness of our approach we apply the method to explore the melting curve in the Lennard-Jones phase diagram.

2016 ◽  
Vol 195 ◽  
pp. 557-568 ◽  
Author(s):  
Pablo M. Piaggi ◽  
Omar Valsson ◽  
Michele Parrinello

We study by computer simulation the nucleation of a supersaturated Lennard-Jones vapor into the liquid phase. The large free energy barriers to transition make the time scale of this process impossible to study by ordinary molecular dynamics simulations. Therefore we use a recently developed enhanced sampling method [Valsson and Parrinello, Phys. Rev. Lett.113, 090601 (2014)] based on the variational determination of a bias potential. We differ from previous applications of this method in that the bias is constructed on the basis of the physical model provided by the classical theory of nucleation. We examine the technical problems associated with this approach. Our results are very satisfactory and will pave the way for calculating the nucleation rates in many systems.


Energies ◽  
2020 ◽  
Vol 14 (1) ◽  
pp. 11
Author(s):  
Shi Yu ◽  
Ruizhi Chu ◽  
Xiao Li ◽  
Guoguang Wu ◽  
Xianliang Meng

Both molecular dynamics (MD) and nonequilibrium molecular dynamics (NEMD) simulations were performed to simulate coal ashes using the Guillot-Sator model in this work. The structural and transport properties of coal ashes at high temperatures have been obtained. Superheating of coal ash system with anorthite crystal structure initial configuration has been observed for MD simulation which explains the discrepancy between previous MD simulation results and FactSage thermochemical calculations. The fluxing effects of both calcium oxide and sodium oxide have been investigated systematically through MD and NEMD simulations. Moreover, the viscosities of coal ash systems have been computed by two methods: (1) Stokes-Einstein equation; (2) NEMD simulations. Estimations of viscosities for various coal ash systems based on Stokes-Einstein equation exhibit a strong temperature dependence of viscosity, which agrees with previous experimental results. On the other hand, NEMD simulation results that showed a strong shear-thinning feature, failed to reproduce this strong temperature dependence of viscosity, possibly due to the short simulation time. Nevertheless, NEMD simulations not only provide us detailed information about atoms dynamics under shear, but also allow us to model the coal ash system far from equilibrium which cannot be accessed by thermodynamics calculation using software like FactSage.


2011 ◽  
Vol 134 (1) ◽  
Author(s):  
John C. Duda ◽  
Timothy S. English ◽  
Donald A. Jordan ◽  
Pamela M. Norris ◽  
William A. Soffa

Many random substitutional solid solutions (alloys) will display a tendency to atomically order given the appropriate kinetic and thermodynamic conditions. Such order–disorder transitions will result in major crystallographic reconfigurations, where the atomic basis, symmetry, and periodicity of the alloy change dramatically. Consequently, phonon behavior in these alloys will vary greatly depending on the type and degree of ordering achieved. To investigate these phenomena, the role of the order–disorder transition on phononic transport properties of Lennard–Jones type binary alloys is explored via nonequilibrium molecular dynamics simulations. Particular attention is paid to regimes in which the alloy is only partially ordered. It is shown that by varying the degree of ordering, the thermal conductivity of a binary alloy of fixed composition can be tuned across an order of magnitude at 10% of the melt temperature, and by a factor of three at 40% of the melt temperature.


2018 ◽  
Vol 20 (34) ◽  
pp. 22241-22248 ◽  
Author(s):  
Sayyed Jalil Mahdizadeh ◽  
Elaheh K. Goharshadi ◽  
Golnoosh Akhlamadi

Herein, the applicability and efficiency of two types of pillared graphene nanostructures, namely, (6,6)@G and (7,7)@G, were investigated as membranes in reverse osmosis seawater desalination using extensive nonequilibrium molecular dynamics simulations.


Author(s):  
Toshihiro Kaneko ◽  
Kenji Yasuoka ◽  
Ayori Mitsutake ◽  
Xiao Cheng Zeng

Multicanonical molecular dynamics simulations are applied, for the first time, to study the liquid-solid and solid-solid transitions in Lennard-Jones (LJ) clusters. The transition temperatures are estimated based on the peak position in the heat capacity versus temperature curve. For LJ31, LJ58 and LJ98, our results on the solid-solid transition temperature are in good agreement with previous ones. For LJ309, the predicted liquid-solid transition temperature is also in agreement with previous result.


Entropy ◽  
2021 ◽  
Vol 23 (2) ◽  
pp. 146
Author(s):  
Alessandro Coretti ◽  
Lamberto Rondoni ◽  
Sara Bonella

We illustrate how, contrary to common belief, transient Fluctuation Relations (FRs) for systems in constant external magnetic field hold without the inversion of the field. Building on previous work providing generalized time-reversal symmetries for systems in parallel external magnetic and electric fields, we observe that the standard proof of these important nonequilibrium properties can be fully reinstated in the presence of net dissipation. This generalizes recent results for the FRs in orthogonal fields—an interesting but less commonly investigated geometry—and enables direct comparison with existing literature. We also present for the first time a numerical demonstration of the validity of the transient FRs with nonzero magnetic field via nonequilibrium molecular dynamics simulations of a realistic model of liquid NaCl.


1993 ◽  
Vol 317 ◽  
Author(s):  
N.A. Marks ◽  
P. Guan ◽  
D.R. Mckenzie ◽  
B.A. PailThorpe

ABSTRACTMolecular dynamics simulations of nickel and carbon have been used to study the phenomena due to ion impact. The nickel and carbon interactions were described using the Lennard-Jones and Stillinger-Weber potentials respectively. The phenomena occurring after the impact of 100 e V to 1 keV ions were studied in the nickel simulations, which were both two and three-dimensional. Supersonic focussed collision sequences (or focusons) were observed, and associated with these focusons were unexpected sonic bow waves, which were a major energy loss mechanism for the focuson. A number of 2D carbon films were grown and the stress in the films as a function of incident ion energy was Measured. With increasing energy the stress changed from tensile to compressive and reached a maximum around 50 eV, in agreement with experiment.


Sign in / Sign up

Export Citation Format

Share Document