Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELASTIC 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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Extending the capillary wave model to include the effect of bending rigidity: X-ray reflectivity and diffuse scattering

The surface roughness of a thin film at a liquid interface exhibits contributions of thermally excited fluctuations. This thermal roughness depends on temperature (𝑇), surface tension (𝛾), and elastic material properties, specifically the bending modulus (𝜅) of the film. A nonzero 𝜅 suppresses the thermal roughness at small length scales compared to an interface with zero 𝜅, as expressed by the power spectral density of the thermal roughness. The description of the x-ray scattering of the standard capillary wave model (CWM), which is valid for zero 𝜅, is extended to include the effect of 𝜅. The extended CWM (eCWM) provides a single analytical form for both the specular x-ray reflectivity (XRR) and the diffuse scattering around the specular reflection, and recovers the expression of the CWM at its zero 𝜅 limit. This theoretical approach enables the use of single-shot grazing incidence x-ray off-specular scattering (GIXOS) measurements for characterizing the structure of thin films on a liquid surface. The eCWM analysis approach decouples the thermal roughness factor from the surface scattering signal, providing direct access to the intrinsic surface-normal structure of the film and its bending modulus. Moreover, the eCWM facilitates the calculation of reflectivity at any desired resolution (pseudo-XRR approach). The transformation into pseudo-XRR provides the benefit of using widely available XRR software to perform GIXOS analysis. The extended range of the vertical scattering vector (𝑄 𝑧 ) available with the GIXOS pseudo-XRR approach allows for a higher spatial resolution than with conventional XRR. Experimental results are presented for various lipid systems, showing strong agreement between conventional specular XRR and pseudo-XRR methods. This agreement validates the proposed approach and highlights its utility for analyzing soft, thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

Here, we analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call “complexity reduction”. In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature T g of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2–6 nm at T g . Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

cancer↗
DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗
DOE OSTI2024

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultrasonic Resonance Techniques for Materials Research

Mechanical resonances are directly related to the physical behavior of a system at the bulk and microscopic levels. In materials science, resonant ultrasound spectroscopy (RUS) has long been a preferred nondestructive method to study mechanical resonances of solids and precisely measure quantitative material properties, namely elasticity. In recent years, advances in computational power and hardware have enabled RUS to be relevant for an increasing range of applications, such as advanced manufacturing. An extension of this technique, nonlinear RUS (NRUS), has been demonstrated to provide unmatched sensitivity to early-stage damage. NRUS was originally developed to probe geologic materials but has become a vital tool in nondestructive evaluation and materials research, offering a powerful means of quantifying and characterizing microstructural nonlinearity in a broad range of materials. This review summarizes recent developments and growth opportunities in RUS and NRUS techniques, modeling, and applications across a wide range of material systems including metals, composites, geomaterials, and explosives.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nanoscale Polarization-Dependent Young’s Modulus of Ferroelectric BaTiO 3 Single Crystals

Flexoelectric fields produced by strong strain gradients at the nanoscale couple to ferroelectric polarization, promoting changes in the mechanical properties of ferroelectric materials as a function of the direction of the ferroelectric polarization. In this work, we calculate the asymmetry in the Young’s modulus found in oppositely polarized out-of-plane domains of BaTiO 3 single crystals by means of contact resonance AFM, and we evaluate its impact on the electromechanical response as measured by piezoresponse force microscopy, both using band excitation modes. Furthermore, we analyze the electromechanical response of the different BaTiO 3 domains using k-means to build up their mean elastic and electromechanical features and quantify the flexoelectrically induced modulation of the Young’s modulus. Finally, we discuss the use of contact Kelvin probe force microscope measurements to decouple electrostatic artifacts from purely piezoelectric ones when flexoelectric coupling appears.

BaTiO 3 single crystals↗
DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗
DOE OSTI2025

Thick graded interfaces increase wear resistance in Ti/TiN nanolayered thin films

Multilayered composites with nanoscale layer thickness incorporating titanium and titanium nitride (Ti/TiN) are used as a model system to study the effects of heterophase interface structure on elastic and plastic deformation, as well as wear behavior. Here, in this work, hardness, modulus, and wear rate under dry reciprocating sliding contact are quantified as a function of Ti-TiN heterophase interfacial nitrogen gradient thickness for Ti/TiN multilayers with 10–80 nm layer thickness. Hardness and modulus are found to be inversely proportional to layer thickness and independent of interface gradient for most specimens. Wear rate is found to be inversely proportional to interface gradient thickness at constant layer thickness, demonstrating that control of nanoscale interface structure is a valid approach to enhancing wear behavior. The materials studied in this work wear comparably or slower than other Ti- and TiN-based composites in the literature, providing a promising avenue for engineering wear-resistant materials for use in industrially relevant applications.

Graded interfaces↗
DOE OSTI2026

Achieving ultrahigh modulus of resilience and enhanced thermal stability in ZnO x /SU-8 interpenetrating network polymer nanocomposite nanopillars

The modulus of resilience, a mechanical property that quantifies the maximum strain energy density a material can store during elastic deformation, is a crucial parameter for materials used in flexible displays, micro/nano-electro-mechanical system (M/NEMS) actuators, and ultra-sensitive pressure sensors. In this study, ZnO x /SU-8 nanocomposite nanopillars with a diameter of 300 nm, fully infiltrated with a uniformly distributed, interpenetrating amorphous ZnO x filler network, were synthesized via vapor-phase infiltration (VPI). In-situ uniaxial nano-compression tests revealed that the modulus of resilience of ZnO x /SU-8 reaches ∼ 12 MJ/m 3 , which is an ultrahigh value among all engineering materials with comparable strength. In addition, the synthesis fidelity, inorganic infiltration depth, and mechanical performance were all significantly improved compared to VPI-synthesized AlO x nanocomposites. Thermal stability, another key requirement for M/NEMS device materials operating under extreme environments, was also notably enhanced. Furthermore, partial crystallization of the amorphous ZnO x fillers during annealing contributed to an additional increase in modulus of resilience, reaching up to ∼ 13.9 MJ/m 3 . This work presents an effective fabrication strategy for producing nanostructured organic–inorganic hybrid nanocomposites with ultrahigh modulus of resilience and superior thermal stability, paving the way for their integration into next-generation flexible displays and high-performance M/NEMS devices working under harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elastic properties of the W 0.75 R e 0.25 alloy at high pressure up to 183 GPa

The high pressure equation of state for the W 0.75 R e 0.25 alloy is experimentally determined up to 183 GPa with synchrotron angle-dispersive powder x-ray diffraction in the diamond-anvil cell and to ∼925 GPa with density-functional theory. W-Re alloys are used in many industrial high-temperature applications and as a confining gasket material in high-pressure diamond-anvil cell research. The inclusion of 25 wt. % Re achieves the highest performance in terms of strength and ductility while also maintaining the body-centered-cubic (bcc) crystal structure, yet to date there has been no investigation into its elastic behavior at high pressure. We present the experimentally and theoretically determined volumetric and elastic pressure response and systematically compare these results to other W-Re alloys, finding that the bulk modulus of W-Re alloys varies nonlinearly with Re content and W 0.75 Re 0.25 becomes less incompressible than W at 85 GPa. Published by the American Physical Society 2025

Alloys↗
DOE OSTI2024

Ductile Glassy Polymer Networks Capable of Large Plastic Deformation and Heat-Induced Elastic Recovery

Many thermoplastic polymers are ductile by combining strength and large deformations. These deformations are irreversible - known as plastic deformation. Elastomers can deform reversibly but have low strength. To this end, we developed glassy and ductile polyamide networks capable of large plastic deformation (>200% strain) and high strength (~50 MPa tensile strength and ~1500 MPa Young’s modulus), similar to those of polyolefins and Nylon-66. We discovered that hydrogen bonding between meta-phthalamide groups was essential to the ductility. Since these polyamide networks are covalently bonded, we demonstrated their unique durability by repeatable elastic recovery at elevated temperatures, exhibiting indifferent tensile properties in each cycle. Furthermore, when we fixed the strain during the elastic recovery, these polyamide networks actuated stresses of 9–18 MPa, among the highest reported in shape-memory polymer actuators. Here, we envision these ductile, glassy polymer networks as promising alternatives to ductile thermoplastics, given the combined benefit of ductility and durability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Increased AGE Cross-Linking Reduces the Mechanical Properties of Osteons

Abstract The osteon is the primary structural component of bone, contributing significantly to its unique toughness and strength. Despite extensive research on osteonal structure, the properties of osteons have not been fully investigated, particularly within the context of bone fragility diseases like type 2 diabetes mellitus (T2DM). This study aims to isolate osteons from bovine bone, simulate the effects of increased advanced glycation end-products (AGEs) in T2DM through ribosylation, and evaluate the mechanical properties of isolated osteons. Osteons extracted from the posterior section of bovine femur mid-diaphysis were processed to achieve a sub-millimeter scale for microscale imaging. Subsequently, synchrotron radiation micro-computed tomography was employed to precisely localize and isolate the osteon internally. While comparable elastic properties were observed between control and ribosylated osteons, the presence of AGEs led to decreased strain to failure. Young’s modulus was quantified (9.9 ± 4.9 GPa and 8.7 ± 3 GPa, respectively), aligning closely with existing literature. This study presents a novel method for the extraction and isolation of osteons from bone and shows the detrimental effect of AGEs at the osteonal level.

Materials Science↗
DOE OSTI2025

Effects cascade debris and helium bubbles on the strength of aged plutonium

The radioactive decay of aging Pu is dominated by α-decay. This persistent α-decay produces crystalline defects in the form of dislocation loops and helium bubbles that evolve with time. Comparable defects are produced in other metallic alloys when subject to neutron irradiation, and these defects are known to modify the plastic deformation of irradiated materials. Models have been developed for these irradiated materials and validated against experimental confirmations of yield strength and the concomitant microstructural evolution. In this paper, we deploy those previously developed models and apply their mechanics to plutonium aging.

Chemical elements↗
DOE OSTI2025

Ultrasonic characterization of material heterogeneities in stainless steel components produced by laser powder bed fusion

We introduce pulse-echo ultrasound as a method for characterizing the impact of powder bed fusion parameters on the properties of additively manufactured stainless-steel components, their material anisotropy, and location-dependent heterogeneity. Our results indicate that accurate characterization requires careful selection of ultrasonic propagation paths, which must consider the direction of additive layering, variations in processing parameters, and the component's geometry. We employed two distinct methods to estimate material properties from ultrasonic data: One assumes isotropy, while the other accounts for anisotropic interactions during the propagation of elastic waves. When applied to samples fabricated with laser energy densities ranging from 24 to 42 J/mm³ , these methods revealed transverse isotropy and weak anisotropy (quantified by small Thomsen parameters, ε = 0.0651 and γ = 0.0092) and less than a ∼6 % change in acoustic impedance. The assumption of isotropy, in this case, leads to small errors (less than 4 % or 1 % for Young's modulus in the build or transverse directions) when estimating orthotropic material properties using ultrasonic data measured along just two orthogonal directions, one of which must align with the build direction. By comparing ultrasonic measurements — which aggregate the spatial variability in material properties along the length of elastic wave propagation into a single value — with localized measurements obtained from surface nanoindentation, we uncovered and spatially profiled significant differences between the surface and interior properties. Specifically, the surface Young's modulus decreased from approximately 210 GPa to 180 GPa within a depth of about 3 mm. We attribute this surface-localized heterogeneity in PBF-fabricated components to distinct thermal histories experienced by the surface and interior regions. Collectively, the results of this study establish a framework for the ultrasonic characterization of material heterogeneity and anisotropy in material properties and demonstrate its application in additively manufactured metal components.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear strength of a refractory high entropy alloy MoNbTaVW under high pressure

Radial X-ray diffraction (R-XRD) was performed in situ using a Panoramic Diamond Anvil Cell on the refractory high entropy alloy MoNbTaVW. The lattice parameters were determined through a Le Bail fit using a body centered cubic (BCC) lattice of symmetry $Im$$\overline{3}$$m$ (international space group number 229). Upper and lower bounds to the shear strength were determined up to 80 GPa nonhydrostatically using copper as a pressure standard. The equation of state was derived at the magic angle ψ = 54.7° and yielded a bulk modulus of K 0 = 220.8 ± 1.65 GPa. The experimental lattice parameters and bulk modulus match closely with corresponding density functional theory (DFT) calculations. The BCC phase remains stable up to the highest pressure of 80 GPa studied and is shown to be elastically anisotropic. The shear strength was found to saturate around 70 GPa with a value of τ = 1.75 GPa, and the shear moduli are presented in different limits of iso-strain and iso-stress.

36 MATERIALS SCIENCE↗