Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient and accurate determination of lattice-vacancy diffusion coefficients via non equilibrium ab initio molecular dynamics
75
0
0.0
(
0
)
Added by
Davide Giuseppe Sangiovanni
Publication date
2015
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We revisit the color-diffusion algorithm [P. C. Aeberhard et al., Phys. Rev. Lett. 108, 095901 (2012)] in nonequilibrium ab initio molecular dynamics (NE-AIMD), and propose a simple efficient approach for the estimation of monovacancy jump rates in crystalline solids at temperatures well below melting. Color-diffusion applied to monovacancy migration entails that one lattice atom (colored-atom) is accelerated toward the neighboring defect-site by an external constant force F. Considering bcc molybdenum between 1000 and 2800 K as a model system, NE-AIMD results show that the colored-atom jump rate k_{NE} increases exponentially with the force intensity F, up to F values far beyond the linear-fitting regime employed previously. Using a simple model, we derive an analytical expression which reproduces the observed k_{NE}(F) dependence on F. Equilibrium rates extrapolated by NE-AIMD results are in excellent agreement with those of unconstrained dynamics. The gain in computational efficiency achieved with our approach increases rapidly with decreasing temperatures, and reaches a factor of four orders of magnitude at the lowest temperature considered in the present study.
rate research
Read More
We have performed quasielastic neutron scattering (QENS) experiments up to 1243 K and ab-initio molecular dynamics (AIMD) simulations to investigate the Na diffusion in various phases of NaAlSiO4 (NASO), namely, low-carnegieite (L-NASO; trigonal), high-carnegieite (H-NASO; cubic) and nepheline (N-NASO; hexagonal) phases. The QENS measurements reveal Na ions localized diffusion behavior in L-NASO and N-NASO, but long-range diffusion behavior in H-NASO. The AIMD simulation supplemented the QENS measurements and showed that excess Na ions in H-NASO enhance the host network flexibility and activate the AlO4/SiO4 tetrahedra rotational modes. These framework modes enable the long-range diffusion of Na across a pathway of interstitial sites. The simulations also show Na diffusion in Na-deficient N-NASO through vacant Na sites along the hexagonal c-axis.
We present a comprehensive ab initio study of structural, electronic, lattice dynamical and electron-phonon coupling properties of the Bi(111) surface within density functional perturbation theory. Relativistic corrections due to spin-orbit coupling are consistently taken into account. As calculations are carried out in a periodic slab geometry, special attention is given to the convergence with respect to the slab thickness. Although the electronic structure of Bi(111) thin films varies significantly with thickness, we found that the lattice dynamics of Bi(111) is quite robust and appears converged already for slabs as thin as 6 bilayers. Changes of interatomic couplings are confined mostly to the first two bilayers, resulting in super-bulk modes with frequencies higher than the optic bulk spectrum, and in an enhanced density of states at lower frequencies for atoms in the first bilayer. Electronic states of the surface band related to the outer part of the hole Fermi surfaces exhibit a moderate electron-phonon coupling of about 0.45, which is larger than the coupling constant of bulk Bi. States at the inner part of the hole surface as well as those forming the electron pocket close to the zone center show much increased couplings due to transitions into bulk projected states near Gamma_bar. For these cases, the state dependent Eliashberg functions exhibit pronounced peaks at low energy and strongly deviate in shape from a Debye-like spectrum, indicating that an extraction of the coupling strength from measured electronic self-energies based on this simple model is likely to fail.
We have performed quasielastic and inelastic neutron scattering (QENS and INS) measurements from 300 K to 1173 K to investigate the Na-diffusion and underlying host dynamics in Na2Ti3O7. The QENS data show that the Na atoms undergo localized jumps up to 1173 K. The ab-initio molecular dynamics (AIMD) simulations supplement the measurements and show 1-d long-ranged diffusion along the a-axis above 1500 K. The simulations indicate that the occupancy of the interstitial site is critical for long-range diffusion. The nudged-elastic-band (NEB) calculation confirmed that the activation energy barrier is lowest for diffusion along the a-axis. In the experimental phonon spectra the peaks at 10 and 14 meV are dominated by Na dynamics that disappear on warming, suggesting low-energy phonons significantly contribute to large Na vibrational amplitude at elevated temperatures that enhances the Na hopping probability. We have also calculated the mode Gruneisen parameters of the phonons and thereby calculated the volume thermal expansion coefficient, which is found to be in excellent agreement with available experimental data.
Available simulation methods, suitable to describe solid-solid phase transitions occurring upon increasing of presssure and/or temperature, are based on empirical interatomic potentials: this restriction reduces the predictive power, and thus the general usefulness of numeric simulations in this very relevant field. We present a new simulation scheme which allows, for the first time, the simulation of these phenomena with the correct quantum-mechanical description of interatomic forces and internal stress, along with the correct statistical mechanics of ionic degrees of freedom. The method is obtained by efficiently combining the Car-Parrinello method for ab- initio molecular dynamics with the Parrinello Rahman method to account for a variable cell shape. Within this scheme phase trasformations may spontaneously take place during the simulation with variation of external pressure and/or temperature. The validity of the method is demonstrated by simulating the metal-insulator transition in Silicon (from diamond structure to simple hexagonal structure) under high pressure.
Despite their rich information content, electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized. We introduce a transferable high-fidelity neural network representation of such data in the form of tight-binding Hamiltonians for crystalline materials. This predictive representation of ab initio electronic structure, combined with machine-learning boosted molecular dynamics, enables efficient and accurate electronic evolution and sampling. When applied to a one-dimension charge-density wave material, carbyne, we are able to compute the spectral function and optical conductivity in the canonical ensemble. The spectral functions evaluated during soliton-antisoliton pair annihilation process reveal significant renormalization of low-energy edge modes due to retarded electron-lattice coupling beyond the Born-Oppenheimer limit. The availability of an efficient and reusable surrogate model for the electronic structure dynamical system will enable calculating many interesting physical properties, paving way to previously inaccessible or challenging avenues in materials modeling.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
