SEARCH · Engineering Papers

Results for “bulk modulus”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Enhancing the accuracy and generality of the Debye–Grüneisen Model: Optimizing the volume dependence for accurate predictions across varied compositions

In this work, we have introduced an optimized Debye-Grüneisen model that revolutionizes the determination of the Debye temperature and Grüneisen parameters. Unlike conventional methods, our model requires only the 0 K energy volume data for a material as input, eliminating the need to determine the bulk modulus and its pressure derivative, which often pose challenges due to numerical uncertainties. This unique feature sets our model apart from existing approaches and streamlines the process, enabling accurate predictions of thermal expansion behavior across various materials. To demonstrate its effectiveness, we showcase its excellent agreement with measured coefficients of thermal expansion (CTE) for the nickel-cobalt-chromium-aluminum-yttrium (Ni-Co-Cr-Al-Y) bond-coating system. Additionally, we apply our approach by conducting a high-throughput search for potential bond-coating materials among 90,000 compositions within the aluminum-cobalt-chromium-iron-nickel (Al-Co-Cr-Fe-Ni) system. From this extensive search, four compositions are synthesized, and the measured CTE values agree very well with theoretical predictions, hence validating our approach. In conclusion, the current optimized Debye-Grüneisen model combined with Density Functional Theory (DFT)-based thermodynamic database enables reliable and efficient high-throughput calculations of CTE of of a material without expensive phonon calculations.

Bond coating materials↗

DOE OSTI2025

Intrinsic mechanical properties and seeding effect of tobermorite synthesized in supercritical water

This publication reports for the first time the physiochemical, the intrinsic mechanical properties and the seed effect of anomalous Al-substituted 11 Å tobermorite synthesized via the innovative supercritical hydrothermal flow process at 400 °C and 25 MPa. This approach allows synthesizing highly crystalline tobermorite fibers in only 8 s, with characteristics very close to the natural tobermorite. The anomalous 11 Å tobermorite exhibits aluminosilicate chains with a high polymerisation degree and a less defective structure compared to materials produced via the conventional hydrothermal method. Furthermore, the intrinsic mechanical properties of Al-tobermorite synthesized in supercritical water are investigated for the first time by High-Pressure XRD. This Al-substituted 11 Å tobermorite is characterised by higher incompressibility along the b-axis and bulk modulus K0 in comparison with what is commonly observed for other synthetic tobermorite. The tobermorite acts as nucleation points to trigger the quick formation of the hydration product in Portland cement paste.

Calcium silicate hydrate↗

DOE OSTI2026

Pressure-Driven Helium Insertion for Structural Stability of CH 3 NH 3 PbBr 3 Hybrid Perovskites

Organic−inorganic metal halide perovskites (MHPs) have garnered significant attention due to their outstanding performance in optoelectronic devices, but they are prone to degradation through a variety of pathways. High-pressure studies are being used to understand the mechanical and structural stability of MHPs and investigate new methods for improving these. In this study, we map the high-pressure, low-temperature structural phase diagram of CH 3 NH 3 PbBr 3 (MAPbBr 3 ) between 15−300 K and up to 1.5 GPa. We first compare the temperature and pressure effects on the global and local structures of MAPbBr 3 using synchrotron X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS), respectively, and find evidence of PbBr 6 octahedral distortion. Our results also suggest that He inserts into the MAPbBr 3 structure. The thermodynamics of He insertion into MAPBHe x are calculated to be plausible at room temperature for x ≤ 1 and P = 1−3 GPa and shown to stabilize higher-symmetry phases. This phenomenon increases the fitted bulk modulus of MAPbBr 3 . Structural stabilization by inert atom insertion is proposed as a method to improve the mechanical robustness of MHPs and other soft hybrid materials.

Diffraction↗

DOE OSTI2023

High-Pressure Structural and Thermodynamic Properties of Cerium Orthosilicates (CeSiO 4 )

Pressure-induced phase transitions from the zircon structure-type (I4 1 /amd) to the scheelite structure type (I4 1 /a) are known for many ternary oxides systems (ABO 4 ). In this work, we present the first high-pressure study on synthetic stetindite (CeSiO 4 ) by a combination of in situ high-pressure synchrotron powder X-ray diffraction up to 36 GPa, implemented with and without dual sided laser heating, and in situ high-pressure Raman spectroscopy up to 43 GPa. Two phase transitions were identified: zircon to a high-pressure low-symmetry (HPLS) phase at 15 GPa and then to a scheelite at 18 GPa. The latter from HPLS scheelite phase was found irreversible; i.e., scheelite is fully quenchable at ambient conditions, as in other zircon-type phases. The bulk moduli (K 0 ) of stetindite, HPLS, and high-pressure scheelite phases were determined, respectively, as 171(5), 105(4), and 221(40) GPa by fitting to a second-order Birch-Murnaghan equation of state. The pressure derivatives of vibrational modes and Gru''neisen parameters of the zircon-structured polymorph are similar to those of other orthosilicate minerals. In conclusion, due to the larger ionic radii of Ce 4+ , with respect to Zr 4+ , stetindite was found to possess a softer bulk modulus and undergo the phase transitions at a lower pressure than zircon (ZrSiO 4 ), such observations are consistent with what were found in coffinite (USiO 4 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Synthesis and Properties of the Helium Clathrate and Defect Perovskite [He 2– x $^\square$ x ][CaNb]F 6

The defect double perovskite [He 2– x $^\square$ x ][CaNb]F 6 , with helium on its A-site, can be prepared by the insertion of helium into ReO 3 -type CaNbF 6 at high pressure. Upon cooling from 300 to 100 K under 0.4 GPa helium, ~60% of the A-sites become occupied. Helium uptake was quantified by both neutron powder diffraction and gas insertion and release measurements. After the conversion of gauge pressure to fugacity, the uptake of helium by CaNbF 6 can be described by a Langmuir isotherm. The enthalpy of absorption for helium in [He 2– x $^\square$ x ][CaNb]F 6 is estimated to be ~+3(1) kJ mol –1 , implying that its formation is entropically favored. Helium is able to diffuse through the material on a time scale of minutes at temperatures down to ~150 K but is trapped at 100 K and below. The insertion of helium into CaNbF 6 reduces the magnitude of its negative thermal expansion, increases the bulk modulus, and modifies its phase behavior. On compressing pristine CaNbF 6 , at 50 and 100 K, a cubic (Fm$\overline{3}$m) to rhombohedral (R$\overline{3}$) phase transition was observed at <0.20 GPa. However, a helium-containing sample remained cubic at 0.4 GPa and 50 K. CaNbF 6 , compressed in helium at room temperature, remained cubic to >3.7 GPa, the limit of our X-ray diffraction measurements, in contrast to prior reports that upon compression in a nonpenetrating medium, a phase transition is detected at ~0.4 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Coupled Chemical and Mechanical Control of Phase Stability in Lanthanide-Substituted BiVO 4

Doping is widely used to enhance the photoelectrochemical performance of BiVO 4 , yet solubility limits and polymorphic stability constrain compositional tuning. Here, in this study, the role of trivalent cation substitution (Ln = La, Nd, Dy, Ho, Y) on pressure-induced phase transformations in Bi 1–x Ln x VO 4 (x ≤ 0.5) is described. Powder X-ray and neutron diffraction reveal that increasing Ln content stabilizes the tetragonal zircon-type polymorph under ambient conditions, while applied pressures of up to ∼5 GPa promote conversion to the monoclinic fergusonite-type polymorph. In-situ neutron diffraction on Bi 0.8 La 0.2 VO 4 shows a reversible monoclinic to tetragonal transition near 2–3 GPa with a bulk modulus of 147 GPa. The extent of conversion depends strongly on dopant identity, concentration, and synthetic route, with mixed-phase solid-state samples converting more efficiently than phase-pure coprecipitated materials. These results demonstrate pressure as a viable pathway to access metastable, doped BiVO 4 compositions beyond conventional solubility limits.

Sypkes, Kathryn I. [University of Sydney, NSW (Aus↗

DOE OSTI2026

Polymorphism of Magnesium Chloride Decahydrate under High Pressures and Low Temperatures Found in Icy Moons

Magnesium chloride hydrates show a rich diversity of structural forms, with different hydration numbers ranging from 1 to 12. Here, we show that this structural versatility is also present within the same hydration numbers, with the formation of a new phase of magnesium chloride decahydrate (MgCl 2 ·10H 2 O-II) stable at the conditions of 0.03–0.40(14) GPa at 234(11) K and up to 2.11(3) GPa at 220(2) K, as determined by in situ single-crystal X-ray diffraction. The hydrogen bonding interactions of MgCl 2 ·10H 2 O-II are identified and the bulk modulus determined at B 0 = 18.9(12) GPa at 220(2) K. We discuss the implication of this high-pressure MgCl 2 hydrate for icy moon compositional investigations and showcase the structural and density changes between the various hydrates in MgCl 2 ·nH 2 O.

high pressure↗

DOE OSTI2024

Critical Role of Framework Flexibility and Disorder in Driving High Ionic Conductivity in LiNbOCl 4

Understanding Li-ion transport is key for the rational design of superionic solid electrolytes with exceptional ionic conductivities. LiNbOCl 4 is reported to be one of the most highly conducting materials in the recently realized new class of soft oxyhalide solid electrolytes, exhibiting an ionic conductivity of ~11 mS·cm -1 . Here, we apply X-ray/neutron diffraction and pair distribution function analysis - coupled with density functional theory/ab-initio molecular dynamics - to determine a structural model that provides a rationale for the high conductivity that we observe experimentally in this nanocrystalline solid. We show that it arises from unusually high framework flexibility at room temperature. This owes to isolated 1-D [NbOCl 4 ] - anionic chains which exhibit energetically favorable orientational disorder that is - in turn - correlated to multiple, disordered and equi-energetic Li + sites in the lattice. As the Li-ions sample the 3-D energy landscape with a fast predicted diffusion coefficient of 5.1 x 10 -7 cm 2 /s at room temperature (σ i calc = 17.4 mS·cm -1 ), the inorganic polymer chains can reorient or vice versa. The activation energy barrier for Li migration through the frustrated energy landscape is especially reduced by the elastic nature of the NbO 2 Cl 4 octahedra evident from very widely dispersed Cl-Nb-Cl bond angles in AIMD snapshots at 300 K. The phonon spectra are predominantly influenced by Cl vibrations in the low energy range, and there is strong overlap between the framework (Cl, Nb) and Li partial pDOS in the region between 1.2 - 4.0 THz. The framework flexibility is also reflected in a relatively low bulk modulus of 22 GPa. In conclusion, our findings pave the way for investigation of future “flex-ion” inorganic solids and open up a new direction for the design of high conductivity, soft solid electrolytes for all-solid-state batteries.

AIMD↗

DOE OSTI2025

Symmetrization of Strong Hydrogen Bond under High Pressure in Bihydroxide-Ion-Containing NaCu 2 (SO 4 ) 2 ·H 3 O 2 Revealed by Experimental Charge Density, Single-Crystal Electron Diffraction, and Neutron Diffraction Studies

In minerals and inorganic compounds, strong hydrogen bonding can lead to the formation of complex ionic species such as the H 3 O 2 – bihydroxide anion and Zundel cation H 5 O 2 + . We studied [NaCu 2 (SO 4 ) 2 ·H 3 O 2 ] natrochalcite, which contains bihydroxide anions and undergoes hydrogen bond symmetrization at the lowest pressure reported so far among inorganic compounds. Hydrogen bond symmetrization leads to changes in the bulk modulus, seismic wave velocities, and proton mobility and plays a primary role in high-temperature superconductivity, but its characteristics are not well understood due to a lack of systematic studies and limitations of experimental methods sensitive to this subtle change. In this work, we applied experimental charge density analysis based on in situ single-crystal X-ray diffraction data, along with the single-crystal neutron and electron diffraction experiments, to probe the behavior of hydrogen atoms during the hydrogen bond symmetrization process under high-pressure conditions. On the way to the symmetrical H-bonding, natrochalcite undergoes a series of complex redistributions of electron density, which we trace with multipole refinement and detailed analysis of changes in the Laplacian of electron density values. Additionally, we deconvoluted the equation of state (volume of the unit cell vs pressure relation) into the atomic equation of states describing dependencies of atomic charges or volumes vs pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗

DOE OSTI2023

The Post‐Perovskite Transition in Fe‐ and Al‐Bearing Bridgmanite: Effects on Seismic Observables

Abstract The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al 3+ , Fe 3+ )‐, (Al 3+ )‐, (Fe 2+ )‐, and (Fe 3+ )‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe 3+ ‐, Al 3+ ‐, or (Al 3+ , Fe 3+ )‐alloying increases the transition pressure, while Fe 2+ ‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO 3 , and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.

Geochemistry & Geophysics↗

DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Vickers hardness prediction from machine learning methods

Abstract The search for new superhard materials is of great interest for extreme industrial applications. However, the theoretical prediction of hardness is still a challenge for the scientific community, given the difficulty of modeling plastic behavior of solids. Different hardness models have been proposed over the years. Still, they are either too complicated to use, inaccurate when extrapolating to a wide variety of solids or require coding knowledge. In this investigation, we built a successful machine learning model that implements Gradient Boosting Regressor (GBR) to predict hardness and uses the mechanical properties of a solid (bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio) as input variables. The model was trained with an experimental Vickers hardness database of 143 materials, assuring various kinds of compounds. The input properties were calculated from the theoretical elastic tensor. The Materials Project’s database was explored to search for new superhard materials, and our results are in good agreement with the experimental data available. Other alternative models to compute hardness from mechanical properties are also discussed in this work. Our results are available in a free-access easy to use online application to be further used in future studies of new materials at www.hardnesscalculator.com .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Phase comparison and equation of state for Ta 2 O 5

Abstract Tantalum pentoxide (Ta 2 O 5 ) is among the most technologically useful heavy transition metal oxides. Unfortunately, its crystal structure is the subject of long-standing and unresolved disagreement. Among other consequences, this uncertainty has made it impossible to formulate a robust high pressure equation of state for Ta 2 O 5 . Here, we report the results of high pressure x-ray diffraction experiments indexed against a comprehensive list of proposed Ta 2 O 5 phases. Five of the proposed phases produce good matches to experimental observations, and the compressibility parameters for these phases are all consistent within uncertainty. This means that regardless of the particular phase chosen, the Ta 2 O 5 equation of state can be described with bulk modulus K 0 = 138 ± 3.68 GPa and pressure derivative K 0 ′ = 1.82 ± 0.45 . Combining these experimental results with new density-functional theory calculations allows us to identify the λ phase as the best structural model of Ta 2 O 5 at ambient conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Hybridization-driven strong anharmonicity in Yb-filled skutterudites

A high-pressure study of the structural and thermal transport properties is carried out on one of the most efficient filled skutterudites, $\mathrm{Yb_{0.3}Co_4Sb_{12}}$, to understand the relatively low thermal conductivity behavior in this family. By combining x-ray diffraction and Raman scattering measurements, we detect a phase transition at around 12.4 GPa. The mode Grüneisen parameters of the observed phonon modes are obtained from the determined bulk modulus and the phonon frequency shifts with pressure. The strong anharmonicity in this material is demonstrated by the obtained large average Grüneisen parameter. We also find the depressed group velocity within the low-frequency range, the flat guest mode avoided crossing with the acoustic-phonon mode, and the significant contribution of optical phonons. The hybridization of the guest atom and host lattice and the related enhanced anharmonicity are suggested to account for the low lattice thermal conductivity in the studied system. So, these findings provide new insight into how phonon-phonon interactions lower lattice thermal conductivity in this important thermoelectric family.

36 MATERIALS SCIENCE↗