Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1986

On the bonding in Fe2(CO)9

The present paper is concerned with an analysis of the bonding in Fe2(CO)9, taking into account the results of self-consistent-field (SCF) calculations. The basis sets used are the same as those employed in the study of Fe(CO)5 conducted by Bauschlicher and Bagus (1985). The Fe basis set starts with the 14s9p5d primitive set optimized by Wachters (1970), which is contracted to 8s4p5d. An SCF wave function, consisting of only the core electrons, was determined, and then electrons were slowly added to the virtual orbitals with the lowest eigenvalues. The conducted calculations show no direct Fe-Fe bonding, but rather bonding resulting from the bridging CO's. The bonding is found to involve both the 5 sigma and 2 pi(asterisk) orbitals of the bridging CO's.

Bauschlicher, C. W., Jr.↗
NASA NTRS2012

Accurate ab initio Quartic Force Fields of Cyclic and Bent HC2N Isomers

Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.

Inostroza, Natalia↗
NASA NTRS1985

Theoretical transition probabilities between the lowest 2S, 2P and 2D states of Na, K, Rb and Cs

Theoretical transition probabilities between the lowest 2S, 2P and 2D states of the alkali atoms Na through Cs have been computed using near Hartree-Fock quality Slater basis sets. The important core-valence correlation effects are incorporated explicitly by a configuration-interaction procedure. For Cs, the calculations were repeated using a Gaussian basis set so that relativistic effects could be incorporated through an effective core potential procedure. The best calculated electric quadrupole Einstein coefficients are Na(196.3/s), K(103.6/s), Rb(72.4/s) and Cs(19.7/s). Core-valence effects become increasingly important down the column, and reduce the quadrupole transition strengths to about the same degree as for the 2P-2S and 2D-2P dipole-allowed transitions. Relativistic effects increase the quadrupole moment of Cs, but less so than in Ba, presumably because the alkali 2D states are more diffuse.

Langhoff, S. R.↗
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↗
NASA NTRS1986

A full CI treatment of Ne atom - A benchmark calculation performed on the NAS CRAY 2

Full CI calculations are performed for Ne atom using Gaussian basis sets of up to triple-zeta plus double polarization quality. The total valence correlation energy through double, triple, quadruple and octuple excitations is compared for eight different basis sets. These results are expected to be an important benchmark for calibrating methods for estimating the importance of higher excitations.

Bauschlicher, C. W., Jr.↗
DOE OSTI2022

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

A correlated ab initio study of the A2 pi <-- X2 sigma+ transition in MgCCH

The A2 pi <-- X2 sigma+ transition in MgCCH was studied with correlation consistent basis sets and single- and multireference correlation methods. The A2 pi excited state was characterized in detail; the x2 sigma+ ground state has been described elsewhere recently. The estimated complete basis set (CBS) limits for valence correlation, including zero-point energy corrections, are 22668, 23191, and 22795 for the RCCSD(T), MRCI, and MRCI + Q methods, respectively. A core-valence correction of +162 cm-1 shifts the RCCSD(T) value to 22830 cm-1, in good agreement with the experimental result of 22807 cm-1.

Non-NASA Center↗
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity (EA) of oxygen is computed to be 1.287 eV, using 2p electron full configuration-interaction (CI) wave functions expanded in a 6s5p3d2f Slater-type orbital basis. The best complete active space self-consistent field - multireference CI (CASSCF-MRCI) result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell correlation increases the computed EA to 1.290 at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. This study clearly establishes the synergistic effect between the higher excitations and basis set completeness on the electron affinity when the 2s electrons are correlated.

Bauschlicher, C. W., Jr.↗
DOE OSTI2022

Evaluating fast methods for static polarizabilities on extended conjugated oligomers

Given the importance of accurate polarizability calculations to many chemical applications, coupled with the need for efficiency when calculating the properties of sets of molecules or large oligomers, we present a benchmark study examining possible calculation methods for polarizable materials. We first investigate the accuracy of the additive model used in GFN2, a highly-efficient semi-empirical tight-binding method, and the D4 dispersion model, comparing its predicted additive polarizabilities to ωB97XD results for a subset of PubChemQC and a compiled benchmark set of molecules spanning polarizabilities from approximately 3 Å 3 to 600 Å 3 , with some compounds in the range of approximately 1200–1400 Å 3 . Although we find additive GFN2 polarizabilities, and thus D4, to have large errors with polarizability calculations on large conjugated oligomers, it would appear an empirical quadratic correction can largely remedy this. We also compare the accuracy of DFT polarizability calculations run using basis sets of varying size and level of augmentation, determining that a non-augmented basis set may be used for large, highly polarizable species in conjunction with a linear correction factor to achieve accuracy extremely close to that of aug-cc-pVTZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Theoretical study of the 7Sigma(u)+ state of N2

Theoretical potentials for the 7Sigma(u)+ state are reported using both extended Slater and Gaussian basis sets. Electron correlation is included using the interacting correlated fragments (ICF), the singles plus doubles configuration-interaction (SDCI), and the coupled-pair functional (CPF) approaches. The best potential, corrected for basis set superposition errors, has a well depth of about 21/cm and an r(e) of 7.52 bohr. The inclusion of the nitrogen 2s correlation significantly increases the well depth. These results are in reasonably good agreement with the empirical potential of Ferrante and Stwalley (1983), and support the contention that spin-polarized atomic nitrogen should behave like a classical solid such as Ne. It is shown that the modified Buckingham potential used by Ferrante and Stwalley is consistent with the present ab initio potentials if smaller values are used for the dispersion coefficients of the nitrogen atom.

Partridge, H.↗
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity of oxygen is computed to be 1.287 eV, at the full CI level using a 6s5p3d 2f Slater-type orbital basis and correlating only the 2p electrons. The best CASSCF-MRCI result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell coorelation increases the computed EA to 1.290 eV at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. The higher excitation contribution to the electron affinity is found to increase substantially with basis set completeness, especially when the 2s electrons are correlated. Relativistic effects are shown to make a small (less than 0.01 eV) change in the EA.

Bauschlicher, C. W., Jr.↗
NASA NTRS1972

Performance of correlation receivers in the presence of impulse noise.

An impulse noise model, which assumes that each noise burst contains a randomly weighted version of a basic waveform, is used to derive the performance equations for a correlation receiver. The expected number of bit errors per noise burst is expressed as a function of the average signal energy, signal-set correlation coefficient, bit time, noise-weighting-factor variance and probability density function, and a time range function which depends on the crosscorrelation of the signal-set basis functions and the noise waveform. Unlike the performance results for additive white Gaussian noise, it is shown that the error performance for impulse noise is affected by the choice of signal-set basis function, and that Orthogonal signaling is not equivalent to On-Off signaling with the same average energy. Furthermore, it is demonstrated that the correlation-receiver error performance can be improved by inserting a properly specified nonlinear device prior to the receiver input.

Moore, J. D.↗
NASA NTRS1987

Theoretical study of the dipole moment function of OH(X 2Pi)

A complete theoretical electric dipole moment function mu(r) is computed for OH(X 2Pi) at the full configuration-interaction (FCI) level for DZP + diffuse Gaussian basis set. The energy and dipole moment are also determined with 2d functions on oxygen at seven r values. The following methods are compared to the FCI for both the oxygen 1d and 2d basis sets: single-reference singles-plus-doubles configuration interaction (SDCI), SDCI with Davidson's correction, the coupled pair functional approach (CPF) and a modified form of CPF, CASSCF, and CASSCF-MRSDCI, and the external contracted CI(CCI) method. The dipole moments are evaluated both as an expectation value and by finite-field methods. The results support the superiority of evaluating the dipole moment as an energy derivative rather than as an expectation value.

Langhoff, Stephen R.↗
NASA NTRS1986

Theoretical dipole moments for the first-row transition metal hydrides

Spectroscopic parameters (D sub e, r sub e, mu) are determined for the first-row transition metal hydrides using better than DZP basis sets at the modified coupled pair functional (MCPF) level. Extensive comparisons between MCPF and complete-active space self-consistent field (CASSCF)/MRCI calculations with natural orbital iterations, and studies with more extensive basis sets, show this level of treatment to supply an accurate and cost-effective treatment of these systems. For the transition metal hydrides, the bonding can arise from either the 3d(n)4s(2) or 3d(n + 1)4s(1) atomic asymptotes, or a mixture of both. Since the dipole moment arising from these two bonding mechanisms is very different, the dipole moment is found to be directly related to the 3d population. Thus, the magnitude of the dipole moments provides a sensitive test of the wave function, and gives insight into the nature of the bonding.

Chong, D. P.↗
NASA NTRS2017

The Failure of Correlation to Describe Out-of-Plane Carbon=Carbon Bending

Carbon-carbon multiply bonded systems are improperly described with standard correlation methods and basis sets. For computations of vibrational modes, the out-of- plane bends can be reported as imaginary at worst or simply too low at best. Utilizing the simplest of aromatic structures (cyclopropenylidene) and various levels of theory, this work diagnoses this known behavior for the first time. A combined 1-particle and n-particle basis set effect conspire to produce these non-physical results. When moving from sp2 to sp3 hybridization in the carbon atoms, the larger number of basis functions overcorrects the energy. This is exacerbated by correlation methods. These allow for occupation of the π and π∗ orbitals in the expanded wave function that combine with the hydrogen s orbitals. As a result, the improperly described space can be further and non-physically stabilized by post-Hartree-Fock correlation. This represents a fundamental problem with at least Hartree-Fock based methods of all flavors in describing carbon. Beyond being a flaw in quantum chemical theory, other repercussions will be present in computations regarding spectroscopy as well as energy and environmental studies where highly-accurate hydrocabon vibrational transitions or thermochemical data are needed.

Fortenberry, Ryan C.↗
NASA NTRS1991

Theoretical study of the low-lying bound states of O2

It is demonstrated that a complete-active-space self-consistent-field (CASSCF) (2p)/MRCI + Q (multireference configuration interaction with a Davidson correction) description in a (13s8p6d 4f2g)/((5s4p3d 2f1g) atomic natural orbits (ANO) basis set supplemented with diffuse functions provides a quantitative description of the six lowest states of O2. The calculated potentials are within 0.05 eV (1.2 kilocal/mol) of accurate experimental results. The importance of substantially expanding the primitive basis set has been investigated, and it is demonstrated that such expansions yield insignificant improvement in the spectroscopic constants. Potential energy curves have also been reported for the weakly bound states of O2. The 5Pi(g) state is estimated to have a D(e) of 0.16 +/- 0.03 eV. The upper bound of D(e) is found to be sufficiently large that the importance of this state as a precursor for the formation of O2 (b 1Sigma(t)(+)) and O(1S) should be reconsidered.

Partridge, Harry↗