Diffusion in Li₂O studied by non-equilibrium molecular dynamics for 873 < T/K < 1603.

The use of non-equilibrium molecular dynamics facilitates the calculation of the cation diffusion constant of Li2O at temperatures too low to be accessible by other methods. Excellent agreement with experimental diffusion coefficients has been obtained over the temperature range 873 < T/K < 1603. Diffusion below 1200 K was shown to be dominated by a concerted nearest-neighbour hopping process, whereas in the high-temperature superionic region an additional mechanism involving a six-coordinate interstitial cation site in the anti-fluorite structure becomes increasingly dominant.

Molecular dynamics investigation of the disordered crystal structure of hexagonal LiBH4.

The crystal structure of the hexagonal phase of solid lithium borohydride (LiBH4) is studied by ab initio molecular dynamics simulations of both the low and high-temperature phases. A temperature range of 200-535 K is simulated with the aim of characterising the disorder in the high-temperature structure in detail. The mechanism and kinetics of the reorientational motion of the borohydride units (BH4(-)) are determined and are consistent with published neutron scattering experiments; it is found that rotational diffusivity increases by an order of magnitude at the phase transition temperature.

Ab initio nonequilibrium molecular dynamics in the solid superionic conductor LiBH4.

The color-diffusion algorithm is applied to ab initio molecular dynamics simulation of hexagonal LiBH(4) to determine the lithium diffusion coefficient and diffusion mechanisms. Even in the best solid lithium ion conductors, the time scale of ion diffusion is too long to be readily accessible by ab initio molecular dynamics at a reasonable computational cost. In our nonequilibrium method, rare events are accelerated by the application of an artificial external field acting on the mobile species; the system response to this perturbation is accurately described in the framework of linear response theory and is directly related to the diffusion coefficient, thus resulting in a controllable approximation.

Hydrogen bonding described using dispersion-corrected density functional theory.

In recent works, dispersion-corrected atom-centered potentials (DCACPs) were developed as a method to account for long-range dispersion forces between molecules in density functional theory calculations within the generalized gradient approximation (GGA). Here, we test the ability of DCACPs to improve the GGA treatment of hydrogen-bonded systems. We assessed both BLYP and dispersion-corrected BLYP with respect to benchmark calculations for the hydrogen bond lengths and binding energies of 20 complexes containing the elements C, H, N, O, and S.

Accurate DFT Descriptions for Weak Interactions of Molecules Containing Sulfur.

Dispersion corrected atom centered potentials (DCACPs) have been shown to significantly improve the density functional theory (DFT) description of weak interactions. In this work, we have calibrated a DCACP for sulfur in combination with the widely used Generalized Gradient Approximation (GGA) BLYP, thereby augmenting the existing library of DCACPs for the first- and second-row elements H, C, N, O, and rare gases. Three weakly bound complexes as well as elemental (orthorhombic) sulfur are used as test cases to evaluate the transferability of the DCACP to different chemical environments.