Molecular Dynamics Simulation of Mass and Charge Transport in Superionic Conductors, and Structural Correlations in Chalcogenide Glasses

1988 ◽  
Vol 135 ◽  
Author(s):  
P. Vashishta ◽  
José P. Rino ◽  
Rajiv K. Kalia

AbstractStructural properties, single-particle dynamics, and the charge transport are studied in superionic conductor Ag2Se using the molecular dynamics (MD) technique. The calculations are based on a model of interionic potentials in which ions interact through Coulomb interaction, steric repulsion and charge-dipole interaction due to the large electronic polarizability of the selenium ions. Structural and dynamics correlations are studied at five temperatures in the superionic phase. Among the structural correlations the results are presented for partial pair correlation function, coordination numbers, bond angle distributions and wave-vector dependence of the Bragg intensities. Detailed comparison with neutron and x-ray single crystal diffraction experiments. The calculated temperature dependence of the self-diffusion constant of silver is in good agreement with the tracer diffusion measurements. The spectra of velocity autocorrelation functions and the frequency dependent ionic conductivity are calculated. The Haven's ratio is also in good agreement with experiments.Effective interatomic potentials consisting of two-body (steric effect, charge transfer and charge-dipole interactions) and three-body covalent forces are proposed for GeSe2. Using these interaction potentials in MD simulations, the nature of short-range and medium-range order is investigated in glassy and molten GeSe2. All the features in the static structure factor, S(q), including the first sharp diffraction peak (FSDP), are in good agreement with experiments. The FSDP arises from Ge-Ge and Ge-Se correlations between 4-8Å, and the anomalous decrease in its height on cooling is due to frustration enhanced by the increased density.

1992 ◽  
Vol 291 ◽  
Author(s):  
Andrew Horsfield ◽  
Paulette Clancy

ABSTRACTThe melting of crystalline silicon and the cooling of liquid silicon are investigated using Molecular Dynamics. Both the Stillinger-Weber (SW) potential and the Tight-Binding Bond Model are used to calculate the forces. The electrical properties are investigated using an empirical pseudopotential method with a plane wave basis. The melting point of the solid is found to be about 2300K. The dependency of this temperature with cell size is investigated. On cooling, there are changes in some of the properties of the liquid: the energy per particle decreases, the diffusion constant decreases, and the low frequency electrical conductivity decreases slightly as the temperature decreases. Between 1180K and 980K the liquid undergoes a transition to a glassy phase. There are large changes in the pair correlation function, the SW three-body energy distribution, the diffusion constant, the density of electron single particle states and the electrical conductivity. All of these changes are consistent with increased tetrahedral bonding.


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.


2006 ◽  
Vol 258-260 ◽  
pp. 522-530 ◽  
Author(s):  
José Pedro Rino ◽  
Paulo S. Branício ◽  
Denílson S. Borges

An effective inter-atomic potential is proposed in order to describe structural and dynamical properties of II-VI and III-V semiconductors. The interaction potential consists of twoand three-body interactions. The two-body term takes into account steric repulsion, charge-induce dipole interaction due to the electronic polarizability of ions, Coulomb interaction due to charge transfer between ions, and dipole-dipole (van der Waals) interactions. The three-body term, which has a modified Stillinger-Weber form, describes bond-bending as well as bond-stretching effects. Here we report the fitting and the application of this interaction potential for InP in the crystalline phase and for CdTe in the crystalline and liquid phases. The structural correlations are discussed through pair distribution, coordination number and bond-angle functions. Vibrational density of states for InP and CdTe as well as the static structure factor for liquid CdTe are in very good agreement with experimental data.


RSC Advances ◽  
2016 ◽  
Vol 6 (11) ◽  
pp. 9096-9105 ◽  
Author(s):  
Gholam Hossien Rounaghi ◽  
Mostafa Gholizadeh ◽  
Fatemeh Moosavi ◽  
Iman Razavipanah ◽  
Hossein Azizi-Toupkanloo ◽  
...  

The variation of molar conductance versus mole ratio for (kryptofix 22DD·La)3+ complex in methanol solution at different temperatures is in accordance with the variation of pair correlation function of oxygen atoms.


2021 ◽  
Author(s):  
Kazushi Fujimoto ◽  
Tetsuro Nagai ◽  
Tsuyoshi Yamaguchi

<div>The position-dependent diffusion coefficient along with free energy profile are important parameters needed to study mass transport in heterogeneous systems such as biological and polymer membranes, and molecular dynamics (MD) calculation is a popular tool to obtain them. Among many methodologies, the Marrink-Berendsen (MB) method is often employed to calculate the position-dependent diffusion coefficient, in which the autocorrelation function of the force on a fixed molecule is related to the friction on the molecule. However, the diffusion coefficient is shown to be affected by the period of the removal of the center-of-mass velocity, which is necessary when performing MD calculations using the Ewald method for Coulombic interaction. We have clarified theoretically in this study how this operation affects the diffusion coefficient calculated by the MB method, and the theoretical predictions are proven by MD calculations. Therefore, we succeeded in providing guidance on how to select an appropriate the period of the removal of the center-of-mass velocity in estimating the position-dependent diffusion coefficient by the MB method. This guideline is applicable also to the Woolf-Roux method.</div>


1993 ◽  
Vol 07 (09n10) ◽  
pp. 1779-1788 ◽  
Author(s):  
JASON A.C. GALLAS ◽  
HANS J. HERRMANN ◽  
STEFAN SOKOLOWSKI

When sand or other granular materials are shaken, poured or sheared many intriguing phenomena can be observed. We will model the granular medium by a packing of elastic spheres and simulate it via Molecular Dynamics. Dissipation of energy and shear friction at collisions are included. The onset of fluidization can be determined and is in good agreement with experiments. On a vibrating plate we observe the formation of convection cells due to walls or amplitude modulations. Density and velocity profiles on conveyor belts are measured and the influence of an obstacle discussed. We mention various types of rheology for flow down an inclined chute or through a pipe and outflowing containers.


2012 ◽  
Vol 26 (20) ◽  
pp. 1250117 ◽  
Author(s):  
L. T. VINH ◽  
N. V. HUY ◽  
P. K. HUNG

Molecular dynamics simulation is carried out for liquid SiO 2 at pressure ranged from zero to 30 GPa and by using BKS, Born–Mayer type and Morse–Stretch potentials. The constructed models reproduce well the experimental data in terms of mean coordination number, bond angle and pair radial distribution function. Furthermore, the density of all samples can be expressed by a linear function of fractions SiO x. It is found that the topology of units SiO x and linkages OSi y is unchanged upon compression although the liquid undergoes substantial change in its network structure. Consequently, the partial bond angle distribution for SiO x and OSi y is identical for all samples constructed by the same potential. This result allows to establishing a simple expression between total bond angle distribution (BAD) and fraction of SiO x and OSi y. The simulation shows a good agreement between the calculation and simulation results for both total O–Si–O and Si–O–Si BADs. This supports a technique to estimate amount of units SiO x and linkages OSi y on base of total Si–O–Si and O–Si–O BADs measured experimentally.


1980 ◽  
Vol 35 (5) ◽  
pp. 493-499 ◽  
Author(s):  
Isao Okada ◽  
Ryuzo Takagi ◽  
Kazutaka Kawamura

Abstract A new transport property, the self-exchange velocity (SEV) of neighbouring unlike ions, has been evaluated from molecular dynamics simulations of molten LiCl, RbCl and LiRbCl2 at 1100 K and the mixture at 750 K. From the increase of the SEV's in the order Rb+ (pure salt) <Li+ (mixture) < Rb+ (mixture) < Li+ (pure salt), it is conjectured that there is a strong correlation between the SEV’s and the internal mobilities. An interpretation of the Chemla effect in its dependence on temperature is given. The pair correlation functions and the self-diffusion coefficients are also calculated and discussed.


Sign in / Sign up

Export Citation Format

Share Document