Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extrapolation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding

The optical spectra of neutral oxygen vacancies (F 0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. Here, to estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet–singlet decay pathway.

embedding↗
DOE OSTI2025

Structure and Superconductivity of Hydrogenated TiZr Alloys Under High Pressure

The effects of pressure on the structure and electrical transport properties, including electrical resistance and critical magnetic field values for superconductivity in pristine and hydrogenated TiZr alloys, were investigated. As pressure increased from ambient to 55 GPa, the structure of the pristine TiZr alloys transformed from hcp to bcc via an intermediate ω phase. The superconducting transition temperature (T c ) significantly increased from 1.76 K at ambient pressure to 15.5 K at 52.2 GPa, accompanied by an extrapolated upper critical field of 24.5 T at T = 0 K. In contrast, the hydrogenated samples exhibited substantially lower T c values compared to the pristine ones, with a T c of 4.6 K at 50.5 GPa and an extrapolated upper critical field of 1.8 T at T = 0 K. Upon decompression, the superconductivity in the hydrogenated sample persisted down to 25 GPa, disappearing when the pressure was reduced to 19 GPa. Synchrotron X-ray diffractions (XRD) analysis revealed the formation of three hydrides: I4/mmm-(Ti, Zr)H 2 , P6 3 /mmc-(Ti, Zr)H 3 , and I4̅3d-(Ti,Zr) 4 H 15 . Finally, by combining XRD and electrical resistance data, we tentatively propose that the P6 3 /mmc-(Ti, Zr)H 3 predominantly induces superconductivity. In contrast to superhydrides, where electron–phonon interactions increase the superconducting temperature at high pressure, hydrogen atoms in hcp-(Ti, Zr)H 3 reduce T c .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shallow Rate-Redox Potential Scaling in Aqueous Molecular Oxygen Reduction Electrocatalysis Across a Family of Iron Macrocycles

Rate-overpotential scaling relationships have been employed widely to understand trends in oxygen reduction reaction (ORR) electrocatalysis by dissolved metal macrocycles in organic electrolytes. Similar scaling relationships remain unknown for surface-adsorbed ORR electrocatalysts in the acidic aqueous environments germane to proton-exchange membrane (PEM) fuel cells. Herein, we examine ORR catalysis in aqueous perchloric acid media for a structurally diverse array of iron macrocycle complexes adsorbed on Vulcan carbon black. The macrocycles encompass Fe– N 4 , Fe–N 2 N' 2 and Fe–N x C 4-x motifs bearing pyrrolic, pyridinic, and N-heterocyclic carbene (NHC) moieties in the primary ligation sphere, giving rise to a 670 mV range in Fe(III/II) redox potentials, E Fe(III/II) . Experimental Tafel data in the micropolarization regime were extrapolated to the E Fe(III/II) to furnish estimated per-site-normalized current density (j per-site ) values that span ~4.6 orders of magnitude across the family of compounds. Despite the structural diversity of this family of compounds, extrapolated j per-site values correlate with the Fe(III/II) redox potentials in a roughly log-linear fashion with a shallow scaling factor of approximately 145 mV/decade. Further, these findings highlight that negative shifts in E Fe(III/II) lead to diminishing returns in catalytic rate promotion and suggest that changes to the primary ligating environment in a macrocycle are insufficient to break fundamental rate-potential scaling relationships in aqueous ORR catalysis. Together these studies motivate the further development of higher-potential iron complexes that employ motifs beyond the equatorial ligation plane to enhance ORR catalysis.

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

A near-explicit mechanistic evaluation of Isoprene photochemical seecondary organic aerosol formation and evolution: simulations of multiple chamber experiments with and without added NOx

Experimentally determined yields of secondary organic aerosol (SOA) from the photochemical oxidation of isoprene in the absence of aqueous acidic aerosol vary substantially, both within a given experiment and across different environmental chamber conditions. The underlying mechanisms driving this variation remain poorly evaluated, leading to significant uncertainty in how to extrapolate laboratory chamber results to the atmosphere. Herein, we compare SOA predictions from a near-explicit gas-phase chemical mechanism of isoprene oxidation by the hydroxyl radical (OH) in the presence and absence of nitrogen oxide radicals (NOx), to multiple chamber experiments on non-aqueous isoprene photochemical SOA (ipSOA) conducted by different groups in different chambers. SOA is predicted by volatility-driven gas-particle partitioning of hundreds of individual reaction products. The mechanism includes simplified descriptions of particle-phase organic chemistry, including organic hydroperoxide photolysis, and organic nitrate hydrolysis and accretion reactions. The model has good skill (mean normalized bias typically within 25%) at predicting the observed formation and evolution of ipSOA across a range of chambers and conditions at low NOx. The model has much less skill at describing the observed non-linear response of ipSOA to elevated NOx. Organic nitrate hydrolysis is unable to explain significant ipSOA at high NOx, whereas particle-phase accretion reactions of tertiary nitrates may play a role. Uncertainties in the chamber radical environment and fate of key organic peroxy radicals (RO2) remain as or even more important than vapor losses to chamber walls in determining how best to extrapolate chamber-based yields to the atmosphere. Implications for likely atmospheric yields of ipSOA and recommendations for future chamber experiments are discussed.

Thornton, Joel A.↗
DOE OSTI2025

On the Relationship Between Methane Production in Anaerobic Incubations of Peat Material and In Situ Methane Emissions

Anaerobic incubations of peat have been widely used to explore soil processes, but this in vitro technique raises many questions as to how well it reproduces in situ conditions. To investigate this, we conducted 60–100 days (+25 days pre‐incubation) anaerobic, temperature‐controlled incubation experiments across a temperature range of 1–26°C on samples from bog and fen habitats, at two different depths (9–19 and 25–35 cm). Here, we observed exponential increases in CO 2 and methane production with temperature in all conditions. We then compared field‐based measurements of methane emission with modeled expectations by extrapolating incubation‐determined methane production rates based on (a) soil temperature profiles, (b) the observed incubation temperature‐methane production relationship, and (c) seasonal thaw depth from each site. The resulting incubation‐extrapolated methane production agreed with measured emission rates within a factor of two at both sites and corresponded to 182 ± 54% and 59 ± 14% of the measured average yearly fluxes from the field for the bog and fen, respectively. The underestimation of fen methane fluxes may be due to the lack of living plant root‐derived dissolved organic carbon inputs in incubations, a key process in fens. Conversely, the overestimation in bogs could be attributed to methane oxidation in the field, which is absent in anaerobic incubation conditions. Nonetheless incubations predicted greenhouse gas emissions from a northern peatland within a factor of two.

Biological and medical sciences,Geosciences,Enviro↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2023

Evidence for charge delocalization crossover in the quantum critical superconductor CeRhIn5

Abstract The nature of charge degrees-of-freedom distinguishes scenarios for interpreting the character of a second order magnetic transition at zero temperature, that is, a magnetic quantum critical point (QCP). Heavy-fermion systems are prototypes of this paradigm, and in those, the relevant question is where, relative to a magnetic QCP, does the Kondo effect delocalize their f -electron degrees-of-freedom. Herein, we use pressure-dependent Hall measurements to identify a finite-temperature scale E loc that signals a crossover from f -localized to f -delocalized character. As a function of pressure, E loc ( P ) extrapolates smoothly to zero temperature at the antiferromagnetic QCP of CeRhIn 5 where its Fermi surface reconstructs, hallmarks of Kondo-breakdown criticality that generates critical magnetic and charge fluctuations. In 4.4% Sn-doped CeRhIn 5 , however, E loc ( P ) extrapolates into its magnetically ordered phase and is decoupled from the pressure-induced magnetic QCP, which implies a spin-density-wave (SDW) type of criticality that produces only critical fluctuations of the SDW order parameter. Our results demonstrate the importance of experimentally determining E loc to characterize quantum criticality and the associated consequences for understanding the pairing mechanism of superconductivity that reaches a maximum T c in both materials at their respective magnetic QCP.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measurement of the proton spin structure at long distances

Measuring the spin structure of protons and neutrons tests our understanding of how they arise from quarks and gluons, the fundamental building blocks of nuclear matter. At long distances, the coupling constant of the strong interaction becomes large, requiring non-perturbative methods to calculate quantum chromodynamics processes, such as lattice gauge theory or effective field theories. In this work, we report proton spin structure measurements from scattering a polarized electron beam off polarized protons. The spin-dependent cross-sections were measured at large distances, corresponding to the region of low momentum transfer squared between 0.012 and 1.0GeV 2 . This kinematic range provides unique tests of chiral effective field theory predictions. Our results show that a complete description of the nucleon spin remains elusive, and call for further theoretical works, for example, in lattice quantum chromodynamics. Finally, our data extrapolated to the photon point agree with the Gerasimov–Drell–Hearn sum rule, a fundamental prediction of quantum field theory that relates the anomalous magnetic moment of the proton to its integrated spin-dependent cross-sections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

J-PLUS: Support vector machine applied to STAR-GALAXY-QSO classification

Context. In modern astronomy, machine learning has proved to be efficient and effective in mining big data from the newest telescopes. Aims. In this study, we construct a supervised machine-learning algorithm to classify the objects in the Javalambre Photometric Local Universe Survey first data release (J-PLUS DR1). Methods. The sample set is featured with 12-waveband photometry and labeled with spectrum-based catalogs, including Sloan Digital Sky Survey spectroscopic data, the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, and VERONCAT – the Veron Catalog of Quasars & AGN. The performance of the classifier is presented with the applications of blind test validations based on RAdial Velocity Extension, the Kepler Input Catalog, the Two Micron All Sky Survey Redshift Survey, and the UV-bright Quasar Survey. A new algorithm was applied to constrain the potential extrapolation that could decrease the performance of the machine-learning classifier. Results. The accuracies of the classifier are 96.5% in the blind test and 97.0% in training cross-validation. The F1-scores for each class are presented to show the balance between the precision and the recall of the classifier. We also discuss different methods to constrain the potential extrapolation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On Alfvénic turbulence of solar wind streams observed by Solar Orbiter during March 2022 perihelion and their source regions

It has been recently accepted that the standard classification of the solar wind solely according to flow speed is outdated, and particular interest has been devoted to the study of the origin and evolution of so-called Alfvénic slow solar wind streams and to what extent such streams resemble or differ from fast wind. In March 2022, Solar Orbiter completed its first nominal phase perihelion passage. During this interval, it observed several Alfvénic streams, allowing for characterization of fluctuations in three slow wind intervals (AS1-AS3) and comparison with a fast wind stream (F) at almost the same heliocentric distance. This work makes use of Solar Orbiter plasma parameters from the Solar Wind Analyzer (SWA) and magnetic field measurements from the magnetometer (MAG). The magnetic connectivity to the solar sources of selected solar wind intervals was reconstructed using a ballistic extrapolation based on measured solar wind speed down to the (spherical) source surface at 2.5 R s below which a potential field extrapolation was used to map back to the Sun. The source regions were identified using SDO/AIA observations. A spectral analysis of in situ measured magnetic field and velocity fluctuations was performed to characterize correlations, Alfvénicity, normalized cross-helicity, and residual energy in the frequency domain as well as intermittency of the fluctuations and spectral energy transfer rate estimated via mixed third-order moments. A machine learning technique was used to separate proton core, proton beam, and alpha particles and to study v − b correlations for the different ion populations in order to evaluate the role played by each population in determining the Alfvénic content of solar wind fluctuations. The comparison between fast wind and Alfvénic slow wind intervals highlights the differences between the two solar wind regimes: The fast wind is characterized by larger amplitude fluctuations, and magnetic and velocity fluctuations are closer to equipartition of energy. In fact the Alfvénic slow wind streams appear to be on a spectrum of wind types, with AS1, originating from open field lines neighboring active regions and displaying similarities with the fast wind in terms of fluctuation amplitude and turbulence characteristics, but not with respect to the alpha particles and proton beams. The other two slow streams differed both in their sources as well as plasma characteristics, with AS2 coming from the expansion of a narrow coronal hole corridor and AS3 from a region straddling a pseudostreamer. The latter displayed the coldest and highest density but the slowest stream with the smallest fluctuation amplitude and greatest magnetic energy excess. It also showed the largest scatter in proton beam speeds and the greatest difference in speed between proton beam and alpha particles. This study shows how the old fast–slow solar wind dichotomy, already called into question by the observations of slower Alfvénic solar wind streams, should further be refined, as the Alfvénic slow wind, originating in different solar wind regions, show significant differences in density, temperature, and proton and alpha-particle properties in the inner heliosphere. The observations presented here provide the starting point for a better understanding of the origin and evolution of different solar wind streams as well as the evolving turbulence contained within.

magnetohydrodynamics (MHD)↗
DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

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

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

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

Role of microtearing mode in DIII-D and future high- β p core plasmas

This paper provides modeling evidence for a new variant of microturbulence instability that regulates the electron thermal fluxes in the transport barrier region of DIII-D "high β p " discharges. This scenario is being investigated as a possible scenario for future steady-state reactors, and thus it is essential to understand the physics processes which govern its performance in order to confidently extrapolate to future devices. Additional analysis examining the extrapolation of this scenario to ITER shows that while the new microtearing mode identified here is unlikely to be unstable in that regime (due to increased temperatures and reduced collisionality), the large shift in magnetic axis (Shafranov shift) that arises in these plasmas may still favor stabilization of common ballooning instabilities which peak on the low field side, but destabilize electrostatic slab modes on the high field side.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Photoelectric absorption cross section of silicon near the bandgap from room temperature to sub-Kelvin temperature

The use of cryogenic silicon as a detector medium for dark matter searches is gaining popularity. Many of these searches are highly dependent on the value of the photoelectric absorption cross section of silicon at low temperatures, particularly near the silicon bandgap energy, where the searches are most sensitive to low mass dark matter candidates. While such cross section data have been lacking from the literature, previous dark matter search experiments have attempted to estimate this parameter by extrapolating it from higher temperature data. However, discrepancies in the high temperature data have led to order-of-magnitude differences in the extrapolations. In this paper, we resolve these discrepancies by using a novel technique to make a direct, low temperature measurement of the photoelectric absorption cross section of silicon at energies near the bandgap (1.2 eV–2.8 eV).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗
DOE OSTI2022

Time-domain self-broadened and air-broadened nitrogen S-branch Raman linewidths at 80–200 K recorded in an underexpanded jet

We report pure-rotational N 2 –N 2 , N 2 –air, and O 2 –air S-branch linewidths for temperatures of 80–200 K by measuring the time-dependent decay of rotational Raman coherences in an isentropic free-jet expansion from a sonic nozzle. We recorded pure-rotational hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) spectra along the axial centerline of the underexpanded jet, within the barrel shock region upstream of the Mach disk. The dephasing of the pure-rotational Raman coherence was monitored using probe-time-delay scans at different axial positions in the jet, corresponding to varying local temperatures and pressures. The local temperature was obtained by fitting CARS spectra acquired at zero probe time delay, where the impact of collisions was minimal. The measured decay of each available Raman transition was fit to a dephasing constant and corrected for the local pressure, which was obtained from the CARS-measured static temperature and thermodynamic relationships for isentropic expansion from the known stagnation state. Nitrogen self-broadened transitions decayed more rapidly than those broadened in air for all temperatures, corresponding to higher Raman linewidths. In general, the measured S-branch linewidths deviated significantly in absolute and relative magnitudes from those predicted by extrapolating the modified exponential gap model to low temperatures. The temperature dependence of the Raman linewidth for each measured rotational state of nitrogen (J ≤ 10) and oxygen (N ≤ 11) was fit to a temperature-dependent power law over the measurable temperature domain (80–200 K) and extrapolated to both higher rotational states and room temperature. Furthermore, the measured and modeled low-temperature linewidth data provided here will aid low temperature gas-phase pressure measurements with fs/ps CARS.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗