SEARCH · Search NASA

Results for “Basis sets”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗

DOE OSTI2025

An Ab Initio Electronic Structure Investigation of the Ground and Excited States of ScH + , YH + , and LaH +

Multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), coupled-cluster singles, doubles, and perturbative triples [CCSD(T)], and frozen-core full configuration interaction (fcFCI) calculations were carried out using large, correlation-consistent basis sets to investigate the excited states of the Sc atom and the spin–free and spin–orbit coupled potential energy profiles, energetics, spectroscopic constants, and electron populations of low-lying states of MH + (M = Sc, Y, La). The core electron correlation effects, complete basis set effects, and spin–orbit coupling effects were also evaluated. The first four electronic states of all MH + are 1 2 Δ, 1 2 Σ + , 1 2 Π, and 2 2 Σ + with 1σ 2 1δ 1 , 1σ 2 2σ 1 , 1σ 2 1π 1 , and 1σ 2 3σ 1 single-reference electron configurations, respectively. These states of MH + can be represented by the M 2+ H – ionic structure. The ground states of ScH + , YH + , and LaH + are 1 2 Δ 3/2 , 1 2 Σ + 1/2 , and 1 2 Δ 3/2 with 55.45, 60.54, and 62.34 kcal/mol bond energies, respectively. The core electron correlation was found to be vital for gaining accurate predictions on the ground and excited state properties of MH + . The spin–orbit coupling effects are minor for ScH + but become substantial moving to YH + and LaH + . Overall, the results of this work are in good agreement with the limited set of experimental findings of MH + available in the literature and will be of use for future investigations. Furthermore, the theoretical approaches, findings, and trends reported here are expected to aid studies of similar species.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Toward a Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

By correlating only the 1-particle states occupied in the reference determinant, the conventional design for the single-reference R12/F12 explicitly correlated methods biases them toward the ground-state description, thereby making the treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP-projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation length scale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground-state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of ground-state coupled-cluster singles and doubles (CCSD) energies, vertical excitation energies, and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.

Hamiltonians↗

DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗

DOE OSTI2024

Revisiting the Optical Spectrum of the Plutonyl Ion (PuO 2 ) 2+ in 1 M HClO 4

The analysis of the solution absorption spectrum of the plutonyl ion in an aqueous environment was given by Eisenstein and Pryce (E&P) in 1968. In 2011 a new spectrum was published of the (PuO 2 ) 2+ ion in 1 M HClO 4 . We have been provided with the original data of this spectrum and have found in the data a previously unreported low-lying transition at 7385 cm –1 which we have assigned as a magnetic dipole transition. We have fit most of the near-infrared and optical transitions with Gaussian fits and tabulated a new energy level list up to 22,000 cm –1 which mostly agrees with the data of E&P. We assumed a crystal field of D ∞h (only axial symmetry) and utilized the intensity calculations published for the isoelectronic (NpO 2 ) 1+ ion using a complete basis set for the 5f 2 problem including the Coulombic, spin–orbit as well as the crystal field Hamiltonian. Our results differ substantially from those of E&P. Subsequently, we used a truncated Hamiltonian to try to establish the effects of assuming the σ antibonding orbitals are at such high energies that we can ignore their contributions to the lower lying φ and δ orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Spatial Signatures of Electron Correlation in Least-Squares Tensor Hypercontraction

Least Squares Tensor Hypercontraction (LS-THC) has received some attention in recent years as an approach to reduce the significant computational costs of wavefunc- tion based methods in quantum chemistry. However, previous work has demonstrated that the LS-THC factorization performs disproportionately worse in the description of wavefunction components (e.g. cluster amplitudes T 2 ) than Hamiltonian compo- nents (e.g. electron repulsion integrals (pq|rs)). This work develops novel theoretical methods to study the source of these errors in the context of the real-space T 2 kernel, and reports, for the first time, the existence of a “correlation feature” in the errors of the LS-THC representation of the “exchange-like” correlation energy EX and T 2 that is remarkably consistent across ten molecular species, three correlated wavefunctions, and four basis sets. This correlation feature portends the existence of a “pair-point kernel” missing in the usual LS-THC representation of the wavefunction, which critically depends upon pairs of grid points situated close to atoms and with inter-pair distances between one and two Bohr radii. These findings point the way for future LS-THC developments to address these shortcomings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗

DOE OSTI2025

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

We introduce the first version of GPU4P Y SCF, a module that provides GPU acceleration of methods in P Y SCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of P Y SCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗

DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ab Initio Structures and Energetics of Hydrated Flat and Terrace-Step Surfaces of Forsterite (Mg 2 SiO 4 )

Forsterite (Mg 2 SiO 4 ), a model divalent metal silicate mineral, has been extensively studied in the context of mineral carbonation. Although dissolution is a key step in this process, the mechanisms by which forsterite dissolves under high CO 2 conditions remain poorly understood. Atomistic simulations could aid in exploring these mechanisms, but it is essential first to understand the structures and energetics of the relevant forsterite surfaces. We present an ab initio study of the structure and surface energy at 0 K of the flat $(010), (110), (001), (111), (021), (101)$ and $(120)$ faces of forsterite using the density functional PBE Hamiltonian and a plane-wave basis set. Dry surfaces became stabilized upon hydration through the formation of bonds between surface Mg and O from water, as well as by the formation of hydrogen bonds. According to surface energy values, the stability order of the hydrated forsterite faces was found to be $(120) < (101) < (021) < (111) < (001) < (110) < (010)$. We also investigated the energetics of the terrace-step $(0\bar{41})$ surface as a model site for forsterite dissolution. Among all the facets, the $(0\bar{41})$ surface is the least stable termination in water. Hydration of Mg atoms on the $(0\bar{41})$ surface increases their susceptibility to dissolution. The presence of a step and its hydration destabilizes the terraces, making step retreat more likely than a dissolution front advancing along the [010] direction. This research will support future simulations to investigate forsterite dissolution in water under CO 2 -rich conditions.

PBE Hamiltonian↗

DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗

DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗