Towards a Linear-Scaling DFT Technique: The Density Matrix Approach

Abstract

A recently proposed linear-scaling scheme for density-functional pseudopotential calculations is described in detail. The method is based on a formulation of density functional theory in which the ground state energy is determined by minimization with respect to the density matrix, subject to the condition that the eigenvalues of the latter lie in the range [0,1]. Linear-scaling behavior is achieved by requiring that the density matrix should vanish when the separation of its arguments exceeds a chosen cutoff. The limitation on the eigenvalue range is imposed by the method of Li, Nunes and Vanderbilt. The scheme is implemented by calculating all terms in the energy on a uniform real-space grid, and minimization is performed using the conjugate-gradient method. Tests on a 512-atom Si system show that the total energy converges rapidly as the range of the density matrix is increased. A discussion of the relation between the present method and other linear-scaling methods is given, and some problems that still require solution are indicated.

REVTeX file, 27 pages with 4 uuencoded postscript figures

References in corpus (7)

• Total energy global optimizations using non orthogonal localized orbitals
• A Self-Consistent First-Principles Technique Having Linear Scaling
• Large Scale Electronic Structure Calculations with Multigrid Acceleration
• Generalization of the density-matrix method to a non-orthogonal basis
• Electronic structure calculations and molecular dynamics simulations with linear system-size scaling
• Real-space approach to calculation of electric polarization and dielectric constants
• Ab initio Molecular Dynamics in Adaptive Coordinates

Cited by in corpus (28)

• Maximally localized Wannier functions: Theory and applications
• Numerical atomic orbitals for linear scaling
• O(N) methods in electronic structure calculations
• Maximally-localized Wannier functions for disordered systems: application to amorphous silicon
• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Calculations on millions of atoms with DFT: Linear scaling shows its potential
• Expansion algorithm for the density matrix
• Large scale and linear scaling DFT with the CONQUEST code
• Basis Functions for Linear-Scaling First-Principles Calculations
• Obtaining Wannier Functions of a Crystalline Insulator within a Hartree-Fock approach: Applications to LiF and LiCl
• Parallel Sparse Matrix Multiplication for Linear Scaling Electronic Structure Calculations
• Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids
• Linear-scaling quantum Monte Carlo with non-orthogonal localized orbitals
• Next generation extended Lagrangian first principles molecular dynamics
• Optimized multi-site local orbitals in the large-scale DFT program CONQUEST
• Efficient Calculations with Multisite Local Orbitals in the large-scale DFT Code CONQUEST
• Spin unrestricted linear scaling electronic structure theory and its application to magnetic carbon doped BN nanotubes
• Density-functional embedding using a plane-wave basis
• Non-parametric Local Pseudopotentials with Machine Learning: a Tin Pseudopotential Built Using Gaussian Process Regression
• Canonical density matrix perturbation theory
• First-principles density-functional calculations using localized spherical-wave basis sets
• Fast method for force computations in electronic structure calculations
• Linear scaling calculation of band edge states and doped semiconductors
• Canonical-ensemble extended Lagrangian Born-Oppenheimer molecular dynamics for the linear scaling density functional theory
• Large-scale DFT methods for calculations of materials with complex structures
• High-accuracy large-scale DFT calculations using localized orbitals in complex electronic systems: The case of graphene-metal interfaces
• A Time-Dependent Random State Approach for Large-scale Density Functional Calculations
• Positivity Preserving Density Matrix Minimization at Finite Temperatures via Square Root