First-principles prediction of doped graphane as a high-temperature electron-phonon superconductor.

We predict by first-principles calculations that p-doped graphane is an electron-phonon superconductor with a critical temperature above the boiling point of liquid nitrogen. The unique strength of the chemical bonds between carbon atoms and the large density of electronic states at the Fermi energy arising from the reduced dimensionality give rise to a giant Kohn anomaly in the optical phonon dispersions and push the superconducting critical temperature above 90 K. As evidence of graphane was recently reported, and doping of related materials such as graphene, diamond, and carbon nanostructures is well established, superconducting graphane may be feasible.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

First-principles energetics of water: 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 solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

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.

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.

First-principles prediction of the magnetism of 3d transition-metal-doped Rocksalt MgO

Using first-principles band structure methods, we have systematically studied the electronic structures, magnetic stabilities, and half-metal properties of 3d transition-metal (TM) doped Rocksalt MgO compounds TMMg3O4 (TM = V, Cr, Mn, Fe, Co, and Ni). The calculations reveal that only CrMg3O4 has a ferromagnetic stability among the six compounds, which is explained by double-exchange mechanism. The magnetic stability is affected by the doping concentration of TM if the top valance band is composed of partially occupied t(2g) states. In addition, CrMg3O4 is a half-metallic ferromagnet. The origins of half-metallic and ferromagnetic properties are explored. The Curie temperature (T-c) of CrMg3O4 is 182 K. And it is hard for CrMg3O4 to deform due to the large bulk modulus and shear modulus, so it is a promising spintronic material. Our calculations provide the first available information on the magnetic properties of 3d TM-doped MgO.

First-principles LDA plus U and GGA plus U study of neptunium dioxide

We have performed a systematic first-principles investigation to calculate the electronic structures, mechanical properties, and phonon-dispersion curves of NpO2. The local-density approximation+U and the generalized gradient approximation+U formalisms have been used to account for the strong on-site Coulomb repulsion among the localized Np 5f electrons. By choosing the Hubbard U parameter around 4 eV, the orbital occupancy characters of Np 5f and O 2p are in good agreement with recent experiments [A. Seibert, T. Gouder, and F. Huber, J. Nucl. Mater. 389, 470 (2009)]. Comparing to our previous study of ThO2, we note that stronger covalency exists in NpO2 due to the more localization behavior of 5f electrons of Np in line with the localization-delocalization trend exhibited by the actinides series.

First-principles study of UC2 and U2C3

The electronic structure and mechanical properties Of UC2 and U2C3 have been systematically investigated using first-principles calculations by the projector-augmented-wave (PAW) method. Furthermore, in order to describe precisely the strong on-site Coulomb repulsion among the localized U 5f electrons, we adopt the generalized gradient approximation +U formalisms for the exchange-correlation term. We show that our calculated structural parameters and electronic properties for UC2 and U2C3 are in good agreement with the experimental data by choosing an appropriate Hubbard U = 3 eV. As for the chemical bonding nature, the contour plot of charge density and total density of states suggest that UC2 and U2C3 are metallic mainly contributed by the 5f electrons, mixed with significant covalent component resulted from the strong C-C bonds. The present results also illustrate that the metal-carbon (U-C) bonding and the carbon-carbon covalent bonding in U2C3 are somewhat weaker than those in UC2, leading to the weaker thermodynamic stability at high temperature as observed by experiments.

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.

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.

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.

Defects in gallium nitride nanowires: First principles calculations

Atomic configurations and formation energies of native defects in an unsaturated GaN nanowire grown along the [001] direction and with (100) lateral facets are studied using large-scale ab initio calculation. Cation and anion vacancies, antisites, and interstitials in the neutral charge state are all considered. The configurations of these defects in the core region and outermost surface region of the nanowire are different. The atomic configurations of the defects in the core region are same as those in the bulk GaN, and the formation energy is large. The defects at the surface show different atomic configurations with low formation energy. Starting from a Ga vacancy at the edge of the side plane of the nanowire, a N-N split interstitial is formed after relaxation. As a N site is replaced by a Ga atom in the suboutermost layer, the Ga atom will be expelled out of the outermost layers and leaves a vacancy at the original N site. The Ga interstitial at the outmost surface will diffuse out by interstitialcy mechanism. For all the tested cases N-N split interstitials are easily formed with low formation energy in the nanowires, indicating N-2 molecular will appear in the GaN nanowire, which agrees well with experimental findings.

First-principles study of p-type transparent conductive oxides CuXO2 (X=Y, Sc, and Al)

Using first-principles methods we have calculated electronic structures, optical properties, and hole conductivities of CuXO2 (X=Y, Sc, and Al). We show that the direct optical band gaps of CuYO2 and CuScO2 are approximately equal to their fundamental band gaps and the conduction bands of them are localized. The direct optical band gaps of CuXO2 (X=Y, Sc, and Al) are 3.3, 3.6, and 3.2 eV, respectively, which are consistent with experimental values of 3.5, 3.7, and 3.5 eV. We find that the hole mobility along long lattice c is higher than that along other directions through calculating effective masses of the three oxides. By analyzing band offset we find that CuScO2 has the highest valence band maximum (VBM) among CuXO2 (X=Y, Sc, and Al). In addition, the approximate transitivity of band offset suggests that CuScO2 has a higher VBM than CuGaO2 and CuInO2 [Phys. Rev. Lett. 88, 066405 (2002)]. We conclude that CuScO2 has a higher p-type doping ability in terms of the doping limit rule. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2991157]