Minimal Basis Set Research Articles

Partial density of states representation for accurate deep neural network predictions of X-ray spectra.

The performance of a machine learning (ML) algorithm for chemistry is highly contingent upon the architect's choice of input representation. This work introduces the partial density of states (p-DOS) descriptor: a novel, quantum-inspired structural representation which encodes relevant electronic information for machine learning models seeking to simulate X-ray spectroscopy. p-DOS uses a minimal basis set in conjunction with a guess (non-optimised) electronic configuration to extract and then discretise the density of states (DOS) of the absorbing atom to form the input vector. We demonstrate that while the electronically-focused p-DOS performs well in isolation, optimal performance is achieved when supplemented with nuclear structural information imparted via a geometric representation. p-DOS provides a description of the key electronic properties of a system which is not only concise and computationally efficient, but also independent of molecular size or choice of basis set. It can be rapidly generated, facilitating its application with large training sets. Its performance is demonstrated using a wide variety of examples at the sulphur K-edge, including the prediction of ultrafast X-ray spectroscopic signal associated with photoexcited 2(5H)-thiophenone. These results highlight the potential for ML models developed using p-DOS to contribute to the interpretation and prediction of experimental results e.g. in operando measurements of batteries and/or catalysts and femtosecond time-resolved studies, especially those made possible by emergent cutting-edge technologies, especially X-ray free electron lasers.

Targeting spectroscopic accuracy for dispersion bound systems from abinitio techniques: Translational eigenstates of Ne@C70.

We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm-1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126-2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice.

Obtaining Robust Density Functional Tight-Binding Parameters for Solids across the Periodic Table.

The density functional tight-binding (DFTB) approach allows electronic structure-based simulations at length and time scales far beyond what is possible with first-principles methods. This is achieved by using minimal basis sets and empirical approximations. Unfortunately, the sparse availability of parameters across the periodic table is a significant barrier to the use of DFTB in many cases. We therefore propose a workflow that allows the robust and consistent parametrization of DFTB across the periodic table. Importantly, our approach requires no element-pairwise parameters so that the parameters can be used for all element combinations and are readily extendable. This is achieved by parametrizing all elements on a consistent set of artificial homoelemental crystals, spanning a wide range of coordination environments. The transferability of the resulting periodic table baseline parameters to multielement systems and unknown structures is explored and the model is extensively benchmarked against previous specialized and general DFTB parametrizations.

Which Options Exist for NISQ-Friendly Linear Response Formulations?

Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136-9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed "proj LRSD" and "all-proj LRSD".

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive.

We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C5H5)2], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.

Hartree-fock roothaan calculations using optimized huzinaga orbitals on small molecules

We present the ground-state solution of some small molecules using the Hartree–Fock Roothaan method with the optimized Huzinaga basis set. Unlike the previously used least-square methods, the contraction coefficients and exponents of Huzinaga-parameterized primitive Gaussian functions for minimal basis sets are energy-optimized at the atomic level for each molecule. Consequently, as an alternative to and in comparison with standard parameterization, the optimized orbitals significantly improve the total energy and the equilibrium bond length with substantial enhancement shown for heavier nuclei. Despite similar computational cost, the application of our scheme leads to much improved minimal-basis-set Hartree–Fock calculations with less required parameters to match the large basis set calculations. Furthermore, the localization of the electrons near the nuclei which is missing with the standard parameterization is observed with the current scheme.

Toward Accurate Post-Born-Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals.

Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born-Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H2 and D2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born-Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory.

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter (α) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively α can be viewed as a regularization parameter, where each value of α represents a valid variant of BW-s2, which we denote as BW-s2(α). In this work, we semi-empirically optimize α for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the abinitio parameterization of BW-s2 and MP2. We demonstrate that the optimal α parameter (α = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2(α) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2(α) is iterative O(N5), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

An atom-in-molecule adaptive polarized valence single-ζ atomic orbital basis for electronic structure calculations.

Many low-cost or semiempirical quantum mechanical-based electronic structure methods suffer from the use of unpolarized minimal atomic orbital (AO) basis sets. In this work, we overcome this limitation by a fully DFT variationally optimized, adaptive minimal basis set consistently available for the elements up to radon (Z = 86). The new key feature is to make the linear coefficients of the primitive Gaussians in a contracted AO dependent on the effective atomic charge of the atom in the molecule, i.e., each symmetry-unique atom obtains its "own" specifically adapted basis functions. In this way, the physically important "breathing" of the AOs in a molecule with (a) atomic charge (expansion/contraction for anionic/cationic states) and (b) the number of close-lying bonded neighbor atoms is accounted for. The required atomic charges are obtained from a specially developed extended Hückel type Hamiltonian and the coordination numbers from the molecule geometry. Proper analytical derivatives of the resulting adaptive basis functions can easily be derived. Moreover, the basis functions are electric field-dependent, thus improving the description of, e.g., dipole moments and polarizabilities. The new basis set termed q-vSZP (charge dependent valence single-ζ, polarized) is thoroughly benchmarked for atomic/molecular and thermochemical properties compared to standard minimal and double-ζ basis sets at the DFT level with the accurate ωB97X-D4 functional. It is shown that q-vSZP is clearly superior to existing minimal basis sets, often reaching double-ζ quality or even better results. We expect it to be the optimal choice in future semiempirical quantum mechanical methods.

Use of Gaussian-Type Functions for Describing Fast Ion-Matter Irradiation with Time-Dependent Density Functional Theory.

The electronic stopping power is an observable property that quantifies the ability of swift ions to penetrate matter to transfer energy to the electron cloud. The recent literature has proven the value of Real-Time Time-Dependent Density Functional Theory to accurately evaluate this property from first-principles, but questions remain regarding the capability of computer codes relying on atom-centered basis functions to capture the physics at play. In this Perspective, we draw attention to the fact that irradiation by swift ions triggers electron emission into the continuum, especially at the Bragg peak. We investigate the ability of Gaussian atomic orbitals (AOC), which were fitted to mimic continuum wave functions, to improve electronic stopping power predictions. AOC are added to standard correlation-consistent basis sets or STO minimal basis sets. Our benchmarks for water irradiation by fast protons clearly advocate for the use of AOC, especially near the Bragg peak. We show that AOC only need to be placed on the molecules struck by the ion. The number of AOC that are added to the usual basis set is relatively small compared to the total number of atomic orbitals, making the use of such a basis set an excellent choice from a computational cost point of view. The optimum basis set combination is applied for the calculation of the stopping power of a proton in water with encouraging agreement with experimental data.

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy.

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

Quantum-inspired tempering for ground state approximation using artificial neural networks

A large body of work has demonstrated that parameterized artificial neural networks (ANNs) can efficiently describe ground states of numerous interesting quantum many-body Hamiltonians. However, the standard variational algorithms used to update or train the ANN parameters can get trapped in local minima, especially for frustrated systems and even if the representation is sufficiently expressive. We propose a parallel tempering method that facilitates escape from such local minima. This methods involves training multiple ANNs independently, with each simulation governed by a Hamiltonian with a different "driver" strength, in analogy to quantum parallel tempering, and it incorporates an update step into the training that allows for the exchange of neighboring ANN configurations. We study instances from two classes of Hamiltonians to demonstrate the utility of our approach using Restricted Boltzmann Machines as our parameterized ANN. The first instance is based on a permutation-invariant Hamiltonian whose landscape stymies the standard training algorithm by drawing it increasingly to a false local minimum. The second instance is four hydrogen atoms arranged in a rectangle, which is an instance of the second quantized electronic structure Hamiltonian discretized using Gaussian basis functions. We study this problem in a minimal basis set, which exhibits false minima that can trap the standard variational algorithm despite the problem’s small size. We show that augmenting the training with quantum parallel tempering becomes useful to finding good approximations to the ground states of these problem instances.

Periodic Bootstrap Embedding.

Bootstrap embedding (BE) is a recently developed electronic structure method that has shown great success at treating electron correlation in molecules. Here, we extend BE to treat surfaces and solids where the wave function is represented in periodic boundary conditions using reciprocal space sums (i.e., k-point sampling). The major benefit of this approach is that the resulting fragment Hamiltonians carry no explicit dependence on the reciprocal space sums, allowing one to apply traditional nonperiodic electronic structure codes to the fragments even though the entire system requires careful consideration of periodic boundary conditions. Using coupled cluster singles and doubles (CCSD) as an example method to solve the fragment Hamiltonians, we present minimal basis set CCSD-in-HF results on 1D conducting polymers. We show that periodic BE-CCSD can typically recover ∼99.9% of the electron correlation energy. We further demonstrate that periodic BE-CCSD is feasible even for complex donor-acceptor polymers of interest to organic solar cells─despite the fact that the monomers are sufficiently large that even a Γ-point periodic CCSD calculation is prohibitive. We conclude that BE is a promising new tool for applying molecular electronic structure tools to solids and interfaces.

Adaptive basis sets for practical quantum computing

AbstractElectronic structure calculations on small systems such as H2, H2O, LiH, and BeH2 with chemical accuracy are still a challenge for the current generation of noisy intermediate‐scale quantum (NISQ) devices. One of the reasons is that due to the device limitations, only minimal basis sets are commonly applied in quantum chemical calculations, which allows one to keep the number of qubits employed in the calculations at a minimum. However, the use of minimal basis sets leads to very large errors in the computed molecular energies as well as potential energy surface shapes. One way to increase the accuracy of electronic structure calculations is through the development of small basis sets better suited for quantum computing. In this work, we show that the use of adaptive basis sets, in which exponents and contraction coefficients depend on molecular structure, provides an easy way to dramatically improve the accuracy of quantum chemical calculations without the need to increase the basis set size and thus the number of qubits utilized in quantum circuits. As a proof of principle, we optimize an adaptive minimal basis set for quantum computing calculations on an H2 molecule, in which exponents and contraction coefficients depend on the HH distance, and apply it to the generation of H2 potential energy surface on IBM‐Q quantum devices. The adaptive minimal basis set reaches the accuracy of the double‐zeta basis sets, thus allowing one to perform double‐zeta quality calculations on quantum devices without the need to utilize twice as many qubits in simulations. This approach can be extended to other molecular systems and larger basis sets in a straightforward manner.

Generalized complete orthonormal sets of Exponential-Type functions

New generalized complete orthonormal sets of ψα∗,μ-exponential type functions are described and applied to LCAO calculations for the total energies of atomic systems within the framework of the minimal basis sets approximation in combined Hartree-Fock-Roothaan method. A new variational parameter μ is proposed in the complete orthonormal sets of ψα∗,μ-exponential type functions where α∗ parameter changes in the range -∞<α∗<3. The combined Hartree-Fock-Roothaan calculations have been performed to determine the performance of new basis sets. The comparison has been made by the use of different generalized exponential type basis sets used in the literature. The optimal values of parameters are determined by minimizing the total energies. According to obtained total energy results, the ψα∗,μ basis sets are more effective than standard definition of ψα∗ basis sets and there is no restriction on α∗ values studied systems in this work.

Local Potentials Reconstructed within Linearly Independent Product Basis Sets of Increasing Size.

Given a matrix representation of a local potential v(r) within a one-electron basis set of functions that form linearly independent products (LIP), it is possible to construct a well-defined local potential that is equivalent to v(r) within that basis set and has the form of an expansion in basis function products. Recently, we showed that for exchange-correlation potentials vXC(r) defined on the infinite-dimensional Hilbert space, the potentials reconstructed from matrices of vXC(r) within minimal LIP basis sets of occupied Kohn-Sham orbitals bear only qualitative resemblance to the originals. Here, we show that if the LIP basis set is enlarged by including low-lying virtual Kohn-Sham orbitals, the agreement between and vXC(r) improves to the extent that the basis function products are appropriate as a basis for vXC(r). These findings validate the LIP technology as a rigorous potential reconstruction method.

Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules.

We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one s-type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra. Benchmarked on organic molecules and compared to standard TDDFT, TDDFT-ris has an average energy error of only 0.06 eV and yields absorption spectra in close agreement with TDDFT. Thus, TDDFT-ris enables simulation of realistic absorption spectra in large molecules that would be inaccessible from standard TDDFT.

Aspherical atom refinements on X-ray data of diverse structures including disordered and covalent organic framework systems: a time-accuracy trade-off.

Aspherical atom refinement is the key to achieving accurate structure models, displacement parameters, hydrogen-bond lengths and analysis of weak interactions, amongst other examples. There are various quantum crystallographic methods to perform aspherical atom refinement, including Hirshfeld atom refinement (HAR) and transferable aspherical atom model (TAAM) refinement. Both HAR and TAAM have their limitations and advantages, the former being more accurate and the latter being faster. With the advent of non-spherical atoms in Olex2 (NoSpherA2), it is now possible to overcome some limitations, like treating disorder, twinning and network structures, in aspherical refinements using HAR, TAAM or both together. TAAM refinement in NoSpherA2 showed significant improvement in refinement statistics compared with independent atom model (IAM) refinements on a diverse set of X-ray diffraction data. The sensitivity of TAAM towards poor data quality and disorder was observed in terms of higher refinement statistics for such structures. A comparison of IAM with TAAM and HAR in NoSpherA2 indicated that the time taken by TAAM refinements was of the same order of magnitude as that taken by IAM, while in HAR the time taken using a minimal basis set was 50 times higher than for IAM and rapidly increased with increasing size of the basis sets used. The displacement parameters for hydrogen and non-hydrogen atoms were very similar in both HAR and TAAM refinements. The hydrogen-bond lengths were slightly closer to neutron reference values in the case of HAR with higher basis sets than in TAAM. To benefit from the advantages of each method, a new hybrid refinement approach has been introduced, allowing a combination of IAM, HAR and TAAM in one structure refinement. Refinement of coordination complexes involving metal-organic compounds and network structures such as covalent organic frameworks and metal-organic frameworks is now possible in a hybrid mode such as IAM-TAAM or HAR-TAAM, where the metal atoms are treated via either the IAM or HAR method and the organic part via TAAM, thus reducing the computational costs without compromising the accuracy. Formal charges on the metal and ligand can also be introduced in hybrid-mode refinement.

Improving the Accuracy of Variational Quantum Eigensolvers with Fewer Qubits Using Orbital Optimization

Near-term quantum computers will be limited in the number of qubits on which they can process information as well as the depth of the circuits that they can coherently carry out. To date, experimental demonstrations of algorithms such as the Variational Quantum Eigensolver (VQE) have been limited to small molecules using minimal basis sets for this reason. In this work we propose incorporating an orbital optimization scheme into quantum eigensolvers wherein a parametrized partial unitary transformation is applied to the basis functions set in order to reduce the number of qubits required for a given problem. The optimal transformation is found by minimizing the ground state energy with respect to this partial unitary matrix. Through numerical simulations of small molecules up to 16 spin orbitals, we demonstrate that this method has the ability to greatly extend the capabilities of near-term quantum computers with regard to the electronic structure problem. We find that VQE paired with orbital optimization consistently achieves lower ground state energies than traditional VQE when using the same number of qubits and even frequently achieves lower ground state energies than VQE methods using more qubits.

A simple quantum model for solving two-electron diatomics approximately

Analytical expressions are obtained for the total energy of two-electron heteronuclear diatomics using a minimal basis set consisting of a 1S-Gaussian orbital centered on each atom and a Heitler-London trial wavefunction. The equations are applied to the recently discovered helium hydride cation (2019, Nature 568, 357–9) to understand its stability. Special attention was given to the fulfillment of the virial theorem and the different limiting cases (homonuclear diatomics, short and long internuclear distances). The correctness of equations and numerical results was confirmed via a brute-force Metropolis Monte Carlo integration. Expressions for the moments of the electron density (including the dipole moments) were obtained. The effect of the spin state on the dipole moment is discussed and the model agrees qualitatively with results from more advanced ab initio calculations. I also applied the model to study the metastable helium dimer dication () and the predicted bond length (0.7359 Å) is in good agreement with the experimental value (0.75 ± 0.02 Å). The advantages and limitations of this simple quantum model are discussed.