Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained density functional theory”

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 145 records · Page 8

Enhanced orbital anisotropy through the proximity to a SrTiO 3 layer in the perovskite iridate superlattices

We have used angle-dependent soft x-ray absorption spectroscopy (XAS) at the O K edge and first-principles calculations to investigate the electronic structures of iridate-based superlattices ( SrIr O 3 ) m / ( SrTi O 3 ) ( m = 1 , 2, 3, and ∞). We focus on the pre-edge Ir 5 d t 2 g – O 2 p orbital hybridization feature in the XAS spectra. By varying the measurement geometry relative to the incident photon polarization, we are able to extract the dichroic contrast and observe the systematic increase in the anisotropy of Ir 5 d orbitals as m decreases. First-principles calculations elucidate the orbital anisotropy coming mainly from the enhanced out-of-plane compression of Ir O 6 octahedra in the SrIr O 3 layers that are adjacent to the inserted SrTi O 3 layers. As m decreases, the increased volume fraction of these interfacial SrIr O 3 layers and their contact with the SrTi O 3 layers within the ( SrIr O 3 ) m / ( SrTi O 3 ) supercell lead to enhanced orbital anisotropy. Furthermore, the tilt and rotation of Ir O 6 octahedra are shown to be essential to understand the subtle orbital anisotropy in these superlattices, and constraining these degrees of freedom will give an incorrect trend. Overall, our results demonstrate that the structural constraint from the inserted SrTi O 3 layers, in addition to other electronic means such as polar interface and charge transfer, can serve as a knob to control the orbital degree of freedom in these iridate-based superlattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From chiral effective field theory to perturbative QCD: A Bayesian model mixing approach to symmetric nuclear matter

Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory (𝜒⁢EFT ) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars, while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with quantified theoretical uncertainties. However, describing the full range of densities in between with a single EOS that has well-quantified uncertainties is a challenging problem. Bayesian multimodel inference from 𝜒⁢EFT and pQCD can help bridge the gap between the two theories. In this work, we introduce a correlated Bayesian model mixing framework that uses a Gaussian process (GP) to assimilate different information into a single QCD EOS for symmetric nuclear matter. The present implementation uses a stationary GP to infer this mixed EOS solely from the EOSs of 𝜒⁢EFT and pQCD while accounting for the truncation errors of each theory. The GP is trained on the pressure as a function of number density in the low- and high-density regions where 𝜒⁢EFT and pQCD are, respectively, valid. We impose priors on the GP kernel hyperparameters to suppress unphysical correlations between these regimes. This, together with the assumption of stationarity, results in smooth 𝜒⁢EFT-to-pQCD curves for both the pressure and the speed of sound. We show that using uncorrelated mixing requires uncontrolled extrapolation of at least one of 𝜒⁢EFT or pQCD into regions where the perturbative series breaks down and leads to an acausal EOS. Here, we also discuss extensions of this framework to nonstationary and less differentiable GP kernels, its future application to neutron-star matter, and the incorporation of additional constraints from nuclear theory, experiment, and multimessenger astronomy.

Bayesian methods↗

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search for beyond-mean-field signatures in heavy-ion fusion reactions

Examination of high-resolution, experimental fusion excitation functions for 16,17,18 O + 12 C reveals a remarkable irregular behavior that is rooted in the structure of both the colliding nuclei and the quasi-molecular composite system. The impact of the l-dependent fusion barriers is assessed using a time-dependent Hartree-Fock model, viewed as a mean-field reference. Barrier penetrabilities, taken directly from a density-constrained calculation, provide a significantly improved description of the experimental data as compared to the standard Hill-Wheeler approach. Furthermore, the remaining deviations between the parameter-free theoretical mean-field predictions and experimental fusion cross sections are exposed and discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Highly accurate and constrained density functional obtained with differentiable programming

Using an end-to-end differentiable implementation of the Kohn-Sham self-consistent field equations, we obtain a highly accurate neural network–based exchange and correlation (XC) functional of the electronic density. The functional is optimized using information on both energy and density while exact constraints are enforced through an appropriate neural network architecture. Here we evaluate our model against different families of XC approximations and show that at the meta-GGA level our functional exhibits unprecedented accuracy for both energy and density predictions. For nonempirical functionals, there is a strong linear correlation between energy and density errors. We use this correlation to define an XC functional quality metric that includes both energy and density errors, leading to an improved way to rank different approximations.

36 MATERIALS SCIENCE↗

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of Doppler back-scattering and charge exchange measurements of E × B plasma rotation in the DIII-D tokamak under varying torque conditions

Measurements of the E × B toroidal angular velocity, $\omega_{E \times B}\,=\,E_\mathrm{r}/RB_\theta$ ($E_\mathrm{r}$ is the radial electric field, B θ is the poloidal magnetic field), are made using the Doppler back-scattering (DBS) and charge-exchange recombination (CER) spectroscopy diagnostics. DBS uses the Doppler shift of wavenumber-resolved density fluctuations while CER uses the Doppler shift of impurity emission lines to independently measure plasma parameters for calculating the local radial electric field. DBS and CER profiles of $\omega_{E \times B}$ as a function of normalized toroidal flux (ρ) are compared at various levels of neutral beam applied torque on the plasma. Under standard neoclassical theory $\omega_{E\times B}$ is a flux surface quantity, making it appropriate to compare across diagnostics. DBS and CER generally show good agreement when comparing $\omega_{E \times B}$ profiles at different levels of neutral beam injection-applied torque. Furthermore, the DBS values have close to the same precision as CER values when averaged over a similar time-scale and effects, such as prompt-torque are considered. DBS is able to observe the rapid ($\lt$10 ms) modification of the $E_\mathrm{r}$ profile by the diagnostic neutral beam 'blips'. This modification is most pronounced when the blip applies a large relative change in torque on the plasma. Overall, these results could have implications on transport analysis and suggests using DBS and CER in conjunction to constrain values of the E × B-shear (sometimes called $\gamma_{E \times B}$).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian uncertainty quantification for nuclear matter incompressibility

Within a Bayesian statistical framework using the standard Skyrme-Hartree-Fock model, the maximum a posteriori (MAP) values and uncertainties of nuclear matter incompressibility and isovector interaction parameters are inferred from the experimental data of giant resonances and neutron-skin thicknesses of typical heavy nuclei. Furthermore, with the uncertainties of the isovector interaction parameters constrained by the data of the isovector giant dipole resonance and the neutron-skin thickness, we have obtained K 0 = 223$^{+7}_{–8}$ MeV at 68% confidence level using the data of the isoscalar giant monopole resonance in 208 Pb measured at the Research Center for Nuclear Physics (RCNP), Japan, and at the Texas A&M University, USA. Although the corresponding 120 Sn data gives a MAP value for K 0 about 5 MeV smaller than the 208 Pb data, there are significant overlaps in their posterior probability distribution functions.

190 ≤ A ≤ 219↗

Highly dense atomic Fe–Ni dual metal sites for efficient CO 2 to CO electrolyzers at industrial current densities

Carbon-supported, atomically dispersed, nitrogen-coordinated metal sites (e.g., Fe and Ni) are arguably the most promising catalysts for the electrochemical reduction of CO 2 to CO due to their unique catalytic properties and the use of earth-abundant elements. However, conventional single metal sites are constrained by their structural simplicity, causing either too weak or too strong absorption/desorption of multiple critical intermediates (e.g., *COOH and *CO). Current catalysts also suffer from ultra-low loadings (<1.0 wt%) of atomic metal active sites in catalysts, leading to inadequate performance for CO 2 -to-CO conversion. Here, we develop dual Ni/Fe metal site catalysts with significantly increased atomically dispersed metal loadings (up to 4.8 wt%). A gas-phase chemical vapor deposition (CVD) approach to introducing single Ni sites was integrated with Fe 2 O 3 /ZIF-8 precursors, followed by an optimal thermal activation. The optimized CVD-Ni/Fe–N–C catalyst exhibited remarkable electrocatalytic performance for the CO 2 reduction to CO in a continuous membrane-electrode-assembly electrolyzer, achieving a maximum CO faradaic efficiency (FE CO ) of 96% at a current density of 700 mA cm −2 in a near-neutral electrolyte. Furthermore, a desirable but challenging acidic flow-cell electrolyzer was designed using this dual metal site catalyst to improve CO 2 utilization, accomplishing a FE CO of up to 95% at a CO partial current density close to 600 mA cm −2 . Density functional theory (DFT) calculations suggest a synergetic effect between Fe–Ni pairs facilitating *COOH intermediate formation and *CO desorption simultaneously during CO 2 to CO conversion. This is key to breaking the linear scaling relationship of conventional single-metal site catalysts during the CO 2 reduction reaction.

36 MATERIALS SCIENCE↗

Learning thermodynamically constrained equations of state with uncertainty

Numerical simulations of high energy-density experiments require equation of state (EOS) models that relate a material’s thermodynamic state variables—specifically pressure, volume/density, energy, and temperature. EOS models are typically constructed using a semi-empirical parametric methodology, which assumes a physics-informed functional form with many tunable parameters calibrated using experimental/simulation data. Since there are inherent uncertainties in the calibration data (parametric uncertainty) and the assumed functional EOS form (model uncertainty), it is essential to perform uncertainty quantification (UQ) to improve confidence in EOS predictions. Model uncertainty is challenging for UQ studies since it requires exploring the space of all possible physically consistent functional forms. Thus, it is often neglected in favor of parametric uncertainty, which is easier to quantify without violating thermodynamic laws. This work presents a data-driven machine learning approach to constructing EOS models that naturally captures model uncertainty while satisfying the necessary thermodynamic consistency and stability constraints. We propose a novel framework based on physics-informed Gaussian process regression (GPR) that automatically captures total uncertainty in the EOS and can be jointly trained on both simulation and experimental data sources. A GPR model for the shock Hugoniot is derived, and its uncertainties are quantified using the proposed framework. We apply the proposed model to learn the EOS for the diamond solid state of carbon using both density functional theory data and experimental shock Hugoniot data to train the model and show that the prediction uncertainty is reduced by considering thermodynamic constraints.

Sharma, Himanshu (ORCID:000900050235718X)↗

Correlations between triaxial shapes and formation history of dark matter haloes

ABSTRACT The shape of dark matter (DM) haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in DM haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the normalized peak-centroid offset. These correlations are nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multidimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-sized haloes. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring and enhancing the accuracy of interior-scaled Perdew–Zunger self-interaction correction

The Perdew–Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange–correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z σ , which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional’s energy density. Here, the resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear and hybrid equations of state in light of the low-mass compact star in HESS J1731-347

Here, we sample over 500 relativistic mean-field theories constrained by chiral effective-field theory and properties of isospin-symmetric nuclear matter and test them against known stellar structure constraints. This includes a recent mass and radius measurement of a compact object in supernova remnant HESS J1731-347, with an unusually low mass of M = 0.77$^{+0.20}_{–0.17}$M ⊙ and a compact radius of R = 10.4$^{+0.86}_{–0.78}$ km. We show that none of the sampled nuclear models meet all constraints at the 68% credibility level, but that hybrid equations of state with a quark-matter inner core and nuclear outer core easily can. This indicates a tension between astrophysical constraints and low-energy nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Synergistic double-atom catalysts of metal-boron anchored on g-C 2 N for electrochemical nitrogen reduction: Mechanistic insight and catalyst screening

The rational design of a novel catalytic center with a sound basis remains both challenging and rewarding for the electrochemical reduction of N 2 (eNRR), which has provided a feasible route for achieving clean and sustainable NH 3 production under ambient conditions. Herein, using density functional theory calculations, we demonstrate that hybrid metal (M)-boron (B) double-atom catalysts (DACs) embedded in g-C 2 N substrate (M-B@C 2 N, M = 3d, 4d and 5d transition metals) can achieve both high catalytic activity and high selectivity in eNRR. The proposed M-B@C 2 N DACs have exhibited impressive feasibility and stability thanks to the resilient and robust C 2 N substrate with abundant pyridinic N atoms distributed among right-sized pore structures. Our results reveal that like the metal center, the embedded B atom can actively involve in Ntriple bondN bond activation via π*-backdonation mechanism concomitant with the substantial charge transfer to adsorbed *N 2 , leading to sizable Nsingle bondN bond elongation. Accordingly, both adsorption energy and Nsingle bondN bond length of *N 2 can be employed as catalytic descriptors for predicting eNRR activity in terms of the limiting potentials (UL). Using high-throughput screening method, we found that six M-B@C 2 N candidates have stood out as the outstanding electrocatalysts for driving eNRR, namely, M = Ti (U L = 0 V), Mo (U L = 0 V), Nb (U L = -0.04 V), W (U L = -0.23 V), Zr (U L = -0.26 V), V (U L = -0.28 V). The underlying origin is attributed to the balanced and constrained N-affinity of M-B dual site working in synergy, which can thus be used as one important guide of catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced electric propulsion and space plasma contactor research

A theory of the plasma contacting process is described and experimental results obtained using three different hollow cathode-based plasma contactors are presented. The existence of a sheath across which the bulk of the voltage drop associated with the contacting process occurs is demonstrated. Test results are shown to agree with a model of a spherical, space-charge-limited double sheath. The concept of ignited mode contactor operation is discussed, which is shown to enhance contactor efficiency when it is collecting electrons. An investigation of the potentials in the plasma plumes downstream of contactors operating at typical conditions is presented. Results of tests performed on hollow cathodes operating at high interelectrode pressures (up to about 1000 Torr) on ammonia are presented and criteria that are necessary to ensure that the cathode will operate properly in this regime are presented. These results suggest that high pressure hollow cathode operation is difficult to achieve and that special care must be taken to assure that the electron emission region remains diffuse and attached to the low work function insert. Experiments conducted to verify results obtained previously using a ring cusp ion source equipped with a moveable anode are described and test results are reported. A theoretical study of hollow cathode operation at high electron emission currents is presented. Preliminary experiments using the constrained sheath optics concept to achieve ion extraction under conditions of high beam current density, low net accelerating voltage and well columniated beamlet formation are discussed.

Wilbur, Paul J.↗

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗

Platform for probing radiation transport properties of hydrogen at conditions found in the deep interiors of red dwarfs

We describe an experimental concept at the National Ignition Facility for specifically tailored spherical implosions to compress hydrogen to extreme densities (up to ∼800× solid density, electron number density ne∼4×1025 cm−3) at moderate temperatures (T∼200 eV), i.e., to conditions, which are relevant to the interiors of red dwarf stars. The dense plasma will be probed by laser-generated x-ray radiation of different photon energy to determine the plasma opacity due to collisional (free–free) absorption and Thomson scattering. The obtained results will benchmark radiation transport models, which in the case for free–free absorption show strong deviations at conditions relevant to red dwarfs. This very first experimental test of free–free opacity models at these extreme states will help to constrain where inside those celestial objects energy transport is dominated by radiation or convection. Moreover, our study will inform models for other important processes in dense plasmas, which are based on electron–ion collisions, e.g., stopping of swift ions or electron–ion temperature relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constraints on the Energy Content of the Universe from a Combination of Galaxy Cluster Observables

We demonstrate that constraints on cosmological parameters from the distribution of clusters as a function of redshift (dN/dz) are complementary to accurate angular diameter distance (D(sub A)) measurements to clusters, and their combination significantly tightens constraints on the energy density content of the Universe. The number counts can be obtained from X-ray and/or SZ (Sunyaev-Ze'dovich effect) surveys, and the angular diameter distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung X-ray emission and the SZ effect. We combine constraints from simulated cluster number counts expected from a 12 deg(sup 2) SZ cluster survey and constraints from simulated angular diameter distance measurements based on the X-ray/SZ method assuming a statistical accuracy of 10% in the angular diameter distance determination of 100 clusters with redshifts less than 1.5. We find that Omega(sub m), can be determined within about 25%, Omega(sub lambda) within 20% and w within 16%. We show that combined dN/dz+(sub lambda) constraints can be used to constrain the different energy densities in the Universe even in the presence of a few percent redshift dependent systematic error in D(sub lambda). We also address the question of how best to select clusters of galaxies for accurate diameter distance determinations. We show that the joint dN/dz+ D(lambda) constraints on cosmological parameters for a fixed target accuracy in the energy density parameters are optimized by selecting clusters with redshift upper cut-offs in the range 0.55 approx. less than 1. Subject headings: cosmological parameters - cosmology: theory - galaxies:clusters: general

Molnar, Sandor M.↗