Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical calculations, Polarization”

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 OSTI2020

Ion Pairing and Redissociaton in Low-Permittivity Electrolytes for Multivalent Battery Applications

Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.

25 ENERGY STORAGE↗
DOE OSTI2022

Magnetic transition behavior in epitaxial Fe 47 Rh 47 Pd 6 films

The properties of Fe–Rh–Pd epitaxial thin films grown on MgO(001) were studied as a function of growth temperature. Films grown above 400°C exhibit a first-order antiferromagnetic to ferromagnetic magnetic phase transition with a transition temperature that decreases as the growth temperature is increased. The chemical order parameter computed from the ratio of intensities of the (001) and (002) diffraction peaks is nearly independent of the growth temperature, while the lattice constants change slightly. A comparison of our structural, magnetic, and electrical transport results with first-principle-based calculations as well as literature results indicates that the transition temperature of Fe–Rh-based alloy films depends sensitively on the lattice parameters and is of electronic origin. The transition temperature and its width can be tuned over a wide range by controlling the crystal structure via growth conditions or postdeposition annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multifunctional Cu 2 TSiS 4 (T = Mn and Fe): Polar Semiconducting Antiferromagnets with Nonlinear Optical Properties

In this work, Cu 2 TSiS 4 (T = Mn and Fe) polycrystalline and single-crystal materials were prepared with high-temperature solid-state and chemical vapor transport methods, respectively. The polar crystal structure (space group Pmn2 1 ) consists of chains of corner-sharing and distorted CuS 4 , Mn/FeS 4 , and SiS 4 tetrahedra, which is confirmed by Rietveld refinement using neutron powder diffraction data, X-ray single-crystal refinement, electron diffraction, energy-dispersive X-ray spectroscopy, and second harmonic generation (SHG) techniques. Magnetic measurements indicate that both compounds order antiferromagnetically at 8 and 14 K, respectively, which is supported by the temperature-dependent (100–2 K) neutron powder diffraction data. Additional magnetic reflections observed at 2 K can be modeled by magnetic propagation vectors k = (1/2,0,1/2) and k = (1/2,1/2,1/2) for Cu 2 MnSiS 4 and Cu 2 FeSiS 4 , respectively. The refined antiferromagnetic structure reveals that the Mn/Fe spins are canted away from the ac plane by about 14°, with the total magnetic moments of Mn and Fe being 4.1(1) and 2.9(1) μ B , respectively. Both compounds exhibit an SHG response with relatively modest second-order nonlinear susceptibilities. Density functional theory calculations are used to describe the electronic band structures.

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

Activation of Low-Valent, Multiply M–M Bonded Group VI Dimers toward Catalytic Olefin Metathesis via Surface Organometallic Chemistry

Olefin metathesis is a broadly employed reaction with applications that range from fine chemicals to materials and petrochemicals. The design and investigation of olefin metathesis catalysts have been ongoing for over half a century, with advancements made in terms of activity, stability, and selectivity. Immobilization of organometallic complexes onto solid supports such as silica or alumina is a promising strategy for catalyst heterogenization, often resulting in increased activity and stability. Consequently, a broad range of early transition metal catalysts bearing alkyl, oxide/alkoxide, and amide ligands have been grafted onto silica and their reactivities investigated. In this paper, we report a series of silica-supported tungsten and molybdenum dimers (X 3 M≡MX 3, where M = W and Mo; X = neopentyl, tert-butoxide, and dimethyl amide) and their reactivities toward catalytic olefin metathesis. Dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (SSNMR), diffuse reflectance infrared Fourier transform (DRIFT), UV resonance Raman, and X-ray absorption (XAS) spectroscopies suggest that upon heterogenization the dimers bind to the surface in a monopodal fashion, with the M≡M triple bond remaining intact. These structural assignments were further corroborated by density functional theory (DFT) calculations. While the homogeneous dimer counterparts are inert, the supported low-valent alkyl W and Mo dimers become active for the disproportionative self-metathesis of propylene to ethylene and butenes and 4-nonene to 4-octene and 5-decene under mild conditions. The lack of activity observed for the free and supported tert-butoxide and dimethyl amide dimers likely suggests that the neopentyl groups are necessary for the formation of a putative alkylidene active species. The difference in reactivity between the free and supported dimers could be explained either by the lowering of the activation barrier of the complex through the electronic effects of the surface or by site isolation of catalytically relevant reactive intermediates.

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

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

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

Structure and thermodynamics of water adsorption in NU-1500-Cr

Abstract Metal-organic frameworks (MOFs) are a class of materials with diverse chemical and structural properties, and have been shown to effectively adsorb various types of guest molecules. The mechanism of water adsorption in NU-1500-Cr, a high-performance atmospheric water harvesting MOF, is investigated using a combination of molecular dynamics simulations and infrared spectroscopy. Calculations of thermodynamic and dynamical properties of water as a function of relative humidity allow for following the adsorption process from the initial hydration stage to complete filling of the MOF pores. Initial hydration begins at the water molecules that saturate the open Cr 3+ sites of the framework, which is then followed by the formation of water chains that extend along the channels connecting the hexagonal pores of the framework. Water present in these channels gradually coalesces and fills the hexagonal pores sequentially after the channels are completely hydrated. The development of hydrogen-bond networks inside the MOF pores as a function of relative humidity is characterized at the molecular level using experimental and computational infrared spectroscopy. A detailed analysis of the OH-stretch vibrational band indicates that the low-frequency tail stems from strongly polarized hydrogen-bonded water molecules, suggesting the presence of some structural disorder in the experimental samples. Strategies for designing efficient water harvesting MOFs are also proposed based on the mechanism of water adsorption in NU-1500-Cr.

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

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

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

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Enhanced Interlayer Coupling and Excitons in Twin-Stacked Two-Dimensional Magnetic CrSBr Bilayers

The degree of electronic coupling between individual layers in van der Waals heterostructures offers a route to engineer their magnetic, electronic, and optical functionalities. Using state-of-the-art first-principles calculations, in this work, we demonstrate that the electronic coupling between two monolayers of CrSBr─an anisotropic two-dimensional magnetic semiconductor─is highly nonlinear and nonmonotonic with respect to their relative twist angle, exhibiting a pronounced maximum at the twin-stacking configuration. The coupling strength scales with both the degree of overlap of Br orbitals adjacent to the van der Waals gap and the cosine of half of the interlayer spin angle. This enhanced interlayer electronic coupling leads to excitons delocalized across the two layers, with a polarization dependence that reflects the interlayer spin alignment. Our results reveal a sensitive interplay among twist angle, magnetism, and excitonic properties in twin-stacked CrSBr bilayers, suggesting twin stacking as an effective means for engineering interlayer coupling.

2D materials↗
DOE OSTI2024

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗
DOE OSTI2020

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

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

Effects of Hydrogen Bonding and Molecular Chain Flexibility of Substituted n-Alkyldimethylsilanes On Impact Ice Adhesion Shear Strength

The effects of hydrogen bonding and molecular flexibility upon ice adhesion shear strength were investigated using aluminum substrates coated with substituted n-alkyldimethylalkoxysilanes. The location of the chemical group substitution was on the opposing end of the linear n-alkyl chain with respect to silicon. Three hydrogen-bonding characteristics were evaluated: 1) non-hydrogen bonding, 2) donor/acceptor, and 3) acceptor. Varying the length of the n-alkyl chain provided an assessment of molecular chain flexibility. Coated and uncoated aluminum surfaces were characterized by receding water contact angle and surface roughness. Ice adhesion shear strength was determined in the Adverse Environment Rotor Test Stand facility from -16 to -8°C that simulated aircraft in-flight icing conditions within the FAR Part 25/29 Appendix C icing envelope. Surface roughness of the coatings was similar allowing for comparison of the test results. An adhesion reduction factor, based on the ice adhesion shear strength data with respect to uncoated aluminum obtained at the same temperature, was calculated to compare the data. The results revealed complex interactions with impacting supercooled water droplets that were interdependent upon ice accretion temperature, surface energy characteristics of water and ice, hydrogen bonding characteristic of the substituent, and length of the n-alkyl chain. To aid in explaining the results, 1) changes in the surface energy component (i.e., non-polar and polar) values that water undergoes during its’ phase change from liquid to solid that arise from the freezing of impacting supercooled water droplets on the surface depended upon the temperature during accretion were taken into account and 2) the physical properties (i.e., water solubility and melting point) of small compounds analogous to the substituted n-alkyldimethylalkoxysilanes used in this study were compared.

Smith, Joseph G.↗
DOE OSTI2021

Revealing the Configuration and Conformation of Surface Organometallic Catalysts with DNP-Enhanced NMR

Although single-site, supported organometallic catalysts were designed with homogeneity in mind, little is strictly known about their atomic-level structure and uniformity. This is in large part due to the inability of conventional characterization tools to provide structural information with adequate resolution and range for this purpose. In this work, we show that dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (NMR) enables the measurement of distances between the surface and individual carbons, making it possible to orient complexes and molecules on surfaces. We use this approach to determine the orientation and configuration of naturally 13 C-abundant acetate and supported organoiridium(III) pincer complex, both supported on γ-Al 2 O 3 . By combining solid-state NMR and extended X-ray absorption fine structure spectroscopy (EXAFS) experiments with periodic density functional theory (DFT) calculations, we are able to determine the three-dimensional arrangement of ligands around the metal center and thereby identify the structure of the supported complex at a resolution reminiscent of that associated with single-crystal diffractometry.

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

First-principles modeling of conductivity at the (001), (110), and (111) SrTiO 3 /LaAlO 3 heterointerfaces

Here, the complex polar oxide heterojunction of SrTiO 3 /LaAlO 3 (STO/LAO) is of great interest due to the emergent physical phenomena observed at the interface. STO and LAO separately are wide band-gap insulators. However, upon joining them at the 001, 110, and 111 crystallographic planes, the interface undergoes a transition to a conductive state. Although first-principles modeling of the 001 plane interface has been widely studied, there is a lack of reports regarding the 110 and 111. This paper expands the theoretical model of the STO/LAO heterointerface to the three crystallographic planes (001, 110, and 111) where the conductivity has been experimentally reported. The calculations showed that whereas at the 001 interface the conductivity appears at a critical thickness of 4 monolayers of LAO, the 110 and 111 planes have no clear critical thickness; these two interfaces were always conductive. Nevertheless, the number of conductive electrons per unit cell increases with the thickness of the LAO layer in the 110 and 111 interfaces. This is related to the energy levels downshifting due to the electrostatic potential buildup (which was in the opposite direction respect to the 001 interface), increasing the number of conductive sates below the Fermi level. Given the absence of a critical thickness and the fact that chemical intermixing and oxygen vacancies at the interface were not considered, the main mechanism responsible for the conductivity in the 110 and 111 planes was attributed to the large structure reconstruction that locally changes the energy levels at the interface causing charge transfer and accumulation at the layers close to the interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of high-dose gamma-ray irradiation on electrical characteristics of N-polar and Ga-polar GaN p–n diodes

We investigate the impact of high-dose gamma-ray irradiation on the electrical performance of Ga-polar and N-polar GaN-based p-n diodes grown by metalorganic chemical vapor deposition. We compare the current density-voltage (J-V), capacitance-voltage (C-V), and circular transfer length method (CTLM) characteristics of the p-n diodes fabricated on Ga-polar and N-polar orientations before and after irradiation. The relative turn-on voltage increases for the Ga-polar diodes with increasing irradiation dose, while it increases initially and then starts to decrease for the N-polar diodes. The p-contact total resistance increases for Ga-polar and decreases for N-polar samples, which we attribute to the formation of point defects and additional Mg activation after irradiation. The J-V characteristics of most of the tested diodes recovered over time, suggesting the changes in the J-V characteristics are temporary and potentially due to metastable occupancy of traps after irradiation. X-ray photoelectron spectroscopy (XPS) and photoluminescence (PL) measurements reveal the existence of different types of initial defects and surface electronic states on Ga-polar and N-polar samples. Gallium vacancies (V Ga ) are dominant defects in Ga-polar samples, while nitrogen vacancies (V N ) are dominant in N-polar samples. The presence of a higher concentration of surface states on Ga-polar surfaces compared to N-polar was confirmed by calculating the band bending and the corresponding screening effect due to opposite polarization bound charge and ionized acceptors at the surface. The difference in surface stoichiometry in these two orientations is responsible for the different behavior in electrical characteristics after gamma-ray interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optical and mass-spectral characterization of mixed-gas flowing atmospheric-pressure afterglow sources

Plasma-based ambient desorption/ionization (ADI) sources for mass spectrometry (MS) often use helium as the primary discharge gas due to the large reaction cross section of excited helium species with atmospheric gases, which ultimately leads to efficient reagent-ion formation. However, some studies have shown that mixed-gas plasmas provide unique advantages. For instance, our group showed that a helium-oxygen flowing atmospheric-pressure afterglow (He:O 2 -FAPA) source yielded enhanced ion signals for small polar analytes, but compounds with aromatic rings underwent chemical modification to produce pyrylium species. Here, it is shown that the addition of H 2 to the helium FAPA gas significantly decreased the reagent-ion signals, but analyte-ion signal was not affected to the same degree. In addition, mass spectra obtained with He:H 2 -FAPA were chemically cleaner with fewer ions stemming from analyte species and produced less oxidation of analytes. Addition of nitrogen increased the abundance of negative reagent ions (e.g., NO 2 - and NO 3 - ) which led to enhanced ion signal for RDX. To better understand these unique advantages of a mixed-gas FAPA, the impact of molecular-gas (O 2 , N 2 , or H 2 ) addition on the optical-emission spectra from the discharge was also measured. Spatially resolved emission was obtained from the anode- and negative-glow regions. In general, addition of small fractions (~0.1%) of molecular gases decreased helium and OH emission intensities. Atomic oxygen and hydrogen emission increased with the addition of O 2 and H 2 , respectively. Furthermore, emission characteristics of N 2 , N 2 + , and NO with respect to molecular gas composition on He-FAPA were measured. Finally, fundamental plasma parameters, such as OH rotational temperature (T rot ), were also calculated for these mixed-gas systems.

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

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗