Molecular Dynamics Simulations of Atoms Diffusion in Solid

2018 ◽  
Vol 15 ◽  
pp. 51-64
Author(s):  
Yu Lu Zhou ◽  
Xiao Ma Tao ◽  
Qing Hou ◽  
Yi Fang Ouyang
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Md Simulations ◽  
Einstein Relation ◽  
Interatomic Potentials ◽  
Body System ◽  
Many Body ◽  
Point Particles ◽  
Bulk Surface ◽  
Dynamics Simulations

Molecular dynamics (MD) simulations, which treat atoms as point particles and trace their individual trajectories, are always employed to investigate the transport properties of a many-body system. The diffusion coefficients of atoms in solid can be obtained by the Einstein relation and the Green-Kubo relation. An overview of the MD simulations of atoms diffusion in the bulk, surface and grain boundary is provided. We also give an example of the diffusion of helium in tungsten to illustrate the procedure, as well as the importance of the choice of interatomic potentials. MD simulations can provide intuitive insights into the atomic mechanisms of diffusion.

scholarly journals Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau
Keyword(s):  
Carbon Dioxide ◽  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Hydrodynamic Radius ◽  
Md Simulations ◽  
Molecular Models ◽  
Experimental Values ◽  
Carbon Dioxide Co2 ◽  
Dynamics Simulations ◽  
Apparent Disagreement

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.

MOLECULAR DYNAMICS APPROACH TO CORRELATION CLUSTERING

2008 ◽  
Vol 19 (09) ◽  
pp. 1349-1358 ◽  
Author(s):  
R. SUMI ◽  
Z. NÉDA
Keyword(s):  
Phase Transition ◽  
Molecular Dynamics ◽  
Mechanical Equilibrium ◽  
Body System ◽  
Many Body ◽  
Competing Interactions ◽  
Correlation Clustering ◽  
Simulated System ◽  
The One ◽  
Dynamics Simulations

A many-body system with co-existing attractive and repulsive interactions is considered on a ring. The competing interactions lead to a frustration similar with the one existing in Correlation Clustering (CC). The optimal mechanical equilibrium of the system is searched by molecular dynamics simulations. As a function of the disorder quenched in the interactions, the system exhibits the phase-transition recently reported in CC. The simulated system can be considered as a continuous and efficient approach to the otherwise discrete, NP hard CC problem.

Molecular Dynamics simulations of displacement cascades: role of the interatomic potentials and of the potential hardening

2000 ◽  
Vol 650 ◽  
Author(s):  
C.S. Becquart ◽  
C. Domain ◽  
A. Legris ◽  
J.C. van Duysen
Keyword(s):  
Molecular Dynamics ◽  
Thermal Characteristics ◽  
Md Simulations ◽  
Interatomic Potentials ◽  
Static Properties ◽  
Displacement Cascades ◽  
Energy Transfers ◽  
Primary Damage ◽  
Dynamics Simulations

ABSTRACTThe role of the interatomic potentials on the primary damage has been investigated by Molecular Dynamics (MD) simulations of displacement cascades with three different interatomic potentials dedicated to α-Fe. The primary damage, caused by the neutron interaction with the matter, has been found to be potential sensitive. We have investigated the equilibrium parts of the potential as well as the “short distance interactions” which appear to have a strong influence on the cascade morphology and defects distribution at the end of the cascade. The static properties as well as dynamical (thermal) characteristics of the potentials have been considered; the kinetic and potential energy transfers during the collisions have also been studied.

Molecular Dynamics Simulations of Ionic Liquids: Influence of Polarization on IL Structure and Ion Transport

2008 ◽  
Vol 1082 ◽  
Author(s):  
Oleg Borodin
Keyword(s):  
Molecular Dynamics ◽  
Ionic Liquid ◽  
Ionic Liquids ◽  
Condensed Phase ◽  
Wide Class ◽  
Md Simulations ◽  
Liquid Structure ◽  
Many Body ◽  
Polarizable Force Field ◽  
Dynamics Simulations

ABSTRACTMany-body polarizable force field has been developed and validated for a wide class of ionic liquids. Classical molecular dynamics (MD) simulations have been performed on 29 ionic liquids. This presentation will focus on ability of developed force fields to predict condensed phase properties and on understanding the influence of many-body polarizable interactions on the ionic liquid structure and transport.

Predicting Methane Diffusivity in Polymeric Membranes by Molecular Dynamics

2015 ◽  
Vol 1119 ◽  
pp. 461-465
Author(s):  
M.K. Hadj-Kali ◽  
A. Bessadok-Jemai ◽  
S. Haider ◽  
Y. Alzeghayer
Keyword(s):  
Molecular Dynamics ◽  
Small Molecules ◽  
Diffusion Coefficients ◽  
Md Simulations ◽  
Polymeric Membranes ◽  
Transport Mechanisms ◽  
Fluid Equation ◽  
Cell Modules ◽  
Dynamics Simulations ◽  
Materials Studio

Diffusion coefficients of methane (CH4) have been obtained by Molecular Dynamics (MD) simulations combined with Einstein fluid equation. Three polymers were considered, namely polyethylene, polypropylene and poly (cis-1,4-butadiene). All calculations were performed by means of Polymer Builder and Amorphous Cell modules within Materials Studio (Accelrys). The obtained diffusivity results are within the range of published results for similar small molecules. Molecular dynamics simulations proved to be a useful tool for understanding the detailed descriptions and transport mechanisms occurring within the material.

Lanthanoids(III) and actinoids(III) in water: Diffusion coefficients and hydration enthalpies from polarizable molecular dynamics simulations

2012 ◽  
Vol 85 (1) ◽  
pp. 237-246 ◽  
Author(s):  
Fausto Martelli ◽  
Sacha Abadie ◽  
Jean-Pierre Simonin ◽  
Rodolphe Vuilleumier ◽  
Riccardo Spezia
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Heavy Ions ◽  
Md Simulations ◽  
Contradictory Result ◽  
Dielectric Friction ◽  
Water Interaction ◽  
Diffusive Regime ◽  
Dynamics Simulations ◽  
Good Agreement

By using polarizable molecular dynamics (MD) simulations of lanthanoid(III) and actinoid(III) ions in water, we obtained ionic diffusion coefficients and hydration enthalpies for both series. These values are in good agreement with experiments. Simulations thus allow us to relate them to microscopic structure. In particular, across the series the diffusion coefficients decrease, reflecting the increase of ion–water interaction. Hydration enthalpies also show that interactions increase from light to heavy ions in agreement with experiment. The apparent contradictory result of the decrease of the diffusion coefficient with decreasing ionic radius is tentatively explained in terms of dielectric friction predominance on Stokes’ diffusive regime.

Predicting XAFS scattering path cumulants and XAFS spectra for metals (Cu, Ni, Fe, Ti, Au) using molecular dynamics simulations

2013 ◽  
Vol 20 (4) ◽  
pp. 555-566 ◽  
Author(s):  
M. A. Karolewski ◽  
R. G. Cavell ◽  
R. A. Gordon ◽  
C. J. Glover ◽  
M. Cheah ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Md Simulations ◽  
Embedded Atom Method ◽  
Interatomic Potentials ◽  
Embedded Atom ◽  
Scaling Parameters ◽  
Order Of Magnitude ◽  
Experimental Errors ◽  
Dynamics Simulations ◽  
X Ray Absorption

The ability of molecular dynamics (MD) simulations to support the analysis of X-ray absorption fine-structure (XAFS) data for metals is evaluated. The low-order cumulants (ΔR, σ2,C3) for XAFS scattering paths are calculated for the metals Cu, Ni, Fe, Ti and Au at 300 K using 28 interatomic potentials of the embedded-atom method type. The MD cumulant predictions were evaluated within a cumulant expansion XAFS fitting model, using global (path-independent) scaling factors. Direct simulations of the corresponding XAFS spectra, χ(R), are also performed using MD configurational data in combination with theFEFFab initiocode. The cumulant scaling parameters compensate for differences between the real and effective scattering path distributions, and for any errors that might exist in the MD predictions and in the experimental data. The fitted value of ΔRis susceptible to experimental errors and inadvertent lattice thermal expansion in the simulation crystallites. The unadjusted predictions of σ2vary in accuracy, but do not show a consistent bias for any metal except Au, for which all potentials overestimate σ2. The unadjustedC3predictions produced by different potentials display only order-of-magnitude consistency. The accuracy of direct simulations of χ(R) for a given metal varies among the different potentials. For each of the metals Cu, Ni, Fe and Ti, one or more of the tested potentials was found to provide a reasonable simulation of χ(R). However, none of the potentials tested for Au was sufficiently accurate for this purpose.

Molecular Dynamics Simulations of Thermal Conductivity of Bi2Te3 Nanowires

2009 ◽  
Author(s):  
Bo Qiu ◽  
Lin Sun ◽  
Xiulin Ruan
Keyword(s):  
Thermal Conductivity ◽  
Molecular Dynamics ◽  
Surface Roughness ◽  
Smooth Surface ◽  
Lattice Thermal Conductivity ◽  
Phonon Scattering ◽  
Md Simulations ◽  
Interatomic Potentials ◽  
Phonon Confinement ◽  
Dynamics Simulations

In this paper, by employing the previously developed two-body interatomic potentials for bismuth telluride, molecular dynamics (MD) simulations are used to describe the thermoelectric properties, namely the lattice thermal conductivity, of Bi2Te3 nanowires. Cylindrical nanowires with both smooth surface and sawtooth surface roughness are studied, aiming at revealing the effects of phonon confinement in 1-D structures, phonon boundary scatterings and surface roughness on the lattice thermal conductivity of Bi2Te3 nanowires. In the end, the influence of various phonon scattering mechanisms on the nanostructures under study are summarized, possible paths to reduce lattice thermal conductivity in nanostructured Bi2Te3, which is favorable for enhancing thermoelectric performance, are pointed out.

scholarly journals Bond order redefinition needed to reduce inherent noise in molecular dynamics simulations

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Ibnu Syuhada ◽  
Nikodemus Umbu Janga Hauwali ◽  
Ahmad Rosikhin ◽  
Euis Sustini ◽  
Fatimah Arofiati Noor ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Bond Order ◽  
Interatomic Potential ◽  
Interatomic Potentials ◽  
Many Body ◽  
Range Interaction ◽  
Dynamics Simulations ◽  
Interaction Approach ◽  
Many Body Interactions

AbstractIn this work, we present the bond order redefinition needed to reduce the inherent noise in order to enhance the accuracy of molecular dynamics simulations. We propose defining the bond order as a fraction of energy distribution. It happens due to the character of the material in nature, which tries to maintain its environment. To show the necessity, we developed a factory empirical interatomic potential (FEIP) for carbon that implements the redefinition with a short-range interaction approach. FEIP has been shown to enhance the accuracy of the calculation of lattice constants, cohesive energy, elastic properties, and phonons compared to experimental data, and can even be compared to other potentials with the long-range interaction approach. The enhancements due to FEIP can reduce the inherent noise, then provide a better prediction of the energy based on the behaviour of the atomic environment. FEIP can also transform simple two-body interactions into many-body interactions, which is useful for enhancing accuracy. Due to implementing the bond order redefinition, FEIP offers faster calculations than other complex interatomic potentials.

scholarly journals Structural, Mechanical, and Dynamical Properties of Amorphous Li2CO3 from Molecular Dynamics Simulations

Crystals ◽  
2018 ◽  
Vol 8 (12) ◽  
pp. 473 ◽  
Author(s):  
Mahsa Ebrahiminia ◽  
Justin Hooper ◽  
Dmitry Bedrov
Keyword(s):  
Molecular Dynamics ◽  
Room Temperature ◽  
Md Simulations ◽  
Mechanical Analysis ◽  
Glassy Matrix ◽  
Many Body ◽  
Anion Diffusion ◽  
Similar Temperature Dependence ◽  
Similar Temperature ◽  
Dynamics Simulations

Structural, mechanical, and transport properties of amorphous Li2CO3 were studied using molecular dynamics (MD) simulations and a hybrid MD-Monte Carlo (MC) scheme. A many-body polarizable force field (APPLE&P) was employed in all simulations. Dynamic and mechanical properties of Dilithium carbonate, Li2CO3, in amorphous liquid and glassy phases were calculated over a wide temperature range. At higher temperatures, both anion and cation diffusion coefficients showed similar temperature dependence. However, below the glass transition temperature (T < 450 K) the anions formed a glassy matrix, while Li+ continued to be mobile, showing decoupling of cation and anion diffusion. The conductivity of Li+ at room temperature was estimated to be on the order of 10−6 S/cm. Mechanical analysis revealed that at room temperature the amorphous phase had a shear modulus of about 8 GPa, which was high enough to suppress Li metal dendrite growth on an electrode surface.

Sign in / Sign up

Export Citation Format

Share Document