Insight on Gaussian basis set truncation errors in weak to intermediate magnetic fields with an approximate Hamiltonian

Abstract

Strong magnetic fields such as those found on white dwarfs have significant effects on the electronic structure of atoms and molecules. However, the vast majority of molecular studies in the literature in such fields are carried out with Gaussian basis sets designed for zero field, leading to large basis set truncation errors [Lehtola et al, Mol. Phys. 2020, 118, e1597989]. In this work, we aim to identify the failures of the Gaussian basis sets in atomic calculations to guide the design of new basis sets for strong magnetic fields. We achieve this by performing fully numerical electronic structure calculations at the complete basis set (CBS) limit for the ground state and low lying excited states of the atoms $1\le Z\le 18$ in weak to intermediate magnetic fields. We also carry out finite-field calculations for a variety of Gaussian basis sets, introducing a real-orbital approximation for the magnetic-field Hamiltonian. Our primary focus is on the aug-cc-pVTZ basis set, which has been used in many works in the literature. A study of the differences in total energies of the fully numerical CBS limit calculations and the approximate Gaussian basis calculations is carried out to provide insight into basis set truncation errors. Examining a variety of states over the range of magnetic field strengths from $B=0$ to $B=0.6B_0$, we observe significant differences for the aug-cc-pVTZ basis set, while much smaller errors are afforded by the benchmark-quality AHGBSP3-9 basis set [Lehtola, J. Chem. Phys. 2020, 152, 134108]. This suggests that there is considerable room to improve Gaussian basis sets for calculations at finite magnetic fields.

arXiv comment27 pages, 18 figures

References in corpus (28)

• The Chronus Quantum (ChronusQ) Software Package
• The helium atom in a strong magnetic field
• Coupled-cluster theory for atoms and molecules in strong magnetic fields
• Hydrogen molecule in a magnetic field: The lowest states of the Pi manifold and the global ground state of the parallel configuration
• An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient
• Nonzero angular momentum states of the helium atom in a strong magnetic field
• A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules
• The hydrogen molecule in magnetic fields: The ground states of the Sigma manifold of the parallel configuration
• Ground states of the atoms H, He,...., Ne and their singly positive ions in strong magnetic fields: The high field regime
• The ground state of the Lithium atom in strong magnetic fields
• Lithium in strong magnetic fields
• The ground state of the carbon atom in strong magnetic fields
• Spectrum of Neutral Helium in Strong Magnetic Fields
• Beryllium in strong magnetic fields
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• Large-scale relativistic complete active space self-consistent field with robust convergence
• Fully numerical electronic structure calculations on diatomic molecules in weak to strong magnetic fields
• Cholesky decomposition of complex two-electron integrals over GIAOs: Efficient MP2 computations for large molecules in strong magnetic fields
• Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets
• Polarized Gaussian basis sets from one-electron ions
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• Non-perturbative calculation of orbital- and spin effects in molecules subject to non-uniform magnetic fields
• Exchange and correlation energies of ground states of atoms and molecules in strong magnetic fields
• The beryllium atom and beryllium positive ion in strong magnetic fields
• A DZ white dwarf with a 30 MG magnetic field
• Accurate Anisotropic Gaussian Type Orbital Basis Sets for Atoms in Strong Magnetic Fields
• Time-Dependent Nuclear-Electronic Orbital Hartree-Fock Theory in a Strong Uniform Magnetic Field
• Finite-Difference Calculations for Atoms and Diatomic Molecules in Strong Magnetic and Static Electric Fields

Cited by in corpus (5)

• OpenOrbitalOptimizer -- a reusable open source library for self-consistent field calculations
• Systematic study of confinement induced effects on atomic electronic structure
• Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program
• Atomic Confinement Potentials and the Generation of Numerical Atomic Orbitals
• Importance profiles. Visualization of atomic basis set requirements