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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Quantum Computing to Accelerate High Fidelity Computational Materials Modeling

In this CIF we worked to develop quantum algorithms for material science simulations based on new ideas recently proposed on plane wave basis sets. Many body simulations are not generally performed in plane wave basis sets on classical hardware, but with newly proposed quantum algorithms, it is possible this will be a highly efficient basis set to run quantum simulations on quantum hardware. Using state of the art classical simulations we ran small test simulations to estimate the resources that will be needed to run such plane wave algorithms on quantum hardware. We demonstrate our approach for a series of atoms and molecular systems.

Norman Tubman↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2024

ASMS 2024 Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision Induced Dissociation

Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision Induced Dissociation Brittany D. M. Hodges, Christopher A. Zarzana, JungSoo Kim, Jonathan Martens, and W. C. M. Berden Introduction (120 words max) Effects of electronic structure on chemical bonding and reactivity play critical roles shaping the chemical bonding and reactivity behaviors of heavy elements. Understanding the role of f electrons in bond formation between the actinide-series elements like uranium with other ligands is critical for solving technical challenges associated with these heavy elements, important to nuclear fuel cycles, efficient separations of rare earth metals, and understanding the chemistry of stored nuclear fuels and related environmental management sites. In this study, we further examine the interactions between uranyl and the perchlorate ion in order to understand the structures of these ions through the use of IRMPD. Here we report the IRMPD spectra of [UO2(ClO4)3]-, [UO3(ClO4)2]-, and a proposed transition state. Methods (120 word max) IRMPD spectra and CID spectra were acquired using a Bruker amaZon QIT/MS instrument at the Free-Electron Lasers for Infrared eXperiments (FELIX) laboratory at Radboud University. The FELIX QIT/MS is modified to allow for the high-intensity tunable IR beam from FELIX to be directed into the ion packet, resulting in multiphoton dissociation that is measured only when the IR frequency is in resonance with an adequately high absorption vibrational mode of the mass-selected complex. DFT geometry optimizations and frequency calculations using the Gaussian suite of programs were performed using B3LYP, TPSSh, and PBE0 level of theory with 6-31+G(d) basis for the O, C, H, and N atoms and the SDD basis set for U. The SDD basis set employs the Stuttgart/Dresden effective core potential. Preliminary Data or Plenary Speakers Abstract (300 words max) Metal ion clusters of uranyl perchlorate were formed via direct electrospray ionization. For each metal ligand complex of interest, the parent ion was isolated and collision induced dissociation fragmentation and Infrared Multiphoton Dissociation (IRMPD) fragmentation spectra were acquired. Results presented here are the first look at the IRMPD spectra of [UO3(ClO4)2]-, [UO2(ClO4)3]-. Structures were examined using Gaussian at different levels of theory B3LYP level of theory, TPPSh and PBE0 levels, to reflect the behaviors of uranium metal ligand complexes most accurately. In these structures, we identified an overlap between each of these uranyl stretches resulting in their largely being obscured by a perchlorate mode. The CID product spectra agree with similar structures reported by Groenewold for uranyl nitrate in 2006 (10.1021/ja058106n).

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

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗
NASA NTRS1995

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗
NASA NTRS1990

Theoretical characterization of the potential energy surface for H + N2 yields HN2. II - Computed points to define a global potential

A previous calculation for H + N2 (Walch et al., 1989) focused on the minimum energy path (MEP) region of the potential energy surface and on estimates of the lifetime of the HN2 species. In this paper, energies computed at geometries selected to permit a global representation of the potential energy surface (PES) are reported. As in the previous work, the calculations were performed using the complete active space self-consistent field/externally contracted configuration interaction method. The surface was characterized using the same basis set as in the previous paper except that an improved contraction of the H s-basis is used. Calculations with a larger basis set were carried out along an approximate MEP obtained with the smaller basis set. The new PES exhibits a sharp curvature, which was not present in the previous calculations, and has a slightly narrower and smaller barrier to dissociation. Saddle points for H atom exchange via collinear and T shaped HN2 complexes are also reported.

Walch, Stephen P.↗
NASA NTRS2022

The reliability of the small-core Lanthanide effective core potentials

The reliability of the small-core Lanthanide effective core potentials (ECP) is tested using MF and MF(3), for M=Eu, Gd, Tb, and Yb and the atomic excitation energies for Pr, Nd, Eu, Gd, Tb, Dy, Ho, Er, Tm, and Yb. In some case the ECP and all-electron (AE) results are in good agreement, while in others there are significant difference. The difference are much larger when the segmented basis set is used in conjunction with the ECP than when the atomic natural orbital (ANO) basis set is used. The study of the atoms suggests that problems for lanthanide containing molecules are associated with poor atomic excitation energies in the ECP treatment and even using the ANO basis set does not completely solve the problem. We note that the problem appears to be more severe for density functional approaches than for traditional correlation methods. We suggest that additional studies and new effective core potentials may be required for the Lanthanide atoms.

effective core potentials, atomic excitation, atom↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

An Accurate ab initio Quartic Force Field and Vibrational Frequencies for CH4 and Isotopomers

A very accurate ab initio quartic force field for CH4 and its isotopomers is presented. The quartic force field was determined with the singles and doubles coupled-cluster procedure that includes a quasiperturbative estimate of the effects of connected triple excitations, CCSD(T), using the correlation consistent polarized valence triple zeta, cc-pVTZ, basis set. Improved quadratic force constants were evaluated with the correlation consistent polarized valence quadruple zeta, cc-pVQZ, basis set. Fundamental vibrational frequencies are determined using second-order perturbation theory anharmonic analyses. All fundamentals of CH4 and isotopomers for which accurate experimental values exist and for which there is not a large Fermi resonance, are predicted to within +/- 6 cm(exp -1). It is thus concluded that our predictions for the harmonic frequencies and the anharmonic constants are the most accurate estimates available. It is also shown that using cubic and quartic force constants determined with the correlation consistent polarized double zeta, cc-pVDZ, basis set in conjunction with the cc-pVQZ quadratic force constants and equilibrium geometry leads to accurate predictions for the fundamental vibrational frequencies of methane, suggesting that this approach may be a viable alternative for larger molecules. Using CCSD(T), core correlation is found to reduce the CH4 r(e), by 0.0015 A. Our best estimate for r, is 1.0862 +/- 0.0005 A.

Lee, Timothy J.↗
NASA NTRS1988

Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1989

Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F + H2 yields HF + H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(CO)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule O4(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI,and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in O4(+) has also been studied using the Brueckner coupled-cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
NASA NTRS1991

Determination of a natural basis function set on thick, structured slabs using Prony's method

The calculation of currents on a thick, structured slab such as a thick slab of honeycomb is discussed. Unfortunately, for certain applications, the slab can be on the order of several wavelengths, so that the straightforward application of the method of moments using subdomain basis functions is too expensive. The authors discuss how to apply Prony's method to the currents calculated for a thin, structured slab to obtain a natural set of basis functions to represent the currents in the interior of a thick slab. Prony's method represents the currents as a series of complex exponential functions. The thick slab problem is then solved by the method of moments using subdomain basis functions near the slab interfaces and one or two of the complex exponentials as basis functions within the slab.

Jorgenson, Roy E.↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
NASA NTRS1972

Electric dipole moment of diatomic molecules by configuration interaction. IV.

The theory of basis set dependence in configuration interaction calculations is discussed, taking into account a perturbation model which is valid for small changes in the self-consistent field orbitals. It is found that basis set corrections are essentially additive through first order. It is shown that an error found in a previously published dipole moment calculation by Green (1972) for the metastable first excited state of CO was indeed due to an inadequate basis set as claimed.

Green, S.↗
NASA NTRS1980

Theoretical studies of the HeCN/+/ and NeCN/+/ molecular ions

Self-consistent field calculations are presented for HeCN(+) using several basis sets; calculations for HeH(+) with the same basis sets are given for comparison. These predict that HeCN(+) is stable with a He-CN(+) dissociation energy between 1.5 and 2.0 eV. The binding is found to be mainly electrostatic in nature. Rotation and hyperfine spectral constants are predicted. Similar calculations indicate that the isocyanide HeCN(+) is not stable. Preliminary calculations on NeCN(+) indicate that it is as stable as HeCN(+). The possibility of observing noble gas molecular ions in laboratory experiments and in interstellar space is discussed.

Wilson, S.↗