SEARCH · Search NASA

Results for “Approximations and Expansions”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Urbanization exacerbates continental- to regional-scale warming

Urbanization is usually ignored when estimating past changes in large-scale climate and for future climate projections since cities historically covered a small fraction of the Earth’s surface. Here, by combining global land surface temperature observations with historical estimates of urban area, we demonstrate that the urban contribution to continental- to regional-scale warming has become non-negligible, especially for rapidly urbanizing regions and countries in Asia. Consequently, expected urban expansion over the next century suggests further increased urban influence on large-scale surface climate in the future (approximately 0.16 K for North America and Europe for high-emission scenario in 2100). Based on these results, also seen for air temperature, we argue that, in line with other forms of land use/land cover change, urbanization should be explicitly included in climate change assessments. This requires incorporation of dynamic urban extent and biophysics in current-generation Earth system models to quantify potential urban feedback on the climate system across scales.

Chakraborty, TC↗

DOE OSTI2020

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗

DOE OSTI2020

Large-d phase transitions in holographic mutual information

In the AdS/CFT correspondence, the entanglement entropy of subregions in the boundary CFT is conjectured to be dual to the area of a bulk extremal surface at leading order in G N in the holographic limit. Under this dictionary, distantly separated regions in the CFT vacuum state have zero mutual information at leading order, and only attain nonzero mutual information at this order when they lie close enough to develop significant classical and quantum correlations. Previously, the separation at which this phase transition occurs for equal-size ball-shaped regions centered at antipodal points on the boundary was known analytically only in 3 spacetime dimensions. Inspired by recent explorations of general relativity at large- d , we compute the separation at which the phase transition occurs analytically in the limit of infinitely many spacetime dimensions, and find that distant regions cannot develop large correlations without collectively occupying the entire volume of the boundary theory. We interpret this result as illustrating the spatial decoupling of holographic correlations in the large- d limit, and provide intuition for this phenomenon using results from quantum information theory. We also compute the phase transition separation numerically for a range of bulk spacetime dimensions from 4 to 21, where analytic results are intractable but numerical results provide insight into the dimension-dependence of holographic correlations. For bulk dimensions above 5, our exact numerical results are well approximated analytically by working to next-to-leading order in the large- d expansion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Approximate symmetries in d = 4 CFTs with an Einstein gravity dual

By applying the stress-tensor-scalar operator product expansion (OPE) twice, we search for algebraic structures in d = 4 conformal field theories (CFTs) with a pure Einstein gravity dual. We find that a rescaled mode operator defined by an integral of the stress tensor T ++ on a d = 2 plane satisfies a Virasoro-like algebra when the dimension of the scalar is large. The structure is enhanced to include a Kac-Moody-type algebra if we incorporate the T -- component. In our scheme, the central terms are finite. It remains challenging to directly compute the stress-tensor sector of d = 4 scalar four-point functions at large central charge, which, based on holography and bootstrap methods, were recently shown to have a Virasoro/ $\mathcal{W}$-algebra vacuum block-like structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Improved opacity expansion at NNLO for medium induced gluon radiation

When an energetic parton propagates in a hot and dense QCD medium it loses energy by elastic scatterings or by medium-induced gluon radiation. The gluon radiation spectrum is suppressed at high frequency due to the LPM effect and encompasses two regimes that are known analytically: at high frequencies ω > ω c = q ^ L 2 , where q ^ is the jet quenching transport coefficient and L the length of the medium, the spectrum is dominated by a single hard scattering, whereas the regime ω < ωc is dominated by multiple low momentum transfers. In this paper, we extend a recent approach (dubbed the Improved Opacity Expansion (IOE)), which allows an analytic (and systematic) treatment beyond the multiple soft scattering approximation, matching this result with the single hard emission spectrum. We calculate in particular the NNLO correction analytically and numerically and show that it is strongly suppressed compared to the NLO indicating a fast convergence of the IOE scheme and thus, we conclude that it is sufficient to truncate the series at NLO. We also propose a prescription to compare the GW and the HTL potentials and relate their parameters for future phenomenological works.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗

DOE OSTI2022

Anharmonic thermodynamic properties and phase boundary across the postperovskite transition in MgSi O 3

To address the effects of lattice anharmonicity across the perovskite to postperovskite transition in MgSiO 3 , we conduct calculations using the phonon quasiparticle (PHQ) approach. The PHQ is based on abinitio molecular dynamics and, in principle, captures full anharmonicity. Free energies in the thermodynamic limit (N → ∞) are computed using temperature-dependent quasiparticle dispersions within the phonon gas model. Systematic results on anharmonic thermodynamic properties and phase boundary are reported. Both the local density approximation and the generalized gradient approximation calculations are performed to provide confident constraints on these properties. Anharmonic effects are demonstrated by comparing results with those obtained using the quasiharmonic approximation (QHA). The inadequacy of the QHA is indicated by its overestimation of thermal expansivity and thermodynamic Grüneisen parameter and its converged isochoric heat capacity in the high-temperature limit. The PHQ phase boundary has a Clapeyron slope (dP/dT) that increases with temperature. This result contrasts with the nearly zero curvature of the QHA phase boundary. Anharmonicity bends the phase boundary to lower temperatures at high pressures. Implications for the double-crossing of the phase boundary by the mantle geotherm are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Coupling of magnetism, crystal lattice, and transport in EuCuP and EuCuAs

EuCuP and EuCuAs are members of a family of materials with candidates for realizing magnetically tuned electronic topology. Here, in this work, the magnetic phase transitions and magnetoelastic coupling of EuCuP and EuCuAs have been investigated. Around the Curie temperature, the thermal expansion coefficient and specific-heat capacity of EuCuP are found to have two closely spaced features, at approximately 30 and 31.3 K, which may indicate a cascading magnetic transition or perhaps an electronic transition closely coupled to the magnetic ordering. The zero-field magnetoelastic coupling appears to be stronger in EuCuP than in EuCuAs, and field-dependent measurements suggest this is due to the dominant ferromagnetic interaction in EuCuP as opposed to antiferromagnetic order in EuCuAs below T N =13.5 K. Thermal expansion is anisotropic around T N in EuCuAs, with a maximum expansion along [001] occurring above T N prior to a stiffening of this c -axis component on cooling through T N ; the expansion of the basal plane has a lambda-like peak centered slightly above T N . Matching these behaviors, the ac susceptibility data for EuCuAs suggest the presence of strong ferromagnetic correlations above T N in the vicinity where the c-axis expansion peaks. Both compounds possess similar electrical and thermoelectric transport behaviors, with short-range magnetic order likely playing an important role above T C or T N . Transport properties suggest these are predominantly hole doped due to intrinsic defects, and narrow-gap semiconducting or semimetallic behavior may be achievable if the underlying defects can be tuned.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Toward more accurate adiabatic connection approach for multireference wavefunctions

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the “fixed-RDM” approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Anharmonic Origin of the Giant Thermal Expansion of NaBr

ll phonons in a single crystal of NaBr are measured by inelastic neutron scattering at temperatures of 10, 300, and 700 K. Even at 300 K, the phonons, especially the longitudinal-optical phonons, show large shifts in frequencies and show large broadenings in energy owing to anharmonicity. Ab initio computations are first performed with the quasiharmonic approximation (QHA) in which the phonon frequencies depend only on V and on T only insofar as it alters V by thermal expansion. This QHA is an unqualified failure for predicting the temperature dependence of phonon frequencies, even 300 K, and the thermal expansion is in error by a factor of 4. Ab initio computations that include both anharmonicity and quasiharmonicity successfully predict both the temperature dependence of phonons and the large thermal expansion of NaBr. The frequencies of longitudinal-optical phonon modes decrease significantly with temperature owing to the real part of the phonon self-energy from explicit anharmonicity originating from the cubic anharmonicity of nearest-neighbor NaBr bonds. Anharmonicity is not a correction to the QHA predictions of thermal expansion and thermal phonon shifts but dominates the behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

A fivefold UO$_2^{2+}$ node is a path to dodecagonal quasicrystal approximants in coordination polymers [plus supplemental information]

Aperiodic formations continue to focus interest in areas ranging from advanced scientific theories to practical everyday applications. Starting from diverse and tightly bonded intermetallic compounds, this world showed an important breakthrough toward the so-called soft systems of meso/macroscale: liquid crystals, thin films, polymers, proteins, etc. This work opens a route for making bulk quasicrystals (QC) in an unprecedented but very common area, with molecular ligands. Since these systems are, to a large extent, transparent, they extend the possible areas of QC application to previously unreachable corners, e.g., photonics. We combined efficient bridging ligands with uranyl pentagonal bonding centers and, unexpectedly, brought the unique attributes of f-element coordination chemistry to an interdisciplinary area of aperiodic formations. Taking advantage of the planar coordination of uranyl ions, we were able to direct the structure expansion solely in two directions with a characteristic snub square tiling, a predicted but previously unobtainable dodecagonal approximant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗