Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ethane Dehydrogenation on Single and Dual Centers of Ga-modified γ-Al 2 O 3

Density-functional theory calculations and microkinetic analysis are used to investigate the efficacy of Ga-modified γ-Al 2 O 3 (110) surfaces for the catalytic dehydrogenation of ethane and elucidate the synergy between Ga and Al sites. The model surfaces are modified by either Ga grafting or doping. We consider and analyze numerous active sites and rank them using microkinetic analysis. The kinetic parameters obtained from microkinetic modelling are compared with experimental values for ethane dehydrogenation over Ga 2 O 3 -Al 2 O 3 mixed oxides prepared by co-precipitation. The dominant reaction pathway proceeds via heterolytic C-H bond dissociation to a surface proton and a metal-carbanion intermediate that undergoes β-hydride elimination. We find that grafted Ga sites are catalytically inactive. In contrast, Ga-doped sites exhibit fivefold enhancement in catalytic activity when compared to the sites on pristine Al 2 O 3 , owed to the synergy between neighboring AlIII and GaIV sites. Furthermore, we model and investigate the effect of surface hydroxylation, demonstrate how surface water interferes with the aforementioned synergy between AlIII and GaIV sites and discuss the implications for the catalytic activity of the modified surfaces. Here, an increase in the partial pressure of H 2 O significantly increases the apparent activation energies of dehydrogenation and interestingly changes the most active site.

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

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗
DOE OSTI2023

Migration of Ga vacancies and interstitials in β – Ga 2 O 3

Pathways and energy barriers for the migration of Ga vacancies (V Ga ) and Ga interstitials (Ga i ) in β–Ga 2 O 3 are explored using hybrid functional calculations and the nudged elastic band method. Considering β–Ga 2 O 3 as primarily being an n-type semiconductor, we focus on defect charge states relevant under such conditions: $V^{3 –}_{Ga}$, Ga$^{3+}_{i}$, and Ga$^{+}_{i}$. Notably, we describe a mechanism by which V Ga can transform between its different split configurations. In all cases, the intermediate state consists of a vacancy split between three Ga sites—a three-split vacancy—which enables passage over a significantly lower energy barrier. This is because it avoids the unfavorable simple vacancy at the tetrahedral Ga site. The proposed mechanism lowers the overall barrier for $V^{3–}_{Ga}$ diffusion along the [001] crystal direction from 1.73 to 0.97 eV, whereas the 2.08 eV barrier for the [100] and [010] directions is unaffected. For Ga$^{3+}_{i}$, we obtain similar overall migration barriers of 0.72, 0.80, and 1.02 eV for the [010], [001], and [100] directions, respectively. Ga$^{+}_{i}$ exhibits a strong preference for diffusion within the large eight-sided channel; the overall migration barrier is 0.92 eV for the [010] direction, and 2.16 eV for the [001] and [100] directions. The limiting step for the two latter directions involves ionization of Ga$^{+}_{i}$ followed by a jump to an adjacent large eight-sided channel as Ga$^{3+}_{i}$. Finally, our results are discussed in light of experimental observations of thermally activated recovery processes in irradiated material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A new metric to control nucleation and grain size distribution in hybrid organic–inorganic perovskites by tuning the dielectric constant of the antisolvent

In perovskite research, there is a widely exploited but poorly explained phenomenon in which the addition of “antisolvents (ATS)” to precursor solutions results in higher-quality films. We explain the mechanism and driving force underlying an antisolvent-driven solvent extraction process. Density functional theory calculations uncover the defining effects of antisolvent choice on the extent of complexation between a lead salt and a methylammonium cation in solution. We experimentally validate the computational results using ultraviolet-visible spectroscopy and 207 Pb nuclear magnetic spectroscopy of methylammonium lead iodide solutions, containing both a processing solvent and an antisolvent. Furthermore, we uncover, and subsequently identify, the appearance of new species in solution as a result of the addition of the antisolvent. Furthermore, we observe that the choice of antisolvent has a substantial effect on the nature of the complexation of the methylammonium lead iodide (MAPbI 3 ) precursor species, whose origin we explain at an atomic level; specifically, the lower the dielectric of the antisolvent, the stronger the intermolecular binding energy between methylammonium cation (MA + ) cation and PbI 3 - plumbate, independent of the solvent or antisolvent interaction with the lead salt. Thin films were characterized using scanning electron microscopy; images of the films show how the addition of an antisolvent influences and, importantly, can be used to alter thin-film grain size. Grain size and distribution in thin films is reflected by the choice of antisolvent, promoting slower nucleation rates, a lower nucleation density, and hence larger final grain size.

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

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Chern numbers of topological phonon band crossing determined with inelastic neutron scattering

Topological invariants in the band structure, such as Chern numbers, are important for the understanding and classification of the topological properties of matter and dictate the occurrence of exotic behaviors, yet their direct spectroscopic determination has been largely limited to electronic bands. Here, we use inelastic neutron scattering in conjunction with ab initio calculations to identify a variety of topological phonon band crossings in MnSi and CoSi single crystals. We find a distinct relation between the Chern numbers of a band-crossing node and the scattering intensity modulation in momentum space around the node. Given sufficiently high resolution, our method can be used to determine arbitrarily large Chern numbers of topological phonon band-crossing nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Tuning the structure and superconductivity of SrNi 2⁢ P 2 by Rh substitution

The compound SrNi 2⁢ P 2 is unique among the ThCr 2 ⁢Si 2 class since it exhibits a temperature-induced transition upon cooling from an uncollapsed tetragonal (ucT) state to a one-third-collapsed orthorhombic (tcO) state where one out of every three P-rows bond across the Sr layers. This compound is also known for exhibiting bulk superconductivity below 1.4 K at ambient pressure. Here, in this paper, we report on the effects of Rh substitution in Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 on the structural and superconducting properties. We studied the variation of the nearest P-P distances as a function of the Rh fraction at room temperature, as well as its temperature dependence for selected compositions. We find that increasing the Rh fraction leads to a decrease in the transition temperature between the ucT and tcO states, until a full suppression of the tcO state for x ≥ 0.166. The superconducting transition first remains nearly insensitive to the Rh fraction, and then it increases to 2.3 K after the tcO state is fully suppressed. These results are summarized in a phase diagram, built upon the characterization by energy dispersive x-ray spectroscopy, x-ray diffraction, resistance, magnetization, and specific heat measurements done on crystalline samples with varying Rh content. The relationship between band structure, crystal structure, and superconductivity is discussed based on previously reported band structure calculations on SrRh 2 ⁢P 2 . Moreover, the effect of Rh fraction on the stress-induced structural transitions is also addressed by means of strain-stress studies done by uniaxial compression of single-crystalline micropillars of Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

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

Energy Migration Processes in Re(I) MLCT Complexes Featuring a Chromophoric Ancillary Ligand

We present herein the synthesis, structural characterization, electronic structure calculations, and ultrafast and supra-nanosecond photophysical properties of a series of five Re(I) bichromophores exhibiting metal to ligand charge transfer (MLCT) excited states based on the general formula fac -[Re(N ^ N)(CO) 3 (PNI-py)]PF 6, where PNI-py is 4-piperidinyl-1,8-naphthalimidepyridine and N ^ N is a diimine ligand ( Re1 – 5 ), along with their corresponding model chromophores where 4-ethylpyridine was substituted for PNI-py ( Mod1 – 5 ). The diimine ligands used include 1,10-phenanthroline (phen, 1 ), 2,9-dimethyl-4,7-diphenyl-1,10-phenanthroline (bcp, 2 ), 4,4'-di- tert -butyl-2,2'-bipyridine (dtbb, 3 ), 4,4'-diethyl ester-2,2'-bipyridine (deeb, 4 ), and 2,2'-biquinoline (biq, 5 ). In these metal–organic bichromophores, structural modification of the diimine ligand resulted in substantial changes to the observed energy transfer efficiencies between the two chromophores as a result of the variation in 3 MLCT excited-state energies. The photophysical properties and energetic pathways of the model chromophores were explored in parallel to accurately track the changes that arose from introduction of the organic chromophore pendant on the ancillary ligand. All relevant photophysical and energy transfer processes were probed and characterized using time-resolved photoluminescence spectroscopy, ultrafast and nanosecond transient absorption spectroscopy, and time-dependent density functional theory calculations. Of the five bichromophores in this study, four ( Re1 – 4 ) exhibited a thermal equilibrium between the 3 PNI-py and the 3 MLCT excited state, drastically extending the lifetimes of the parent model chromophores.

14 SOLAR ENERGY↗
DOE OSTI2021

Unraveling the Intermediate Reaction Complexes and Critical Role of Support-Derived Oxygen Atoms in CO Oxidation on Single-Atom Pt/CeO 2

CeO 2 -supported Pt single-atom catalysts have been extensively studied due to their relevance in automobile emission control and for the fundamental understanding of CeO 2 -based catalysts. Though CeO 2 -supported Pt nanoparticles are often more active than their single-atom counterparts, the former could easily redisperse to Pt single atom under oxidizing diesel conditions. Therefore, to maximize the reactivity of every Pt atom, it is important to fully understand the reaction mechanism of CeO 2 -supported Pt single atoms. Here, we report a CO oxidation study on a Pt/CeO 2 single-atom catalyst, where we can account for all of the neighbors using in situ and operando spectroscopy techniques and microcalorimetric measurements. Coupled with density functional theory calculations, we present a comprehensive picture of the dynamics of the surface species, the role of surface intermediates, and explain the observed reaction kinetics. We started with a catalyst containing exclusively single atoms and used in situ/operando spectroscopy to provide evidence for their stability during the reaction and to identify the Pt 1 complexes before and during the reaction and their binding to CO. The results reveal that in the precatalyst, Pt is present as Pt(O) 4 on the CeO 2 (111) step edge sites, but during CO oxidation, we find that two Pt 1 complexes coexist, representing two states of the same active site in the reaction cycle. The dominant state/complex remains Pt(O) 4 , which adsorbs CO very weakly as shown by CO microcalorimetry. The second, minority state/complex, Pt(CO)(O) 3 is generated through the reaction of Pt(O) 4 with CO, and CO is bound strongly to Pt 1 . Labile oxygen adatoms from the CeO 2 surface play a major role in the regeneration of Pt(O) 4 either directly from Pt(O) 3 or by reaction with the strongly adsorbed CO in Pt(CO)(O) 3 . We show that the formation of an oxygen vacancy and generation of a labile O* are not barrierless, which explains the long lifetime of Pt(CO)(O) 3 and its detectability despite being a minority complex. The results help to develop a comprehensive view of the dynamic evolution of Pt 1 complexes along the reaction cycle and provide mechanistic insights to guide the design of Pt-based single-atom catalysts.

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

Comparison of the canonical transformation and energy functional formalisms for ab initio calculations of self-localized polarons

In materials with strong electron-phonon (e-ph) interactions, charge carriers can distort the surrounding lattice and become trapped, forming self-localized (small) polarons. We recently developed an ab initio approach based on canonical transformations to efficiently compute the formation and energetics of small polarons A different approach based on a Landau-Pekar energy functional has been proposed in the recent literature. Here in this work, we analyze and compare these two methods in detail. We show that the small polaron energy is identical in the two formalisms when using the same polaron wave function. We also show that our canonical transformation formalism can predict polaron band structures and can properly treat zero- and finite-temperature lattice vibration effects, although at present using a fixed polaron wave function. Conversely, the energy functional approach can compute the polaron wave function, but as we show here, it neglects lattice vibrations and cannot address polaron self-localization and thermal band narrowing. Taken together, this work relates two different methods developed recently to study polarons from first-principles, highlighting their merits and shortcomings and discussing them both in a unified formalism.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin dynamics in the itinerant antiferromagnet SrCr 2 ⁢As 2

SrCr 2 ⁢As 2 is an itinerant antiferromagnet in the same structural family as the SrFe 2 ⁢As 2 high-temperature superconductors. Here, we report our calculations of exchange-coupling parameters 𝐽 𝑖⁢𝑗 for SrCr 2 ⁢As 2 using a static linear-response method based on first-principles electronic-structure calculations. We find that the dominant nearest-neighbor exchange coupling 𝐽 1 >0 is antiferromagnetic whereas the next-nearest-neighbor exchange coupling 𝐽 2 <0 is ferromagnetic with 𝐽 2 /𝐽 1 = −0.68, reinforcing the checkerboard in-plane magnetic structure. Thus, unlike other transition-metal arsenides based on Mn, Fe, or Co, we find no competing magnetic interactions in SrCr 2 ⁢As 2 , which aligns with experimental findings. Moreover, the orbital resolution of exchange interactions shows that 𝐽 1 and 𝐽 2 are dominated by direct exchange mediated by the Cr 𝑑 orbitals. To validate the calculations we conduct inelastic neutron-scattering measurements on powder samples that show steeply dispersive magnetic excitations arising from the magnetic Γ points and persisting up to energies of at least 175 meV. The spin-wave spectra are then modeled using the Heisenberg Hamiltonian with the theoretically calculated exchange couplings. In conclusion, the calculated neutron-scattering spectra are in good agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modification of LiMn 2 O 4 surfaces by controlling the Acid–Base surface chemistry of atomic layer deposition

In this combined theoretical and experimental study, we report the effects of Al precursor selection on the growth chemistry, interfacial structure, and electrochemistry of Al 2 O 3 coatings on spinel LiMn 2 O 4 (LMO) surfaces by atomic layer deposition (ALD). Five Al precursors, exhibiting a range of Lewis acid-base properties, were used to establish trends in the ALD Al 2 O 3 growth chemistry on LMO powders in comparison to redox-inactive planar silicon substrates. We show that the Lewis acid-base properties of the Al precursor ligands can be used to tune both the Mn oxidation states and the Al 2 O 3 coverage on the LMO surface. Density functional theory calculations are used to examine the reaction mechanism of each precursor on LMO surfaces, elucidating a correlation between the Lewis acidity of the ligands and the decomposition thermochemistry. In-depth X-ray photoelectron spectroscopy and in situ Fourier transform infrared spectroscopy measurements support these theoretical predictions and further reveal how different Al precursors modify the atomic and electronic structure near the LMO surface during ALD. While the Mn oxidation state is strongly influenced by the Lewis acidity of the precursor ligand, the surface coverage and thickness of the Al 2 O 3 coating are a more representative descriptor of the electrochemical performance measured in coin cell experiments. Finally, we show how the ligand acid-base properties of ALD precursors can be used to rationally tailor atomic layer growth mechanisms, which may enable atomic level control over the structure of functionalized interfaces for various applications related to catalysis, semiconductors, and energy storage.

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

Towards cubic symmetry for Ir 4 + : Structure and magnetism of the antifluorite K 2 IrBr 6

Here, crystal structure, electronic state of Ir 4 + , and magnetic properties of the antifluorite compound K 2 IrBr 6 are studied using high-resolution synchrotron x-ray diffraction, resonant inelastic x-ray scattering (RIXS), thermodynamic and transport measurements, and ab initio calculations. The crystal symmetry is reduced from cubic at room temperature to tetragonal below 170 K and eventually to monoclinic below 122 K. These changes are tracked by the evolution of the noncubic crystal-field splitting Δ measured by RIXS. Nonmonotonic changes in Δ are ascribed to the competing effects of the tilt, rotation, and deformation of the IrBr 6 octahedra as well as tetragonal strain on the electronic levels of Ir 4 + . The Néel temperature of T N = 11.9 K exceeds that of the isostructural K 2 IrCl 6 , and the magnitude of frustration on the fcc spin lattice decreases. We argue that the replacement of Cl by Br weakens electronic correlations and enhances magnetic couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complex structure due to As bonding and interplay with electronic structure in superconducting BaNi 2 As 2

BaNi 2 As 2 is a superconductor chemically related to the Fe-based superconductors, with a complex and poorly understood structural phase transition. We show, based on first-principles calculations, that, in fact, there are two distinct competing structures. These structures are different from the electronic, transport, and bonding points of view but are close in energy. These arise due to complex As bonding patterns and drive distortions of the Ni layers. This is supported by photoemission experiments. This is very distinct from views of the distortion as being primarily driven by electrons at the Fermi surface as in a classical charge density wave, or by correlated electron physics associated with the Ni d electrons and the Ni d orbitals. The structural distortion, although connected with As does lead to an interplay of electronic and structural behavior including induced anisotropic electronic transport. Here, the local bonding nature of the instabilities and the competition between distortions is discussed in the context of the complex behavior observed in BaNi 2 As 2 samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Polarization-resolved core exciton dynamics in LiF using attosecond transient absorption spectroscopy

The ability to control absorption by modifying the polarization of light presents an exciting opportunity to experimentally determine the orbital alignment of absorption features. Here, attosecond extreme ultraviolet (XUV) transient absorption spectroscopy is used to investigate the polarization dependence of core exciton dynamics in LiF thin films at the Li + K edge. XUV pulses excite electrons from the Li 1s core level into the conduction band, allowing for the formation of a p-orbital-like core exciton, aligned along the XUV light polarization axis. A sub-5 fs near-infrared (NIR) probe pulse then arrives at variable time delays, perturbing the XUV-excited states and allowing the coherence decay of the core exciton to be mapped. The coherence lifetimes are found to be ≈ 2.4 ± 0.4 fs, which is attributed to a phonon-mediated dephasing mechanism as in previous core exciton studies. The differential absorption features are also shown to be sensitive to the relative polarization of the XUV and NIR fields. The parallel NIR probe induces couplings between the initial XUV-excited p-like bright exciton and s-like dark excitons. When crossed pump and probe polarizations are used, the coupling between the bright and dark states is no longer dipole-allowed, and the transient absorption signal associated with the coupling is suppressed by approximately 90%. This interpretation is supported by simulations of a few-level model system, as well as analysis of the calculated band structure. Here, the results indicate that laser polarization can serve as a powerful experimental tool for exploring the orbital alignment of core excitonic states in solid-state materials.

Attosecond laser spectroscopy↗
DOE OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗