Method for determining the thermodynamic properties of the quantum electron gas.
Effective pair potential obtained for quantum electron gas to determine thermodynamic properties over large temperature and density range
SEARCH · Engineering Papers
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.
Effective pair potential obtained for quantum electron gas to determine thermodynamic properties over large temperature and density range
Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.
Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.
In this work, we investigate magnetic monolayers of the form A i A ii B 4 X 8 based on the well-known intrinsic topological magnetic van der Waals (vdW) material MnBi 2 Te 4 (MBT) using first-principles calculations and machine learning techniques. We select an initial subset of structures to calculate the thermodynamic properties, electronic properties, such as the band gap, and magnetic properties, such as the magnetic moment and magnetic order using density functional theory (DFT). Data analytics approaches are used to gain insight into the microscopic origin of materials’ properties. The dependence of materials’ properties on chemical composition is also explored. For example, we find that the formation energy and magnetic moment depend largely on A and B sites whereas the band gap depends on all three sites. Finally, we employ machine learning tools to accelerate the search for novel vdW magnets in the MBT family with optimized properties. Finally, this study creates avenues for rapidly predicting novel materials with desirable properties that could enable applications in spintronics, optoelectronics, and quantum computing.
Abstract The electronic properties and thermal stabilities of MAlO 2 and M 5 AlO 4 (M = Li, Na, K) are investigated by density functional theory and lattice phonon dynamics. Based on the calculated electronic and lattice thermodynamic properties, their abilities to capture CO 2 as solid sorbents are analyzed. The calculated electronic structural properties of MAlO 2 and M 5 AlO 4 indicate that all these alkali aluminates are semiconductors with a bandgap range of 2.4 ~ 6.4 eV. The 1st valence bands of these alkali aluminates are located 0 ~ − 6 eV under Fermi levels and are mainly contributed by p orbitals of O, s and p orbitals of Al and M. The phonon vibrational frequencies of M 5 AlO 4 spread at a lower frequency range compared to their MAlO 2 phases. With increasing temperature, the calculated phonon free energies of M 5 AlO 4 decrease faster than their corresponding MAlO 2 while their entropies have opposite trends. The reaction 2MAlO 2 + CO 2 = M 2 CO 3 + Al 2 O 3 has higher reaction heat and Gibbs free energy change than those of corresponding reaction 2 / 5 M 5 AlO 4 + CO 2 = M 2 CO 3 + 1 / 5 Al 2 O 3 , which shows the former reaction possesses lower turnover temperature. Among the alkali aluminates studied, the β-NaAlO 2 , lt-KAlO 2 , and γ-LiAlO 2 are better candidates that could be applied for CO 2 capture technologies. Graphical Abstract
Here, this work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate electron relaxation time and its temperature dependence. The model is applied to Ca 3 Ru 2 O 7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity.
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.
Kagome metals are a unique class of quantum materials characterized by their distinct atomic lattice arrangement, featuring interlocking triangles and expansive hexagonal voids. These lattice structures impart exotic properties, including superconductivity, interaction-driven topological many-body phenomena, and magnetism, among others. The kagome metal CsM 3 Te 5 (where M = Ti, Zr, or Hf) exhibits both superconductivity and nontrivial topological electronic properties, offering a promising platform for exploring topological superconductivity. This study employs first-principles density functional theory calculations to systematically analyze the elastic, mechanical, vibrational, thermodynamic, and electronic properties of CsM 3 Te 5 (M = Ti, Zr, Hf). Our calculations reveal that the studied compounds—CsTi 3 Te 5 , CsZr 3 Te 5 , and CsHf 3 Te 5 —are ductile metals with elastic properties akin to the hexagonal Bi and Sb, with average elastic constants, including a bulk modulus of 27 GPa, a shear modulus of 11 GPa, and Young's modulus of 29 GPa. We observe peculiar dispersionless, flat, phonon branches in the vibrational spectra of these metals. Additionally, we thoroughly analyze the symmetries of the zone-center phonon eigenvectors and predict vibrational fingerprints of the Raman- and infrared-active phonon modes. The analysis of thermodynamic properties reveals the Einstein temperature for CsTi 3 Te 5 , CsZr 3 Te 5 , and CsHf 3 Te 5 to be 66, 54, and 53 K, respectively. Our orbital-decomposed electronic structure calculations reveal significant in-plane steric interactions and multiple Dirac band crossings near the Fermi level. We further investigate the role of spin-orbit coupling effect on the studied properties. Furthermore, this theoretical investigation sheds light on the intriguing quantum behavior of kagome metals.
Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.
Double transition metal (DTM) nitride MXenes offer enhanced electrical conductivity and tunable catalytic properties compared to conventional carbide-based MXenes. In this study, we employed first-principles density functional theory (DFT) calculations to discover and computationally validate a novel DTM nitride MXene, Nb 2 TiN 2 , derived from its MAX phase precursor and investigated its potential as an anchoring material (AM) for Li–Se batteries. This newly proposed MXene expands the compositional landscape of DTM nitrides and opens new avenues for functional material design. We performed a comprehensive analysis of the thermodynamic and electronic properties of Nb 2 TiAlN 2 , and the MAX phase precursor to Nb 2 TiN 2 to assess its structural stability and exfoliation potential. Exfoliation energy calculations confirmed the feasibility of synthesizing Nb 2 TiN 2 from Nb 2 TiAlN 2 . We then explored the functionalized form, Nb 2 TiN 2 S 2 , evaluating its capability to serve as an effective anchoring material (AM) in Li–Se batteries by analyzing the reaction mechanisms and kinetics of the selenium reduction reaction (SeRR). Our results indicate that Nb 2 TiN 2 S 2 exhibits a strong binding affinity for lithium polyselenides (Li 2 Se n ), effectively suppressing the shuttle effect. Gibbs free energy calculations for the rate-limiting step of the SeRR reveal favorable kinetics and reduced reaction barriers. Overall, this study provides a detailed evaluation of the structural and electronic properties of a newly proposed DTM nitride MXene and its S-functionalized derivative and the catalyzing effect of Nb 2 TiN 2 S 2 in accelerating the reaction kinetics in Li–Se batteries. These findings underscore the potential importance of the further exploration of MXenes to address current challenges in high-performance Li–Se batteries.
In this work, we have investigated the thermodynamic stability of disordered rare-earth phases SmX 2 and Sm 10 X 21 (X=Ni, Pd) using machine-learning based analytical descriptor and first-principles density functional theory methods. The absence of Laves phase compounds in R-Pd binary systems is a longstanding problem of rare earth science: even though Ni and Pd belong to the same group of the periodic table and have similar electronic structure, the Pd compound crystallizes in a monoclinic (C2/m) phase with 10:21 stoichiometry, i.e., Sm 10 Pd 21 , while the Ni compound adopts a cubic Laves phase (MgCu 2 ) structure. To understand this contrasting phase stability, we performed thermodynamic convex hull analysis of Sm x Ni 1-x and Sm x Pd 1-x binary systems, which is experimentally validated using powder X-ray diffraction (PXRD) analyzes of polycrystalline Sm(Ni x Pd 1-x ) 2 samples with x=0, 0.5, and 1. A detailed electronic-structure (band-structure, charge density, and Fermi-surface) analysis of the differences between SmNi 2 /SmPd 2 and Sm 10 Ni 21 /Sm 10 Pd 21 compounds provides the quantum mechanical origin of the unfavorable mixing of Pd with Ni in cubic Laves phase. We show that the stability of Sm-Pd in 10:21 stoichiometry arises from improved intra-/inter-layer 5d-4d bonding compared to the 1:2 stoichiometry. Our work emphasizes the importance of ab-initio methods and computationally inexpensive analytical descriptors for the detailed analysis of thermodynamic and electronic properties of hard-to-prepare rare-earth compounds.
The Fermi surface plays an important role in controlling the electronic, transport and thermodynamic properties of materials. As the Fermi surface consists of closed contours in the momentum space for well-defined energy bands, disconnected sections known as Fermi arcs can be signatures of unusual electronic states, such as a pseudogap. Another way to obtain Fermi arcs is to break either the time-reversal symmetry or the inversion symmetry of a three-dimensional Dirac semimetal, which results in formation of pairs of Weyl nodes that have opposite chirality, and their projections are connected by Fermi arcs at the bulk boundary. Here, we present experimental evidence that pairs of hole- and electron-like Fermi arcs emerge below the Neel temperature (TN) in the antiferromagnetic state of cubic NdBi due to a new magnetic splitting effect. The observed magnetic splitting is unusual, as it creates bands of opposing curvature, which change with temperature and follow the antiferromagnetic order parameter. This is different from previous theoretically considered and experimentally reported cases of magnetic splitting, such as traditional Zeeman and Rashba, in which the curvature of the bands is preserved. Therefore, our findings demonstrate a type of magnetic band splitting in the presence of a long-range antiferromagnetic order that is not readily explained by existing theoretical ideas.
First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.
The interplay of dynamical correlations and electronic ordering is pivotal in shaping phase diagrams of correlated quantum materials. In magic-angle twisted bilayer graphene, transport, thermodynamic, and spectroscopic experiments pinpoint a competition between distinct low-energy states with and without electronic order, as well as between localized and delocalized charge carriers. In this study, we utilize dynamical mean-field theory on the topological heavy fermion model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We contrast moment formation, Kondo screening, and ordering on a temperature basis and explain the nature of emergent correlated states based on three central phenomena: (i) the formation of local spin and valley isospin moments around 100 K, (ii) the ordering of the local isospin moments around 10 K preempting Kondo screening, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low-energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a “bad metal” above the ordering temperature, where scattering off the disordered local moments suppresses electronic coherence, and a “good metal” in the ordered states with coherence of quasiparticles facilitated by isospin order. This finding reveals coherence from order as the microscopic mechanism behind the Pomeranchuk effect observed experimentally by Rozen [] and by Saito []. Upon doping, there is a periodic charge reshuffling between localized and delocalized electronic orbitals leading to cascades of doping-induced Lifshitz transitions, local spectral weight redistributions, and periodic variations of the electronic compressibility ranging from nearly incompressible to negative. Our findings highlight the essential role of charge transfer, hybridization, and ordering in shaping the electronic excitations and thermodynamic properties in twisted bilayer graphene and provide a unified understanding of the most puzzling aspects of scanning tunneling spectroscopy, transport, and compressibility experiments. Published by the American Physical Society 2024
This final technical report summarizes the research accomplishments and research highlights at the end of the funding period. This project aimed to leverage existing and new computational data produced from first-principles and molecular dynamics simulations to understand the thermodynamic, mechanical, and electronic properties of icosahedral boron compounds. The goal is to achieve targeted material properties by controlling the synthesis routes of these boron-rich compounds.
Owing to its novel electronic and magnetic properties, two-dimensional CrI 3 has great potential in the application of spintronic devices. However, as an inevitable line defect, the properties of the edges of CrI 3 remain elusive. Here, via first-principles calculations with spin–orbit coupling, we investigated the thermodynamic stabilities, electronic and magnetic properties of thirteen CrI 3 edges with different structures. We showed that zigzag edges are more stable than armchair edges, and a CrI 3 nanoribbon can be either metallic or insulating depending on its chemical growth conditions. Here, the edge stability and associated electronic properties can be understood in terms of the octahedron ligand field and electron counting model. In most cases, both the magnetic moment and Curie temperature can be enhanced by edges, which are in startle contrast to the surfaces of three-dimensional ferromagnetic materials, where a magnetic dead layer is often observed.
Metals with partially filled core atomic shells can form quasiparticles at a low temperature arising from the hybridization of the core level and conduction electrons. The thermodynamic and spectroscopic properties of these metals can be understood as those of a simple metal, but with a significant mass enhancement over the free electron mass—commonly referred to as heavy fermions. In most heavy-fermion materials, the hybridization is approximately isotropic in position and momentum space. However, a combination of low dimensionality and symmetry properties of the core-level wavefunctions can give rise to highly anisotropic electronic interactions with the conduction electrons. Here, in this study, we demonstrate anisotropic hybridization that vanishes along specific directions in momentum space—referred to as nodes—in a lanthanide-based two-dimensional van der Waals heavy-fermion compound, CeSiI. Quasiparticle interference measurements reveal a set of discrete hotspots with high spectral intensity on the Fermi surface. Theoretical modelling and comparison with the quasiparticle interference pattern of the non-heavy-fermion isostructural analogue LaSiI suggest that these features arise from an unconventional electron interaction involving hybridization nodes unique to CeSiI. As a result, the effective mass of the quasiparticles varies by orders of magnitude depending on their direction in momentum space.
Abstract We present a benchmark test suite and an automated machine learning procedure for evaluating supervised machine learning (ML) models for predicting properties of inorganic bulk materials. The test suite, Matbench, is a set of 13 ML tasks that range in size from 312 to 132k samples and contain data from 10 density functional theory-derived and experimental sources. Tasks include predicting optical, thermal, electronic, thermodynamic, tensile, and elastic properties given a material’s composition and/or crystal structure. The reference algorithm, Automatminer, is a highly-extensible, fully automated ML pipeline for predicting materials properties from materials primitives (such as composition and crystal structure) without user intervention or hyperparameter tuning. We test Automatminer on the Matbench test suite and compare its predictive power with state-of-the-art crystal graph neural networks and a traditional descriptor-based Random Forest model. We find Automatminer achieves the best performance on 8 of 13 tasks in the benchmark. We also show our test suite is capable of exposing predictive advantages of each algorithm—namely, that crystal graph methods appear to outperform traditional machine learning methods given ~10 4 or greater data points. We encourage evaluating materials ML algorithms on the Matbench benchmark and comparing them against the latest version of Automatminer.