ABSTRACT
Molecular dynamics (MD) simulations have been carried out to study isotope effects in lithium hydride and lithium deuteride crystals. Quantum effects on nuclear motion have been included through a quantum thermal bath (QTB). The interatomic forces were described either within the density functional theory (DFT) in the generalized gradient approximation (GGA) or by the phenomenological approach using the shell model. For both models, the isotopic shift in the lattice parameter can be successfully predicted by QTB-MD simulations. The slope of the experimental isotopic shift in pressure is satisfactorily reproduced by QTB-MD within DFT-GGA, in contrast to both density functional perturbation theory and QTB-MD with the shell model. We have analyzed the reasons for these discrepancies through the vibrational densities of states and the isotopic shifts in bulk modulus. The results illustrate the importance of anharmonic contributions to vibrations and to the isotopic pressure shift between LiH and LiD.
Subject(s)
Deuterium/chemistry , Lithium/chemistry , Molecular Dynamics Simulation , Acoustics , Computer Simulation , Hot Temperature , Molecular Structure , Physics/methods , Pressure , Quantum Theory , TemperatureABSTRACT
The Schottky defects at both the flat MgO(001) surface and the monatomic-step edge have been investigated by equilibrium molecular dynamics simulations. The formation enthalpy as a function of the distance between the Mg and O vacancies that form a Schottky defect have been calculated and discussed. We conclude that a step edge is a stable location for a vacancy. The migration mechanism has been elucidated and an intermediate state has been identified. The associated activation enthalpies have been determined in the 700-1100 K temperature range. Both magnesium and oxygen vacancies at the surface are very mobile and can play a role during the crystal growth.
