SEARCH · Engineering Papers

Results for “Correlation Functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

BH-DFTB/DFT calculations for iron clusters

We present a study on the structural, electronic, and magnetic properties of Fe{sub n}(n = 2 − 20) clusters by performing density functional tight binding (DFTB) calculations within a basin hopping (BH) global optimization search followed by density functional theory (DFT) investigations. The structures, total energies and total spin magnetic moments are calculated and compared with previously reported theoretical and experimental results. Two basis sets SDD with ECP and 6-31G** are employed in the DFT calculations together with BLYP GGA exchange-correlation functional. The results indicate that the offered BH-DFTB/DFT strategy collects all the global minima of which different minima have been reported in the previous studies by different groups. Small Fe clusters have three kinds of packing; icosahedral (Fe{sub 9−13}), centered hexagonal antiprism (Fe{sub 14−17}, Fe{sub 20}), and truncated decahedral (Fe{sub 17(2)}, Fe{sub 18−19}). It is obtained in a qualitative agreement with the time of flight mass spectra that the magic numbers for the small Fe clusters are 7, 13, 15, and 19 and with the collision induced dissociation experiments that the sizes 6, 7, 13, 15, and 19 are thermodynamically more stable than their neighboring sizes. The spin magnetic moment per atom of Fe{sub n}(n = 2 − 20) clusters is between 2.4 and 3.6 μ{sub B} for the most of the sizes. The antiferromagnetic coupling between the central and the surface atoms of the Fe{sub 13} icosahedron, which have already been reported by experimental and theoretical studies, is verified by our calculations as well. The quantitative disagreements between the calculations and measurements of the magnetic moments of the individual sizes are still to be resolved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Comparison of first principles and semi-empirical models of the structural and electronic properties of Ge{sub 1−x}Sn{sub x} alloys

We present and compare three distinct atomistic models—based on first principles and semi-empirical approaches—of the structural and electronic properties of Ge{sub 1−x}Sn{sub x} alloys. Density functional theory calculations incorporating Heyd–Scuseria–Ernzerhof (HSE), local density approximation (LDA) and modified Becke–Johnson (mBJ) exchange-correlation functionals are used to perform structural relaxation and electronic structure calculations for a series of Ge{sub 1−x}Sn{sub x} alloy supercells. Based on HSE calculations, a semi-empirical valence force field (VFF) potential and sp{sup 3}s{sup ∗} tight-binding (TB) Hamiltonian are parametrised. Comparing the HSE, LDA+mBJ and VFF+TB models, and using the HSE results as a benchmark, we demonstrate that: (1) LDA+mBJ calculations provide an accurate first principles description of the electronic structure at reduced computational cost, (2) the VFF potential is sufficiently accurate to circumvent the requirement to perform first principles structural relaxation, and (3) VFF+TB calculations provide a good quantitative description of the alloy electronic structure in the vicinity of the band edges. Our results also emphasise the importance of Sn-induced band mixing in determining the nature of the conduction band structure of Ge{sub 1−x}Sn{sub x} alloys. The theoretical models and benchmark calculations we present inform and enable predictive, computationally efficient and scalable atomistic calculations for disordered alloys and nanostructures. This provides a suitable platform to underpin further theoretical investigations of the properties of this emerging semiconductor alloy.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2020

Unusual electrical conductivity driven by localized stoichiometry modification at vertical epitaxial interfaces

Precise control of lattice mismatch accommodation and interfacial cation diffusion is critical to modulate correlated functionalities in epitaxial heterostructures, particularly when the interface composition is positioned near a compositional phase transition boundary. Here we select La 1-x Sr x MnO 3 (LSMO) as a prototype phase transition oxide and establish vertical epitaxial interfaces with NiO for exploring the strong interplay between strain accommodation, stoichiometry modification, and localized electron transport across the interface. It is found that localized stoichiometry modification overcomes the dead layer problem plaguing LSMO and leads to strongly directional conductivity, as manifested by a more than three orders of magnitude difference between out-of-plane and in-plane conductivity. Comprehensive atomic scale structural characterization and transport measurements reveal that this emerging behavior is created by a compositional change produced by preferential cation diffusion that pushes the LSMO phase transition from insulating into metallic within an ultrathin interface region. This study explores the nature of unusual electrical conductivity at vertical epitaxial interfaces and highlights the beneficial role of controllable cation diffusion that enables emerging functionalities for a broad range of potential applications such as memristors, spintronic devices, and novel nanoelectronic devices using strongly correlated materials.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Observations of non-linear plasmon damping in dense plasmas

Here, we present simulations using finite-temperature density-functional-theory molecular-dynamics to calculate dynamic dielectric properties in warm dense aluminum. The comparison between exchange-correlation functionals in the Perdew, Burke, Ernzerhof approximation, Strongly Constrained and Appropriately Normed Semilocal Density Functional, and Heyd, Scuseria, Ernzerhof (HSE) approximation indicates evident differences in the electron transition energies, dc conductivity, and Lorenz number. The HSE calculations show excellent agreement with x-ray scattering data as well as dc conductivity and absorption measurements. These findings demonstrate non-Drude behavior of the dynamic conductivity above the Cooper minimum that needs to be taken into account to determine optical properties in the warm dense matter regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Curvature-slope correlation of nuclear symmetry energy and its imprints on the crust-core transition, radius, and tidal deformability of canonical neutron stars

Background: The nuclear symmetry energy E sym (ρ) encodes information about the energy necessary to make nuclear systems more neutron-rich. While its slope parameter L at the saturation density ρ 0 of nuclear matter has been relatively well constrained by recent astrophysical observations and terrestrial nuclear experiments, its curvature K sym characterizing the E sym (ρ) around 2ρ 0 remains largely unconstrained. Over 520 calculations for E sym (ρ) using various nuclear theories and interactions in the literature have predicted several significantly different K sym –L correlations. Purpose: If a unique K sym –L correlation of E sym (ρ) can be firmly established, it will enable us to progressively better constrain the high-density behavior of E sym (ρ) using the available constraints on its slope parameter L. Here, we investigate if and by how much the different K sym –L correlations may affect neutron star observables. Method: A meta-model of nuclear Equation of States (EOSs) with three representative K sym –L correlation functions is used to generate multiple EOSs for neutron stars. We then examine effects of the K sym –L correlation on the crust-core transition density and pressure as well as the radius and tidal deformation of canonical neutron stars. Results: The K sym –L correlation affects significantly both the crust-core transition density and pressure. It also has strong imprints on the radius and tidal deformability of canonical neutron stars especially at small L values. The available data from LIGO/VIRGO and NICER set some useful limits for the slope L but can not distinguish the three representative K sym –L correlations considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Charge localization and ordering in A 2 Mn 8 O 16 hollandite group oxides: Impact of density functional theory approaches

The phases of A 2 Mn 8 O 16 hollandite group oxides emerge from the competition between ionic interactions, Jahn-Teller effects, charge ordering, and magnetic interactions. Their balanced treatment with feasible computational approaches can be challenging for commonly used approximations in density functional theory. Three examples (A = Ag, Li, and K) are studied with a sequence of different approximate exchange-correlation functionals. Starting from a generalized gradient approximation (GGA), an extension to include van der Waals interactions and a recently proposed meta-GGA are considered. Then local Coulomb interactions for the Mn 3d electrons are more explicitly considered with the DFT + U approach. Finally, selected results from a hybrid functional approach provide a reference. Results for the binding energy of the A species in the parent oxide highlight the role of van der Waals interactions. Relatively accurate results for insertion energies can be achieved with a low-U and a high-U approach. In the low-U case, the materials are described as band metals with a high-symmetry, tetragonal crystal structure. In the high-U case, the electrons donated by A result in formation of local Mn 3+ centers and corresponding Jahn-Teller distortions characterized by a local order parameter. The resulting degree of monoclinic distortion depends on charge ordering and magnetic interactions in the phase formed. The reference hybrid functional results show charge localization and ordering. Comparison to low-temperature experiments of related compounds suggests that charge localization is the physically correct result for the hollandite group oxides studied here. Lastly, while competing effects in the local magnetic coupling are subtle, the fully anisotropic implementation of DFT + U gives the best overall agreement with results from the hybrid functional.

meta-GGA↗

DOE OSTI2018

Structure and dynamics of water confined in nanoporous carbon

We study the structure and dynamics of water confined in a nanoporous graphitic carbon structure using molecular-dynamics (MD) simulations. The carbon structure is generated by a reactive MD simulation of oxidation of a silicon carbide nanoparticle. We embed water molecules in the nanopores and study structural and dynamical properties of nanoconfined water as a function of temperature. MD simulation results indicate the presence of high-density water (HDW) and low-density water. Radial distribution functions and spatial density functions indicate that the second solvation shell of the HDW is broken. We calculate the self-diffusion coefficient of confined water molecules as a function of temperature and find a significant decrease in the diffusion of water molecules around T = 190K. The cage correlation function c(t) of confined water molecules at T = 200K exhibits stretched exponential decay, c(t) = exp[–(t/τ) β ], with β = 0.43, which matches exactly with the theoretical prediction β = 3/7. Finally, the self-intermediate scattering function at T = 200K indicates differences in small-scale and large-scale dynamics of water molecules.

36 MATERIALS SCIENCE↗

NASA NTRS1971

Test-particle method in kinetic theory of a plasma.

The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.

Matsuda, K.↗

NASA NTRS1994

Voids and constraints on nonlinear clustering of galaxies

Void statistics of the galaxy distribution in the Center for Astrophysics Redshift Survey provide strong constraints on galaxy clustering in the nonlinear regime, i.e., on scales R equal to or less than 10/h Mpc. Computation of high-order moments of the galaxy distribution requires a sample that (1) densely traces the large-scale structure and (2) covers sufficient volume to obtain good statistics. The CfA redshift survey densely samples structure on scales equal to or less than 10/h Mpc and has sufficient depth and angular coverage to approach a fair sample on these scales. In the nonlinear regime, the void probability function (VPF) for CfA samples exhibits apparent agreement with hierarchical scaling (such scaling implies that the N-point correlation functions for N greater than 2 depend only on pairwise products of the two-point function xi(r)) However, simulations of cosmological models show that this scaling in redshift space does not necessarily imply such scaling in real space, even in the nonlinear regime; peculiar velocities cause distortions which can yield erroneous agreement with hierarchical scaling. The underdensity probability measures the frequency of 'voids' with density rho less than 0.2 -/rho. This statistic reveals a paucity of very bright galaxies (L greater than L asterisk) in the 'voids.' Underdensities are equal to or greater than 2 sigma more frequent in bright galaxy samples than in samples that include fainter galaxies. Comparison of void statistics of CfA samples with simulations of a range of cosmological models favors models with Gaussian primordial fluctuations and Cold Dark Matter (CDM)-like initial power spectra. Biased models tend to produce voids that are too empty. We also compare these data with three specific models of the Cold Dark Matter cosmogony: an unbiased, open universe CDM model (omega = 0.4, h = 0.5) provides a good match to the VPF of the CfA samples. Biasing of the galaxy distribution in the 'standard' CDM model (omega = 1, b = 1.5; see below for definitions) and nonzero cosmological constant CDM model (omega = 0.4, h = 0.6 lambda(sub 0) = 0.6, b = 1.3) produce voids that are too empty. All three simulations match the observed VPF and underdensity probability for samples of very bright (M less than M asterisk = -19.2) galaxies, but produce voids that are too empty when compared with samples that include fainter galaxies.

Vogeley, Michael S.↗

DOE OSTI2022

Fuzzy dark matter and the Dark Energy Survey Year 1 data

ABSTRACT Gravitational weak lensing by dark matter haloes leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass m ∼ 10−22 eV, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimized halo model HMCode, accounting for additional suppression of the mass function and halo concentration in FDM as observed in N-body simulations. We combine Dark Energy Survey Year 1 (DES-Y1) data with the Planck cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold dark matter model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalizing over uncertainties related to the non-linear model of FDM, we find a new, independent 95 per cent C.L. lower limit log10m > −23 combining Planck and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal {O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Goldstone bosons on celestial sphere and conformal soft theorems

In this paper, we study celestial amplitudes of Goldstone bosons and conformal soft theorems. Motivated by the success of soft bootstrap in momentum space and the important role of the soft limit behavior of tree-level amplitudes, our goal is to extend some of the methods to the celestial sphere. The crucial ingredient of the calculation is the Mellin transformation, which transforms four-dimensional scattering amplitudes to correlation functions of primary operators in the celestial CFT. The soft behavior of the amplitude is then translated to the singularities of the correlator. Only for amplitudes in “UV completed theories” (with sufficiently good high energy behavior) the Mellin integration can be properly performed. In all other cases, the celestial amplitude is only defined in a distributional sense with delta functions. We provide many examples of celestial amplitudes in UV-completed models, including linear sigma models and Z-theory, which is a certain completion of the SU(N) non-linear sigma model. We also comment on the BCFW-like and soft recursion relations for celestial amplitudes and the extension of soft bootstrap ideas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Critical phenomena of the layered ferrimagnet Mn 3 Si 2 Te 6 following proton irradiation

The critical phenomena and magnetic entropy of the quasi-2D ferrimagnetic crystal, Mn 3 Si 2 Te 6 (MST), is analyzed along the easy axis (H || ab) as a function of proton irradiance. The critical exponents β and γ do not fall into any particular universality class upon proton irradiation. However, for pristine and irradiated samples, the critical exponents lie closer to mean field-like interactions; therefore, long-range interactions are presumed to be sustained in MST. The effective spatial dimensionality reveals that MST remains at d =3 under proton irradiation, whereas spin dimensionality transitions from an initial n =1 to n =2 and n =3 for 1 × 10 15 and 5 × 10 15 H + /cm 2 , indicating XY and Heisenberg interactions, respectively. The spin correlation function reveals an increase in magnetic correlations at 5 × 10 15 H + /cm 2 . Maximum change in magnetic entropy at 3 T is the largest for 5 × 10 15 H + /cm 2 at 2.45 J/kg K, in comparison to 1.60 J/kg K for pristine MST. These results intriguingly align with previous findings on MST where magnetization increased by ~50% at 5 × 10 15 H + /cm 2 , in comparison to its pristine counterpart [Martinez et al., Appl. Phys. Lett. 116, 172404 (2020)]. Magnetic entropy derived from heat capacity analysis shows no large deviations across the proton irradiated samples suggesting that the antiferromagnetic (AFM) coupling between the Mn sites is stable even after proton irradiation. This implies that magnetization is enhanced through a strengthening of the super-exchange interaction between Mn atoms mediated through Te rather than a weakening of the AFM component.

2D materials↗

DOE OSTI2017

Chaos and complexity by design

We study the relationship between quantum chaos and pseudorandomness by developing probes of unitary design. A natural probe of randomness is the “frame poten-tial,” which is minimized by unitary k-designs and measures the 2-norm distance between the Haar random unitary ensemble and another ensemble. A natural probe of quantum chaos is out-of-time-order (OTO) four-point correlation functions. We also show that the norm squared of a generalization of out-of-time-order 2k-point correlators is proportional to the kth frame potential, providing a quantitative connection between chaos and pseudorandomness. In addition, we prove that these 2k-point correlators for Pauli operators completely determine the k-fold channel of an ensemble of unitary operators. Finally, we use a counting argument to obtain a lower bound on the quantum circuit complexity in terms of the frame potential. This provides a direct link between chaos, complexity, and randomness.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2016

A bound on chaos

We conjecture a sharp bound on the rate of growth of chaos in thermal quantum systems with a large number of degrees of freedom. Chaos can be diagnosed using an out-of-time-order correlation function closely related to the commutator of operators separated in time. We conjecture that the influence of chaos on this correlator can develop no faster than exponentially, with Lyapunov exponent λ L ≤ 2πk B T/ℏ. We give a precise mathematical argument, based on plausible physical assumptions, establishing this conjecture.

1/N Expansion↗

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

Azimuthal dependence of two-particle transverse momentum current correlations

Two-particle transverse momentum correlation functions are a powerful technique for understanding the dynamics of relativistic heavy-ion collisions. Among these, the transverse momentum correlator G 2 (Δη,Δφ) is of particular interest for its potential sensitivity to the shear viscosity per unit of entropy density η/s of the quark-gluon plasma formed in heavy-ion collisions. We use the UrQMD, AMPT, and EPOS models for Au–Au at √ s NN = 200 GeV and Pb–Pb at √ s NN = 2760 GeV to investigate the long range azimuthal dependence of G 2 (Δη,Δφ), and explore its utility to constrain η/s based on charged particle correlations. We find that the three models yield quantitatively distinct transverse momentum Fourier harmonics coefficients $a$ $^{pT}_{n}$. We also observe these coefficients exhibit a significant dependence on η/s in the context of the AMPT model. These observations suggest that exhaustive measurements of the dependence of G 2 (Δφ) with collision energy, system size, collision centrality, in particular, offer the potential to distinguish between different theoretical models and their underlying assumptions. Exhaustive analyses of G 2 (Δφ) obtained in large and small systems should also be instrumental in establishing new constraints for precise extraction of η/s.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗