Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density correlation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Reproducibility of protein x-ray diffuse scattering and potential utility for modeling atomic displacement parameters

Protein structure and dynamics can be probed using x-ray crystallography. Whereas the Bragg peaks are only sensitive to the average unit-cell electron density, the signal between the Bragg peaks—diffuse scattering—is sensitive to spatial correlations in electron-density variations. Although diffuse scattering contains valuable information about protein dynamics, the diffuse signal is more difficult to isolate from the background compared to the Bragg signal, and the reproducibility of diffuse signal is not yet well understood. We present a systematic study of the reproducibility of diffuse scattering from isocyanide hydratase in three different protein forms. Both replicate diffuse datasets and datasets obtained from different mutants were similar in pairwise comparisons (Pearson correlation coefficient ≥0.8). The data were processed in a manner inspired by previously published methods using custom software with modular design, enabling us to perform an analysis of various data processing choices to determine how to obtain the highest quality data as assessed using unbiased measures of symmetry and reproducibility. The diffuse data were then used to characterize atomic mobility using a liquid-like motions (LLM) model. This characterization was able to discriminate between distinct anisotropic atomic displacement parameter (ADP) models arising from different anisotropic scaling choices that agreed comparably with the Bragg data. Our results emphasize the importance of data reproducibility as a model-free measure of diffuse data quality, illustrate the ability of LLM analysis of diffuse scattering to select among alternative ADP models, and offer insights into the design of successful diffuse scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

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

Impact of the equation of state on f - and p - mode oscillations of neutron stars

We investigate the impact of the neutron-star matter equation of state on the f- and p1-mode oscillations of neutron stars obtained within the Cowling approximation and linearized general relativity. The f- and p 1 -mode oscillation frequencies, and their damping times are calculated using representative sets of Skyrme Hartree-Fock and relativistic mean-field models, all of which reproduce nuclear systematics and support 2M ⊙ neutron stars. Our study shows strong correlations between the frequencies of f- and p 1 -modes and their damping times with the pressure of β-equilibrated matter at densities equal to or slightly higher than the nuclear saturation density ρ 0 . Such correlations are found to be almost independent of the composition of the stars. The frequency of the p 1 -mode of 1.4M ⊙ star is strongly correlated with the slope of the symmetry energy L 0 and β-equilibrated pressure at density ρ 0 . Compared to GR calculations, the error in the Cowling approximation for the f-mode is about 30% for neutron stars of low mass, whereas it decreases with increasing mass. Here, the accuracy of the p 1 -mode is better than 15% for neutron stars of maximum mass, and improves for lower masses and higher number of radial nodes. I.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An analytic hybrid halo + perturbation theory model for small-scale correlators: baryons, halos, and galaxies

In this work, we update Halo Zeldovich Perturbation Theory (HZPT, [1]), an analytic model for the two-point statistics of dark matter, to describe halo and galaxy clustering, and galaxy-matter cross-correlation on nonlinear scales. The model correcting Zeldovich has an analytic Fourier transform, and therefore is valid in both configuration space and Fourier space. The model is accurate at the 2%-level or less for P_mm (k< 1 h/Mpc), P_hm (k< 1 h/Mpc), P_hh (k< 2 h/Mpc), P gm (k< 1 h/Mpc), P gg (k< 1 h/Mpc), ξ mm (r> 1 Mpc/h), ξ hm (r> 2 Mpc/h), ξ hh (r> 2 Mpc/h), ξ gm (r> 1 Mpc/h), ξ gg (r> 2 Mpc/h), for LRG-like mock galaxies. We show that the HZPT model for matter correlators can account for the effects of a wide range of baryonic feedback models and provide two extended dark matter models which are of 1% (3%) accuracy for k < 10 (8) h/Mpc. We explicitly model the non-perturbative features of halo exclusion for the halo-halo and galaxy-galaxy correlators, as well as the presence of satellites for galaxy-matter and galaxy-galaxy correlation functions. We perform density estimation using N-body simulations and a wide range of HOD galaxy mocks to obtain correlations of model parameters with the cosmological parameters Ω m and σ 8 . HZPT can provide a fast, interpretable, and analytic model for combined-probe analyses of redshift surveys using scales well into the non-linear regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

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

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances, Weak Lensing, and Galaxy Correlations

We present the first joint analysis of cluster abundances and auto or cross-correlations of three cosmic tracer fields: galaxy density, weak gravitational lensing shear, and cluster density split by optical richness. From a joint analysis (4×2pt+N) of cluster abundances, three cluster cross-correlations, and the auto correlations of the galaxy density measured from the first year data of the Dark Energy Survey, we obtain Ωm=0.305-0.038+0.055 and σ8=0.783-0.054+0.064. This result is consistent with constraints from the DES-Y1 galaxy clustering and weak lensing two-point correlation functions for the flat νΛCDM model. Consequently, we combine cluster abundances and all two-point correlations from across all three cosmic tracer fields (6×2pt+N) and find improved constraints on cosmological parameters as well as on the cluster observable-mass scaling relation. This analysis is an important advance in both optical cluster cosmology and multiprobe analyses of upcoming wide imaging surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
DOE OSTI2020

A benchmark of Gutzwiller conjugate gradient minimization method in ground state energy calculations of dimers

Herein we present numerical results of ground-state energies of 9 molecules in the well-established G2 molecule set given by the Gutzwiller conjugate gradient minimization (GCGM) method. The method, beyond the commonly used Gutzwiller approximation, was recently developed based on Gutzwiller variational wave functions. We find that compared to benchmark data given by full configuration interaction, GCGM total energies are reasonably well reproduced with the minimum basis set. To include the dynamical correlation beyond the minimal basis calculations, we adopt the local density approximation for the dynamical correlation energy $E_c$. By comparing the results with benchmark data given by experiments and large-basis configuration interaction, the GCGM total energies with $E_c$ are in general better reproduced, but discrepancies are still observed for some dimers.

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

Cosmological constraints from density-split clustering in the BOSS CMASS galaxy sample

We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1 \, h^{-1}\, \text{Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-ΛCDM framework, including redshift-space and Alcock–Paczynski distortions. Our base-ΛCDM analysis finds ω cdm = 0.1201 ± 0.0022, σ 8 = 0.792 ± 0.034, and n s = 0.970 ± 0.018, corresponding to fσ 8 = 0.462 ± 0.020 at z ≈ 0.525, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-ΛCDM, varying the running of the spectral index, the dark energy equation of state, and the density of mass-less relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy–halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: cosmology from galaxy×CMB lensing and shear×CMB lensing

ABSTRACT We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence 〈δgκCMB〉, and source galaxy weak lensing shear with CMB lensing convergence 〈γκCMB〉. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of 〈δgκCMB〉 + 〈γκCMB〉 changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from 〈δgκCMB〉 and 〈γκCMB〉 are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in 〈δgκCMB〉 becomes subdominant to cosmic variance, preventing 〈δgκCMB〉 to further improve the constraints. This implies that to maximize the cosmological constraints from future 〈δgκCMB〉 + 〈γκCMB〉 analyses, we should focus more on the requirements on 〈γκCMB〉 instead of 〈δgκCMB〉. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2024

Quantifying carrier density in monolayer MoS 2 by optical spectroscopy

The successful design and device integration of nanoscale heterointerfaces hinges upon precise manipulation of both ground- and excited-state charge carrier (electron and hole) densities. However, it is particularly challenging to quantify these charge carrier densities in nanoscale materials, leading to uncertainties in the mechanisms of many carrier density-dependent properties and processes. Here, we demonstrate a method that utilizes steady-state and transient absorption spectroscopies to correlate monolayer MoS 2 electron density with the easily measured metric of excitonic optical absorption quenching in a variety of mixed-dimensionality s-SWCNT/MoS 2 heterostructures. By employing a 2D phase-space filling model, the resulting correlation elucidates the relationship between charge density, local dielectric environment, and concomitant excitonic properties. The phase-space filling model is also able to describe existing trends from the literature on transistor-based measurements on MoS 2 , WS 2 , and MoSe 2 monolayers that were not previously compared to a physical model, providing additional support for our method and results. The findings provide a pathway to the community for estimating both ground- and excited-state carrier densities in a wide range of TMDC-based systems.

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

Pathways to High-Power-Density Redox Flow Batteries

Redox flow batteries (RFBs) promise to fill a crucial missing link in the energy transition: inexpensive and widely deployable grid and industrial-scale energy storage for intermittent renewable electricity. While numerous lab-scale and demonstration-scale RFBs have been delivered, widespread commercial deployment is still limited by high electrolyte, stack, and balance of plant capital costs. Increasing the power density of RFBs is correlated with lower stack costs, primarily because the area needed for expensive electrode and membrane components to reach a target power density is reduced. Here, in the present contribution, we summarize the areal power densities reported for lab-scale RFBs, critically evaluate major pathways employed for power optimization, and identify opportunities for developing yet-higher power density systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗