Calculations of two- and three-body dispersion coefficients for ions in crystals

Environmental modifications of atomic properties: The ground and 1s2p excited states of compressed helium

Atoms remaining as recognizably distinct constituents of bulk condensed phases can have properties modified from those of the isolated species. Dense helium bubbles at high pressures are a common form of radiation damage degrading the mechanical and electrical properties of host materials. Detailed knowledge is critical for predicting their long term performance. Modifications of the ground and first singlet excited states of confined compressed helium are investigated using an entirely non-empirical theory based on the results of ab initio self-consistent field calculations with corrections for the effects of electron correlation. For finite sized portions representing bulk condensed fcc and bcc phases of helium atoms, Hartree-Fock wavefunctions, energies, and charge distributions were computed as a function of different atomic densities using two models. The first model for the first excited state localizes the excitation on the central atom; in the second model, this is partially delocalized over the closest atomic neighbors. Total energies for the finite size portions are derived by adding the inter-atomic dispersive attractions and a density functional description of the short-range inter-atomic correlation energy. The experimental energy of the first allowed electronic transition increases with density being larger than in an isolated atom. The intra-atomic correlation energy does not contribute to this energy shift. The calculated energy shifts agree well with experiment for both bulk solid and liquid helium. The 2p orbital is increasingly compressed by density enhancement, thus generating the energy shifts. Consequently, calculations of the inelastic electron scattering cross sections are substantially incorrect if the compression of the final 1s2p state is not included. The character of the excitations is examined, and it is argued that these are of Frenkel rather than the Wannier type.

Analysis of the relative stability of lithium halide crystal structures: Density functional theory and classical models

All lithium halides exist in the rock salt crystal structure under ambient conditions. In contrast, common lithium halide classical force fields more often predict wurtzite as the stable structure. This failure of classical models severely limits their range of application in molecular simulations of crystal nucleation and growth. Employing high accuracy density functional theory (DFT) together with classical models, we examine the relative stability of seven candidate crystal structures for lithium halides. We give a detailed examination of the influence of DFT inputs, including the exchange-correlation functional, basis set, and dispersion correction. We show that a high-accuracy basis set, along with an accurate description of dispersion, is necessary to ensure prediction of the correct rock salt structure, with lattice energies in good agreement with the experiment. We also find excellent agreement between the DFT-calculated rock salt lattice parameters and experiment when using the TMTPSS-rVV10 exchange-correlation functional and a large basis set. Detailed analysis shows that dispersion interactions play a key role in the stability of rock salt over closely competing structures. Hartree-Fock calculations, where dispersion interactions are absent, predict the rock salt structure only for LiF, while LiCl, LiBr, and LiI are more stable as wurtzite crystals, consistent with radius ratio rules. Anion-anion second shell dispersion interactions overcome the radius ratio rules to tip the structural balance to rock salt. We show that classical models can be made qualitatively correct in their structural predictions by simply scaling up the pairwise additive dispersion terms, indicating a pathway toward better lithium halide force fields.

First-principle modelling of forsterite surface properties: Accuracy of methods and basis sets

The seven main crystal surfaces of forsterite (Mg2 SiO4 ) were modeled using various Gaussian-type basis sets, and several formulations for the exchange-correlation functional within the density functional theory (DFT). The recently developed pob-TZVP basis set provides the best results for all properties that are strongly dependent on the accuracy of the wavefunction. Convergence on the structure and on the basis set superposition error-corrected surface energy can be reached also with poorer basis sets. The effect of adopting different DFT functionals was assessed. All functionals give the same stability order for the various surfaces. Surfaces do not exhibit any major structural differences when optimized with different functionals, except for higher energy orientations where major rearrangements occur around the Mg sites at the surface or subsurface. When dispersions are not accounted for, all functionals provide similar surface energies. The inclusion of empirical dispersions raises the energy of all surfaces by a nearly systematic value proportional to the scaling factor s of the dispersion formulation. An estimation for the surface energy is provided through adopting C6 coefficients that are more suitable than the standard ones to describe O-O interactions in minerals. A 2 × 2 supercell of the most stable surface (010) was optimized. No surface reconstruction was observed. The resulting structure and surface energy show no difference with respect to those obtained when using the primitive cell. This result validates the (010) surface model here adopted, that will serve as a reference for future studies on adsorption and reactivity of water and carbon dioxide at this interface.

The effect of dispersion interactions on the properties of LiF in condensed phases

Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap repulsion and polarization terms of the interaction potential are derived on a purely non-empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parametrizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained from direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the densities of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.

Van der Waals interactions in ionic and semiconductor solids

van der Waals (vdW) energy corrected density-functional theory [Phys. Rev. Lett. 102, 073005 (2009)] is applied to study the cohesive properties of ionic and semiconductor solids (C, Si, Ge, GaAs, NaCl, and MgO). The required polarizability and dispersion coefficients are calculated using the dielectric function obtained from time-dependent density-functional theory. Coefficients for "atoms in the solid" are then calculated from the Hirshfeld partitioning of the electron density. It is shown that the Clausius-Mossotti equation that relates the polarizability and the dielectric function is accurate even for covalently-bonded semiconductors. We find an overall improvement in the cohesive properties of Si, Ge, GaAs, NaCl, and MgO, when vdW interactions are included on top of the Perdew-Burke-Ernzerhof or Heyd-Scuseria-Ernzerhof functionals. The relevance of our findings for other solids is discussed.

System-dependent dispersion coefficients for the DFT-D3 treatment of adsorption processes on ionic surfaces

Dispersion-corrected density functional theory calculations (DFT-D3) were performed for the adsorption of CO on MgO and C(2) H(2) on NaCl surfaces. An extension of our non-empirical scheme for the computation of atom-in-molecules dispersion coefficients is proposed. It is based on electrostatically embedded M(4)X(4) (M=Na, Mg) clusters that are used in TDDFT calculations of dynamic dipole polarizabilities. We find that the C(MM)(6) dispersion coefficients for bulk NaCl and MgO are reduced by factors of about 100 and 35 for Na and Mg, respectively, compared to the values of the free atoms. These are used in periodic DFT calculations with the revPBE semi-local density functional. As demonstrated by calculations of adsorption potential energy curves, the new C(6) coefficients lead to much more accurate energies (E(ads)) and molecule-surface distances than with previous DFT-D schemes. For NaCl/C(2) H(2) we obtained at the revPBE-D3(BJ) level a value of E(ads) =-7.4 kcal mol(-1) in good agreement with experimental data (-5.7 to -7.1 kcal mol(-1)). Dispersion-uncorrected DFT yields an unbound surface state. For the MgO/CO system, the computed revPBE-D3(BJ) value of E(ads) =-4.1 kcal mol(-1) is also in reasonable agreement with experimental results (-3.0 kcal mol(-1)) when thermal corrections are taken into account. Our new dispersion correction also improves computed lattice constants of the bulk systems significantly compared to plain DFT or previous DFT-D results. The extended DFT-D3 scheme also provides accurate non-covalent interactions for ionic systems without empirical adjustments and is suggested as a general tool in surface science.

Ab initio DFT+U study of He atom incorporation into UO(2) crystals

We present and discuss results of the density functional theory (DFT) for perfect UO(2) crystals with He atoms in octahedral interstitial positions therein. We have calculated basic bulk crystal properties and He incorporation energies into the low temperature anti-ferromagnetic UO(2) phase using several exchange-correlation functionals within the spin-polarized local density (LDA) and generalized gradient (GGA) approximations. In all DFT calculations we included the on-site correlation corrections using the Hubbard model (DFT+U approach). We analysed a potential crystalline symmetry reduction from tetragonal down to orthorhombic structure and confirmed the presence of the Jahn-Teller effect in a perfect UO(2). We discuss also the problem of a conducting electronic state arising when He is placed into a tetragonal antiferromagnetic phase of UO(2) commonly used in defect modelling. Consequently, we found a specific monoclinic lattice distortion which allowed us to restore the semiconducting state and properly estimate He incorporation energies. Unlike the bulk properties, the He incorporation energy strongly depends on several factors, including the supercell size, the use of spin polarization, the exchange-correlation functionals and on-site correlation corrections. We compare our results for the He incorporation with the previous shell model and ab initio DFT calculations.

Direct observation of H2 binding to a metal oxide surface

Inelastic neutron scattering is used to probe the dynamical response of H2 films adsorbed on MgO(100) as a function of film thickness. Concomitant diffraction measurements and a reduced-dimensionality quantum dynamical model provide insight into the molecule-surface interaction potential. At monolayer thickness, the rotational motion is strongly influenced by the surface, so that the molecules behave like quasiplanar rotors. These findings have a direct impact on understanding how molecular hydrogen binds to the surface of materials used in catalytic and storage applications.

The stability of fission products in uranium dioxide

A review of experimental data concerning the behaviour of fission products in nuclear fuels is used to illustrate the significant variation in solubility exhibited by the different species. To understand the reasons for this variation, it is necessary to obtain a reliable estimate of the solution energies and thus to determine the most stable solution site. This we suggest will be critical in predicting the behaviour of nuclear fuels in both accident and normal operating conditions. We have therefore used the M ott-Littleton simulation technique to calculate solution energies for the fission products Br, Kr, Rb, Sr, Y, Zr, Te, I, Xe, Cs, Ba, La and Ce in UO 2. We considered solution at both uranium and oxygen vacancies, the interstitial site and at the di-, tri- and tetra-vacancy complexes. Non-stoichiometry and variable charge state are important components of the model. From these results we conclude that the solubility is significantly affected by non-stoichiometry. In UO 2 and UO 2-x ., products such as Cs, Rb and Ba are thermodynamically more stable as binary oxide precipitates. Conversely, Y, La and Sr are soluble in UO 2 and UO 2+x , while Cs, Rb, Sr and Ba are only soluble in UO 2+x .The behaviour of I, Br and Te is complicated by the fact that these species are most stable as anions in UO 2 and UO 2-x but as cations in UO 2+x . In our model, Zr and the inert gas species Xe and Kr are always predicted to be insoluble, while CeO 2 will form a solid solution with UO 2.