Accurate screened exchange band structures for the transition metal monoxides MnO, FeO, CoO and NiO.

We report calculations of the band structures and density of states of the four transition metal monoxides MnO, FeO, CoO and NiO using the hybrid density functional sX-LDA ('screened exchange local density approximation'). Late transition metal oxides are prototypical examples of strongly correlated materials, which pose challenges for electronic structure methods. We compare our results with available experimental data and show that our calculations generally yield accurate predictions for the fundamental band gaps and valence bands, in favourable agreement with previously reported theoretical studies. For MnO, the band gaps are still underestimated, suggesting additional many-body effects that are not captured by our screened hybrid functional approach.

Nature of the electronic band gap in lanthanide oxides

Accurate electronic structures of the technologically important lanthanide/rare-earth sesquioxides (Ln2O3, with Ln=La, ⋯,Lu) and CeO2 have been calculated using hybrid density functionals HSE03, HSE06, and screened exchange (sX-LDA). We find that these density functional methods describe the strongly correlated Ln f electrons as well as the recent G0W0@LDA+U results, generally yielding the correct band gaps and trends across the Ln period. For HSE, the band gap between O 2p states and lanthanide 5d states is nearly independent of the lanthanide, while the minimum gap varies as filled or empty Ln 4f states come into this gap. sX-LDA predicts the unoccupied 4f levels at higher energies, which leads to a better agreement with experiments for Sm2O 3, Eu2O3, and Yb2O3. © 2013 American Physical Society.

Communication: energy benchmarking with quantum Monte Carlo for water nano-droplets and bulk liquid water.

We show the feasibility of using quantum Monte Carlo (QMC) to compute benchmark energies for configuration samples of thermal-equilibrium water clusters and the bulk liquid containing up to 64 molecules. Evidence that the accuracy of these benchmarks approaches that of basis-set converged coupled-cluster calculations is noted. We illustrate the usefulness of the benchmarks by using them to analyze the errors of the popular BLYP approximation of density functional theory (DFT). The results indicate the possibility of using QMC as a routine tool for analyzing DFT errors for non-covalent bonding in many types of condensed-phase molecular system.

Between scylla and charybdis: Hydrophobic graphene-guided water diffusion on hydrophilic substrates

The structure of water confined in nanometer-sized cavities is important because, at this scale, a large fraction of hydrogen bonds can be perturbed by interaction with the confining walls. Unusual fluidity properties can thus be expected in the narrow pores, leading to new phenomena like the enhanced fluidity reported in carbon nanotubes. Crystalline mica and amorphous silicon dioxide are hydrophilic substrates that strongly adsorb water. Graphene, on the other hand, interacts weakly with water. This presents the question as to what determines the structure and diffusivity of water when intercalated between hydrophilic substrates and hydrophobic graphene. Using atomic force microscopy, we have found that while the hydrophilic substrates determine the structure of water near its surface, graphene guides its diffusion, favouring growth of intercalated water domains along the C-C bond zigzag direction. Molecular dynamics and density functional calculations are provided to help understand the highly anisotropic water stripe patterns observed.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Stabilization of a complex perovskite superstructure under ambient conditions: Influence of cation composition and ordering, and evaluation as an SOFC cathode

Ba1.6Ca2.3Y1.1Fe5O13 is an Fe3+ oxide adopting a complex perovskite superstructure, which is an ordered intergrowth between the Ca2Fe2O5 and YBa2Fe3O8 structures featuring octahedral, square pyramidal, and tetrahedral B sites and three distinct A site environments. The distribution of A site cations was evaluated by combined neutron and X-ray powder diffraction. Consistent with the Fe3+ charge state, the material is an antiferromagnetic insulator with a Néel temperature of 480-485 °C and has a relatively low d.c. conductivity of 2.06 S cm-1 at 700 °C. The observed area specific resistance in symmetrical cell cathodes with the samarium-doped ceria electrolyte is 0.87 Ω cm2 at 700 °C, consistent with the square pyramidal Fe3+ layer favoring oxide ion formation and mobility in the oxygen reduction reaction. Density functional theory calculations reveal factors favoring the observed cation ordering and its influence on the electronic structure, in particular the frontier occupied and unoccupied electronic states. © 2010 American Chemical Society.

First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

Accuracy and transferability of GAP models for tungsten

We introduce interatomic potentials for tungsten in the bcc crystal phase and its defects within the Gaussian Approximation Potential (GAP) framework, fitted to a database of first principles density functional theory (DFT) calculations. We investigate the performance of a sequence of models based on databases of increasing coverage in configuration space and showcase our strategy of choosing representative small unit cells to train models that predict properties only observable using thousands of atoms. The most comprehensive model is then used to calculate properties of the screw dislocation, including its structure, the Peierls barrier and the energetics of the vacancy-dislocation interaction. All software and raw data are available at www.libatoms.org.

Quantum Confinement Effects and Electronic Properties of SnO2 Quantum Wires and Dots

On the basis of the density functional theory (DFT) within local density approximations (LDA) approach, we calculate the band gaps for different size SnO2 quantum wire (QWs) and quantum dots (QDs). A model is proposed to passivate the surface atoms of SnO2 QWs and QDs. We find that the band gap increases between QWs and bulk evolve as Delta E-g(wire) = 1.74/d(1.20) as the effective diameter d decreases, while being Delta E-g(dot) = 2.84/d(1.26) for the QDs. Though the similar to d(1.2) scale is significantly different from similar to d(2) of the effective mass result, the ratio of band gap increases between SnO2 QWs and QDs is 0.609, very close to the effective mass prediction. We also confirm, although the LDS calculations underestimate the band gap, that they give the trend of band gap shift as much as that obtained by the hybrid functional (PBE0) with a rational mixing of 25% Fock exchange and 75% of the conventional Perdew-Burke-Ernzerhof (PBE) exchange functional for the SnO2 QWs and QDs. The relative deviation of the LDA calculated band gap difference Lambda E-g compared with the corresponding PBE0 results is only within 5%. Additionally, it is found the states of valence band maximum (VBM) and conduction band minimum (CBM) of SnO2 QWs or QDs have a mostly p- and s-like envelope function symmetry, respectively, from both LDA and PBE0 calculations.

First-principles study on the structural and electronic properties of ultrathin ZnO nanofilms

First-principles calculations; ZnO nanofilms; Electronic properties; Quantum effects； NANOBELTS; NANORINGS; WURTZITE; ENERGY Abstract: Using first-principles density-functional calculations, we have studied the structural and electronic properties Of Ultrathin ZnO {0001} nanofilms. The structural parameters, the charge densities, band structures and density of states have been investigated. The results show that there are remarkable charge transfers from Zn to O atoms in the ZOO nanofilms. All the ZOO nanofilms exhibit direct wide band gaps compared with bulk counterpart, and the gap decreases with increased thickness of the nanofilms. The decreased band gap is associated with the weaker ionic bonding within layers and the less localization of electrons in thicker films. A staircase-like density of states occurs at the bottom of conduction band, indicating the two-dimensional quantum effects in ZnO nanofilms.

Mechanical and chemical bonding properties of ground state BeH2

The crystal structure, mechanical properties and electronic structure of ground state BeH2 are calculated employing the first-principles methods based on the density functional theory. Our calculated structural parameters at equilibrium volume are well consistent with experimental results. Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. The bulk modulus B, Shear modulus G, Young's modulus E, and Poisson's ratio upsilon are deduced from the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics in band structure, density of states, and charge distribution. The ionicity in the Be-H bond is mainly featured by charge transfer from Be 2s to H 1s atomic orbitals while its covalency is dominated by the hybridization of H 1s and Be 2p states. The Bader analysis of BeH2 and MgH2 are performed to describe the ionic/covalent character quantitatively and we find that about 1.61 (1.6) electrons transfer from each Be (Mg) atom to H atoms.

First-principles study of ground-state properties and high pressure behavior of ThO2

The mechanical properties, electronic structure and phonon dispersion of ground state ThO2 as well as the structure behavior up to 240 GPa are studied using first-principles density-functional theory. Our calculated elastic constants indicate that both the ground-state fluorite structure and high pressure cotunnite structure of ThO2 are mechanically stable. The bulk modulus, shear modulus, and Young's modulus of cotunnite ThO2 are all smaller by approximately 25% compared with those of fluorite ThO2. The Poisson's ratios of both structures are approximately equal to 0.3 and the hardness of fluorite ThO2 is 22.4 GPa. The electronic structure and bonding nature of fluorite ThO2 are fully analyzed, and show that the Th-O bond displays a mixed ionic/covalent character. The phase transition from the fluorite to cotunnite structure is calculated to occur at the pressure of 26.5 GPa, consistent with recent experimental measurement by ldiri et al. [1]. For the cotunnite phase it is further predicted that an isostructural transition takes place in the pressure region of 80-130 GPa.

Optical properties of UO2 and PuO2

We perform first-principles calculations of electronic structure and optical properties for UO2 and PuO2 based on the density functional theory using the generalized gradient approximation (GGA) + U scheme. The main features in orbital-resolved partial density of states for occupied f and p orbitals, unoccupied d orbitals, and related gaps are well reproduced compared to experimental observations. Based on the satisfactory ground-state electronic structure calculations, the dynamical dielectric function and related optical spectra, i.e., the reflectivity, adsorption coefficient, energy-loss, and refractive index spectrum, are obtained. These results are consistent with the available experiments.

Structural, mechanical, thermodynamic, and electronic properties of thorium hydrides from first-principles

We perform first-principles calculations of the structural, electronic, mechanical, and thermodynamic properties of thorium hydrides (ThH2 and Th4H15) based on the density functional theory with generalized gradient approximation. The equilibrium geometries, the total and partial densities of states, charge density, elastic constants, elastic moduli, Poisson's ratio, and phonon dispersion curves for these materials are systematically investigated and analyzed in comparison with experiments and previous calculations. These results show that our calculated equilibrium structural parameters are well consistent with experiments. The Th-H bonds in all thorium hydrides exhibit weak covalent character, but the ionic properties for ThH2 and Th4H15 are different due to their different hydrogen concentration. It is found that while in ThH2 about 1.5 electrons transfer from each Th atom to H, in Th4H15 the charge transfer from each Th atom is around 2.1 electrons. Our calculated phonon spectrum for the stable body-centered tetragonal phase of ThH2 accords well with experiments. In addition we show that ThH2 in the fluorite phase is mechanically and dynamically unstable.

Au@Si-n: Growth behavior, stability and electronic structure

The configurations, stability, and electronic structure of AuSin (n = 1-16) clusters have been investigated within the framework of the density functional theory at the B3PW91/LanL2DZ and PW91/DNP levels. The results show that the Au atom begins to occupy the interior site for cages as small as Si-11 and for Si-12 the Au atom completely falls into the interior site forming Au@Si-12 cage. A relatively large embedding energy and small HOMO-LUMO gap are also found for this Au@Si-12 structure indicating enhanced chemical activity and good electronic transfer properties. All these make Au@Si-12 attractive for cluster-assembled materials.