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 469 records · Page 26

Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ~4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. In conclusion, these results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

Meteoroid Bulk Density and Ceplecha Types

The determination of asteroid bulk density is an important aspect of Near Earth Object (NEO) characterization. A fraction of meteoroids originate from asteroids (including some NEOs), thus in lieu of mutual perturbations, satellites, or expensive spacecraft missions, a study of meteoroid bulk densities can potentially provide useful insights into the densities of NEOs and PHOs (Potentially Hazardous Objects). Meteoroid bulk density is still inherently difficult to measure, and is most often determined by modeling the ablation of the meteoroid. One approach towards determining a meteoroid density distribution entails using a more easily measured proxy for the densities, then calibrating the proxy with known densities from meteorite falls, ablation modelling, and other sources. An obvious proxy choice is the Ceplecha type, KB (Ceplecha, 1958), which is thought to indicate the strength of a meteoroid and often correlated to different bulk densities in literature. KB is calculated using the air density at the beginning height of the meteor, the initial velocity, and the zenith angle of the radiant; quantities more readily determined than meteoroid bulk density itself. Numerical values of K(sub B) are sorted into groups (A, B, C, etc.), which have been matched to meteorite falls or meteor showers with known composition such as the porous Draconids. An extensive survey was conducted to establish the strength of the relationship between bulk density and K(sub B), specifically looking at those that additionally determined K(sub B) for the meteors. In examining the modeling of high-resolution meteor data from Kikwaya et al. (2011), the correlation between K(sub B) and bulk density was not as strong as hoped. However, a distinct split by dynamical type was seen with Jovian Tisserand parameter (T(sub J)), with meteoroids from Halley Type comets (T(sub J) < 2) exhibiting much lower bulk densities than those originating from Jupiter Family comets and asteroids (T(sub J) > 2). Therefore, this work indicates that the dynamical classification of a meteoroid is a better indicator of the density than the strength proxy, a somewhat surprising result.

Blaauw, R. C.↗

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↗

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↗

On the geometric nature of high energy nucleus-nucleus reaction cross sections

Within the context of a high energy double-folding optical potential approximation to the exact nucleus-nucleus multiple-scattering series, eikonal scattering theory is used to investigate the validity of geometric reaction cross sections in relativistic heavy ion collisions. The potential used includes a finite range interaction and nuclear single-particle densities extracted from nuclear charge distributions by unfolding the finite proton charge distribution. Pauli correlation effects are also included in an approximate way. The sensitivity of the predictions to the assumed interaction, Pauli correlation approximation, and nuclear density distributions is investigated. These results are in agreement with early predictions concerning the geometric nature of relativistic heavy ion collisions and in disagreement with a recent analysis, utilizing the zero range approximation, which suggested otherwise. Reasons for the lack of agreement between the analyses are also presented. Finally, approximate applicability limits for geometric reaction cross sections are determined.

Townsend, L. W.↗

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↗

Atmospheric attenuation for correlated satellite communication ground sites

Link margin estimates are calculated for 30/20 GHz satellite communication systems employing closely-spaced (4 to 100 km) 'dual diversity' switched ground sites. The link margin estimates are based on a new analysis in which the bivariate rain attenuation density function for two correlated ground sites is modelled by an exponential density function. The results of the exponential density function analysis enable derivation of a direct relation between rain margin estimates and probability of exceedance (link availability). Margins typically in the range 2 to 12 dB are calculated for various ground site separations and summarized for seven city locations and five satellite orbit positions over the U.S. The results can be easily extended to other EHF satellite frequencies and other locations.

Christopher, P.↗

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↗

Ultraviolet observations of galaxies with the FAUST experiment

We have used the set of point sources detected by the Far Ultraviolet Space Telescope (FAUST) instrument to identify galaxies and study the total galaxy flux in a 250 A wide band peaking at 1650 A. A sample of 144 galaxies has been obtained after cross-reference with the RC3 catalog, elimination of objects confused with stars and various corrections for the photometry. The UV-B color dispersion is found to increase while the galaxies get redder from late to early types. The irregular galaxies appear on average redder and the Sbc galaxies bluer than indicated by the spectral energy distributions currently used for the calculations of K-corrections. Various arguments lead us to make the assumption of a constant dust extinction within each galaxy. The UV flux per unit area decreases on average from late to early type spirals. We find a weak correlation between the UV and far infra-red emission while the infra-red to UV flux ratio gets lower when galaxies get bluer (as measured by the UV to B flux ratio). The UV flux per unit area correlates with the HI gas surface density and the total gas surface density when this quantity is available. The correlation with the molecular gas alone is weak. In the Virgo cluster, the UV flux per unit area does not decrease in direct proportion to the HI deficiency. Galaxy counts per square degree and per magnitude interval have been obtained at high-galactic latitudes. Combined with data at fainter magnitudes, they show a variation as a function of magnitude with a near-euclidean slope over a range of 8 magnitudes.

Deharveng, J.-M.↗

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↗

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↗

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↗

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↗

A hydrodynamic approach to cosmology - Texture-seeded cold dark matter and hot dark matter cosmogonies

Hydrodynamical simulations of galaxy formation in a texture-seeded cosmology are presented, with attention given to Omega = 1 galaxies dominated by both hot dark matter (HDM) and cold dark matter (CDM). The simulations include both gravitational and hydrodynamical physics with a detailed treatment of collisional and radiative thermal processes, and use a cooling criterion to estimate galaxy formation. Background radiation fields and Zel'dovich-Sunyaev fluctuations are explicitly computed. The derived galaxy mass function is well fitted by the observed Schechter luminosity function for a baryonic M/L of 3 and total M/L of 60 in galaxies. In both HDM and CDM texture scenarios, the 'galaxies' and 'clusters' are significantly more strongly correlated than the dark matter due to physical bias processes. The slope of the correlation function in both cases is consistent with observations. In contrast to Gaussian models, peaks in the dark matter density distributrion are less correlated than average.

Cen, R. Y.↗

Neutral hydrogen and star formation in irregular galaxies

The Very Large Array and WSTR H I synthesis observations of seven irregular galaxies are presented. The total H I images of four Local Group dwarf irregular galaxies and three larger more distant irregular galaxies are constructed at the identical resolution of 500 pc. When compared to H II region distributions derived from H alpha images, all galaxies studied show an excellent correlation between the H I surface density and the presence of H II regions. This correlation is most easily interpreted in terms of a requisite threshold H I surface density for massive star formation. This threshold is 1 x 10 to the 21st power H I atoms/sq cm for a resolution of 500 pc. Giant extragalactic H II regions are only found near H I surface densities of a factor of 3 to 5 times this threshold level. The observed threshold implies a Jeans length of 150 pc, which is the same as the size scale at which the structure in the H I complexes correlates well with the H II region distribution. This, combined with the fact that in none of the galaxies observed is there H I above the threshold level with concomitant H II regions, implies an exclusively gravitational origin for the star formation events. That is, there is no need to involve a trigger as in the SSPSF theory (Seiden 1983) or feedback as in Dopita (1985).

Skillman, Evan D.↗

Observation of forbidden C II 158 micron emission from the diffuse interstellar medium at high Galactic latitude

We report the first detection of forbidden C II 158 micron line emission from the diffuse interstellar medium at high Galactic latitude. We have measured the integrated line intensity in a 36 arcmin field of view along a triangular scan path in a 5 x 20 deg region in Ursa Major using a rocket-borne, liquid-helium-cooled spectrophotometer. The scan included high-latitude infrared cirrus, molecular clouds, a bright external galaxy, M82, and the H I Hole, which is a region of uniquely low neutral hydrogen column density. Emission from forbidden C II is observed in all regions, and, in the absence of appreciable CO emission, it is well correlated with neutral hydrogen column density. We observe a forbidden C II gas cooling rate which varies from (1.18 +/- 0.4 to 3.25 +/- 0.8) x 10 exp -26 ergs/s/H atom. Regions with CO emission have enhanced forbidden C II line emission over that expected from the correlation with neutral hydrogen column density. We measure a line-to-continuum ratio which varies from 0.002 to 0.008 in comparison with the all-sky average of 0.0082 reported by FIRAS, which is heavily weighted toward the Galactic plane.

Bock, J. J.↗