Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solid solution”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Exploring the deformation behavior of nanotwinned Al–Zr alloy via in situ compression

Nanotwinned metals have demonstrated the capacity for concomitant high strength and ductility. However, metals with high stacking fault energies, such as aluminum (Al), have a low propensity for twin formation. Here, we show the fabrication of supersaturated solid-solution Al–Zr alloys with a high density of growth twins. Incoherent twin boundaries (ITBs) are strong barriers to dislocation motion, while mobile partial dislocations promote plasticity. These deformable nanotwinned Al–Zr alloys reach a flow stress of ∼1 GPa, as demonstrated using in situ micropillar compression tests. Density functional theory calculations uncover the role Zr solute plays in the formation and deformation of the nanotwinned microstructure. This study features a strategy for incorporating ITBs and 9R phase into Al alloys for simultaneous benefits to strength and deformability.

Richter, N. A. (ORCID:0000000332732276)↗
DOE OSTI2022

Optimization of gold germanium (Au 0.17 Ge 0.83 ) thin films for high sensitivity resistance thermometry

Gold–germanium (Au x Ge 1–x ) solid solutions have been demonstrated as highly sensitive thin film thermometers for cryogenic applications. However, little is known regarding the performance of the films for thicknesses less than 100 nm. In response, here we report on the resistivity and temperature coefficient of resistance (TCR) for sputtered films with thicknesses ranging from 10 to 100 nm and annealed at temperatures from 22 to 200°C. The analysis is focused upon composition x=0.17, which demonstrates a strong temperature sensitivity over a broad range. The thinnest films are found to provide an enhancement in TCR, which approaches 20% K –1 at 10 K. Furthermore, reduced anneal temperatures are required to crystallize the Ge matrix and achieve a maximum TCR for films of reduced thickness. These features favor the application of ultra-thin films as high-sensitivity, on-device thermometers in micro- and nanolectromechanical systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Pyrochlore-type lanthanide titanates and zirconates: Synthesis, structural peculiarities, and properties

This contribution provides a thorough examination of the structural characteristics of pyrochlore-type lanthanide titanates and zirconates Ln2Ti2O7 and Ln2Zr2O7, across various length scales. This paper also examines their processing, interesting physical properties (electrical, magnetic, and thermal characteristics), and responses to high pressure and ion irradiation. Brief sections on the elemental oxides' crystal chemistry, pertinent phase diagrams, and energetics of defect formation are also provided. Pyrochlore-type Ln2Ti2O7 and Ln2Zr2O7 stand out as truly multifunctional materials. Moreover, they have emerged as fascinating materials due to magnetic geometrical frustration, arising from the ordering of magnetic Ln3+ and non-magnetic Ti4+ (or Zr4+) cations into separate, interpenetrating lattices of corner-sharing tetrahedra. This results in a diverse array of exotic magnetic ground states, such as spin-ice (e.g., Dy2Ti2O7 or Ho2Ti2O7) or quantum spin ice (e.g., Tb2Ti2O7), observed at both low and room temperatures. They also exhibit varied electrical and electrochemical characteristics. Some members such as Gd2Zr2O7, function as fast ion conductors with a conductivity (σ) of ≈10−2 S·cm−1 at 800 °C and activation energy (Ea) ranging from 0.85 to 1.52 eV, depending on the degree of structural disorder. Others, such as Gd2TiMoO7, are mixed ionic-electronic conductors with σ ≈ 25 S·cm−1 at 1000 °C, making them promising candidate materials for applications in energy conversion and storage devices and oxygen separation membranes. Their exceptionally low thermal conductivity (e.g., κ ∼ 1.1–1.7 W·m−1·K−1 between 700 and 1200 °C for Ln2Zr2O7), close to the glass-like lower limit of highly disordered solids, positions them as valuable materials for thermal barrier coatings. They can also effectively accommodate actinides (e.g., Pu, Np, Cm, Am) in solid solutions and sustain prolonged exposure to radiation due to alpha-decay events, while preserving the integrity of the periodic atomic structure. Proposed as major components in actinide-bearing ceramics, they contribute to the long-term immobilization and disposal of long-lived waste radionuclides from nuclear programs. Some of these properties are displayed simultaneously, opening avenues for new applications. Despite the wealth of data available in the literature, this review highlights the need for a better understanding of order/disorder processes in pyrochlore-type materials and the influence of the structural length scale on their physical and chemical properties. Recent experimental evidence has revealed that pyrochlore short-range structure is far more complex than originally thought. Moreover, pyrochlore local structure is now believed to include short-range, lower symmetry, ordered domains, such as the orthorhombic weberite-type of structure. Notably, short- and long-range structures appear decoupled across different length scales and temperature regimes, and these differences persist even in well-ordered samples. We believe that the pyrochlore structure offers a unique opportunity for examining the interplay between chemical composition, defect chemistry, and properties. In Memoriam: Rodney C. Ewing, Fondly Remembered.

Physics↗
DOE OSTI2024

Effective optimization of atomic decoration in giant and superstructurally ordered crystals with machine learning

Crystals with complicated geometry are often observed with mixed chemical occupancy among Wyckoff sites, presenting a unique challenge for accurate atomic modeling. Similar systems possessing exact occupancy on all the sites can exhibit superstructural ordering, dramatically inflating the unit cell size. In this work, a crystal graph convolutional neural network (CGCNN) is used to predict optimal atomic decorations on fixed crystalline geometries. This is achieved with a site permutation search (SPS) optimization algorithm based on Monte Carlo moves combined with simulated annealing and basin-hopping techniques. Our approach relies on the evidence that, for a given chemical composition, a CGCNN estimates the correct energetic ordering of different atomic decorations, as predicted by electronic structure calculations. This provides a suitable energy landscape that can be optimized according to site occupation, allowing the prediction of chemical decoration in crystals exhibiting mixed or disordered occupancy, or superstructural ordering. Verification of the procedure is carried out on several known compounds, including the superstructurally ordered clathrate compound Rb8Ga27Sb16 and vacancy-ordered perovskite Cs2SnI6, neither of which was previously seen during the neural network training. In addition, the critical temperature of an order–disorder phase transition in solid solution CuZn is probed with our SPS routines by sampling site configuration trajectories in the canonical ensemble. This strategy provides an accurate method for determining favorable decoration in complex crystals and analyzing site occupation at unprecedented speed and scale.

Chemistry↗
DOE OSTI2024

Effects of electron beam irradiation on CrMnV and CrMnTiV high entropy alloys: Nano-mechanical, structural, and thermodynamic perspectives

Beam exit windows are crucial components of any particle accelerator as they provide an interface between the beamline vacuum and target material at atmospheric media. For high beam power machines, special materials and designs are required to withstand high radiation and mechanical loads, while minimizing energy loss during transition and maximizing window lifetime. This research investigates the impact of electron beam exposure to bulk CrMnV and CrMnTiV high entropy alloys (HEAs) with the primary goal of identifying suitable candidate materials for the design of robust and durable exit window settings. The selection criteria include intrinsic characteristics, power dissipation, and mechanical responses. According to the thermodynamic calculations, both equiatomic CrMnV and the addition of 7% of Ti with equiatomic CrMnV yield solid-solutions phases. The structural and mechanical properties of CrMnV and CrMnTiV samples were tested using field emission scanning electron microscopy, atomic force microscopy, scanning electron microcopy with energy dispersive x-ray spectroscopy, x-ray diffraction, and nanoindentation before and after exposure to a dose of ~66 kGy from a 10 MeV e-beam accelerator. Despite exhibiting beam transmission characteristics comparable to Cr and V, the indentation hardness of HEAs exceeded that of the Cr and V samples by five to six times. The examination of the CrMnTiV irradiated samples revealed organized deformation patterns depicting new features, which we suspect twinning and twin boundaries due to the addition of Ti to CrMnV. Ti, a hexagonal-close-packed crystal structure, is commonly known for deformation twinning behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High pressure structural and lattice dynamics study of α-In2Se3

Layered α-In2Se3 has been studied using a combined in situ synchrotron angle-dispersive powder x-ray diffraction and Raman spectroscopy study in a diamond anvil cell up to 60+ GPa, at room temperature. Helium, which remains fairly hydrostatic up to the highest pressure in this study, was used as the pressure-transmitting medium. The results from both experimental methods reveal a pressure-induced structural phase transition from α-In2Se3 to a monoclinic β'-In2Se3 structure at ≈1 GPa, in agreement with previous studies. Based on our detailed measurements using both experimental techniques and the F-f formalism, the β'-In2Se3 structure remains stable up to 45 GPa, without a clear indication of a phase transition toward the previously reported β-In2Se3 phase. Above this pressure, In2Se3 adopts a disordered solid-solution-like orthorhombic structure, phase IV. The results are discussed in comparison with the relevant previous studies of α-In2Se3 under pressure.

Feng, Shiyu↗
DOE OSTI2023

Effect of local chemical order on the irradiation-induced defect evolution in CrCoNi medium-entropy alloy

High- (and medium-) entropy alloys have emerged as potentially suitable structural materials for nuclear applications, particularly as they appear to show promising irradiation resistance. Recent studies have provided evidence of the presence of local chemical order (LCO) as a salient feature of these complex concentrated solid-solution alloys. However, the influence of such LCO on their irradiation response has remained uncertain thus far. Here, in this work, we combine ion irradiation experiments with large-scale atomistic simulations to reveal that the presence of chemical short-range order, developed as an early stage of LCO, slows down the formation and evolution of point defects in the equiatomic medium-entropy alloy CrCoNi during irradiation. In particular, the irradiation-induced vacancies and interstitials exhibit a smaller difference in their mobility, arising from a stronger effect of LCO in localizing interstitial diffusion. This effect promotes their recombination as the LCO serves to tune the migration energy barriers of these point defects, thereby delaying the initiation of damage. These findings imply that local chemical ordering may provide a variable in the design space to enhance the resistance of multi-principal element alloys to irradiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Designing Nuclear Fuels with a Multi-Principal Element Alloying Approach

Previous research has shown that multi-principal element alloys (MPEAs) using chromium, molybdenum, niobium, tantalum, titanium, vanadium, and zirconium can form stable body-centered-cubic (BCC) structures across a large temperature region (25°C to 1000°C). This is the same crystal structure as γ-uranium (U), which has shown desirable thermal and irradiation behavior in previous alloy fuel research. It is hypothesized then that the MPEA alloying approach can be used to produce a stable BCC uranium-bearing alloy and to retain its stability throughout anticipated operating regimes of power-producing reactors. Candidate elements were assessed using Monte Carlo N-Particle (MCNP) analysis to determine uranium densities necessary to make the alloy an economically viable fuel compared to conventional fuel forms. Following neutronic considerations, materials property databases and empirical predictors were used to determine the compositions with a high potential to form a BCC solid solution alloy. The final four alloys were MoNbTaU 2 , MoNbTiU 2 , NbTaTiU 2 , and NbTaVU 2 , which were cast using arc melting of raw elemental foils and chunks. Characterization of the fabricated alloys included scanning electron microscopy, X-ray diffraction, energy dispersive X-ray spectroscopy, and transmission electron microscopy. The results showed a two-phase system with a U-rich matrix phase surrounding the refractory precipitates. The U phase was found to contain varying concentrations of the alloying elements and was a BCC γ-U phase. These results warrant further research to identify ideal compositions for use as an advanced alloy fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Metal–insulator transition in composition-tuned nickel oxide films

Abstract Thin films of the solid solution Nd 1 − x La x NiO 3 are grown in order to study the expected 0 K phase transitions at a specific composition. We experimentally map out the structural, electronic and magnetic properties as a function of x and a discontinuous, possibly first order, insulator–metal transition is observed at low temperature when x = 0.2. Raman spectroscopy and scanning transmission electron microscopy show that this is not associated with a correspondingly discontinuous global structural change. On the other hand, results from density functional theory (DFT) and combined DFT and dynamical mean field theory calculations produce a 0 K first order transition at around this composition. We further estimate the temperature-dependence of the transition from thermodynamic considerations and find that a discontinuous insulator–metal transition can be reproduced theoretically and implies a narrow insulator–metal phase coexistence with x . Finally, muon spin rotation ( µ SR) measurements suggest that there are non-static magnetic moments in the system that may be understood in the context of the first order nature of the 0 K transition and its associated phase coexistence regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Effective pair potential for random fcc CoCrFeMnNi alloys

The single-phase equiatomic CoCrFeMnNi alloy is a random solid solution of five elements on the face-centered cubic lattice, whose pure constituents crystallize in very different structures and exhibit diverse magnetic properties. Due to the randomness of the alloy, 80% of nearest neighbor bonds are between unlike elements and thus the details of bonding in pure structures are less important. The elastic moduli of this alloy give rise to small Cauchy pressure C 12 - C 44 , which suggests that the dominant part of bonding may be described by a simple pair potential. We test this hypothesis by developing a long-range Lennard-Jones potential in which the equilibrium crystal structures of pure constituents are taken as reference. The standard mixing rules for regular solutions are then adopted to obtain parameters for bonds between unlike elements in the quinary system. The transferability of this potential to quaternary CoCrFeNi, ternary CoCrNi, and binary FeNi alloys is investigated and the predictions compared with experiments and density functional theory calculations. By sampling over a large number of random configurations, we investigate the effect of compositional randomness on misfit volumes, energies of point defects and stacking faults, and the dislocation friction stresses experienced by moving edge and screw dislocations.

Materials Science↗
DOE OSTI2021

New insight into strain and composition of BaZrO 3 nanorods in REBCO superconductor

We report on strain and composition effects associated with growth of self-assembled BaZrO 3 (BZO) nanorods in REBa 2 Cu 3 O 7–δ (REBCO) superconductors (RE = rare earth = Y and Gd), which have a profound effect on flux pinning and in-field critical current performance. The a-b plane mismatch between BZO and REBCO is never fully coherently accommodated. Instead, the nanorods always assume a size at least one unit cell smaller than the corresponding 'hole' in the REBCO matrix, thus providing deep minima of in-plane mismatch strain. Next, we show that the nominal BZO nanorods are in fact solid solution Ba 2+ (Zr 4+ 1–z RE 3+ z )O 3–δ perovskite, thus strongly affecting the stoichiometry and relative amounts of REBCO and BZO. We demonstrate that by varying only the Ba content in the nominal composition of 15 mol.% BZO + REBCO, the unit cell density of BZO can be tuned from 5% to 23%, and the linear density of RE 2 O 3 (REO) precipitates from 18 to 1 μm –1 . Furthermore, the results explain the wide range of pinning performances observed for the same nominal amount of Zr addition and provide insight into the mechanisms behind the complex phenomenon of growth of nanorods by self-assembly in REBCO superconductors.

2G-HTS↗
DOE OSTI2020

Laser metal deposition of refractory high-entropy alloys for high-throughput synthesis and structure-property characterization

Progress in materials development is often paced by the time required to produce and evaluate a large number of alloys with different chemical compositions. This applies especially to refractory high-entropy alloys (RHEAs), which are difficult to synthesize and process by conventional methods. To evaluate a possible way to accelerate the process, high-throughput laser metal deposition was used in this work to prepare a quinary RHEA, TiZrNbHfTa, as well as its quaternary and ternary subsystems by in-situ alloying of elemental powders. Compositionally graded variants of the quinary RHEA were also analyzed. Our results show that the influence of various parameters such as powder shape and purity, alloy composition, and especially the solidification range, on the processability, microstructure, porosity, and mechanical properties can be investigated rapidly. The strength of these alloys was mainly affected by the oxygen and nitrogen contents of the starting powders, while substitutional solid solution strengthening played a minor role.

additive manufacturing↗
DOE OSTI2023

Will high-entropy carbides and borides be enabling materials for extreme environments?

Abstract The concept of multi-principal component has created promising opportunities for the development of novel high-entropy ceramics for extreme environments encountered in advanced turbine engines, nuclear reactors, and hypersonic vehicles, as it expands the compositional space of ceramic materials with tailored properties within a single-phase solid solution. The unique physical properties of some high-entropy carbides and borides, such as higher hardness, high-temperature strength, lower thermal conductivity, and improved irradiation resistance than the constitute ceramics, have been observed. These promising properties may be attributed to the compositional complexity, atomic-level disorder, lattice distortion, and other fundamental processes related to defect formation and phonon scattering. This manuscript serves as a critical review of the recent progress in high-entropy carbides and borides, focusing on synthesis and evaluations of their performance in extreme high-temperature, irradiation, and gaseous environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TEM Approaches for Microstructure-Informed Prediction of Mechanical Properties in Structural Alloys

Predicting the mechanical performance of structural alloys from their evolving microstructure remains a major challenge in materials science, particularly for nuclear structural materials, where irradiation-induced defects span multiple types and length scales and interact through complex mechanisms. The dispersed barrier hardening (DBH) [1] and Friedel–Kroupa–Hirsch (FKH) [2,3] models have been widely used to evaluate the hardening contributions of individual obstacles and to estimate tensile strength from quantified microstructures; however, when multiple size-dependent obstacles coexist and evolve, predicting temperature-dependent tensile strength becomes significantly more complex, and a fully consistent hardening model is still lacking. Transmission electron microscopy (TEM) plays a central role in refining hardening models and enabling predictive assessments of tensile strength evolution by providing quantitative characterization of dislocations, irradiation-induced defects (e.g., dislocation loops and cavities), precipitates, and grain structure (Fig. 1.). These experimentally measured defect densities are incorporated into physically based hardening models with size- and shape- dependent obstacle strengths [4], using root-sum-square superposition for obstacles of comparable strength and linear superposition for dissimilar ones [5]. In addition, recent advances in TEM [6-8], including high-resolution imaging, 4D-STEM strain mapping, EDS/EELS elemental analysis, and flash-polishing-based TEM specimen preparation and extraction-replica methods (Fig. 2), further improve the accuracy of microstructural quantification. By comparison with prior studies as well as our own results, we show that when TEM-derived microstructural information is carefully integrated with physically grounded hardening models, yield strength (or irradiation-induced hardening) measured at room temperature can be predicted with good quantitative agreement across multiple alloy classes. In-situ TEM combined with high-temperature mechanical testing represents an important next step for refining hardening models by directly probing dislocation–obstacle interactions across varying irradiation doses and temperatures [9]. Because the barrier strength factor (α) depends on both temperature and obstacle size, it should not be treated as a constant fitting parameter; rather, it must be explicitly evaluated to achieve physically meaningful predictions of mechanical behaviour at operating temperatures. This presentation therefore discusses why all strengthening contributions (e.g., Peierls stress, solid-solution strengthening, voids, bubbles, dislocation loops, dislocation lines, and grain boundaries) must be considered collectively, why appropriate superposition methods are essential when obstacles possess different barrier strength factors, how hardness measurements can be meaningfully related to tensile properties, and how TEM-derived microstructural information can be systematically incorporated into hardening models. More broadly, it outlines a pathway toward microstructure-informed prediction of mechanical properties and supports the goal of establishing science-based tools for evaluating structural materials in extreme environments [10].

Lin, Yan-Ru [ORNL] (ORCID:0000000339991473)↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
DOE OSTI2020

Averaged cluster approach to including chemical short-range order in KKR-CPA

The single-site Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA) ignores short-range ordering present in disordered metallic systems. Here in this paper, we establish a technique to fix this shortcoming by embedding an averaged cluster that displays chemical short-range order (SRO). The degree of SRO can be tuned by externally defined order parameters. This averaged cluster can be embedded in the single-site CPA medium or a self-consistently obtained effective medium that contains SRO information. The validity of this method is demonstrated by applying it to two alloy systems—the CuZn body-centered cubic (BCC) solid solution, and AlCrTiV, a four-element BCC high entropy alloy. A comparison between the non-self-consistent and self-consistent modes is also provided for the two above-mentioned systems. We make the code available on the internet. Planned extensions to this paper are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phonon dispersion throughout the iron spin crossover in ferropericlase

Ferropericlase (Fp), (Mg 1– x Fex)O, is the second most abundant phase in the Earth's lower mantle. At relevant pressure-temperature conditions, iron in Fp undergoes a high spin (HS), $\textit{S}$ = 2, to low spin (LS), $\textit{S}$ = 0, state change. The nature of this phenomenon is quite well understood now, but there are still basic questions regarding the structural stability and the existence of soft phonon modes during this iron state change. General theories exist to explain the volume reduction, the significant thermoelastic anomalies, and the broad nature of this HS-LS crossover. These theories make extensive use of the quasiharmonic approximation. Therefore, dynamical and structural stability is essential to their validity. Here, we investigate the vibrational spectrum of Fp throughout this spin crossover using ab initio density-functional theory $+U_{sc}$ calculations. We address vibrational modes associated with isolated and (second-)nearest-neighbor iron ions undergoing the HS-LS state change. As expected, acoustic modes of this solid solution are resilient, while optical modes are the most affected. In this work, we show that there are no soft phonon modes across this HS-LS crossover, and Fp is dynamically stable at all relevant pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗