Search NASA⌕ Search

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.

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

At least 163 records · Page 9

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.↗

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prospects for detecting UF6 hydrolysis intermediates by mass spectrometry, resonance Raman, and NQR spectroscopy

Simulations are performed to consider several spectrometric and spectroscopic candidates for elucidating the mechanism of the hydrolysis of uranium hexafluoride (UF6). This study is among the first to benchmark the def-mTZVP basis sets for actinide-containing molecules, and it is shown to be a suitable basis set for surveying geometrical structures and vibrational spectra when used in conjunction with density functional theory. An experiment is proposed coupling mass spectrometric ion selection with vibrational spectroscopy, and supporting infrared and Raman spectral simulations demonstrate that ionization blue-shifts bands and can change their qualitative features. Ultraviolet resonance Raman is shown to have good prospects for discriminating U–O–U bridged intermediates, evidence for which was observed recently [L. E. McNamara et al., J. Phys. Chem. A 130, 775–786 (2026)]. As a more speculative approach, we also consider the prospect of addressing 233U or 235U nuclei by quadrupole resonance (NQR) spectroscopy for assigning early stage intermediate complexes. In doing so, we provide a first order-of-magnitude estimate for the collision-induced NQR signal for the UF6 dimer, which is of interest for the interpretation of a fast decoherence time observed in liquid-phase nuclear magnetic resonance. Having provided critical insights into the formation and stability of UF6 hydrolysis intermediates in previous studies, this computational spectroscopy survey is expected to help guide and expedite future laboratory campaigns.

Lutz, Jesse J. [Center for Computing Research, San↗

Calculation of Mg(+)-ligand relative binding energies

The calculated relative binding energies of 16 organic molecules to Mg(+) are compared with experimental results where available. The geometries of the ligands and the Mg(+)-ligand complexes arc optimized at the self-consistent field level using a 6-31G* basis set. The Mg(+) binding energies are evaluated using second-order perturbation theory and basis sets of triple-sigma quality augmented with two sets of polarization functions. This level of theory is calibrated against higher levels of theory for selected systems. The computed binding energies are accurate to about 2 kcal/mol.

Partridge, Harry↗

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗

The 2S(+) - 2P separation in KO

The 2S(+) - 2P separation in KO is investigated using large basis sets and high levels of correlation treatment. Relativistic effects are included at the Dirac-Fock level and reduce the separation only slightly. The basis set superposition error is considered in detail. On the basis of these calculations, our best estimate places the 2p sub 3/2 state about 200 cm(exp -1) above the ground 2 sigma(+) state in agreement with our previous estimate.

Bauschlicher, Charles W., Jr.↗

The 2Sigma(+)-2Pi separation in KO

The 2Sigma(+)-2Pi separation in KO is investigated using large basis sets and high levels of correlation treatment. Relativistic effects are included at the Dirac-Fock level and reduce the separation only slightly. The basis set superposition error is considered in detail. On the basis of these calculations, our best estimate places the 2Pi3/2 state about 200/cm above the ground 2Sigma(+) state, in agreement with our previous estimate.

Bauschlicher, Charles W., Jr.↗

Theoretical dissociation energies for ionic molecules

Ab initio calculations at the self-consistent-field and singles plus doubles configuration-interaction level are used to determine accurate spectroscopic parameters for most of the alkali and alkaline-earth fluorides, chlorides, oxides, sulfides, hydroxides, and isocyanides. Numerical Hartree-Fock (NHF) calculations are performed on selected systems to ensure that the extended Slater basis sets employed for the diatomic systems are near the Hartree-Fock limit. Extended Gaussian basis sets of at least triple-zeta plus double polarization equality are employed for the triatomic system. With this model, correlation effects are relatively small, but invariably increase the theoretical dissociation energies. The importance of correlating the electrons on both the anion and the metal is discussed. The theoretical dissociation energies are critically compared with the literature to rule out disparate experimental values. Theoretical (sup 2)Pi - (sup 2)Sigma (sup +) energy separations are presented for the alkali oxides and sulfides.

Langhoff, S. R.↗

Calculations of H2O microwave line broadening in collisions with He atoms - Sensitivity to potential energy surfaces

Theoretical computations of broadening parameters are reported for three microwave lines of H2O in a bath of He atoms. The potential-energy surfaces employed are corrected for basis-set superposition error, and their reliability is checked by repeating the calculations with a different basis set for orbital expansion. The results are presented in extensive tables and discussed in detail. The corrections applied are shown to have a significant impact on the accuracy of the room-temperature broadenings determined: 8.9 sq A for the 22.2-GHz line, 11.8 sq A for the 183,3-GHz line, and 10.0 sq A for the 380.2-GHz line, in good agreement with published experimental data. The importance of collisional broadening for the atmospheric transmission of radiation and for remote-sensing applications is indicated.

Green, Sheldon↗

Theoretical study of the hyperfine parameters of OH

In the present study of the hyperfine parameters of O-17H as a function of the one- and n-particle spaces, all of the parameters except oxygen's spin density, b sub F(O), are sufficiently easily tractable to allow concentration on the computational requirements for accurate determination of b sub F(O). Full configuration-interaction (FCI) calculations in six Gaussian basis sets yield unambiguous results for (1) the effect of uncontracting the O s and p basis sets; (2) that of adding diffuse s and p functions; and (3) that of adding polarization functions to O. The size-extensive modified coupled-pair functional method yields b sub F values which are in fair agreement with FCI results.

Chong, Delano P.↗

The Calculation of Accurate Harmonic Frequencies of Large Molecules: The Polycyclic Aromatic Hydrocarbons, a Case Study

The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.

Bauschlicher, Charles W., Jr.↗

Nonlinear Dynamics of Cantilever-Sample Interactions in Atomic Force Microscopy

The interaction of the cantilever tip of an atomic force microscope (AFM) with the sample surface is obtained by treating the cantilever and sample as independent systems coupled by a nonlinear force acting between the cantilever tip and a volume element of the sample surface. The volume element is subjected to a restoring force from the remainder of the sample that provides dynamical equilibrium for the combined systems. The model accounts for the positions on the cantilever of the cantilever tip, laser probe, and excitation force (if any) via a basis set of set of orthogonal functions that may be generalized to account for arbitrary cantilever shapes. The basis set is extended to include nonlinear cantilever modes. The model leads to a pair of coupled nonlinear differential equations that are solved analytically using a matrix iteration procedure. The effects of oscillatory excitation forces applied either to the cantilever or to the sample surface (or to both) are obtained from the solution set and applied to the to the assessment of phase and amplitude signals generated by various acoustic-atomic force microscope (A-AFM) modalities. The influence of bistable cantilever modes of on AFM signal generation is discussed. The effects on the cantilever-sample surface dynamics of subsurface features embedded in the sample that are perturbed by surface-generated oscillatory excitation forces and carried to the cantilever via wave propagation are accounted by the Bolef-Miller propagating wave model. Expressions pertaining to signal generation and image contrast in A-AFM are obtained and applied to amplitude modulation (intermittent contact) atomic force microscopy and resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM). The influence of phase accumulation in A-AFM on image contrast is discussed, as is the effect of hard contact and maximum nonlinearity regimes of A-AFM operation.

Cantrell, John H.↗

The vibrational frequencies of difluoroethyne

Ab initio coupled-cluster calculations with single and double excitations and with a perturbational treatment of connected triple excitations are reported for difluoroethyne using large basis sets. The results for the transbending mode nu-4 are extremely sensitive to electron correlation and basis set effects. The best theoretical and experimental estimates for the fundamental vibrational frequencies are in excellent agreement.

Breidung, Juergen↗

On the role of the seasonal cycle

The definition of the seasonal cycle as the projection of an atmospheric time series onto a suitably defined subset of orthogonal basis functions guarantees that any atmospheric covariance may be expressed as the sum of a seasonal cycle part and a transient part (where transient refers to departures from the seasonal cycle rather than the time mean). Results are presented for the seasonal cycle contribution to the zonally averaged fluxes of momentum and heat, and for the zonally averaged height variance, in data for: the winter season, the winter season with the Fourier basis set appropriate for the seven-year time series, and the seven-year mean. The Fourier basis set calculation indicates that the interannual variability of the momentum flux is dominated by the interactions between the seasonal cycle and the meteorologically low frequency flow.

Straus, D. M.↗

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller–Plesset perturbation method

Here, using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller–Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ~0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗