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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Quantum Chemical Evaluation of the Astrochemical Significance of Reactions between S Atom and Acetylene or Ethylene

Addition-elimination reactions of S atom in its P-3 ground state with acetylene (C2H2) and ethylene (C2H4) were characterized with both molecular orbital and density functional theory calculations employing correlation consistent basis sets in order to assess the likelihood either reaction might play a general role in astrochemistry or a specific role in the formation of S2 (X (sup 3 SIGMA (sub g) (sup -)) via a mechanism proposed by Saxena and Misra (Mon. Not. R. Astron. Soc. 1995, 272, 89). The acetylene and ethylene reactions proceed through C2H2S ((sup 3)A")) and C2H4S ((sup 3)A")) intermediates, respectively, to yield HCCS ((sup 2)II)) and C2H3S ((sup 2)A')). Substantial barriers were found in the exit channels for every combination of method and basis set considered in this work, which effectively precludes hydrogen elimination pathways for both S + C2H2 and S + C2H4 in the ultracold interstellar medium where only very modest barriers can be surmounted and processes without barriers tend to predominate. However, if one or both intermediates is formed and stabilized efficiently under cometary or dense interstellar cloud conditions, they could serve as temporary reservoirs for S atom and participate in reactions such as S + C2H2S (right arrow) S2 = C2H2 or S + C2H4S (right arrow) S2 + C2H4. For formation and stabilization to be efficient, the reaction must possess a barrier height small enough to be surmountable at low temperatures yet large enough to prevent redissociation to reactants. Barrier heights computed with B3LYP and large basis sets are very low, but more rigorous QCISD(T) and RCCSD(T) results indicate that the barrier heights are closer to 3-4 kcal/mol. The calculations therefore indicate that S + C2H2 or S + C2H4 could contribute to the formation of S2 in comets and may serve as a means to gauge coma temperature. The energetics of the ethylene reaction are more favorable.

Woon, David E.↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
NASA NTRS1987

On the bond distance in methane

The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Bowen-Jenkins, Philippa↗
NASA NTRS1988

On the bond distance in methane

The equilibrium bond distance in methane has been optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valance correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Bowen-Jenkins, Philippa↗
NASA NTRS1993

Structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structures and energetics of Cr(CO)6 and Cr(CO)5 are determined at the modified coupled-pair functional, single and double excitation coupled-cluster (CCSD), and CCSD(T) levels of theory. For Cr(CO)6, the structure and force constants for the totally symmetric representation are in good agreement with experimental data once basis set constants are taken into account. In the largest basis set at the CCSD(T) level of theory, the total binding energy of CR(CO)6 is estimated at around 140 kcal/mol, or about 86 percent of the experimental value. In contrast, the first bond energy of Cr(CO)6 is very well described at the CCSD(T) level of theory, with the best estimated value of 38 kcal/mol being within the experimental uncertainty.

Barnes, Leslie A.↗
DOE OSTI2021

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

GDB-9-Ex_EOM-CCSD: Dataset containing Equation of Motion Coupled Cluster (EOM-CCSD) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 80,593 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Equation of Motion Coupled Cluster (EOM-CCSD) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The similarity-transformed EOM-CCSD method that used domain-based local pair natural orbitals (DLPNO) approximation which constitutes the STEOM-DLPNO-CCSD method was used. This method is based on the STEOM approach and was found to make accurate predictions of transition energies for organic molecules. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2024

GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

Building the Case for SNAP: Creation of Multi-Band, Simulated Images With Shapelets

Dark energy has simultaneously been the most elusive and most important phenomenon in the shaping of the universe. A case for a proposed space-telescope called SNAP (SuperNova Acceleration Probe) is being built, a crucial component of which is image simulations. One method for this is "Shapelets," developed at Caltech. Shapelets form an orthonormal basis and are uniquely able to represent realistic space images and create new images based on real ones. Previously, simulations were created using the Hubble Deep Field (HDF) as a basis Set in one band. In this project, image simulations are created.using the 4 bands of the Hubble Ultra Deep Field (UDF) as a basis set. This provides a better basis for simulations because (1) the survey is deeper, (2) they have a higher resolution, and (3) this is a step closer to simulating the 9 bands of SNAP. Image simulations are achieved by detecting sources in the UDF, decomposing them into shapelets, tweaking their parameters in realistic ways, and recomposing them into new images. Morphological tests were also run to verify the realism of the simulations. They have a wide variety of uses, including the ability to create weak gravitational lensing simulations.

Ferry, Matthew A.↗
DOE OSTI2025

Method-independent cusps for atomic orbitals in quantum Monte Carlo

Here, we present an approach for augmenting Gaussian atomic orbitals with correct nuclear cusps. Like the atomic orbital basis set itself and unlike previous cusp corrections, this approach is independent of the many-body method used to prepare wave functions for quantum Monte Carlo. Once the basis set and molecular geometry are specified, the cusp-corrected atomic orbitals are uniquely specified, regardless of which density functionals, quantum chemistry methods, or subsequent variational Monte Carlo optimizations are employed. We analyze the statistical improvement offered by these cusps in a number of molecules and find them to offer similar advantages as molecular-orbital-based approaches while remaining independent of the choice of many-body method.

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

Accurate quantum chemical calculations

An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

An Accurate Quartic Force Field and Vibrational Frequencies for HNO and DNO

An accurate ab initio quartic force field for HNO has been determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with the correlation consistent polarized valence triple zeta (cc-pVTZ) basis set. Improved harmonic frequencies were determined with the cc-pVQZ basis set. Fundamental vibrational frequencies were determined using a second-order perturbation theory analysis and also using variational calculations. The N-0 stretch and bending fundamentals are determined well from both vibrational analyses. The H-N stretch, however, is shown to have an unusually large anharmonic correction, and is not well determined using second-order perturbation theory. The H-N fundamental is well determined from the variational calculations, demonstrating the quality of the ab initio quartic force field. The zero-point energy of HNO that should be used in isodesmic reactions is also discussed.

Dateo, Christopher E.↗
NASA NTRS1994

Accurate Energetics for the Unimolecular Decomposition of HN2

We report internally contracted configuration interaction (ICCI) calculations along the constrained energy minimum (CEM) path for the addition of H to N2 with both CC-pVQZ and augmented CC-pVQZ basis sets. The results agree to within 0.5 kcal/mol with the larger basis set results of Walch, which were used to calibrate the HN2 potential energy surface used by Koizumi, Schatz, and Walch in coupled channel calculations to determine the lifetimes for unimolecular decay of the lowest vibrational levels of HN2. These small changes in the potential energy surface should have no significant effect on the computed lifetimes.

Walch, Stephen P.↗
DOE OSTI2024

Convergent Concordant Mode Approach for Molecular Vibrations: CMA-2

The concordant mode approach (CMA) is a promising new scheme for dramatically increasing the system size and level of theory achievable in quantum chemical computations of molecular vibrational frequencies. Here, we achieve advances in the CMA hierarchy by computations targeting CCSD(T)/cc-pVTZ (coupled cluster singles and doubles with perturbative triples using a correlation-consistent polarized-valence triple-ζ basis set) benchmarks within the G2 molecular test set, executing a statistical analysis for 1501 frequencies from 111 compounds and then separately solving the refractory case of pyridine. First, MP2/cc-pVTZ (second-order Møller–Plesset perturbation theory with the same basis set) proves to be an excellent and preferred choice for generating the underlying (Level B) normal modes of the CMA scheme. Utilizing this Level B within the CMA-0A method reproduces the 1501 benchmark frequencies with a mean absolute error (MAE) of only 0.11 cm –1 and an attendant standard deviation of 0.49 cm –1 . Second, a convergent CMA-2 method is constituted that allows efficient computation of higher level (Level A) frequencies to any reasonable accuracy threshold by using only Hartree–Fock (HF) and MP2 or density functional theory (DFT) data to generate ξ parameters, which select the sparse off-diagonal force field elements for explicit evaluation at Level A. When Level B = MP2/cc-pVTZ, a cutoff of ξ = 0.02 provides an average maximum absolute error per molecule of only 0.17 cm –1 by incurring merely a 33% increase in average cost over CMA-0A. This CMA-2 method also eradicates the 4 problematic CMA-0A outliers of pyridine with even less effort (ξ = 0.04, 22% increase). Finally, the newly developed CMA procedures are shown to be highly successful when applied to 1-(1H-pyrrol-3-yl)ethanol, a new test molecule with diverse types of vibration.

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

The Heats of Formation of GaCl3 and its Fragments

The heats of formation of GaC13 and its fragments are computed. The geometries and frequencies are obtained at the B3LYP level. The CCSD(T) approach is used to solve the correlation problem. The effect of Ga 3d correlation is studied, and found to affect the bond energies by up to 1 kcal/mol. Both basis set extrapolation and bond functions are considered as ways to approach the basis set limit. Spin-orbit and scalar relativistic effects are also considered.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2022

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗