Search NASA⌕ Search

SEARCH · Search NASA

Results for “lattice constants”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Instability and evolution of the magnetic ground state in metallic perovskites GdRh 3 C 1-x B x

In this paper, we report investigations of the structural, magnetic, electrical transport and thermal properties of five compositions of the metallic perovskite GdRh 3 C 1–x B x (0.00 ≤ x ≤ 1.00). Our results show that all the five compositions undergo magnetic ordering at low temperatures, but the nature of the ordered state is significantly different in the carbon- and the boron-rich compositions, where the former shows signatures of an amplitude-modulated magnetic structure and the latter exhibits evidences of an equal-moment incommensurate antiferromagnetic ordering. We also observe a remarkable field-dependent evolution of conduction carrier polarization in the compositionally disordered compounds. The outcomes indicate that this system is energetically situated in proximity to a magnetic instability where small variations in the control parameter(s), such as lattice constant and/or electron density, lead to considerably different ground states.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Long-range hydrogen-binding effects of carbide interfaces in iron

A micromechanics model was developed to evaluate the elastic binding energy between carbide precipitates and hydrogen interstitials using Eshelby's equivalent inclusion method. In this work, density functional theory (DFT) simulations were performed to obtain the material-specific quantities, e.g., lattice constants and the elastic constants, for the continuum model. Using this model, we find that for coherent carbide precipitates, hydrogen atoms are more likely to bind on the broad surfaces of the disk-like precipitates, which is consistent with experimental observations. For semicoherent and incoherent precipitates, our model suggests that it is possible for semicoherent precipitates to have significant hydrogen binding capability while there is no hydrogen-binding capability of incoherent precipitates, which also agrees with experimental findings. In addition, several factors that influence the binding energies between hydrogen atoms and carbide precipitates were quantitatively analyzed, including the precipitate size, morphology, orientation, and interface. These collective results include both the position and the value of the strongest hydrogen-binding interaction for a wide range of carbide stoichiometries, which contributes to our understanding of hydrogen trapping in steel-based materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solid state reduction of nickelate thin films

The square-planar nickelates are a class of superconductors analogous to the cuprates and promise to provide insight into the pairing mechanism in high-temperature superconducting oxides. The parent phase of doped superconducting films is NdNiO 2 , which is prepared by reducing 3 + Ni in NdNiO 3 films. In this paper, we develop an ultrahigh vacuum reduction method using aluminum deposited on top of the 3 + nickelates and monitor the reduction process in real time through in situ crystal truncation rod measurements and diffraction x-ray absorption near edge spectroscopy measurements across the Ni K edge. Further, we establish a relation between Ni valence and the lattice constant of NdNiO 3-x and show that the process can precisely control the oxygen content in the films. Finally, we extend the reduction process to quintuple square-planar nickelates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Random forest prediction of crystal structure from electron diffraction patterns incorporating multiple scattering

Diffraction is the most common method to solve for unknown or partially known crystal structures. However, it remains a challenge to determine the crystal structure of a new material that may have nanoscale size or heterogeneities. Here, in this study, we train an architecture of hierarchical random forest models capable of predicting the crystal system, space group, and lattice parameters from one or more unknown two-dimensional electron diffraction patterns. Our initial model correctly identifies the crystal system of a simulated electron diffraction pattern from a 20-nm-thick specimen of arbitrary orientation 67% of the time. We achieve a topline accuracy of 79% when aggregating predictions from ten patterns of the same material but different zone axes. The space group and lattice predictions range from 70% to 90% accuracy and median errors of 0.01-0.5Å, respectively, for cubic, hexagonal, trigonal, and tetragonal crystal systems while being less reliable on orthorhombic and monoclinic systems. We apply this architecture to a four-dimensional scanning transmission electron microscopy scan of gold nanoparticles, where it accurately predicts the crystal structure and lattice constants. These random forest models can be used to significantly accelerate the analysis of electron diffraction patterns, particularly in the case of unknown crystal structures. Additionally, due to the speed of inference, these models could be integrated into live transmission electron microscopy experiments, allowing real-Time labeling of a specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tunable polar distortions and magnetism in Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb epitaxial films

Hexagonal 𝐴⁢𝐵⁢𝐶 intermetallics are predicted to have tunable ferroelectric, topological, and magnetic properties as a function of the polar buckling of 𝐵⁢𝐶 atomic planes. Here, we report the impact of isovalent lanthanide substitution on the buckling, structural phase transitions, and electronic and magnetic properties of Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films grown by molecular beam epitaxy (MBE) on 𝑐 plane sapphire substrates. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films form a solid solution from 𝑥=0 to 𝑥=1 and retain the polar hexagonal structure (𝑃⁢6 3 ⁢𝑚⁢𝑐) out to 𝑥 ≤ 0.95. With increasing 𝑥, the PtSb buckling increases and the out-of-plane lattice constant 𝑐 decreases due to the lanthanide contraction. While hexagonal LaPtSb is a highly conductive polar metal, the carrier density decreases with 𝑥 until an abrupt phase transition to a zero band overlap semimetal is found for cubic GdPtSb at 𝑥=1. The magnetic susceptibility peaks at small but finite 𝑥, which we attribute to Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling between localized 4⁢𝑓 moments, whose concentration increases with 𝑥, and free carriers that decrease with 𝑥. Samples with 𝑥≥0.3 show antiferromagnetic Curie-Weiss behavior and a Neel temperature that increases with 𝑥. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb system provides opportunities to dramatically alter the polar buckling and concentration of local 4⁢𝑓 moments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Straintronics with single-layer MoS 2 : A quantum Monte Carlo study

Using state-of-the-art quantum Monte Carlo (QMC) methods, we study straintronic properties of a single MoS 2 monolayer. 2D MoS 2 is a quintessential straintronic material for which many experiments have been performed. First, we determine the equilibrium atomic structure which is not known experimentally and is strictly needed to correctly determine the straintronic properties. That enables us to precisely analyze the quasiparticle band gaps for any applied biaxial strain, which we describe by a bivariate paraboloid function of lattice constant and internal structural parameter. Using the fixed-node QMC calculations fitted by analytical formulas, we localize the following excited state crossings between the direct, K → K , and indirect Γ → X and K → K /2 excitations. Based on this highly accurate many-body treatment, we predict a gauge factor of 136 meV/% for the K → K transition and a fairly narrow window of ≈ 2.8 % from compressive to mildly tensile strains, accounting for only ≈ 0.3 eV band gap change maintaining the direct character of the gap. Consequently, we suggest that, compared to other 2D materials, such as phosphorene, there is only a limited straintronic tuneability in this material often studied for its straintronic properties. QMC results are compared to results of standard DFT modeling, which reveal insights into the corresponding inaccuracies and therefore open a window for educated use of rapid DFT approaches. Published by the American Physical Society 2024

Huang, Y.↗
DOE OSTI2021

Metamorphic Tunnel Junctions Grown Via MOCVD Designed for GaAs 0.75 P 0.25 /Si Tandem Solar Cells

A high-performance metamorphic Al 0.2 Ga 0.8 As 0.75 P 0.25 /GaAs 0.75 P 0.25 heterojunction tunnel junction structure was developed for application to monolithic epitaxial GaAs 0.75 P 0.25 /Si 1.7 eV/1.1 eV bandgap tandem solar cells produced via metal-organic chemical vapor deposition. Doping optimization focused on both bulk doping concentrations and mitigating the transient incorporation effects that otherwise causes non-abrupt doping profiles in the Te doped layer. Furthermore, these efforts resulted in a fully-relaxed metamorphic tunnel diode at the target lattice constant having a peak tunneling current density (JP) of 279.1 A·cm –2 and a zero-bias resistance area product (RA) of 3.0×10 –4 Ω·cm 2 . After post-growth annealing to emulate the thermal load of the GaAs 0.75 P 0.25 top cell growth and elimination of the performance enhancing carbon activation anneal, the device achieved JP = 13.1 A·cm –2 and RA = 1.5×10 –3 Ω·cm 2 . These worst-case values are very promising as they enable the operation of GaAs 0.75 P 0.25 /Si tandems under one-sun AM1.5G illumination with negligible series resistance and at AM1.5D concentrations up to 240 suns with only 0.1% absolute efficiency loss.

14 SOLAR ENERGY↗
DOE OSTI2021

Fabrication of Quasi-Random Photonic Crystals for Ultrathin Solar Cells

Quasi-random structures for light management present an intermediate ordering between a completely random rough surface and a perfect ordered photonic crystal. They are attractive for broadband applications due to their power spectral density. This makes them excellent candidates for photovoltaic applications. Here we present preliminary fabrication results creating quasi-random structures using a self-assembly approach based on polymer blends and wet etching. We are able to create patterns with average lattice constant of five to two microns. This structure is proposed to be used as a rear patterning layer for ultrathin GaAs cells for better sustainability and radiation- hard devices.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2022

Stranger than metals

In traditional metals, the temperature ( T ) dependence of electrical resistivity vanishes at low or high T , albeit for different reasons. Here, we review a class of materials, known as “strange” metals, that can violate both of these principles. In strange metals, the change in slope of the resistivity as the mean free path drops below the lattice constant, or as T → 0, can be imperceptible, suggesting continuity between the charge carriers at low and high T . We focus on transport and spectroscopic data on candidate strange metals in an effort to isolate and identify a unifying physical principle. Special attention is paid to quantum criticality, Planckian dissipation, Mottness, and whether a new gauge principle is needed to account for the nonlocal transport seen in these materials.

Science & Technology - Other Topics↗
DOE OSTI2022

LAVA

We introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion.

Dang, Khanh↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Metalens formed by structured arrays of atomic emitters

Abstract Arrays of atomic emitters have proven to be a promising platform to manipulate and engineer optical properties, due to their efficient cooperative response to near‐resonant light. Here, we theoretically investigate their use as an efficient metalens. We show that, by spatially tailoring the (subwavelength) lattice constants of three consecutive two‐dimensional arrays of identical atomic emitters, one can realize a large transmission coefficient with arbitrary position‐dependent phase shift, whose robustness against losses is enhanced by the collective response. To characterize the efficiency of this atomic metalens, we perform large‐scale numerical simulations involving a substantial number of atoms (N∼ 5 × 10 5 ) that is considerably larger than comparable works. Our results suggest that low‐loss, robust optical devices with complex functionalities, ranging from metasurfaces to computer‐generated holograms, could be potentially assembled from properly engineered arrays of atomic emitters.

Materials Science↗
DOE OSTI2023

Temperature Control of Crystal Optics for Ultrahigh-Resolution Applications

The temperature control of crystal optics is critical for ultrahigh resolution applications such as those used in meV-resolved Inelastic Scattering. Due to the low count rate and long acquisition time of these experiments, for 1-meV energy resolution, the absolute temperature stability of the crystal optics must be maintained below 4 mK to ensure the required stability of lattice constant, thereby ensuring the energy stability of the optics. Furthermore, the temperature control with sub-mK precision enables setting the absolute temperature of individual crystal, making it possible to align the reflection energy of each crystal’s rocking curve in sub-meV resolution thereby maximizing the combined efficiency of the crystal optics. In this contribution, we report the details of an EPICS control system using PT1000 sensors, Keithley 3706A 7.5 digits sensor scanner, and Wiener MPOD LV power supply for the analyzer crystals of the Inelastic X-ray Scattering (IXS) beamline 10-ID at NSLS-II**. We were able to achieve absolute temperature stability below 1 mK and sub-meV energy alignment for several asymmetrically cut analyzer crystals. The EPICS ePID record was used for the control of the power supplies based on the PT1000 sensor input that was read with 7.5 digits accuracy from the Keithley 3706A scanner. The system enhances the performance of the meV-resolved IXS spectrometer with currently a 1.4 meV total energy resolution and unprecedented spectral sharpness for studies of atomic dynamics in a broad range of materials.

Gofron, Kaz↗
DOE OSTI2023

Atomistic insight into the ferroelastic post-stishovite transition by high-pressure single-crystal X-ray diffraction

Abstract The post-stishovite transition is a classic pseudo-proper typed ferroelastic transition with a symmetry-breaking spontaneous strain. This transition has been studied using high-pressure spontaneous strains, optic modes, and elastic moduli (Cij) based on the Landau modeling, but its atomistic information and structural distortion remain poorly understood. Here we have conducted synchrotron single-crystal X-ray diffraction measurements on stishovite crystals up to 75.3 GPa in a diamond-anvil cell. Analysis of the data reveals atomic positions, bond lengths, bond angles, and variations of SiO6 octahedra across the transition at high pressure. Our results show that the O coordinates split at ~51.4 GPa, where the apical and equatorial Si-O bond lengths cross over, the SiO6 octahedral distortion vanishes, and the SiO6 octahedra start to rotate about the c axis. Moreover, distortion mode analysis shows that an in-plane stretching distortion (GM1+ mode) occurs in the stishovite structure at high pressure while a rotational distortion (GM2+ mode) becomes dominant in the post-stishovite structure. These results are used to correlate with elastic moduli and Landau parameters (symmetry-breaking strain e1–e2 and order parameter Q) to provide atomistic insight into the ferroelastic transition. When the bond lengths of two Si-O bonds are equal due to the contribution from the GM1+ stretching mode, C11 converges with C12, and the shear wave VS1[110] polarizing along [110] and propagating along [110] vanishes. Values of e1–e2 and Q are proportional to the SiO6 rotation angle from the occurrence of the GM1+ rotational mode in the post-stishovite structure. Our results on the pseudo-proper type transition are also compared with that for the proper type in albite and improper type in CaSiO3 perovskite. The symmetry-breaking strain, in all these types of transitions, arises as the primary effect from the structural angle (such as SiO6 rotation or lattice constant angle) and its relevant distortion mode in the low-symmetry ferroelastic phase.

Geochemistry & Geophysics↗