Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 415 records · Page 23

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles study and experimental characterization of metal incorporation in germanium telluride

Germanium telluride is a well-known phase change material (PCM) used in non-volatile memory cells and radio frequency switches. Controlling the properties of GeTe for improved PCM device performance has sometimes been achieved by doping and/or alloying with metals, often at concentrations greater than 10 at. % and using non-equilibrium methods. Here, since switching PCMs between the low-resistance crystalline and high-resistance amorphous states requires a heating cycle, the stability of metal-incorporated GeTe (Ge 0.5–x M x Te 0.5 ) films is also critical to practical implementation of these materials in electronic and optoelectronic devices. In this work, we use both density-functional theory and experimental characterization methods to probe the solubility and critical properties of Ge 0.5–x M x Te 0.5 films. Using first-principles calculations, we determine the enthalpy of formation for GeTe with 2.08, 4.17, and 6.25 at. % of Cu, Fe, Mn, Mo, and Ti and show trends between the stability of the Ge 0.5–x M x Te 0.5 systems and the atomic position, composition, and distribution of the metal atoms in the GeTe matrix.

36 MATERIALS SCIENCE↗

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

Here, we use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100–200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

36 MATERIALS SCIENCE↗

Organic Photovoltaic Efficiency Predictor: Data-Driven Models for Non-Fullerene Acceptor Organic Solar Cells

In the design of organic solar cells, there has been a need for materials with high power conversion efficiencies. Scharber’s model is commonly used to predict efficiency; however, it exhibits poor performance with new non-fullerene acceptor (NFA) devices, since it was designed for fullerene-based devices. In this work, an empirical model is proposed that can be a more accurate alternative for NFA organic solar cells. Additionally, many screening studies use computationally expensive methods. A model based on using semiempirical simplified time-dependent density functional theory (sTD-DFT) as an alternative method can accelerate the calculations and yield a similar accuracy. The models presented in this paper, termed organic photovoltaic efficiency predictor (OPEP) models, have shown significantly lower errors than previous models, with OPEP/B3LYP yielding errors of 1.53% and OPEP/sTD- DFT of 1.55%. As a result, the proposed computational models can be used for the fast and accurate screening of new high-efficiency NFAs/donor pairs.

14 SOLAR ENERGY↗

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear level density and γ -decay strength of Sr 93

This work presents the first experimentally determined nuclear level density and γ-ray strength function of the short-lived fission product 93 Sr, accomplished using the β-Oslo method. Direct measurement of the 92 Sr(n, γ) 93 Sr cross section is not currently possible, as the half-life of 2.66 hours is too short; instead, 93 Sr was formed through β decay of 93 Rb to excitation energies around the neutron separation energy. The γ-ray spectra were measured using a total absorption spectrometer at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (MSU). The statistical properties of the 93 Sr nucleus were experimentally determined, including the γ-ray strength function and nuclear level density. At low energies, the γ-ray strength function exhibits a constant γ-decay strength, rather than a slightly increasing strength with decreasing γ-ray energy as had been previously observed for several nuclei in this mid-mass region. Finally, these statistical properties were then implemented in the reaction code TALYS1.95 to calculate the 92 Sr(n, γ) 93 Sr cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design and fabrication of TiO 2 -based dye sensitized solar cells using plant-derived organic dyes

The utilization of solar energy, an abundant and renewable resource, shows great promise. Various technologies have been employed to capture this vast energy potential. In this study, the effectiveness of natural dye-sensitized solar cells (DSSCs) was explored through experimental and computational methods. Organic dyes derived from specific plant species were examined, with a solvent mixture of ethanol, methanol, and tetrahydrofuran used for extraction. A UV–vis spectrophotometer was utilized to measure the sensitizers’ optical characteristics. The study employed density functional theory (DFT), implemented in Gaussian 09 W software, to perform computational calculations. The B3LYP method and a basis set of 6-31G++(d,p) were chosen for optimizing the geometries and energetics of dye molecules. The DFT results indicated that the dye molecules had a bandgap in the range of 2.16–2.38 eV. The photochemical performance of all four fabricated DSSC-based devices was evaluated to be between 0.33% and 1.04%, and the Euryops pectinatus sensitizer demonstrated the highest efficiency of 1.04% among the studied sensitizers.

36 MATERIALS SCIENCE↗

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hierarchical screening for Li-based solid electrolytes using fast, interpretable machine-learned potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates. Finally, we apply the hierarchical workflow to screen for ionic conductivity across a database of Li-containing compounds.

Materials discovery↗

Hierarchical Screening for Li-Based Solid Electrolytes Using Fast, Interpretable Machine-Learned Potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates.

Materials discovery↗

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural Fluctuations and Thermophysical Properties of Molten II-VI Compounds

The objectives of the project are to conduct ground-based experimental and theoretical research on the structural fluctuations and thermophysical properties of molten II-VI compounds to enhance the basic understanding of the existing flight experiments in microgravity materials science programs as well as to study the fundamental heterophase fluctuation phenomena in these melts by: 1) conducting neutron scattering analysis and measuring quantitatively the relevant thermophysical properties of the II-VI melts (such as viscosity, electrical conductivity, thermal diffusivity and density) as well as the relaxation characteristics of these properties to advance the understanding of the structural properties and the relaxation phenomena in these melts and 2) performing theoretical analyses on the melt systems to interpret the experimental results. All the facilities required for the experimental measurements have been procured, installed and tested. Thermal diffusivity of molten tellurium has been measured by a laser flash method in the temperature range of 500 C to 900 C. The measured diffusivity as a function of temperature agrees fairly well with published data. However, a relaxation phenomenon, which shows a slow drift of the measured thermal conductivity toward the equilibrium value after cooling of the melt, was observed for the first time. An apparatus based on the transient torque induced by a rotating magnetic field has been developed to determine the viscosity and electrical conductivity of semiconducting liquids. Viscosity measurements on molten tellurium showed a similar relaxation behavior to the measured diffusivity. The density and volume expansion coefficients for pure Te and HgTe melts were measured as a function of temperature using a pycnometric method. A density maximum was found for both melts but no relaxation behavior was observed. Neutron scattering experiments were performed on the HgTe and HgZnTe melts and the results on pair distribution showed better resolution than previously reported. A simple kinetic theory, which contains the formation reactions of Te polymerization that increases the molecular weight, was proposed to interpret the relaxation behavior of HgZnTe viscosity data.

Su, Ching-Hua↗

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗