Large Basis Sets Research Articles

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications.

The increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH3) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

Automatic Orbital Pair Selection for Multilevel Local Coupled-Cluster Based on Orbital Maps.

We present an automatic, orbital-map based orbital-pair selection scheme for multilevel local coupled-cluster approaches that exploits the locality of chemical reactions by focusing on the part of the molecule directly involved in the reaction. The previously introduced pair-selected multilevel extension to domain-based local pair natural orbital coupled-cluster with singles, doubles, and semicanonical perturbative triples [DLPNO-CCSD(T0)] partitions the orbital pairs according to relative changes in pair correlation energies [Bensberg, M.; Neugebauer, J. Orbital pair selection for relative energies in the domain-based local pair natural orbital coupled-cluster method. J. Chem. Phys. 2022, 157, 064102. 10.1063/5.0100010]. To this end, maps between localized orbitals are required which in turn require maps between the atoms of structures along reaction paths. So far, these atom maps have been manually determined, which can be a (human) time-consuming procedure. Here, we present an automatic atom mapping algorithm based on the principle of minimum chemical distance that incorporates orientation dependence through dihedral angles. A similar strategy is then introduced to obtain orbital maps, which proves advantageous over the previously used direct orbital selection. Along with a modified orbital pair prescreening, this results in an improved variant of the pair-selected multilevel DLPNO-CCSD(T0) method. The performance of this approach is demonstrated for various reaction types showing a significant efficiency gain and accurate results due to beneficial, systematic error cancellation. The presented method operates in a black-box manner due to its fully automatized algorithms with only the need to specify a single target-accuracy parameter. Additionally, we demonstrate that basis set extrapolation techniques can be applied. In this context, the approach shows deficiencies for the use of large basis sets, especially with diffuse basis functions, which can be traced back to the semicanonical triples correction.

A static quantum embedding scheme based on coupled cluster theory.

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

Effect of Fluorine Substitution on Magnetizabilities: Insights from Density Functional Theory and Benchmark Coupled-Cluster Calculations.

The quantitative calculation of the magnetizability tensor of fluorine-containing molecules has been a longstanding challenge for quantum chemistry. We report a benchmark study on the effect of fluorine substitution on the magnetizability (both isotropic and anisotropic) of 24 small closed-shell molecules using coupled-cluster (CCSD(T)) theory, large basis sets, and London atomic orbitals. By extrapolation to a complete basis set (CBS) result, we establish the CCSD(T)/CBS limit and take the opportunity to assess the performance of various density functional theory approximations. Correcting for zero-point vibration further allows us to directly compare theory with (gas-phase) experimental data. We revisit Flygare's hypothesis on molecular magnetizabilities, which relates the successive replacement of hydrogen by fluorine atoms to changes in the magnetizability anisotropy. Some exceptions to this hypothesis are documented, and we rationalize these observations on the basis of hybridization on carbon.

Gas phase proton affinities of proline-containing peptides. 1: ProGly, ProAla, ProVal, ProLeu, ProIle, and ProPro

The gas-phase proton affinities (PA) for a series of proline-containing dipeptides have been measured in an ESI triple quadrupole instrument using the extended kinetic method. Proton affinities for ProGly (1), ProAla (2), ProVal (3), ProLeu (4), ProIle (5), and ProPro (6) were determined to be 969.6 ± 7.8, 990.4 ± 7.7, 987.6 ± 7.9, 982.8 ± 8.0, 988.8 ± 10.1, and 996.5 ± 12.2 kJ/mol, respectively. Predictions for the proton affinities for 1–6 were also obtained through isodesmic calculations at the B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d) level of theory. The predicted proton affinities for 1 and 6 of 966.9 and 991.0 kJ/mol are in agreement with the experimental values. However, the predicted proton affinities for 2–5 of 973.5, 975.9, 975.7, and 975.9 are between 8 and 15 kJ/mol lower than the experimental values. Additional calculations with a larger basis set (B3LYP/6-311++G(2df,2p), inclusion of dispersion (B3LYP-D3/6-311++G(d,p)), switching to second order perturbation theory (MP2/6-31++G(d,p) and MP2/6-311++G(2df,2p), or switching density functional (M06-2x/6-311++G(d,p) and M06-2x/6-311++G(2df,2p) show only modest changes in derived thermochemistry lending support to the original calculations. We recommend using the experimental proton affinities for ProGly and ProPro and using the calculated values for ProAla, ProVal, ProLeu, and ProIle with the experimental proton affinities as upper limits.

Molecular line lists for the singlet states in gaseous YN at high temperature

ABSTRACT Our aim in this study is to investigate the electronic structure of the gaseous YN molecule and provide an accurate molecular line list of spectroscopic transitions between the five lowest singlet states. The calculations were done by using large basis set for yttrium aug-cc-pVQZ-PP with relativistic effective core potentials at the spin-free level. We first used the method of complete active space self-consistent field, which was followed by the multireference singles and doubles configuration interaction. Potential energy curves and permanent and transition dipole moments were calculated at the internuclear distance range of 1.3–2.96 Å. For each state, we calculated the spectroscopic constants (Te, ωe, ωexe, ωeye, Be, αe, Re). The calculated values of the spectroscopic constants are in excellent agreement with the experimental constants available, with a relative difference of less than 5 cm−1. We further calculated the transition intensities of allowed transitions, and we provided the absorption spectra up to a temperature of 4000 K. The computed molecular line list contains 2.6 million transitions calculated between almost 18 186 vibro-rotational energy levels. It covers the frequency range between 0 and 30 000 cm−1. The effect of temperature on the spectra of hot YN was investigated with partition functions calculated between 298 and 4000 K. This study should help in identifying the singlet state spectrum of the hot YN molecule, which is possibly found in the atmospheres of cool stars and hot rocky super-Earths.

High-Resolution Spectroscopy and Structure of Heavy Carbon Subchalcogenides: Tricarbon Selenide, C3Se.

Linear tricarbon selenide, C3Se, has been studied spectroscopically for the first time using a combination of high-resolution infrared and microwave techniques. Probing laser ablation products from carbon-selenium targets in a free jet expansion with He, initial spectroscopic detection was accomplished in the infrared at a wavelength of 5 μm in search of the ν1 vibrational fundamental. Along with the band of the most abundant isotopic species C380Se found centered at about 2039 cm-1, the corresponding bands of the C382Se, C378Se, C376Se, and C377Se isotopologues were also detected. Pure rotational spectra of the five C3Se isotopologues in the 8-18 GHz frequency range were observed in a supersonic jet expansion using chirped-pulse microwave spectroscopy of the discharge products of selenophene, c-C4H4Se. Spectroscopic analyses were guided by results from high-level quantum-chemical calculations carried out at the coupled-cluster level of theory using large correlation-consistent basis sets. Using the experimentally derived ground-state rotational constants of five isotopologues and calculated zero-point vibrational corrections, an accurate semiexperimental equilibrium carbon-selenium bond length is derived.

Jump discontinuities of finite-basis-set exchange–correlation potentials at atomic nuclei

The kinetic energy density of electrons and the gradient of the electron density have pronounced jump discontinuities at the positions of the atomic nuclei in molecules. Certain exact relations then imply that molecular Kohn–Sham exchange–correlation potentials may also be discontinuous at atomic nuclei. Here, we confirm that exchange–correlation potentials derived from Hartree–Fock and correlated wavefunctions within Slater-type basis sets do exhibit such discontinuities. Despite their persistence even in large basis sets, these discontinuities are almost certainly artifacts of basis set finiteness and are expected to disappear in the basis-set limit. The findings imply that imposing electron–nucleus cusp conditions in spherically averaged form may not always be appropriate.

Molecular structure and spectra of ytterbium dihalides from first principles

The molecular structure, vibrational frequencies, infrared intensities and electric dipole moments of ytterbium dihalides YbX2 (X = F, Cl, Br, I) in their ground electronic state are studied at the coupled-cluster singles, doubles and perturbative triples CCSD(T) level using a series of large all-electron basis sets together with the complete basis set (CBS) extrapolation procedure, and with both scalar-relativistic and spin–orbit coupling effects taken into account. Excitation energies of low-lying electronic states are calculated at the multireference configuration interaction, equation-of-motion CCSD and single-reference CCSD(T) levels of theory. For all of the molecules under study, the CCSD(T)/CBS equilibrium geometries are found to be non-linear (C2v symmetry). There is a rapid decrease in heights of the barriers to linearity on passing through the YbX2 molecular series from fluoride to iodide: 1587 → 261 → 110 → 7 (all in cm−1) for YbF2 → YbCl2 → YbBr2 → YbI2, respectively. The uncertainties in the molecular structure data published to date are discussed and resolved in the light of the present results.

Compactification of determinant expansions via transcorrelation.

Although selected configuration interaction (SCI) algorithms can tackle much larger Hilbert spaces than the conventional full CI method, the scaling of their computational cost with respect to the system size remains inherently exponential. In addition, inaccuracies in describing the correlation hole at small interelectronic distances lead to the slow convergence of the electronic energy relative to the size of the one-electron basis set. To alleviate these effects, we show that the non-Hermitian, transcorrelated (TC) version of SCI significantly compactifies the determinant space, allowing us to reach a given accuracy with a much smaller number of determinants. Furthermore, we note a significant acceleration in the convergence of the TC-SCI energy as the basis set size increases. The extent of this compression and the energy convergence rate are closely linked to the accuracy of the correlation factor used for the similarity transformation of the Coulombic Hamiltonian. Our systematic investigation of small molecular systems in increasingly large basis sets illustrates the magnitude of these effects.

MP2-Based Composite Extrapolation Schemes Can Predict Core-Ionization Energies for First-Row Elements with Coupled-Cluster Level Accuracy.

X-ray photoelectron spectroscopy (XPS) measures core-electron binding energies (CEBEs) to reveal element-specific insights into the chemical environment and bonding. Accurate theoretical CEBE prediction aids XPS interpretation but requires proper modeling of orbital relaxation and electron correlation upon core-ionization. This work systematically investigates basis set selection for extrapolation to the complete basis set limit of CEBEs from ΔMP2 and ΔCC energies across 94 K-edges in diverse organic molecules. We demonstrate that an alternative composite scheme using ΔMP2 in a large basis corrected by ΔCC-ΔMP2 difference in a small basis can quantitatively recover optimally extrapolated ΔCC CEBEs within 0.02 eV. Unlike ΔCC, MP2 calculations do not suffer from convergence issues and are computationally cheaper, and thus, the composite ΔMP2/ΔCC scheme balances accuracy and cost, overcoming limitations of solely using either method. We conclude by providing a comprehensive analysis of the choice of small and large basis sets for the composite schemes and provide practical recommendations for highly accurate (within 0.10-0.15 eV MAE) ab initio prediction of XPS data.

A molecular ground electronic state with an occupied 5g spinor-The superheavy (E125)F molecule.

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac-Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s25g16f3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125+, [Og]8s1/225g7/216f5/23 with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og]8s1/225g7/216f5/227d3/218p1/21 ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. The non-relativistic ground state of the E125 atom was calculated to have a 5g5 ground state with an IE of just 3.4eV. The net effect of relativity on (E125)F is to stabilize its bonding.

Simulation of Vibrational Circular Dichroism Spectra Using Second-Order Møller-Plesset Perturbation Theory and Configuration Interaction Doubles.

We present the first single-reference calculations of the atomic axial tensors (AATs) with wave function-based methods including dynamic electron correlation effects using second-order Møller-Plesset perturbation theory (MP2) and configuration interaction doubles (CID). Our implementation involves computing the overlap of numerical derivatives of the correlated wave functions with respect to both nuclear displacement coordinates and the external magnetic field. Out test set included three small molecules, including the axially chiral hydrogen molecule dimer and (P)-hydrogen peroxide, and the achiral H2O. For our molecular test set, we observed deviations of the AATs for MP2 and CID from that of the Hartree-Fock (HF) method upward of 49%, varying with the choice of basis set. For (P)-hydrogen peroxide, electron correlation effects on the vibrational circular dichroism (VCD) rotatory strengths and corresponding spectra were particularly significant, with maximum deviations of the rotatory strengths of 62 and 49% for MP2 and CID, respectively, using our largest basis set. The inclusion of dynamic electron correlation to the computation of the AATs can have a significant impact on the resulting rotatory strengths and VCD spectra.

Evaluating Noncovalent Interactions in Halogenated Molecules with Double-Hybrid Functionals and a Dedicated Small Basis Set.

We present here an extension of our recently developed PBE-QIDH/DH-SVPD basis set to halogen atoms, with the aim of obtaining, for weakly interacting halogenated molecules, interaction energies close to those provided by a large basis set (def2-TZVPP) coupled to empirical dispersion potential. The core of our approach is the split-valence basis set, DH-SVPD, that has been developed for F, Cl, Br, and I atoms using a self-consistent formula, containing only energy terms computed for dimers and the corresponding monomers at the same level of theory. The basis set developed considering four systems, one for each halogen atoms, has been then tested on the X40, X4 × 10 benchmarks as well as on other two, less standard, data sets. Finally, a large system (380 atoms) has been also considered as a "crash" test. Our results show that the simple and nonempirical PBE-QIDH/DH-SVPD approach is able to provide accurate results for interaction energies of all the considered systems and can thus be considered as a cheaper alternative to DH functionals paired with empirical dispersion corrections and a large basis set of triple-ζ quality.

Reduced Density Matrix Formulation of Quantum Linear Response.

The prediction of spectral properties via linear response (LR) theory is an important tool in quantum chemistry for understanding photoinduced processes in molecular systems. With the advances of quantum computing, we recently adapted this method for near-term quantum hardware using a truncated active space approximation with orbital rotation, named quantum linear response (qLR). In an effort to reduce the classic cost of this hybrid approach, we here derive and implement a reduced density matrix (RDM) driven approach of qLR. This allows for the calculation of spectral properties of moderately sized molecules with much larger basis sets than so far possible. We report qLR results for benzene and R-methyloxirane with a cc-pVTZ basis set and study the effect of shot noise on the valence and oxygen K-edge absorption spectra of H2O in the cc-pVTZ basis.

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

SummaryComputational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

First Principles Simulations of Optical Rotation of Chiral Molecular Crystals.

In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT-PBC). Calculations are compared with experimental measurements and show semi-quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X-ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.

Revisiting the Most Stable Structures of the Benzene Dimer.

The benzene dimer (BD) is an archetypal model of π∙∙∙π and C-H∙∙∙π noncovalent interactions as they occur in its cofacial and perpendicular arrangements, respectively. The enthalpic stabilization of the related BD structures has been debated for a long time and is revisited here. The revisit is based on results of computations that apply the coupled-cluster theory with singles, doubles and perturbative triples [CCSD(T)] together with large basis sets and extrapolate results to the complete basis set (CBS) limit in order to accurately characterize the three most important stationary points of the intermolecular interaction energy (ΔE) surface of the BD, which correspond to the tilted T-shaped (TT), fully symmetric T-shaped (FT) and slipped-parallel (SP) structures. In the optimal geometries obtained by searching extensive sets of the CCSD(T)/CBS ΔE data of the TT, FT and SP arrangements, the resulting ΔE values were -11.84, -11.34 and -11.21 kJ/mol, respectively. The intrinsic strength of the intermolecular bonding in these configurations was evaluated by analyzing the distance dependence of the CCSD(T)/CBS ΔE data over wide ranges of intermonomer separations. In this way, regions of the relative distances that favor BD structures with either π∙∙∙π or C-H∙∙∙π interactions were found and discussed in a broader context.

Can Aromaticity Be Evaluated Using Atomic Partitions Based on the Hilbert-Space?

Aromaticity is a fundamental concept in chemistry that explains the stability and reactivity of many compounds by identifying atoms within a molecule that form an aromatic ring. Reliable aromaticity indices focus on electron delocalization and depend on atomic partitions, which give rise to the concept of an atom-in-the-molecule (AIM). Real-space atomic partitions present two important drawbacks: a high computational cost and numerical errors, limiting some aromaticity measures to medium-sized molecules with rings up to 12 atoms. This restriction hinders the study of large conjugated systems like porphyrins and nanorings. On the other hand, traditional Hilbert-space schemes are free of the latter limitations but can be unreliable for the large basis sets required in modern computational chemistry. This paper explores AIMs based on three robust Hilbert-space partitions - meta-Löwdin, Natural Atomic Orbitals (NAO), and Intrinsic Atomic Orbitals (IAO) - which combine the advantages of real-space partitions without their disadvantages. These partitions can effectively replace real-space AIMs for evaluating the aromatic character. For the first time, we report multicenter index (MCI) and Iring values for large rings and introduce ESIpy, an open-source Python code for aromaticity analysis in large conjugated rings.

Theoretical Modeling of Sr(q+) He (q = 0, 1, 2) van der Waals Systems Including Spin-Orbit Coupling.

Using an ab initio methodology that incorporates pseudopotential technique in conjunction with pair potential approaches, core polarization potentials (CPP), large basis sets of Gaussian type, and full configuration interaction calculations, we investigate interaction of neutral and charged Srq+(q = 0,1,2) with helium atom. In this context, the core-core interaction of Sr2+-He is included using an accurately performed potential for the ground state at CCSD(T) level of calculation. Also, the potential energy curves and permanent and transition dipole moments of the ground state and numerous excited states have been performed respectively for Sr+He and SrHe systems. Subsequently, the spin-orbit effect is considered by utilizing a semiempirical method for states dissociating into Sr+(5p) + He, Sr+(6p) + He, Sr+(4d) + He, Sr+(5d) + He, Sr(5s5p) + He, and Sr(5s4d) + He. The spectroscopic constants of the Srq+(q = 0, 1, 2) He states, with and without spin-orbit interaction, are derived and assessed in comparison to the existing theoretical and experimental studies. Such comparison has revealed good agreement, especially, for the Sr+He ionic system. Additionally, the spin-orbit effect is considered for the X2Σ+ → 22Π1/2,3/2 and X2Σ+ → 32Σ1/2 + transition dipole moments for Sr+He.