Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators

DOI10.1137/17M1140236zbMATH Open1409.65118arXiv1707.07346OpenAlexW3102895150WikidataQ128512597 ScholiaQ128512597MaRDI QIDQ4627445FDOQ4627445

Lin Lin, Yingzhou Li

Publication date: 11 March 2019

Published in: Multiscale Modeling & Simulation (Search for Journal in Brave)

Abstract: Spectral projectors of second order differential operators play an important role in quantum physics and other scientific and engineering applications. In order to resolve local features and to obtain converged results, typically the number of degrees of freedom needed is much larger than the rank of the spectral projector. This leads to significant cost in terms of both computation and storage. In this paper, we develop a method to construct a basis set that is adaptive to the given differential operator. The basis set is systematically improvable, and the local features of the projector is built into the basis set. As a result the required number of degrees of freedom is only a small constant times the rank of the projector. The construction of the basis set uses a randomized procedure, and only requires applying the differential operator to a small number of vectors on the global domain, while each basis function itself is supported on strictly local domains and is discontinuous across the global domain. The spectral projector on the global domain is systematically approximated from such a basis set using the discontinuous Galerkin (DG) method. The global construction procedure is very flexible, and allows a local basis set to be consistently constructed even if the operator contains a nonlocal potential term. We verify the effectiveness of the globally constructed adaptive local basis set using one-, two- and three-dimensional linear problems with local potentials, as well as a one dimensional nonlinear problem with nonlocal potentials resembling the Hartree-Fock problem in quantum physics.

Full work available at URL: https://arxiv.org/abs/1707.07346

zbMATH Keywords

random samplingdifferential operatordiscontinuous Galerkinquantum physicsspectral projectoradaptive local basis

Mathematics Subject Classification ID

Applications to the sciences (65Z05) Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems (65M99)

Cites Work

• Electronic structure calculations of solids using the WIEN2k package for material sciences
• Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method
• Chebyshev-filtered subspace iteration method free of sparse diagonalization for solving the Kohn-Sham equation
• Title not available (Why is that?)
• Finding structure with randomness: Probabilistic algorithms for constructing approximate matrix decompositions
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• Iterative Procedures for Nonlinear Integral Equations
• An Interior Penalty Finite Element Method with Discontinuous Elements
• Title not available (Why is that?)
• THE PARTITION OF UNITY METHOD
• Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
• Barycentric Lagrange Interpolation
• Solution of Sparse Indefinite Systems of Linear Equations
• Spectral Analysis of the Discrete Helmholtz Operator Preconditioned with a Shifted Laplacian
• Electronic Structure
• Ab initio molecular simulations with numeric atom-centered orbitals
• Self-consistent-field calculations using Chebyshev-filtered subspace iteration
• Fast algorithm for extracting the diagonal of the inverse matrix with application to the electronic structure analysis of metallic systems
• Classical and enriched finite element formulations for Bloch-periodic boundary conditions
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. I: Total energy calculation
• Some Recent Advances in Density Matrix Theory
• Optimized local basis set for Kohn-Sham density functional theory
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework. II: force, vibration, and molecular dynamics calculations
• Zolotarev Quadrature Rules and Load Balancing for the FEAST Eigensolver
• Elliptic preconditioner for accelerating the self-consistent field iteration in Kohn-Sham density functional theory
• A posteriorierror estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems
• A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. I: Second order linear PDE.

Cited In (1)

• Numerical methods for Kohn–Sham density functional theory

Uses Software

• ONETEP
• FEAST
• WIEN2k
• SIESTA
• Quantum Espresso
• CheFSI
• lobpcg.m
• DGDFT

This page was built for publication: Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4627445)