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 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 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.

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]

First-principles study of defects in CuGaO2

Using the first-principles methods, we study the electronic structure, intrinsic and extrinsic defects doping in transparent conducting oxides CuGaO2. Intrinsic defects, acceptor-type and donor-type extrinsic defects in their relevant charge state are considered. The calculation result show that copper vacancy and oxygen interstitial are the relevant defects in CuGaO2. In addition, copper vacancy is the most efficient acceptor. Substituting Be for Ga is the prominent acceptor, and substituting Ca for Cu is the prominent donors in CuGaO2. Our calculation results are expected to be a guide for preparing n-type and p-type materials in CuGaO2.

Origin of the doping bottleneck in semiconductor quantum dots: A first-principles study

Doping difficulty in semiconductor nanocrystals has been observed and its origin is currently under debate. It is not clear whether this phenomenon is energetic or depends on the growth kinetics. Using first-principles method, we show that the transition energies and defect formation energies of the donor and acceptor defects always increase as the quantum dot sizes decrease. However, for isovalent impurities, the changes of the defect formation energies are rather small. The origin of the calculated trends is explained using simple band-energy-level models.

First-Principles Study of Magnetic Properties of 3d Transition Metals Doped in ZnO Nanowires

The defect formation energies of transition metals (Cr, Fe, and Ni) doped in the pseudo-H passivated ZnO nanowires and bulk are systematically investigated using first-principles methods. The general chemical trends of the nanowires are similar to those of the bulk. We also show that the formation energy increases as the diameter of the nanowire decreases, indicating that the doping of magnetic ions in the ZnO nanowire becomes more difficult with decreasing diameter. We also systematically calculate the ferromagnetic properties of transition metals doped in the ZnO nanowire and bulk, and find that Cr ions of the nanowire favor ferromagnetic state, which is consistent with the experimental results. We also find that the ferromagnetic coupling state of Cr is more stable in the nanowire than in the bulk, which may lead to a higher T (c) useful for the nano-materials design of spintronics.

Controlling electronic structures by irradiation in single-walled SiC nanotubes: a first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, the displacement threshold energy and defect configurations are determined in SiC nanotubes. The simulation results reveal that a rich variety of defect structures (vacancies, Stone-Wales defects and antisite defects) are formed with threshold energies from 11 to 64 eV. The threshold energy shows an anisotropic behavior and exhibits a dramatic decrease with decreasing tube diameter. The electronic structure can be altered by the defects formed by irradiation, which suggests that the electron irradiation may be a way to use defect engineering to tailor electronic properties of SiC nanotubes.

Band-gap bowing and p-type doping of (Zn, Mg, Be)O wide-gap semiconductor alloys: a first-principles study

Using a first-principles band-structure method and a special quasirandom structure (SQS) approach, we systematically calculate the band gap bowing parameters and p-type doping properties of (Zn, Mg, Be)O related random ternary and quaternary alloys. We show that the bowing parameters for ZnBeO and MgBeO alloys are large and dependent on composition. This is due to the size difference and chemical mismatch between Be and Zn(Mg) atoms. We also demonstrate that adding a small amount of Be into MgO reduces the band gap indicating that the bowing parameter is larger than the band-gap difference. We select an ideal N atom with lower p atomic energy level as dopant to perform p-type doping of ZnBeO and ZnMgBeO alloys. For N doped in ZnBeO alloy, we show that the acceptor transition energies become shallower as the number of the nearest neighbor Be atoms increases. This is thought to be because of the reduction of p-d repulsion. The N-O acceptor transition energies are deep in the ZnMgBeO quaternary alloy lattice-matched to GaN substrate due to the lower valence band maximum. These decrease slightly as there are more nearest neighbor Mg atoms surrounding the N dopant. The important natural valence band alignment between ZnO, MgO, BeO, ZnBeO, and ZnMgBeO quaternary alloy is also investigated.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

Electronic and magnetic properties of YBa2Fe3O8 from a first-principles study

The electronic and magnetic properties of YBa2Fe3O8 have been systematically investigated within the framework of density-functional theory using the standard generalized gradient approximation (GGA) as well as the GGA plus Hubbard U(GGA + U) method. The GGA results show that the G-type antiferromagnetic (AFM) state is preferred among the considered magnetic configurations. The striking ionic character is shown for Y and Ba atoms while very strong hybridization is found between Fe 3d and O 2p orbitals

Structural stability and electronic properties of Ba2MIrO6 (M = La, Y) by first principles

We investigate the structural stability and electronic properties of ordered perovskite-type compounds Ba2MIrO6 (M = La, Y) by use of density functional theory. Cubic (Fm-3m), rhombohedral (R-3) and monoclinic (P2(1)/n) phases are considered for each compound. It was found that the most energetically stable phase for Ba2YIrO6 and Ba2LaIrO6 is P2(1)/n andR-3, respectively. It is also interesting to find that Ba2YIrO6 in R-3 phase, which was not reported in experiment, has a slightly lower energy than experimentally observed cubic Fm-3m phase.

Elastic anisotropy of OsB2 and RuB2 from first-principles study

The elastic anisotropy of the potential low compressible and hard materials OsB2 and RuB2 were studied by first-principles investigation within density functional theory. The structure, elastic constants, bulk modulus, shear modulus, Poisson's ratio and Debye temperature have been calculated within both local density approximation (LDA) and generalized gradient approximation (GGA). The results indicated that the calculated bulk modulus and shear modulus were in good agreement with the experimental and previous theoretical studies.

Structural stability and magnetic coupling in CaCu3Co4O12 from first principles

The structural, electronic and magnetic properties of CaCu3Co4O12 were studied by use of the full-potential linearized augmented plane wave method. The calculated results indicate that CaCu3Co4O12 is stable both thermodynamically and mechanically. Both GGA (generalized gradient approximation) and GGA + U methods predict that CaCu3Co4O12 is metallic. The ferromagnetic configuration is only slightly more stable in energy compared with the non-magnetic configuration (3.7 meV), suggesting that they are competitive for being the ground state. Co is in the low spin state (S = 1/2).