Linear-scaling density functional theory using the projector augmented wave method

Quantum mechanical simulation of realistic models of nanostructured systems, such as nanocrystals and crystalline interfaces, demands computational methods combining high-accuracy with low-order scaling with system size. Blöchl's projector augmented wave (PAW) approach enables all-electron (AE) calculations with the efficiency and systematic accuracy of plane-wave pseudopotential calculations. Meanwhile, linear-scaling (LS) approaches to density functional theory (DFT) allow for simulation of thousands of atoms in feasible computational effort. This article describes an adaptation of PAW for use in the LS-DFT framework provided by the ONETEP LS-DFT package. ONETEP uses optimisation of the density matrix through in situ-optimised local orbitals rather than the direct calculation of eigenstates as in traditional PAW approaches. The method is shown to be comparably accurate to both PAW and AE approaches and to exhibit improved convergence properties compared to norm-conserving pseudopotential methods

Van der Waals interactions in DFT made easy by Wannier functions

Ubiquitous Van der Waals interactions between atoms and molecules are important for many molecular and solid structures. These systems are often studied from first principles using the Density Functional Theory (DFT). However, the commonly used DFT functionals fail to capture the essence of Van der Waals effects. Many attempts to correct for this problem have been proposed, which are not completely satisfactory because they are either very complex and computationally expensive or have a basic semiempirical character. We here describe a novel approach, based on the use of the Maximally-Localized Wannier functions, that appears to be promising, being simple, efficient, accurate, and transferable (charge polarization effects are naturally included). The results of test applications are presented.Comment: submitted to Phys. Rev. Let

Anisotropic charge screening and supercell size convergence of defect formation energies

One of the main sources of error associated with the calculation of defect formation energies using plane-wave density functional theory (DFT) is finite size error resulting from the use of relatively small simulation cells and periodic boundary conditions. Most widely used methods for correcting this error, such as that of Makov and Payne, assume that the dielectric response of the material is isotropic and can be described using a scalar dielectric constant . However, this is strictly only valid for cubic crystals, and cannot work in highly anisotropic cases. Here we introduce a variation of the technique of extrapolation based on the Madelung potential that allows the calculation of well-converged dilute limit defect formation energies in noncubic systems with highly anisotropic dielectric properties. As an example of the implementation of this technique we study a selection of defects in the ceramic oxide Li2TiO3 which is currently being considered as a lithium battery material and a breeder material for fusion reactors

Factors influencing the distribution of charge in polar nanocrystals

We perform first-principles calculations of wurtzite GaAs nanorods to explore the factors determining charge distributions in polar nanostructures. We show that both the direction and magnitude of the dipole moment $\mathbf{d}$ of a nanorod, and its electic field, depend sensitively on how its surfaces are terminated and do not depend strongly on the spontaneous polarization of the underlying lattice. We identify two physical mechanisms by which $\mathbf{d}$ is controlled by the surface termination, and we show that the excess charge on the nanorod ends is not strongly localized. We discuss the implications of these results for tuning nanocrystal properties, and for their growth and assembly.Comment: Accepted for publication in Phys. Rev. B Rapid Communication

Energy landscape and band-structure tuning in realistic MoS2/MoSe2 heterostructures

While monolayer forms of two-dimensional materials are well characterized both experimentally and theoretically, properties of bilayer heterostructures are not nearly so well known. We employ high-accuracy linear-scaling density functional theory calculations utilizing nonlocal van der Waals functionals to explore the possible constructions of the MoS2/MoSe2 interface. Utilizing large supercells, we vary rotation, translation, and separation of the layers without introducing unrealistic strain. The energy landscape shows very low variations under rotation, with no strongly preferred alignments. By unfolding the spectral function into the primitive cells, we show that the monolayers are more independent than in homo-bilayers and that the electronic band structure of each layer is tunable through rotation, thus influencing hole effective masses.The authors acknowledge the support of the Winton Programme for the Physics of Sustainability. Computing resources were provided by the Darwin Supercomputer of the University of Cambridge High Performance Computing Service and the Argonne Leadership Computing Facility at Argonne National Laboratory, supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-06CH11357. G.C.C. acknowledges the support of the Cambridge Home and EU Scholarship Scheme.This is the accepted manuscript of a paper published in Physical Review B (Constantinescu GC, Hine NDM, Physical Review B, 2015, 91, 195416, doi:10.1103/PhysRevB.91.195416). The final version is available at http://dx.doi.org/10.1103/PhysRevB.91.19541

Tracing potential energy surfaces of electronic excitations via their transition origins: application to Oxirane

We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a widerange of molecular geometries. This is achieved by locating the switching of transition origins of adiabatic potential surfaces as the geometry changes. The transition vectors for analysing transition origins are provided by linear response time-dependent density functional theory (TDDFT) calculations under the Tamm-Dancoff approximation. We study the photochemical CO ring opening of oxirane as an example and show that the results corroborate the traditional Gomer-Noyes mechanism derived experimentally. The knowledge of specific states for the reaction also agrees well with that given by previous theoretical work using TDDFT surface-hopping dynamics that was validated by high-quality quantum Monte Carlo calculations. We also show that QNTO can be useful for considerably larger and more complex systems: by projecting the excitations to those of a reference oxirane molecule, the approach is able to identify and analyse specific excitations of a trans-2,3-diphenyloxirane molecule.Comment: 14 pages, 12 figure

Projector self-consistent DFT+U using non-orthogonal generalized Wannier functions

We present a formulation of the density-functional theory + Hubbard model (DFT+U) method that is self-consistent over the choice of Hubbard projectors used to define the correlated subspaces. In order to overcome the arbitrariness in this choice, we propose the use of non-orthogonal generalized Wannier functions (NGWFs) as projectors for the DFT+U correction. We iteratively refine these NGWF projectors and, hence, the DFT+U functional, such that the correlated subspaces are fully self-consistent with the DFT+U ground-state. We discuss the convergence characteristics of this algorithm and compare ground-state properties thus computed with those calculated using hydrogenic projectors. Our approach is implemented within, but not restricted to, a linear-scaling DFT framework, opening the path to DFT+U calculations on systems of unprecedented size.Comment: 4 pages, 3 figures. This version (v2) matches that accepted for Physical Review B Rapid Communications on 26th July 201

Accurate ionic forces and geometry optimization in linear-scaling density-functional theory with local orbitals

Linear scaling methods for density-functional theory (DFT) simulations are formulated in terms of localized orbitals in real space, rather than the delocalized eigenstates of conventional approaches. In local-orbital methods, relative to conventional DFT, desirable properties can be lost to some extent, such as the translational invariance of the total energy of a system with respect to small displacements and the smoothness of the potential-energy surface. This has repercussions for calculating accurate ionic forces and geometries. In this work we present results from onetep, our linear scaling method based on localized orbitals in real space. The use of psinc functions for the underlying basis set and on-the-fly optimization of the localized orbitals results in smooth potential-energy surfaces that are consistent with ionic forces calculated using the Hellmann-Feynman theorem. This enables accurate geometry optimization to be performed. Results for surface reconstructions in silicon are presented, along with three example systems demonstrating the performance of a quasi-Newton geometry optimization algorithm: an organic zwitterion, a point defect in an ionic crystal, and a semiconductor nanostructure.<br/