Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

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

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2023

Measuring the Neutron Spin Asymmetry A1n in the Valence Quark Region in Hall C at Jefferson Lab

The quest to understand how the nucleon spin is decomposed into its constituent quark and gluon spin and orbital angular momentum (OAM) components has been at the forefront of nuclear physics for decades. Due to the non-perturbative nature of Quantum Chromodynamics (QCD) - the theory describing how quarks and gluons bind together to form protons and neutrons - making absolute predictions of nucleon spin structure is generally difficult, especially as a function of its quark and gluon longitudinal momentum fraction x. Measurements involving nucleon spin structure serve as a sensitive test for QCD, including ab-initio lattice QCD calculations due to the advent of the quasi-PDF formalism, and various predictions that diverge at large-x. The neutron spin asymmetry An 1 at high?x is a key observable for probing nucleon spin structure. In the valence domain (x > 0.5), sea effects are expected to be negligible, and so the total nucleon spin is considered to be carried by the valence quarks. The valence region can therefore enable us to study the role of quark OAM and other non-perturbative effects of the strong force. An 1 was measured in the deep inelastic scattering region of 0.40 < x < 0.75 and 6 < Q2 < 10 GeV2 in Hall C at Jefferson Lab using a 10.4 GeV longitudinally polarized electron beam, upgraded polarized 3He target, and the High Momentum Spectrometer (HMS) and Super High Momentum Spectrometer (SHMS). E12- 06-110 provides the first precision data in the valence quark region above x = 0.60, and its preliminary results proved consistent with earlier data disqualifying a pQCD model that excluded quark OAM. Combined with previous world proton data, the ratio of the polarized-to-unpolarized up quark momentum distribution (?u + ?u)/(u + u) remained positive at large-x, and the down quark (?d + ?d)/(d + d) remained negative.

Rehfuss, Melanie↗
DOE OSTI2020

Probing the Local Site Disorder and Distortion in Pyrochlore High-Entropy Oxides

High-entropy oxides (HEOs) have attracted great interest in diverse fields because of their inherent opportunities to tailor and combine materials functionalities. The control of local order/disorder in the class is by extension a grand challenge toward realizing their vast potential. In this paper, we report the first examples of pyrochlore HEOs with five M-site cations, for Nd 2 M 2 O 7 , in which the local structure has been investigated by neutron diffraction and pair distribution function (PDF) analysis. The average structure of the pyrochlores is found to be orthorhombic Imma, in agreement with radius-ratio rules governing the structural archetype. The computed PDFs from density functional theory relaxed special quasirandom structure models are compared with real space PDFs in this work to evaluate M-site order/disorder. Reverse Monte Carlo combined with ab initio molecular dynamics and Metropolis Monte Carlo simulations demonstrates that Nd 2 (Ta 0.2 Sc 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7 is synthesized with its M-site local to nanoscale order highly randomized/disordered, while Nd 2 (Ti 0.2 Nb 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7+x exhibits a strong distortion of the TiO 6 octahedron and small degree of Ti chemical short-range order (SRO) on the subnanometer scale. Calculations suggest that this may be intrinsic, energetically favored SRO rather than due to sample processing. These results offer an important demonstration that the engineered variation of participating ions in HEOs, even among those with very similar radii, provides richly diverse opportunities to control local order/disorder motifs—and therefore materials properties for future designs. This work also hints at the exquisite level of detail that may be needed in computational and experimental data analysis to guide structure–property tuning in the emerging HEO materials class.

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

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H{sub 2}O){sub 128} significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, S{sub OO}(Q), and corresponding oxygen-oxygen radial distribution function, g{sub OO}(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, P{sub OOO}(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.

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

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2600+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6,000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with great emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with a suite of experimental data obtained from x-ray diffraction, infrared spectroscopy, spectroscopic ellipsometry and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Neutrino Flavor Model Building and the Origins of Flavor and CP Violation

The neutrino sector offers one of the most sensitive probes of new physics beyond the Standard Model of Particle Physics (SM). The mechanism of neutrino mass generation is still unknown. The observed suppression of neutrino masses hints at a large scale, conceivably of the order of the scale of a rand unified theory (GUT), which is a unique feature of neutrinos that is not shared by the charged fermions. The origin of neutrino masses and mixing is part of the outstanding puzzle of fermion masses and mixings, which is not explained ab initio in the SM. Flavor model building for both quark and lepton sectors is important in order to gain a better understanding of the origin of the structure of mass hierarchy and flavor mixing, which constitute the dominant fraction of the SM parameters. Recent activities in neutrino flavor model building based on non-Abelian discrete flavor symmetries and modular flavor symmetries have been shown to be a promising direction to explore. The emerging models provide a framework that has a significantly reduced number of undetermined parameters in the flavor sector. In addition, such a framework affords a novel origin of CP violation from group theory due to the intimate connection between physical CP transformation and group theoretical properties of non-Abelian discrete groups. Model building based on non-Abelian discrete flavor symmetries and their modular variants enables the particle physics community to interpret the current and anticipated upcoming data from neutrino experiments. Non-Abelian discrete flavor symmetries and their modular variants can result from compactification of a higher-dimensional theory. Pursuit of flavor model building based on such frameworks thus also provides the connection to possible UV completions: in particular, to string theory. We emphasize the importance of constructing models in which the uncertainties of theoretical predictions are smaller than, or at most compatible with, the error bars of measurements in neutrino experiments. While there exist proof-of-principle versions of bottom-up models in which the theoretical uncertainties are under control, it is remarkable that the key ingredients of such constructions were discovered first in top-down model building. We outline how a successful unification of bottom-up and top-down ideas and techniques may guide us towards a new era of precision flavor model building in which future experimental results can give us crucial insights into the UV completion of the SM.

Almumin, Yahya (ORCID:0000000196396772)↗
DOE OSTI2019

Ab initio folding potentials for nucleon-nucleus scattering based on no-core shell-model one-body densities

Background: Calculating microscopic optical potentials for elastic nucleon-nucleus scattering has already led to large body of work in the past. For folding first-order calculations the nucleon-nucleon ( NN ) interaction and the one-body density of the nucleus were taken as input to rigorous calculations in a spectator expansion of the multiple scattering series. Purpose: Based on the Watson expansion of the multiple scattering series we employ a nonlocal translationally invariant nuclear density derived from a chiral next-to-next-to-leading order (NNLO) and the very same interaction for consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for light nuclei. Methods: The first order effective (optical) folding potential is computed by integrating over the nonlocal, translationally invariant NCSM one-body density and the off-shell Wolfenstein amplitudes A and C. The resulting nonlocal potential serves as input for a momentum-space Lippmann-Schwinger equation, whose solutions are summed to obtain the nucleon-nucleus scattering observables. Results: We calculate scattering observables, such as total, reaction, and differential cross sections as well as the analyzing power and the spin-rotation parameter, for elastic scattering of protons and neutrons from 4 He, 6 He, 12 C, and 16 O, in the energy regime between 100 and 200 MeV projectile kinetic energy, and compare to available data. Conclusions: Our calculations show that the effective nucleon-nucleus potential obtained from the first-order term in the spectator expansion of the multiple scattering expansion describes experiments very well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the NN amplitudes and the one-body nuclear density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Synthesis and thermophysical property determination of NaCl-PuCl3 salts

Currently, a knowledge gap exists in the available data and understanding of thermophysical properties relating to fresh fuel salts, especially those containing plutonium. These data are necessary for designing, constructing, and licensing future commercial molten-salt reactors. Thermophysical properties, such as melting temperature, salt stability, density, and heat capacity were ascertained using NaCl-PuCl3 (36 mol% PuCl3) and a more sodium rich composition containing 25 mol% PuCl3. The NaCl-PuCl3 salt mixture was synthesized for this study and contained 63.4 mol% NaCl, 36.3 mol% PuCl3 and was 99.7% pure. Upon heating, the NaCl–PuCl3 eutectic was stable at temperatures up to 800 °C. The onset of melting occurred at 451 ± 3 °C, and the enthalpy of fusion was determined to be 23.0 ± 1.4 kJ/mol. Heat capacity measurements in the liquid phase ranged from 107.7 to 91.3 J/mol.K, with an average value of 104.6 ± 11.4 J/mol.K between 500 and 650 °C. Three independent trials of the molten NaCl-PuCl3 salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C). In addition, ab initio molecular dynamic simulations to calculate density and heat capacity values are included.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Toward a W4-F12 approach: Can explicitly correlated and orbital-based ab initio CCSD(T) limits be reconciled?

In the context of high-accuracy computational thermochemistry, the valence coupled cluster with all singles and doubles (CCSD) correlation component of molecular atomization energies presents the most severe basis set convergence problem, followed by the (T) component. In the present paper, we make a detailed comparison, for an expanded version of the W4-11 thermochemistry benchmark, between, on the one hand, orbital-based CCSD/AV{5,6}Z + d and CCSD/ACV{5,6}Z extrapolation, and on the other hand CCSD-F12b calculations with cc-pVQZ-F12 and cc-pV5Z-F12 basis sets. This latter basis set, now available for H–He, B–Ne, and Al–Ar, is shown to be very close to the basis set limit. Apparent differences (which can reach 0.35 kcal/mol for systems like CCl{sub 4}) between orbital-based and CCSD-F12b basis set limits disappear if basis sets with additional radial flexibility, such as ACV{5,6}Z, are used for the orbital calculation. Counterpoise calculations reveal that, while total atomization energies with V5Z-F12 basis sets are nearly free of BSSE, orbital calculations have significant BSSE even with AV(6 + d)Z basis sets, leading to non-negligible differences between raw and counterpoise-corrected extrapolated limits. This latter problem is greatly reduced by switching to ACV{5,6}Z core-valence basis sets, or simply adding an additional zeta to just the valence orbitals. Previous reports that all-electron approaches like HEAT (high-accuracy extrapolated ab-initio thermochemistry) lead to different CCSD(T) limits than “valence limit + CV correction” approaches like Feller-Peterson-Dixon and Weizmann-4 (W4) theory can be rationalized in terms of the greater radial flexibility of core-valence basis sets. For (T) corrections, conventional CCSD(T)/AV{Q,5}Z + d calculations are found to be superior to scaled or extrapolated CCSD(T)-F12b calculations of similar cost. For a W4-F12 protocol, we recommend obtaining the Hartree-Fock and valence CCSD components from CCSD-F12b/cc-pV{Q,5}Z-F12 calculations, but the (T) component from conventional CCSD(T)/aug’-cc-pV{Q,5}Z + d calculations using Schwenke’s extrapolation; post-CCSD(T), core-valence, and relativistic corrections are to be obtained as in the original W4 theory. W4-F12 is found to agree slightly better than W4 with ATcT (active thermochemical tables) data, at a substantial saving in computation time and especially I/O overhead. A W4-F12 calculation on benzene is presented as a proof of concept.

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

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2013

High-Accuracy Quartic Force Field Calculations for the Spectroscopic Constants and Vibrational Frequencies of 1(exp 1)A' l-C3H(-): A Possible Link to Lines Observed in the Horsehead Nebula PDR

It has been shown that rotational lines observed in the Horsehead nebula photon-dominated-region (PDR) are probably not caused by l-C3H+, as was originally suggested. In the search for viable alternative candidate carriers, quartic force fields are employed here to provide highly accurate rotational constants, as well as fundamental vibrational frequencies, for another candidate carrier: 1 (sup 1)A' C3H(−). The ab initio computed spectroscopic constants provided in this work are, compared to those necessary to define the observed lines, as accurate as the computed spectroscopic constants for many of the known interstellar anions. Additionally, the computed D-eff for C3H(−) is three times closer to the D deduced from the observed Horsehead nebula lines relative to l-C3H(+). As a result, 1 (sup 1)A' C3H(−). is a more viable candidate for these observed rotational transitions and would be the seventh confirmed interstellar anion detected within the past decade and the first C(sub n)H(−) molecular anion with an odd n.

Horsehead nebula↗
DOE OSTI2016

Experimental and molecular dynamics simulation study of structure of liquid and amorphous Ni 62 Nb 38 alloy

The state-of-the-art experimental and atomistic simulation techniques were utilized to study the structure of the liquid and amorphous Ni 62 Nb 38 alloy. First, the ab initio molecular dynamics (AIMD) simulation was performed at rather high temperature where the time limitations of the AIMD do not prevent to reach the equilibrium liquid structure. A semi-empirical potential of the Finnis-Sinclair (FS) type was developed to almost exactly reproduce the AIMD partial pair correlation functions (PPCFs) in a classical molecular dynamics simulation. This simulation also showed that the FS potential well reproduces the bond angle distributions. The FS potential was then employed to elongate the AIMD PPCFs and determine the total structure factor (TSF) which was found to be in excellent agreement with X-ray TSF obtained within the present study demonstrating the reliability of the AIMD for the simulation of the structure of the liquid Ni–Nb alloys as well as the reliability of the developed FS potential. The glass structure obtained with the developed potential was also found to be in excellent agreement with the X-ray data. We report the analysis of the structure revealed that a network of the icosahedra clusters centered on Ni atoms is forming during cooling the liquid alloy down to T g and the Nb Z14, Z15, and Z16 clusters are attached to this network. This network is the main feature of the Ni 62 Nb 38 alloy and further investigations of the properties of this alloy should be based on study of the behavior of this network.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Ab initio many-body calculations of nucleon- 4 He scattering with three-nucleon forces

We extend the ab initio no-core shell model/resonating-group method to include three-nucleon (3N) interactions for the description of nucleon-nucleus collisions. We outline the formalism, give algebraic expressions for the 3N-force integration kernels, and discuss computational aspects of two alternative implementations. The extended theoretical framework is then applied to nucleon- 4 He elastic scattering using similarity-renormalization-group (SRG)-evolved nucleon-nucleon plus 3N potentials derived from chiral effective field theory. We analyze the convergence properties of the calculated phase shifts and explore their dependence upon the SRG evolution parameter. We include up to six excited states of the 4 He target and find significant effects from the inclusion of the chiral 3N force, e.g., it enhances the spin-orbit splitting between the 3/2 – and 1/2 – resonances and leads to an improved agreement with the phase shifts obtained from an accurate R-matrix analysis of the five-nucleon experimental data. As a result, we find remarkably good agreement with measured differential cross sections at various energies below the d+ 3 H threshold, while analyzing powers manifest larger deviations from experiment for certain energies and angles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Competing Pathways for Nucleation of the Double Perovskite Structure in the Epitaxial Synthesis of La 2 MnNiO 6

Over the past decades a confluence of advanced synthesis techniques, data-driven characterization, and rapid increases in computing power has sparked a renaissance in materials engineering. The emerging “materials-by-design” approach, while radically transforming the development of multicomponent systems, has tended to overlook the complex kinetic pathways that define materials synthesis. Although we are able to envision almost limitless materials combinations, we are unable to synthesize all of them in practice, since existing characterization and modeling approaches often fail to capture the inherent complexity of such systems. There is currently a disconnect between highly local structural characterization and macroscale properties measurements, resulting in oversimplified or incomplete structure-property models. Here we describe a multi-technique approach that combines aberration-corrected transmission electron microscopy with emerging oxide atom probe tomography to measure chemical ordering and extended defects in the model complex oxide, La2MnNiO6. We visualize cation ordering, as well as a three-dimensional network of secondary phases, which we describe in terms of ab initio structure calculations. We propose a defect model in which these phases alter their surrounding octahedral environment, severely disrupting cation superexchange. Through this array of experimental and theoretical techniques, we uncover fundamental structure-property relationships and illustrate a new approach to engineer complex, multicomponent systems.

atom probe tomography↗
DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗