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

Abstract

Linear scaling methods for density-functional theory (DFT) simulations are formulated in terms of localised orbitals in real-space, rather than the delocalised 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 \textsc{onetep}, our linear scaling method based on localised orbitals in real-space. The use of psinc functions for the underlying basis set and on-the-fly optimisation of the localised orbitals results in smooth potential energy surfaces that are consistent with ionic forces calculated using the Hellmann-Feynman theorem. This enables accurate geometry optimisation to be performed. Results for surface reconstructions in silicon are presented, along with three example systems demonstrating the performance of a quasi-Newton geometry optimisation algorithm: an organic zwitterion, a point defect in an ionic crystal, and a semiconductor nanostructure.

11 pages, 9 figures

References in corpus (14)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Selfconsistent order-N density-functional calculations for very large systems
• Nearsightedness of Electronic Matter
• A Self-Consistent First-Principles Technique Having Linear Scaling
• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Exponential decay properties of Wannier functions and related quantities
• Linear Algebraic Calculation of Green's function for Large-Scale Electronic Structure Theory
• Using ONETEP for accurate and efficient O(N) density functional calculations
• Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids
• Model Hessian for accelerating first-principles structure optimizations
• Linear scaling relaxation of the atomic positions in nano structures
• Fast method for force computations in electronic structure calculations
• Geometry Optimization of Crystals by the Quasi-Independent Curvilinear Coordinate Approximation

Cited by in corpus (20)

• O(N) methods in electronic structure calculations
• Challenges in Large Scale Quantum Mechanical Calculations
• The role of spin in the calculation of Hubbard $U$ and Hund's $J$ parameters from first principles
• Electrostatic Interactions in Finite Systems treated with Periodic Boundary Conditions: Application to Linear-Scaling Density Functional Theory
• Linear-scaling time-dependent density-functional theory in the linear response formalism
• Calculating optical absorption spectra for large systems using linear-scaling density-functional theory
• Linear-scaling time-dependent density-functional theory (TDDFT) beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals
• Linear-scaling DFT+U with full local orbital optimization
• Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U
• Factors influencing the distribution of charge in polar nanocrystals
• Fermi level pinning can determine polarity in semiconductor nanorods
• TDDFT+$U$: Hubbard corrected approximate density-functional theory in the excited-state regime
• Low-cost orbital-based linear-scaling \emph{ab initio} molecular dynamics for weakly-interacting systems
• Simulations of Nanocrystals Under Pressure: Combining Electronic Enthalpy and Linear-Scaling Density-Functional Theory
• Bromophenyl functionalization of carbon nanotubes : an ab initio study
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Wannier-function-based constrained DFT with nonorthogonality-correcting Pulay forces in application to the reorganization effects in graphene-adsorbed pentacene
• Band parameters and hybridization in 2D semiconductor heterostructures from photoemission spectroscopy
• Single indium atoms and few-atom indium clusters anchored onto graphene via silicon heteroatoms
• A Unified Heterogeneous Implementation of Numerical Atomic Orbitals-Based Real-Time TDDFT within the ABACUS Package