Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. II: force, vibration, and molecular dynamics calculations
From MaRDI portal
Publication:1685442
DOI10.1016/J.JCP.2016.12.052zbMATH Open1375.82015arXiv1510.06489OpenAlexW2962946439MaRDI QIDQ1685442FDOQ1685442
Authors: Gaigong Zhang, Lin Lin, Wei Hu, Chao Yang, John E. Pask 
Publication date: 14 December 2017
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: Recently, we have proposed the adaptive local basis set for electronic structure calculations based on Kohn-Sham density functional theory in a pseudopotential framework. The adaptive local basis set is efficient and systematically improvable for total energy calculations. In this paper, we present the calculation of atomic forces, which can be used for a range of applications such as geometry optimization and molecular dynamics simulation. We demonstrate that, under mild assumptions, the computation of atomic forces can scale nearly linearly with the number of atoms in the system using the adaptive local basis set. We quantify the accuracy of the Hellmann-Feynman forces for a range of physical systems, benchmarked against converged planewave calculations, and find that the adaptive local basis set is efficient for both force and energy calculations, requiring at most a few tens of basis functions per atom to attain accuracy required in practice. Since the adaptive local basis set has implicit dependence on atomic positions, Pulay forces are in general nonzero. However, we find that the Pulay force is numerically small and systematically decreasing with increasing basis completeness, so that the Hellmann-Feynman force is sufficient for basis sizes of a few tens of basis functions per atom. We verify the accuracy of the computed forces in static calculations of quasi-1D and 3D disordered Si systems, vibration calculation of a quasi-1D Si system, and molecular dynamics calculations of H and liquid Al-Si alloy systems, where we find excellent agreement with independent benchmark results in literature.
Full work available at URL: https://arxiv.org/abs/1510.06489
Recommendations
- Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. I: Total energy calculation
- Optimized local basis set for Kohn-Sham density functional theory
- On the adaptive finite element analysis of the Kohn-Sham equations: methods, algorithms, and implementation
- A posteriori error estimator for adaptive local basis functions to solve Kohn-Sham density functional theory
- Pulay forces in density functional theory for periodic and molecular systems
discontinuous Galerkinmolecular dynamicselectronic structureKohn-Sham density functional theoryPulay forceadaptive local basis setHellmann-Feynman force
Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Quantum equilibrium statistical mechanics (general) (82B10)
Cites Work
- An Interior Penalty Finite Element Method with Discontinuous Elements
- Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
- Nonconforming Elements in the Finite Element Method with Penalty
- Forces in Molecules
- Electronic Structure
- Ab initio molecular simulations with numeric atom-centered orbitals
- Adaptive Finite Element Approximations for Kohn--Sham Models
- Fast algorithm for extracting the diagonal of the inverse matrix with application to the electronic structure analysis of metallic systems
- Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. I: Total energy calculation
- Numerical solution of the Kohn-Sham equation by finite element methods with an adaptive mesh redistribution technique
- Optimized local basis set for Kohn-Sham density functional theory
- A posteriori error estimator for adaptive local basis functions to solve Kohn-Sham density functional theory
- Linear scaling discontinuous Galerkin density matrix minimization method with local orbital enriched finite element basis: 1-D lattice model system
Cited In (17)
- Truncated spherical-wave basis set for first-principles pseudopotential calculations
- A discontinuous Galerkin scheme for full-potential electronic structure calculations
- Numerical methods for Kohn–Sham density functional theory
- A Convergence Analysis of a Structure-Preserving Gradient Flow Method for the All-Electron KohnSham Model
- Optimized local basis set for Kohn-Sham density functional theory
- NURBS-based non-periodic finite element framework for Kohn-Sham density functional theory calculations
- A novel proof of the DFT formula for the interatomic force field of molecular dynamics
- Gaussian basis implementation of the charge patching method
- All-electron formalism for total energy strain derivatives and stress tensor components for numeric atom-centered orbitals
- A posteriori error estimator for adaptive local basis functions to solve Kohn-Sham density functional theory
- Globally constructed adaptive local basis set for spectral projectors of second order differential operators
- Variationally consistent Hellmann-Feynman forces in the finite element formulation of Kohn-Sham density functional theory
- Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. I: Total energy calculation
- Accurate stress calculations based on numerical atomic orbital bases: implementation and benchmarks
- An asymptotics-based adaptive finite element method for Kohn-Sham equation
- SQDFT: spectral quadrature method for large-scale parallel \(\mathcal{O}(N)\) Kohn-Sham calculations at high temperature
- Orbital-enriched flat-top partition of unity method for the Schrödinger eigenproblem
Uses Software
This page was built for publication: Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. II: force, vibration, and molecular dynamics calculations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1685442)
