Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “low frequency”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Reverse phonon thermal flux from an applied electric field

In the phenomenon known as electron drag, a phonon thermal flux is established in a conducting crystal by an electric field applied under isothermal conditions through directed transfer of quasimomentum from the electronic charge current to the phonon subsystem. Prior understanding of this phenomenon involves a thermal current composed of low-frequency acoustic phonons that is in the same direction as the charge current. This results in an increase of the Peltier thermopower. Here, we show that it is also possible to establish a phonon thermal current that is in the opposite direction to the charge current, which we refer to as a reverse phonon thermal flux. We demonstrate this behavior both through a simple qualitative model and through first-principles calculations performed for three materials: p-type 𝜃-TaN, n-type BAs, and n-type Si. The reverse phonon flux is shown to arise through intervalley electron-phonon scattering processes involving high-frequency acoustic phonons. Unlike n-BAs and n-Si, the combined features in the band structure and phonon dispersions of p-type 𝜃-TaN promote a particularly large reverse phonon thermal flux comparable in magnitude to the oppositely directed thermal flux from low-frequency phonons. Finally, this work highlights a previously unrecognized behavior in the coupled electron-phonon system and advances our understanding of the rich physics of transport in solids.

lattice dynamics↗
DOE OSTI2025

Breakdown of the Static Dielectric Screening Approximation of Coulomb Interactions in Atomically Thin Semiconductors

Coulomb interactions in atomically thin materials are remarkably sensitive to variations in the dielectric screening of the environment, which can be used to control exotic quantum many-body phases and engineer exciton potential landscapes. For decades, static or frequency-independent approximations of the dielectric response, where increased dielectric screening is predicted to cause an energy redshift of the exciton resonance, have been sufficient. These approximations were first applied to quantum wells and were more recently extended with initial success to layered transition metal dichalcogenides (TMDs). Here, we use charge-tunable exciton resonances to investigate screening effects in TMD monolayers embedded in materials with low-frequency dielectric constants ranging from 4 to more than 1000, a range of 2 orders of magnitude larger than in previous studies. In contrast to the redshift predicted by static models, we observe a blueshift of the exciton resonance exceeding 30 meV in higher dielectric constant environments. We explain our observations by introducing a dynamical screening model based on a solution to the Bethe-Salpeter equation (BSE). When dynamical effects are strong, we find that the exciton binding energy remains mostly controlled by the low-frequency dielectric response, while the exciton self-energy is dominated by the high-frequency one. Our results supplant the understanding of screening in layered materials and their heterostructures, introduce a knob to tune selected many-body effects, and reshape the framework for detecting and controlling correlated quantum many-body states and designing optoelectronic and quantum devices.

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

Model validation and error attribution for a drifting qubit

Qubit performance is often reported in terms of a variety of single-value metrics, each providing a facet of the underlying noise mechanism limiting performance. However, the value of these metrics may drift over long timescales, and reporting a single number for qubit performance fails to account for the low-frequency noise processes that give rise to this drift. Here, in this work, we demonstrate how we can use the distribution of these values to validate or invalidate candidate noise models. We focus on the case of randomized benchmarking (RB), where typically a single error rate is reported but this error rate can drift over time when multiple passes of RB are performed. We show that using a statistical test as simple as the Kolmogorov-Smirnov statistic on the distribution of RB error rates can be used to rule out noise models, assuming the experiment is performed over a long enough time interval to capture relevant low frequency noise. With confidence in a noise model, we show how care must be exercised when performing error attribution using the distribution of drifting RB error rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First constraints on causal sources of primordial gravitational waves from BICEP/Keck, SPTpol, SPT-3G, Planck and WMAP $B$-mode data

Non-inflationary sources of gravitational waves in the early Universe generically predict causality-limited tensor power spectra at low frequencies. We report the first-ever constraints on such sources based on cosmic microwave background (CMB) $B$-mode polarization measurements. Using data from BICEP/Keck, SPTpol, SPT-3G, Planck, and WMAP, we constrain the amplitude of an early causal tensor (ECT) power spectrum parameterized by $r_{ect}$, the ratio of causal tensor power to total scalar power at $k~=~0.01$ Mpc$^{-1}$, and obtain a 95% CL upper limit of $r_{ect}<$ 0.0077. Since $r_{ect}$ can easily be related to the parameters of a given theory, our bound robustly constrains a broad class of well-motivated gravitational wave sources in the early universe, including first-order cosmological phase transitions, enhanced small-scale density perturbations, and various topological defects. Finally, we translate our limit into a bound on the present-day energy density in gravitational waves at ultra-low frequencies otherwise inaccessible to traditional gravitational wave detection strategies, including pulsar timing arrays, interferometers, and resonant cavities.

Zebrowski, Jessica A. [Fermilab; Chicago U., KICP]↗
DOE OSTI2025

Fermilab PIP II Machine Protection System Digitized Data Noise Elimination Scheme and Its FPGA Implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can con-taminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is re-quired. Several real-time functional blocks were simu-lated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid base-line by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.

Wu, Jinyuan [Fermilab] (ORCID:0000000344329521)↗
DOE OSTI2026

Fermilab PIP II machine protection system digitized data noise elimination scheme and its FPGA implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can contaminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is required. Several real-time functional blocks were simulated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid baseline by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.

Wu, J. [Fermilab] (ORCID:0000000344329521)↗
DOE OSTI2025

Identifying Challenges in Safeguards for Metallic Fuel Fabrication Facilities

As new advanced reactors gain popularity, there is an increasing interest in metallic fuel fabrication for fast reactors. While metallic fuels themselves are not a new idea, as many of the first reactors employed metallic fuels, new designs, compositions, and fabrication methods are appearing throughout the nuclear community. As the interest grows and facilities are constructed, both domestic and international safeguards will need to be heavily involved to support safeguards-by-design (SBD) measures from the start. This work compiles a review of historical and modern fuel types and fabrication methods, fabrication processes, safeguards gaps, and potential safeguards solutions. Metallic nuclear fuel types have been around for many decades and were included in some of the first reactors including the Experimental Breeder Reactor (EBR)-I and -II, the Fermi 1 reactor, the Integral Fast Reactor (IFR), and the Dounreay Fast Reactor (DFR). These reactors used various compositions including pure uranium (U) metal, U-zirconium (Zr) alloys, plutonium (Pu)-aluminum (Al) alloys, U-fissium (Fs) alloys, U-Pu-Zr alloys, and U-molybdenum (Mo) alloys [1, 2, 3, 4, 5]. These small alloying additions are included to improve the material properties of the pure U metal. The alpha-phase U (stable below 661C) suffers elongation in one direction causing grain boundary cracking and increasing creep rate due to irradiation growth, thermal cycling, and preferential crystal orientation. It is ideal to utilize the gamma-phase U (typically stable above 769C) by adding small amounts of alloying elements such as Zr or Mo to stabilize this phase down to room temperature [3]. Additionally, some research has been focused on U with transuranic (TRU) elements present, typically coming from the used fuel recycling process. Including these elements in fast reactor fuel can aid in the reduction of nuclear waste by burning minor long-lived actinides. However, the additions of TRU elements can cause concerns to arise when trying to fabrication or safeguard metallic fuels. A typical metallic fuel element is shown in Figure 1. Sodium is added into the cladding to create a thermal bond between the fuel slug and cladding wall. The fuel slug is then inserted and the end plug is welded on to the top of the fuel element. A gas plenum is left to create a headspace for gaseous fission products to escape rather than continue to build in the fuel itself [1, 5]. Other fuel element geometries exist as well, such as the Lightbridge twisted cruciform geometry shown in Figure 2 [6]. This design allows for better cooling performance and provides room for fuel rod swelling without impacting the fuel rod diameter. There are many different fabrication methods for metallic fuels, which is one of the many benefits of these fuel types. Many of these fabrication methods are relatively easy and cost-efficient. The most popular fabrication method is injection casting, sometimes called vacuum induction melting (VIM), shown in Figure 3 [4, 8, 9, 7, 10]. This method was largely used for EBR-II fuel fabrication. The injection casting system is contained inside of a vessel consisting of a Y2O3-coated graphite crucible surrounded by an induction coil with ZrO2-coated quartz molds suspended above the crucible. The fuel feedstock is placed inside of the graphite crucible and melted using the induction furnace. The induction furnace utilizes a dual frequency with the high frequency melting the feedstock and the low frequency causing stirring of the melted feedstock to form a homogeneous mixture. The mixture is heated to approximately 1600C in an argon environment. The vessel is evacuated and then the quartz molds are lowered into the graphite crucible containing the molten metal and the vessel is repressurized to inject the metal fuel upwards into the molds. The molds are removed and then shattered to release the fuel slugs. This fabrication method was used to fabricate 39,000 metallic fuel pins for EBR-II. While injection casting has been the most common metallic fuel fabrication method throughout the decades, many other methods have been explored including low-pressure gravity casting, microwave casting, continuous casting, centrifugal casting, coextrusion, and many others [11, 12, 8, 13, 14, 15]. Some of these methods aim to mitigate challenges that arise with americium (Am) volatilization during the casting process for TRU-containing fuel feedstocks, an issue with injection casting. Coextrusion is one of the methods explored at the Idaho National Laboratory (INL) and has been utilized for the initial fabrication tests of Lightbridge's unique fuels, as well as other metallic fuels with cladding coextruded. In this process, large billets are formed and machined and then inserted into a molten salt bath for approximately 30 minutes. The billets are then loaded into the extrusion press and extruded. This process can be seen in Figure 4 [15].

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

A generalized effective potential for differentially rotating plasmas

Global stability of differentially rotating plasma is investigated using a generalized effective potential. We first, for a current-free system, obtain a general form of an effective potential in terms of the free energies of global curvature and gradients of rotation for non-axisymmetric disturbances. We then examine the stability of differentially rotating disks for several rotation profiles and present the associated effective potential for the onset of these instabilities in the MHD regime. In particular, results for global axisymmetric magnetorotational instability as well as local and global non-axisymmetric modes are presented. The latter constitute two distinct non-axisymmetric modes, a high frequency local MRI and a global low-frequency non-axisymmetric mode (the magneto-curvature mode, introduced in Ebrahimi and Pharr [Astrophys. J. 936, 145 (2022)]), confined either between two Alfvénic resonances or an Alfvénic resonance and a boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Variable Effects of Dispersed Nanoparticles on Triboelectric Nanogenerators

Technology has recently seen a drastic physical downsizing. Wearable and small devices with lower power demands have become the norm and continue to be more prominent in daily life. With modern devices growing smaller and requiring less electricity, a power source will always be needed. Contemporary batteries are the most common means to power small electronics. However, reliance on conventional batteries may prove insufficient due to the non renewable resources (Li, Ni, Co) required to power the growing number of individual devices each person may own. Additionally, the infrastructure required to harvest and recycle the sheer number of batteries produced presents a further logistic issue to be addressed. A promising alternative to batteries is the usage of triboelectric nanogenerators (TENGs). TENGs are a class of energy harvesting devices that utilize triboelectric generation to convert mechanical/kinetic energy into electrical energy and have exhibited efficiencies up to 85 % at low frequencies. TENGs exhibit a high voltage but low current. Even with the high voltage, the low current output proves to be a significant factor preventing undoped TENGs from being commercially viable. This review will investigate factors that increase the total current produced by TENGs when nanoparticles are utilized in TENGs. Factors such as increasing porosity, surface area, surface charge density, charge storage, deep trap formation, and dielectric constant can be altered to affect the total current by impregnating nanoparticles into the polymer material will be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Variable Effects of Dispersed Nanoparticles on Triboelectric Nanogenerators

Technology has recently seen a drastic physical downsizing. Wearable and small devices with lower power demands have become the norm and continue to be more prominent in daily life. With modern devices growing smaller and requiring less electricity, a power source will always be needed. Contemporary batteries are the most common means to power small electronics. However, reliance on conventional batteries may prove insufficient due to the non renewable resources (Li, Ni, Co) required to power the growing number of individual devices each person may own. Additionally, the infrastructure required to harvest and recycle the sheer number of batteries produced presents a further logistic issue to be addressed. A promising alternative to batteries is the usage of triboelectric nanogenerators (TENGs). TENGs are a class of energy harvesting devices that utilize triboelectric generation to convert mechanical/kinetic energy into electrical energy and have exhibited efficiencies up to 85 % at low frequencies. TENGs exhibit a high voltage but low current. Even with the high voltage, the low current output proves to be a significant factor preventing undoped TENGs from being commercially viable. This review will investigate factors that increase the total current produced by TENGs when nanoparticles are utilized in TENGs. Factors such as increasing porosity, surface area, surface charge density, charge storage, deep trap formation, and dielectric constant can be altered to affect the total current by impregnating nanoparticles into the polymer material will be explored.

42 ENGINEERING↗
DOE OSTI2026

Coming Closer to High Frequency Gravitational Wave Detection With MAGO

In the last years, low frequency gravitational waves (GWs) have been consistently measured by the LIGO&#x2;Virgo collaboration, but little to no attention has been paid to higher frequencies GWs in the range of 10 kHz to 100 MHz, at which confirmation for current theories or even new physics could be hidden. The MAGO 2.0 project aims at filling this gap in the parameter space using superconducting radio-frequency (SRF) cavities. Exploiting the excellent Q-factors of these resonators, we plan to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. We present the results of the first cold tests ran at DESY and FNAL using the cavity prototype built 20 years ago at the end of the MAGO collaboration, characterizing the RF spectrum, Q-factor and surface resistance. Additionally we introduce the mechanical vibration spectrum characterization and the RF response of the cavity with the injection of a “fake GW” signal using piezoelectric actuators.

Marconato, Giovanni [Hamburg U.]↗
DOE OSTI2026

Decadal Spiciness Variability in the Subtropical-Tropical Pacific in the CESM2 Large Ensemble

Tropical Pacific decadal variations impact weather and climate around the world and are also connected to variations in the global warming trend. The mechanisms driving these long-term modulations, particularly the role of subsurface ocean dynamics, are still debated. Here, we investigate the dynamics of spiciness (density-compensated temperature and salinity) anomalies in the tropical and subtropical Pacific, which are hypothesized as a possible driving mechanism of decadal climate variability. Based on the analysis of 100 realizations from the Community Earth System Model Version 2 Large Ensemble (CESM2-LE), we demonstrate a coupling between the subtropics and the equatorial Pacific by propagating spiciness anomalies at decadal time scales. The CESM2-LE simulates spiciness variability along a subduction path from the subtropics to the equator with frequency spectra that show the highest power at low frequencies and a power decay proportional to a −4 slope for frequencies greater than 0.01 cycles per months, corresponding to periods smaller than ∼8.5 years. Signals that originate in the Southern Hemisphere (SH) dominate and arrive with a larger magnitude at the equator compared to spiciness anomalies from the Northern Hemisphere (NH). Spiciness anomalies from the SH have shorter propagation times and are strengthened along their pathway as stochastic wind stress curl forcing generates anomalous baroclinic ocean pressure gradients. These pressure gradients generate spiciness anomalies via anomalous advection across climatological spiciness gradients in the SH. We conclude that the observed spiciness variance at decadal time scales is consistent with a forcing by stochastic wind variations that are low-pass filtered by ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

OC7 phase I: Toward practical sea-state-dependent modeling of hydrodynamic viscous drag and damping

Here, this article presents a collaborative research campaign under the OC7 project on refining the engineering modeling approach for hydrodynamic viscous drag and damping on floating wind platforms, focusing on the adjustment of hydrodynamic drag and damping coefficients for different sea states. The participant simulation results show significant improvements over the previous OC6 project in predicting the low-frequency resonance motion under nonoperational conditions. The improvements are mainly due to enhanced modeling, including the adoption of wave stretching, and directly tuning the coefficients to measured platform motion in waves instead of free decay. For accurate predictions of mean- and slow-drift motion, the better performing models use a decreasing column splash zone drag coefficient and increasing surge damping/drag with increasing wave height. The model tuning for heave and pitch resonance shows less consistency. Generally, both quadratic drag and additional heave or pitch damping are needed for accurate predictions. Alternatively, a quadratic drag formulation with velocity filtering for the rectangular pontoons leads to improved predictions without additional damping. This model is also potentially more predictive, requiring minimal adjustment to its parameters for different conditions.

17 WIND ENERGY↗
DOE OSTI2026

Predicting critical heat flux using localized sensing at invisible vapor-liquid interfaces

Predicting critical heat flux (CHF) in two-phase electronics cooling systems remains a significant challenge due to the sudden onset of boiling crisis and the difficulty in directly visualizing vapor-liquid interfaces. Existing sensing methods rely on lagging temperature measurements, optically accessible systems, or spatially averaged signals that cannot pinpoint CHF initiation at localized high-heat-flux regions. Here, we report a planar capacitive sensing approach that enables real-time, localized detection of vapor-liquid interface dynamics for CHF prediction in boiling heat transfer. The capacitive sensor exploits the dielectric constant difference between liquid and vapor phases to capture bubble nucleation, growth, and departure dynamics with a temporal resolution down to 2 ms. The capacitive sensing reveals distinct signals across boiling regimes: from high-frequency fluctuations during strong nucleate boiling to low-frequency fluctuations with increased amplitudes when approaching CHF. The multi-sensor array experiments demonstrate real-time localized sensing, where each sensor responds exclusively to boiling in its immediate vicinity without crosstalk from neighboring regions. This non-intrusive sensing approach provides predictive rather than lagging sensing signals of CHF occurrence, offering predictive diagnosis of two-phase liquid cooling for the thermal management of high-power-density electronics.

CHF↗
DOE OSTI2025

Theory of tunneling between two-dimensional electron layers driven by spin pumping: Adiabatic regime and beyond

Tunneling spectroscopy between parallel two-dimensional (2D) electronic systems provides a powerful method to probe the underlying electronic properties by measuring tunneling conductance. Here, in this work, we present a theoretical framework for spin transport in 2D-to-2D tunneling systems, driven by spin pumping. This theory applies to a vertical heterostructure where two layers of metallic 2D electron systems are separated by an insulating barrier, with one layer exchange coupled to a magnetic layer driven at resonance. Utilizing a nonperturbative Floquet-Keldysh formalism, we derive general expressions for the tunneling spin and charge currents across a broad range of driving frequencies, extending beyond the traditional adiabatic pumping regime. At low frequencies, we obtain analytical results that recover the known behaviors in the adiabatic regime. However, at higher frequencies, our numerical findings reveal significant deviations in the dependence of spin and charge currents on both frequency and precession angle. This work offers fresh insights into the role of magnetization dynamics in tunneling transport, opening up new avenues for exploring nonadiabatic spin pumping phenomena.

Green's function methods↗
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↗