Subcritical and Supercritical Nanodroplet Evaporation: A Molecular Dynamics Investigation

Author(s):  
S. Mikkilineni ◽  
E. S. Landry ◽  
A. J. H. McGaughey

Molecular dynamics simulations are used to investigate the subcritical and supercritical evaporation of a Lennard-Jones (LJ) argon nanodroplet in its own vapor. Using a new technique to control both the ambient temperature and pressure, a range of conditions are considered to define a transition line between subcritical and supercritical evaporation. The evaporation is considered to be supercritical if the surface temperature of the droplet reaches the LJ argon critical temperature during its lifetime. Between ambient temperatures of 300 K and 800 K, the transition from subcritical to supercritical evaporation is observed to occur at an ambient pressure 1.4 times greater than the LJ argon critical pressure. For subcritical conditions, the droplet lifetimes obtained from the simulations are compared to independently predicted lifetimes from the D2 law.

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.


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.


2013 ◽  
Vol 423-426 ◽  
pp. 935-938 ◽  
Author(s):  
Ji Feng Li ◽  
Xiao Ping Zhao ◽  
Jian Liu

Molecular dynamics simulations were performed to calculate the melting points of perfect crystalline aluminum to high pressures. Under ambientpressure, there exhibits about 20% superheating before melting compared to the experimental melting point. Under high pressures, thecalculated melting temperature increases with the pressure but at a decreasing rate, which agrees well with the Simon's melting equation. Porosity effect was also studied for aluminum crystals with various initial porosity at ambient pressure, which shows that the equilibrium melting point decreases with the initial porosity as experiments expect.


2021 ◽  
Vol 55 (6) ◽  
Author(s):  
M. Gokhan Günay ◽  
Ubade Kemerli

A novel nano-scale pump that can transport atoms or small molecules with a peristaltic motion is designed. It is proven by molecular-dynamics simulations that the introduced nano-pump design works properly. The designed nano-pump consists of one main carbon nanotube named the flow tube and two rotors where multi-walled carbon nanotubes are attached. The pumping of helium atoms by the designed peristaltic carbon nano-pump is investigated by molecular-dynamics simulations. For varying rotor speeds and blade counts, time-averaged velocity, temperature, and pressure results of pumped helium atoms are calculated, and relationships between them are modeled as polynomial surfaces. The results showed that rotor frequency increases the velocity of helium linearly and the temperature and pressure of helium non-linearly. Furthermore, the blade count of the proposed mechanism does not substantially affect the velocity as per the previous studies in the literature.


2021 ◽  
Author(s):  
Martin P. Lautenschlaeger ◽  
Hans Hasse

It was shown recently that using the two-gradient method, thermal, caloric, and transport properties of fluids under quasi-equilibrium conditions can be determined simultaneously from nonequilibrium molecular dynamics simulations. It is shown here that the influence of shear stresses on these properties can also be studied using the same method. The studied fluid is described by the Lennard-Jones truncated and shifted potential with the cut-off radius r*c = 2.5σ. For a given temperature T and density ρ, the influence of the shear rate on the following fluid properties is determined: pressure p, internal energy u, enthalpy h, isobaric heat capacity cp, thermal expansion coefficient αp, shear viscosity η, and self-diffusion coefficient D. Data for 27 state points in the range of T ∈ [0.7, 8.0] and ρ ∈ [0.3, 1.0] are reported for five different shear rates (γ ̇ ∈ [0.1,1.0]). Correlations for all properties are provided and compared with literature data. An influence of the shear stress on the fluid properties was found only for states with low temperature and high density. The shear-rate dependence is caused by changes in the local structure of the fluid which were also investigated in the present work. A criterion for identifying the regions in which a given shear stress has an influence on the fluid properties was developed. It is based on information on the local structure of the fluid. For the self-diffusivity, shear-induced anisotropic effects were observed and are discussed.


2020 ◽  
Author(s):  
Nicholas Francia ◽  
Louise S. Price ◽  
Jonas Nyman ◽  
Sarah (Sally) Price ◽  
Matteo Salvalaglio

<p>Crystal structure prediction methods are prone to overestimate the number of potential polymorphs of organic molecules. In this work, we aim to reduce the overprediction by systematically applying molecular dynamics simulations and biased sampling methods to cluster subsets of structures that can easily interconvert at finite temperature and pressure. Following this approach, we rationally reduce the number of predicted putative polymorphs in CSP-generated crystal energy landscapes. This uses an unsupervised clustering approach to analyze independent finite-temperature molecular dynamics trajectories and hence identify a representative structure of each cluster of distinct lattice energy minima that are effectively equivalent at finite temperature and pressure. Biased simulations are used to reduce the impact of limited sampling time and to estimate the work associated with polymorphic transformations. We demonstrate the proposed systematic approach by studying the polymorphs of urea and succinic acid, reducing an initial set of over 100 energetically plausible CSP structures to 12 and 27 respectively, including the experimentally known polymorphs. The simulations also indicate the types of disorder and stacking errors that may occur in real structures.<br></p>


2013 ◽  
Vol 139 (16) ◽  
pp. 164106 ◽  
Author(s):  
Ross A. Lippert ◽  
Cristian Predescu ◽  
Douglas J. Ierardi ◽  
Kenneth M. Mackenzie ◽  
Michael P. Eastwood ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document