Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “couple cluster approximation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Large field-of-view event-mode camera for high-precision epithermal neutron resonance imaging

A large-area event-mode camera system coupled with a 6 LiF-ZnS:Ag scintillator is applied for neutron resonance imaging (NRI) on the energy-resolved neutron imaging (ERNI) flight path, also known as Flight Path 5 (FP5), at the Los Alamos Neutron Science Center (LANSCE). This novel neutron imaging system, featuring a 120 x 120 mm 2 field of view, efficiently captures resonance information across the entire image in a single acquisition, significantly reducing beam time requirements compared to conventional energy-resolved neutron imaging systems. High-quality neutron radiographs with enhanced spatial resolution are achieved through the reconstruction of neutron events based on observations of individual photons emitted from the scintillator. The system demonstrates reduced background through neutron/gamma discrimination capabilities while maintaining sharpness across a large fields of view. In the measurements presented here, a spatial resolution of approximately 340 μm was achieved using center-of-gravity photon cluster centroiding. We demonstrate the system’s capability for quantitatively determining isotopic distributions in various thin samples, as well as automatically reconstructing complex scenes with overlapping resonances from diverse samples. These results are obtained using standard data analysis tools, despite the relatively slow 6 LiF-ZnS:Ag scintillator, which may not be optimal for absorption resonance detection. The capabilities demonstrated here offer a valuable, versatile, and cost-effective solution for high spatial and temporal resolution, large field-of-view energy-resolved neutron imaging, with potential applications across various scientific and industrial domains.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal–organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule–MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.

Chemistry↗
NASA NTRS1994

Far-ultraviolet photometry of the globular cluster omega Cen

We present far-ultraviolet images of the globular cluster omega Centauri obtained with the Ultraviolet Imaging Telescope (UIT) during the 1990 December Astro-1 mission. A total of 1957 sources are detected at 1620 A to a limiting ultraviolet (UV) magnitude of 16.4 in the central 24 min diameter region of the field and a limit of 15.6 over the remainder of the 40 min diameter field. Over 1400 of these sources are matched with stars on a Stroemgren u band charge coupled devices (CCD) frame obtained with the Cerro Tololo Inter-American Observatory (CTIO) 0.9 m telescope to produce a (far-UV, u) color-magnitude diagram (CMD). Completeness of the sample and error estimates are determined by photometry of artificial stars added to the images. The horizontal branch (HB) of the CMD is heavily populated hotter than 9000 K. A large number of 'extreme HB' stars are found hotter than a conspicuous break in the HB at T(sub e) approximately 16000 K. There is also a significant population of stars above the HB, the brightest of which is 4 mag brighter than the HB. Most of the hotter of these appear to be 'AGB-manque' or 'Post-Early Asymptotic Giant Branch' stars. We compare the observations to recent theoretical evolutionary tracks for the zero-age HB and subsequent phases. The tracks match the data well, with the exception of the hotter HB stars, many of which fall below the zero-age horizontal branch. It is unclear as yet whether these are a special population or an artifact of errors in the models or photometry. We identify 33 stars with T(sub e) greater than or approximately = 50000 K, which are hotter than zero-age HB stars with envelope masses of 0.003 solar mass.

Whitney, Jonathan H.↗
NASA NTRS1997

The High Throughput X-ray Spectroscopy (HTXS) Mission

The HTXS mission concept combines large effective area (approximately 15,000 sq cm at 1 keV), high spectral resolution (E/Delta(E) approximately 300-3000), and broad energy bandpass (0.25-40 keV and possibly up to 100 keV) by using replicated optics together with a complement of spectroscopic instrumentation including reflection gratings readout by charge-coupled device detectors (CCDs), quantum micro-calorimeters, and cadmium zinc telluride (CZT) or comparable high energy detectors. An essential feature of this concept involves minimization of cost (approximately $350M for development and approximately $500-600M including launches) and risk by building six identical modest satellites to achieve the large area. Current mission and technology studies are targeted towards a new start in the 2002 timeframe, with first launch around 2005-2006. The HTXS mission represents a major advance, providing as much as a factor of 100 increase in sensitivity over currently planned high resolution X ray spectroscopy missions. HTXS will mark the start of a new era when high quality X ray spectra will be obtained for all classes of X ray sources, over a wide range of luminosity and distance. With its increased capabilities, HTXS will address many fundamental astrophysics questions such as the origin and distribution of the elements from carbon to zinc, the formation and evolution of clusters of galaxies, the validity of general relativity in the strong gravity limit, the evolution of supermassive black holes in active galactic nuclei, the details of supernova explosions and their aftermath, and the mechanisms involved in the heating of stellar coronae and driving of stellar winds.

White, N. E.↗
DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

UBV stellar photometry of the 30 Doradus region of the large Magellanic Cloud with the Hubble Space Telescope

We report on Planetary Camera observations of the central region of 30 Doradus in the Large Magellanic Cloud (LMC). These images of 30 Doradus are the first `deep' Hubble Space Telescope (HST) exposures that have appropriate photometric calibration. The B band (F439W) image, which shows R136a at the center of the PC6 charge coupled device (CCD) chip, reveals over 200 stars within 3 sec of the center of R13a, and over 800 stars in a 35 sec x 35 sec area. We used Malumuth et al.'s (1991) Point Spread Function (PSF)-fitting method to measure the magnitudes of all stars on the PC6 chip. These new B magnitudes, along with U and V magnitudes from archival PC images, yield a luminosity function, mass density profile, and initial mass function of the 30 Doradus ionizing cluster. The mass distribution is well fit by a King model with a core radius, R(sub c) = 0.96 sec (0.24 pc), a tidal radius, R(sub t) = 110 sec (28 pc), and a total mass, Mass = 16,800 solar mass. Both the luminosity function and initial mass function show evidence for mass segregation, in the sense that the central region has a higher fraction of massive stars than the outer regions. This is the first observational evidence for mass segregation in a very young cluster (age approximately 3 million years). The observations admit the hypothesis that the mass segregation occurred in the process of star formation and/or that the mass segregation is the result of dynamical evolution.

Malumuth, Eliot M.↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2020

Size- and temperature-dependent magnetization of iron nanoclusters

In this paper, the magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics simulations coupled to molecular dynamics, using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M (T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M (T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering. Analytical approximations to the magnetization as a function of temperature and size are proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Radial Segregation between Red Giant Populations in the Globular Cluster NGC 3201: Tentative Role of Varying Mass among the Progenitors

In the globular cluster NGC 3201, we study subpopulations (SPs) of red giants (RGs) distinguished by their characteristics. In its central part, we make use of multicolor Hubble Space Telescope (HST) photometry coupled with spectroscopy of RGs spread out over a larger cluster face. In the color index CI = (m{sub F275W}−m{sub F336W}) −(m{sub F336W}−m{sub F439W}), RGs commonly referred to as first generation, SPI, are separated from the counterparts of the second generation, SPII. The latter are more concentrated than the former at a confidence level of P = 93%, within the field covered by the HST photometry. Their minor portion (∼28%, SPIIm) located in the middle of the δCI range spanned across the split RG branch relative to the splitting edge is more centrally concentrated at a higher confidence level of 96%. We found, from the [Na/Fe]−δCI plot for RGs in common between the photometry and spectroscopy, that SPI and SPII are mostly overlapped within 0.0 < [Na/Fe] < 0.3. RGs with [Na/Fe] < 0.0 (mostly SPI) keep their lower concentration to the cluster center over larger radial distance as compared to SPII RGs with [Na/Fe] > 0.3 at a confidence level of P = 88%. It is increased up to 95% if a restriction on [O/Fe] is imposed. The obtained result, particularly in the central part, is a clue that SPIIm RGs are the descendants of more massive progenitors. Of these RGs approximately one-third probably originate from recent (within the past 1.5 Gyr) blue stragglers with M {sub BS} > 1.1M {sub ☉}, and the other two-thirds could be accounted for by older collision products of primordial MS stars.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Open-Shell Tensor Hypercontraction

The extension of least-squares tensor hypercontracted second- and third-order Møller–Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller–Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Herein we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

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

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
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 OSTI2025

Insight Into the High-Potential Branch of the Alternate Nad+-Dependent NADPH:Ferredoxin Oxidoreductase II (NfnII) from Pyrococcus Furiosus. Evidence for the Gating Step in Electron Bifurcation

We have investigated the rapid-reaction kinetics of the NAD+-dependent NADPH:ferredoxin oxidoreductase II (NfnII) from Pyrococcus furiosus, permitting a comparison with recent work done with the paralog NfnI from the same organism. The half-potentials of the electron-bifurcating L-FAD are highly crossed in both NfnI, meaning the potential for the quinone/semiquinone couple is significantly lower than that for the semiquinone/hydroquinone couple so that the semiquinone oxidation state is thermodynamically unstable. The same appears to be the case with NfnII on the basis of its similar behavior in transient absorption spectroscopy experiments and the absence of any evidence for FAD*- accumulation in the course of reductive titrations (which would be manifested as a transient increase in absorbance at ~380 nm)1. Reductive titrations with both the one-electron donor sodium dithionite and the obligate two-electron donor NADPH demonstrate that, in contrast to NfnI, little FADH* accumulates in the high-potential pathway of NfnII in the course of reduction, as reflected in the absence of a transient increase in absorbance in the 500-600 nm region. This indicates that the half-potentials of the S-FAD are crossed in NfnII by a minimum of 120 mV. Rapid-reaction experiments mixing oxidized NfnII with NADPH also show no evidence of S-FADH* accumulation. Furthermore, the enzyme is only partially reduced at the end of the reaction with NADPH, indicating that there is little electron transfer into the high-potential pathway of NfnII. When the reaction is carried out in the presence of the one-electron carrier ferredoxin, only approximately one equivalent of ferredoxin becomes reduced, in contrast to the maximum of three equivalents seen with NfnI. This is consistent with only a single electron bifurcation event taking place under these conditions with NfnII, with only a single high-potential electron passing into the [2Fe-2S] cluster of the high-potential pathway of NfnII but not progressing further to the S-FAD. This accounts for the observation that NfnII has greatly reduced bifurcating activity compared to NfnI. The lack of electron transfer into the S-FAD of NfnII due to its crossed half-potentials prevents successive bifurcation activity and indicates that electron transfer within the high-potential branch gates bifurcation activity in NfnI.

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

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗