Robust Pipek-Mezey Orbital Localization in Periodic Solids

Abstract

We describe a robust method for determining Pipek-Mezey (PM) Wannier functions (WF), recently introduced by Jónsson et al. (J. Chem. Theor. Chem. 2017, 13, 460), which provide some formal advantages over the more common Boys (also known as maximally-localized) Wannier functions. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) based PMWF solver is demonstrated to yield dramatically faster convergence compared to the alternatives (steepest ascent and conjugate gradient) in a variety of 1-, 2-, and 3-dimensional solids (including some with vanishing gaps), and can be used to obtain Wannier functions robustly in supercells with thousands of atoms. Evaluation of the PM functional and its gradient in periodic LCAO representation used a particularly simple definition of atomic charges obtained by Moore-Penrose pseudoinverse projection onto the minimal atomic orbital basis. An automated "Canonicalize Phase then Randomize" (CPR) method for generating the initial guess for WFs contributes significantly to the robustness of the solver.

arXiv comment24 pages, 1 figure

References in corpus (7)

• Maximally-localized generalized Wannier functions for composite energy bands
• Maximally localized Wannier functions: Theory and applications
• Wannier90 as a community code: new features and applications
• Intrinsic atomic orbitals: An unbiased bridge between quantum theory and chemical concepts
• Automated construction of maximally localized Wannier functions: Optimized projection functions method
• Theory and applications of generalized Pipek--Mezey Wannier functions
• Efficient Evaluation of Exact Exchange for Periodic Systems via Concentric Atomic Density Fitting