Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “highly correlated density”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structure-based mechanism and inhibition of cholesteryl ester transfer protein

Purpose of Review: Cholesteryl ester transfer proteins (CETP) regulate plasma cholesterol levels by transferring cholesteryl esters (CEs) among lipoproteins. Lipoprotein cholesterol levels correlate with the risk factors for atherosclerotic cardiovascular disease (ASCVD). This article reviews recent research on CETP structure, lipid transfer mechanism, and its inhibition. Recent Findings: Genetic deficiency in CETP is associated with a low plasma level of low-density lipoprotein cholesterol (LDL-C) and a profoundly elevated plasma level of high-density lipoprotein cholesterol (HDL-C), which correlates with a lower risk of atherosclerotic cardiovascular disease (ASCVD). However, a very high concentration of HDL-C also correlates with increased ASCVD mortality. Considering that the elevated CETP activity is a major determinant of the atherogenic dyslipidemia, i.e., pro-atherogenic reductions in HDL and LDL particle size, inhibition of CETP emerged as a promising pharmacological target during the past two decades. CETP inhibitors, including torcetrapib, dalcetrapib, evacetrapib, anacetrapib and obicetrapib, were designed and evaluated in phase III clinical trials for the treatment of ASCVD or dyslipidemia. Although these inhibitors increase in plasma HDL-C levels and/or reduce LDL-C levels, the poor efficacy against ASCVD ended interest in CETP as an anti-ASCVD target. Nevertheless, interest in CETP and the molecular mechanism by which it inhibits CE transfer among lipoproteins persisted. Insights into the structural-based CETP-lipoprotein interactions can unravel CETP inhibition machinery, which can hopefully guide the design of more effective CETP inhibitors that combat ASCVD. Summary: Individual-molecule 3D structures of CETP bound to lipoproteins provide a model for understanding the mechanism by which CETP mediates lipid transfer and which in turn, guide the rational design of new anti-ASCVD therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

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

Data Files for P. Mai et al., "Fluctuating charge-density-wave correlations in the three-band Hubbard model" (2024)

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

cuprates↗
DOE OSTI2020

Diagnostic development for parallel wave-number measurement of lower hybrid waves in EAST

In this study, an eight-channel magnetic probe diagnostic system has been designed and installed adjacent to the 4.6 GHz lower hybrid (LH) grill antenna in the low-field side of the Experimental Advanced Superconducting Tokamak (EAST) in order to study the n ∥ evolution of LH waves in the first pass from the launcher to the core plasma. The magnetic probes are separated by 6.6 mm, which allows measurement of the dominant parallel refractive index n ∥ up to n ∥ = 5 for 4.6 GHz LH waves. The magnetic probes are designed to be sensitive to the magnetic field component perpendicular to the background magnetic field with a slit on the casing that encloses the probe. The intermediate frequency stage, which consists of two mixing stages, down-coverts the frequency of the measured wave signals at 4.6 GHz to 20 MHz. A bench test demonstrates the phase stability of the magnetic probe diagnostic system. By evaluating the phase variation of the measured signals along the background magnetic field, the dominant n ∥ of the LH wave in the scrape-off layer has been deduced during the 2019 experimental campaign. In the low density plasma, the measured dominant n ∥ of the LH waves is about 2.1, corresponding to the main peak 2.04 of the launched n ∥ spectrum. n ∥ deduced by the least-squares linear fit method remains near this value in the low density plasma with a high spatial correlation magnitude of 0.9. With an eight-channel probe system, a wave-number spectrum has also been deduced, which has a peak near to the measured dominant n ∥ .

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Diagnostic development for parallel wave-number measurement of lower hybrid waves in EAST

An eight-channel magnetic probe diagnostic system has been designed and installed adjacent to the 4.6 GHz lower hybrid (LH) grill antenna in the low-field side of the EAST tokamak in order to study the n|| evolution of lower hybrid waves in the first pass from the launcher to the core plasma. The magnetic probes are separated by 6.6 mm, which allows measurement of the dominant parallel refractive index n|| up to n|| =5 for 4.6GHz LH waves. The magnetic probes are designed to be sensitive to the magnetic field component perpendicular to the background magnetic field with a slit on the casing that encloses the probe. The intermediate frequency (IF) stage, which consists of two mixing stages, down-coverts the frequency of the measured wave signals at 4.6 GHz to 20 MHz. A bench test demonstrates the phase stability of the magnetic probe diagnostic system. By evaluating the phase variation of the measured signals along the background magnetic field, the dominant n|| of the LH wave in the scrape-off layer (SOL) has been deduced during the 2019 experimental campaign. In the low density plasma, the measured dominant n|| of the lower hybrid waves is about 2.1, corresponding to the main peak 2.04 of the launched n|| spectrum. The n|| deduced by the least square linear fit method remains near this value in the low density plasma with a high spatial correlation magnitude of 0.9. With an eight-channel probe system a wave-number spectrum has also been deduced, which has a peak near to the measured dominant n||.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physical mechanisms for the transition from type-III to large ELMs induced by impurity injection on EAST

Transition from type-III to large-amplitude ELMs induced by neon injection has been observed in the EAST tokamak at overlapping q 95 space between large and small ELMs. With neon injection, pedestal density gradient shows a remarkable increase accompanied by some decrease of pedestal electron temperature, and consequently the pressure gradient increases moderately and edge bootstrap current has minimal change. Further experiment demonstrates that the occurrence of large ELMs after neon injection is highly correlated with the change in edge density. Linear peeling-ballooning stability analysis indicates that the large ELM case is more unstable than the type-III ELM case during the ELM transition. A scan of pedestal density gradient in linear stability analysis shows that the direct destabilizing effect of steep pedestal density gradient on peeling-ballooning instabilities via two-fluid effects could also facilitate the transition to large ELMs. In conclusion, these results could provide more insight into the role of pedestal density gradient on pedestal stability and ELM behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of correlations between intrinsic ductility and electronic density of states in refractory alloys

High entropy alloys (HEAs) correspond to a new and emerging class of materials that allows us to explore a large composition space to tune mechanical strength and thermal stability. Therefore, to design better alloys, it is important to scan the high-dimensional space of chemistry, composition and temperature. Here, to facilitate this search, we present a method to screen intrinsically ductile body centered cubic (BCC) refractory alloys from electronic structure calculations by using the density of states (DOS) at the Fermi level, g(μ F ). This correlation between intrinsic ductility and g(μ F ) is tested by analyzing group V (V, Nb, Ta) and VI (Mo, W) refractory metals, binary alloys, such as W-Nb, W-V, Mo-Nb and Mo-V, and refractory alloys for which experimental stress-strain measurements are available. In addition, we perform a high-throughput exploration of the entire composition space of a recently proposed alloy system, CrMoNbV, and identify compositions that exhibit high intrinsic ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Resident and transient coyotes exhibit differential patterns of movement behavior across heterogeneous landscapes in the southeastern United States

Coyotes (Canis latrans) are a highly adaptable canid species whose behavioral plasticity has allowed them to persist in a wide array of habitats throughout North America. As generalists, coyotes can alter movement patterns and change territorial strategies between residency (high site fidelity) and transiency (low site fidelity) to maximize fitness. Uncertainty remains about resident and transient coyote movement patterns and habitat use because research has reached conflicting conclusions regarding patterns of habitat use by both groups. We quantified effects of habitat on resident and transient coyote movement behavior using first passage time (FPT) analysis, which assesses recursive movement along an individual's movement path to delineate where they exhibit area-restricted search (ARS) behaviors relative to habitat attributes. We quantified monthly movement rates for 171 coyotes (76 residents and 53 transients) and then used estimated FPT values in generalized linear mixed models to quantify monthly habitat use for resident and transient coyotes. Transients had greater movement rates than residents across all months except January. Resident FPT values were positively correlated with agricultural land cover during fall and winter, but negatively correlated with agriculture during spring. Resident FPT values were also negatively correlated with developed habitats during May–August, deciduous land cover during June–August, and wetlands during September–January except November. FPT values of transient coyotes were positively correlated with developed areas throughout much of the year and near wetlands during July–September. Transient FPT values were negatively correlated with agriculture during all months except June and July. High FPT values (ARS behavior) of residents and transients were generally correlated with greater densities of edge habitat. Although we observed high individual variation in space use, our study found substantive differences in habitat use between residents and transients, providing further evidence that complexity and plasticity of coyote habitat use is influenced by territorial strategy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Density weighted angular redshift fluctuations: a new cosmological observable

We propose the use of angular fluctuations in the galaxy redshift field as a new way to extract cosmological information in the Universe. This new probe $\delta z (\hat{\mathbf {n}})$ consists of the statistics of sky maps built by projecting redshifts under a Gaussian window of width σz centred upon a redshift zobs, and weighted by the galaxy density field. We compute the angular power spectrum of the $\delta z (\hat{\mathbf {n}})$ field in both numerical simulations and in linear perturbation theory. From these, we find that the $\delta z (\hat{\mathbf {n}})$ field (i) is sensitive to the underlying density and peculiar velocity fields; (ii) is highly correlated, at the $\gtrsim 60\, {{\ \rm per\ cent}}$ level, to the line-of-sight projected peculiar velocity field; (iii) for narrow windows (σz < 0.03), it is almost completely uncorrelated to the projected galaxy angular density field under the same redshift window; and (iv) it is largely unaffected by multiplicative and additive systematic errors on the observed number of galaxies that are redshift-independent over ~σz. We conclude that $\delta z (\hat{\mathbf {n}})$ is a simple and robust tomographic measure of the cosmic density and velocity fields, complementary to angular clustering, that will contribute to more complete exploitations of current and upcoming galaxy redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Association‐based activity coefficient model for electrolyte solutions

Abstract This work presents an association‐based activity coefficient model that explicitly considers the solution nonideality due to associations among ions and solvent species. Built upon the electrolyte nonrandom two‐liquid (eNRTL) model, the model greatly improves the accuracy of eNRTL model for strongly associating electrolyte solutions due to presence of ionic species with high surface charge density. The model successfully correlates mean ionic activity coefficients of 46 aqueous single‐salt systems from 10 cations and 5 anions at 298.15 K up to their solubility limits. With the ion‐specific association parameters identified, the model accurately predicts activity and osmotic coefficients for aqueous mixed‐salt systems at 298.15 K. The temperature dependence of the model results has also been examined at 273–373 K. With superior accuracy over a wide range of concentration and temperature, the model represents a major advancement over eNRTL model and has a great potential to be a next‐generation model for electrolyte solutions.

Lin, Yu‐Jeng↗
DOE OSTI2020

Sublimation Kinetics for Individual Graphite and Graphene Nanoparticles (NPs): NP-to-NP Variations and Evolving Structure-Kinetics and Structure-Emissivity Relationships

A single nanoparticle (NP) mass spectrometry method was used to measure sublimation rates as a function of nanoparticle temperature (T NP ) for sets of individual graphite and graphene NPs. Initially, the NP sublimation rates were ~400 times faster than that for bulk graphite, and there were large NP-to-NP variations. Over time, the rates slowed substantially, though remaining well above the bulk rate. The initial activation energies (E a values) were correspondingly low and doubled as a few monolayer’s worth of material were sublimed from the surfaces. The high initial rates and low E a values are attributed to large numbers of edge, defect, and other low coordination sites on the NP surfaces, and the changes are attributed to atomic-scale “smoothing” of the surface by preferential sublimation of the less stable sites. The emissivity of the NPs also changed after heating, most frequently increasing. Here, the emissivity and sublimation rates were anti-correlated, leading to the conclusion that high densities of low-coordination sites on the NP surfaces enhances sublimation but suppresses emissivity.

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

Experimental study of high-k turbulence during an energy confinement degradation phase in EAST ohmic plasmas

In this paper, we present experimental studies of both high-k(r) ( cm(-1), and , respectively; ) and high- ( cm(-1), and , respectively; ) turbulence behavior during an ohmic energy confinement degradation phase in experimental advanced superconducting tokamak (EAST). High-k(r) turbulence from density fluctuation at and high- turbulence from density fluctuation at -0.97 were measured by tangential and poloidal CO2 laser collective scattering diagnostics, respectively. Note that k(r), , and are radial wavenumber, poloidal wavenumber, perpendicular wavenumber and ion gyroradius at electron temperature, respectively. Both high-k(r)/ turbulence power and energy confinement time are found to be temporally correlated to line-averaged electron density n(e) in the plasma current flat-top phase (I-p = 0.4 MA): when the n(e) shows continuous increase/decrease, the high-k(r)/ turbulence power increases/decreases correspondingly and the shows corresponding decrease/increase; the stable n(e) is related to stable both and high-k(r)/ turbulence power. Statistical results of high-k(r)/ turbulence power versus further imply that high-k(r)/ turbulence shows a strong correlation with plasma energy confinement degradation in high line-averaged n(e) condition, but high-k(r) and high- turbulence have relatively weak and no dependence on the transition between linear range in low-n(e) condition and energy confinement degradation in the high-n(e) condition, respectively. Moreover, profiles of electron temperature T-e and n(e) as well as their normalized gradients in a part of high- density fluctuation measurement region also have been given to qualitatively explain the enhancement of turbulence power with the increase of line-averaged n(e).

EAST tokamak↗
DOE OSTI2025

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2024

Report Summarizing the Mechanical Properties of a Large ODS Ferritic Alloy Ingot by Forging at High Temperatures for Future Mother Tube Production

Oxide dispersion strengthened (ODS) ferritic alloys are considered the benchmark fuel cladding and core structural material in advanced nuclear energy reactors that require high-temperature strength and creep properties and resistance to radiation damage. The ODS ferritic alloys including 14YWT are produced by mechanical alloying (MA), which is time-consuming and is associated with high manufacturing costs that are not beneficial to being used in advanced nuclear energy reactors. This report summarizes the accomplishments in FY24 in the INM program for producing an ODS ferritic alloy by high-deformation, high-temperature processing of reactive and ferritic alloy powders for advanced reactor fuel cladding applications. Following four hot forging experiments, it was concluded that several significant challenges were encountered that were difficult to overcome. These challenges were related to using high annealing temperatures for pressure assisted sintering processes to produce dense microstructures, but the high annealing temperatures result in long range diffusion of Ti through the bcc Fe lattice and react with the YIG particles that are distributed on the prior surfaces of the ferritic alloy powders. In addition, it was determined that the high deformations induced by forging at high temperatures were not very effective for fracturing the reactive YIG particles into smaller particles and distributing them into the interior of the ferritic alloy powders. Spark plasma sintering was attempted for shortening the time at high temperatures, which led to full densification but the microstructure characterization results still showed that Ti atoms diffused over long distances in the bcc Fe lattice and reacted with the YIG particles. Finally, a short 30 minute high intensity ball milling experiment was conducted on blended 14WT and YIG powders for inducing severe deformations that caused the initially spherical powders of 14WT decorated with the YIG particles on the surfaces to transform to flakes and fracturing of YIG particles into smaller particles that were incorporated into the bcc Fe lattice. Two forgings were performed with the high intensity milled powders: first annealing for 20 minutes at 850ºC followed by forging and second annealing for 20 minutes at 1,100ºC followed by forging. At 850ºC, Ti atoms are effectively immobile, thus allowing for densification albeit incomplete. The deformation by forging after annealing for 20 minutes at 1,100ºC favored dynamic recrystallization processes that led to nano-size grains with very high stored energy due to high dislocation density. The corresponding VH data showed significant hardening that correlated with estimated strengthening of ~1650 MPa. This hybrid processing approach combining high intensity ball milling of powder with annealing and forging showed the most promise for producing an oxide dispersing strengthened ferritic alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Correlating dislocation structures to basal and prismatic slip bands near grain boundaries in tensile-tested Mg–4Al using a multiscale electron microscopy approach

Dislocation structures where basal and prismatic slip bands meet grain boundaries in tensile-tested Mg with 4 wt% Al (Mg–4Al) were studied using a multiscale electron microscopy approach to explore the assumptions made for the dislocation pileup theory: a) slip bands consist of dislocation arrays piling up near grain boundaries, and b) stress concentration in the adjacent grain is solely caused by dislocation pileup. After post-testing electron backscatter diffraction (EBSD) scans, focused-ion-beam (FIB) lift-out specimens were prepared from regions of interest so that the specimen plane is parallel to the bulk sample surface to allow for direct correlation between mesoscale and microscale characterization. High dislocation density and the dislocation arrays in the slip bands corroborated with the first assumption. Evidence of plastic deformation in all grains, however, showed that the second assumption was false. Additionally, dislocation structures including crystal rotation boundaries and dislocations dissociating to the basal plane provide additional stress relief mechanisms near grain boundaries that are not accounted for in many crystal plasticity models, which requires auxiliary models to more accurately predict the stress and strain distribution in a microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward improved property prediction of 2D materials using many-body quantum Monte Carlo methods

The field of 2D materials has grown dramatically in the past two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computing applications. These 2D structures, which are often exfoliated from layered van der Waals materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. To overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. Here, this approach has been applied to 2D systems of wide interest, including graphene, phosphorene, MoS 2 , CrI 3 , VSe 2 , GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.

2D materials↗
DOE OSTI2025

High spatial resolution neutron imaging of lithium-ion batteries: Correlating microstructure and lithium transport

Thick electrodes for lithium-ion batteries can increase the overall energy density, but increasing the electrode thickness introduces charge transport limitations. These limitations may be mitigated through proper electrode structuring. Here, high spatial resolution neutron imaging was used to understand the correlation between microstructure and lithium transport in lithium-ion anodes. Batteries with distinct graphite anode microstructures were produced and studied with high spatial resolution in operando neutron radiography to observe the effects of structure on transport. High spatial resolution neutron computed tomography was performed following in operando neutron radiography. X-ray computed tomography and scanning electron microscopy were used to observe the finer scale anode structure to complement neutron imaging. Solvent-free anodes containing a tightly-packed layered structure confined lithium movement close to the separator. This structure limited capacity, but supported better rate capability. Conversely, a more open pore structure in the wet cast anodes yielded higher capacity with reduced rate capability. Together, these results show that lithium distributions can be controlled by the macroscopic structure of the electrodes, the microstructural pore network, and the microscale active areas that support electrochemical reactions. Furthermore, multimodal imaging applying the complementary strengths of neutron and X-ray methods is shown as a tool for advancing battery design.

25 ENERGY STORAGE↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗