scholarly journals Molecular Dynamics Simulation of Tolman Length and Interfacial Tension of Symmetric Binary Lennard–Jones Liquid

Symmetry ◽  
2021 ◽  
Vol 13 (8) ◽  
pp. 1376
Author(s):  
Hideki Kanda ◽  
Wahyudiono ◽  
Motonobu Goto

The Tolman length and interfacial tension of partially miscible symmetric binary Lennard–Jones (LJ) fluids (A, B) was revealed by performing a large-scale molecular dynamics (MD) simulation with a sufficient interfacial area and cutting distance. A unique phenomenon was observed in symmetric binary LJ fluids, where two surfaces of tension existed on both sides of an equimolar dividing surface. The range of interaction εAB between the different liquids and the temperature in which the two LJ fluids partially mixed was clarified, and the Tolman length exceeded 3 σ when εAB was strong at higher temperatures. The results show that as the temperature or εAB increases, the Tolman length increases and the interfacial tension decreases. This very long Tolman length indicates that one should be very careful when applying the concept of the liquid–liquid interface in the usual continuum approximation to nanoscale droplets and capillary phase separation in nanopores.

2008 ◽  
Vol 19 (09) ◽  
pp. 1315-1319 ◽  
Author(s):  
TIMOTHY C. GERMANN ◽  
KAI KADAU

By utilizing the molecular dynamics code SPaSM on Livermore's BlueGene/L architecture, consisting of 212 992 IBM PowerPC440 700 MHz processors, a molecular dynamics simulation was run with one trillion atoms. To demonstrate the practicality and future potential of such ultra large-scale simulations, the onset of the mechanical shear instability occurring in a system of Lennard-Jones particles arranged in a simple cubic lattice was simulated. The evolution of the instability was analyzed on-the-fly using the in-house developed massively parallel graphical object-rendering code MD_render.


2004 ◽  
Vol 821 ◽  
Author(s):  
Markus J. Buehler ◽  
Alexander Hartmaier ◽  
Mark Duchaineau ◽  
Farid F. Abraham ◽  
Huajian Gao

AbstractWe analyze a large-scale molecular dynamics simulation of work hardening in a ductile model material comprising of 500 million atoms interacting with a Lennard-Jones pair potential within a classical molecular dynamics scheme. With tensile loading, we observe emission of thousands of dislocations from two sharp cracks. The dislocations interact in a complex way, revealing three fundamental mechanisms of work-hardening. These are (1) dislocation cutting processes, jog formation and generation of point defects; (2) activation of secondary slip systems by cross-slip; and (3) formation of sessile Lomer-Cottrell locks. The dislocations self-organize into a complex sessile defect topology. Our analysis illustrates mechanisms formerly only known from textbooks and observed indirectly in experiment. It is the first time that such a rich set of fundamental phenomena has been seen in a single computer simulation.


2019 ◽  
Author(s):  
Liqun Cao ◽  
Jinzhe Zeng ◽  
Mingyuan Xu ◽  
Chih-Hao Chin ◽  
Tong Zhu ◽  
...  

Combustion is a kind of important reaction that affects people's daily lives and the development of aerospace. Exploring the reaction mechanism contributes to the understanding of combustion and the more efficient use of fuels. Ab initio quantum mechanical (QM) calculation is precise but limited by its computational time for large-scale systems. In order to carry out reactive molecular dynamics (MD) simulation for combustion accurately and quickly, we develop the MFCC-combustion method in this study, which calculates the interaction between atoms using QM method at the level of MN15/6-31G(d). Each molecule in systems is treated as a fragment, and when the distance between any two atoms in different molecules is greater than 3.5 Å, a new fragment involved two molecules is produced in order to consider the two-body interaction. The deviations of MFCC-combustion from full system calculations are within a few kcal/mol, and the result clearly shows that the calculated energies of the different systems using MFCC-combustion are close to converging after the distance thresholds are larger than 3.5 Å for the two-body QM interactions. The methane combustion was studied with the MFCC-combustion method to explore the combustion mechanism of the methane-oxygen system.


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.


2018 ◽  
Vol 4 (4) ◽  
Author(s):  
Qiang Zhao ◽  
Yang Li ◽  
Zheng Zhang ◽  
Xiaoping Ouyang

The sputtering of graphite due to the bombardment of hydrogen isotopes is crucial to successfully using graphite in the fusion environment. In this work, we use molecular dynamics to simulate the sputtering using the large-scale atomic/molecular massively parallel simulator (lammps). The calculation results show that the peak values of the sputtering yield are between 25 eV and 50 eV. When the incident energy is greater than the energy corresponding to the peak value, a lower carbon sputtering yield is obtained. The temperature that is most likely to sputter is approximately 800 K for hydrogen, deuterium, and tritium. Below the 800 K, the sputtering yields increase with temperature. By contrast, above the 800 K, the yields decrease with increasing temperature. Under the same temperature and incident energy, the sputtering rate of tritium is greater than that of deuterium, which in turn is greater than that of hydrogen. When the incident energy is 25 eV, the sputtering yield at 300 K increases below an incident angle at 30 deg and remains steady after that.


2014 ◽  
Vol 1700 ◽  
pp. 61-66
Author(s):  
Guttormur Arnar Ingvason ◽  
Virginie Rollin

ABSTRACTAdding single walled carbon nanotubes (SWCNT) to a polymer matrix can improve the delamination properties of the composite. Due to the complexity of polymer molecules and the curing process, few 3-D Molecular Dynamics (MD) simulations of a polymer-SWCNT composite have been run. Our model runs on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), with a COMPASS (Condensed phase Optimized Molecular Potential for Atomistic Simulations Studies) potential. This potential includes non-bonded interactions, as well as bonds, angles and dihedrals to create a MD model for a SWCNT and EPON 862/DETDA (Diethyltoluenediamine) polymer matrix. Two simulations were performed in order to test the implementation of the COMPASS parameters. The first one was a tensile test on a SWCNT, leading to a Young’s modulus of 1.4 TPa at 300K. The second one was a pull-out test of a SWCNT from an originally uncured EPON 862/DETDA matrix.


Sign in / Sign up

Export Citation Format

Share Document