Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic compression”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

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

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2025

Symmetry progression and band vs Mott character of CdPS 3 under pressure

Complex chalcogenides are renowned for their tunable electronic, magnetic, and optical properties under external stimuli. The MPX 3 family (M = Mn, Ni, Co, V; X = S, Se) is a platform for many exciting discoveries—especially under compression—although CdPS3 is thought to be different because the Cd center possesses a filled 4d shell, which precludes Mottness. Here, we combine synchrotron-based infrared absorbance and Raman scattering spectroscopies with diamond anvil cell techniques, complementary lattice dynamics calculations, and an analysis of the energy landscape to reveal a series of structural phase transitions in CdPS 3 . We find four distinct pressure-driven transitions, with low frequency modes detectable over the full 35 GPa range of our investigation. A group–subgroup analysis along with our first-principles calculations allows us to partially unravel the space group sequence. For instance, the first critical pressure is a monoclinic C2/m to trigonal $\overline{R}$3 transition at 10 GPa. Despite the softness and overall sensitivity to pressure, we do not locate an insulator-to-metal transition in this pressure range, indicating that the energy scale for gap closure is significantly higher than expected. We discuss these findings in terms of force-induced color change and Mott vs band character in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multicycle large-eddy simulations of a direct-injection hydrogen-fueled optical engine

Hydrogen (H 2 ) is a carbon-free chemical energy carrier and one promising solution for achieving effective decarbonization of the transportation sector, particularly for internal combustion engines (ICEs). With a focus on ICEs, and compared to port-fuel injection, direct injection (DI) of gaseous H 2 during the compression stroke offers potential advantages, which include backfire avoidance and reduction of preignition occurrence. In these last two decades, much research, experimental and numerical, has been devoted to understanding H 2 's mixing and combustion processes in ICEs. Computational fluid dynamics modeling efforts commonly rely on unsteady Reynolds-averaged Navier Stokes (URANS) turbulence frameworks, mostly due to their computational affordability. However, many authors have pointed out the opportunity to perform large-eddy simulations (LESs) to investigate the cyclic variability of H 2 engines and assess potential advantages of using LES in place of URANS, especially for lean operation. This study addresses this knowledge gap and presents a computational fluid dynamics (CFD) study of the H 2 DI process in an optical engine operating at relatively low tumble conditions, using multicycle LESs. In conclusion, the manuscript presents a thorough validation of the results against experimental data available from the literature as well as direct comparison with URANS, demonstrating the feasibility of multicycle LESs for CFD modeling of DI H 2 -fueled ICEs.

Direct injection↗
DOE OSTI2024

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultrafast structural dynamics of UV photoexcited cis , cis -1,3-cyclooctadiene observed with time-resolved electron diffraction

Conjugated diene molecules are highly reactive upon photoexcitation and can relax through multiple reaction channels that depend on the position of the double bonds and the degree of molecular rigidity. Understanding the photoinduced dynamics of these molecules is crucial for establishing general rules governing the relaxation and product formation. Here, in this study, we investigate the femtosecond time-resolved photoinduced excited-state structural dynamics of cis,cis-1,3-cyclooctadiene, a large-flexible cyclic conjugated diene molecule, upon excitation with 200 nm using mega-electron-volt ultrafast electron diffraction and trajectory surface hopping dynamics simulations. We tracked the photoinduced structural changes from the Franck–Condon region through the conical intersection seam to the ground state. Our findings revealed a novel primary reaction coordinate involving ring distortion, where the ring stretches along one axis and compresses along the perpendicular axis. The nuclear wavepacket remains compact along this reaction coordinate until it reaches the conical intersection seam, and it rapidly spreads as it approaches the ground state, where multiple products are formed.

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

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Vibrational Dynamics and Phase Transitions of Hydrazine to 50 GPa

The high pressure behavior of hydrazine, N 2 H 4 , has been investigated to 50 GPa at room temperature using infrared and Raman spectroscopy to explore pressure induced phase transitions and changes in hydrogen bonding. Three solid–solid phase transitions were detected at 11, 21, and 32 GPa on room temperature compression through dramatic changes in the lattice vibration and N–H stretching regions with increasing pressure in both measurement techniques. The transition to phase IV, which appears at 32 GPa, exhibits increased hydrogen bonding with significant hysteresis, persisting to 9 GPa on decompression. This work presents a detailed analysis of the pressure dependence of mode shifts and calculations of mode Grüneisen parameters as well as a determination of an approximate thermodynamic Grüneisen parameter. We compare these results to the behavior of other small molecular materials such as ammonia and water and explore the evolution of hydrogen bonding in hydrazine toward the symmetrically hydrogen bonded state, which has previously been suggested by theoretical computations.

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

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗
DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hybrid algorithm for the time-dependent Hartree–Fock method using the Yang–Baxter equation on quantum computers *

Abstract The time-dependent Hartree–Fock (TDHF) method is an approach to simulate the mean field dynamics of electrons within the assumption that the electrons move independently in their self-consistent average field and within the space of single Slater determinants. One of the major advantages of performing time dynamics within Hartree–Fock theory is the free fermionic nature of the problem, which makes TDHF classically simulatable in polynomial time. Here, we present a hybrid TDHF implementation for quantum computers. This quantum circuit grows with time; but with our recent work on circuit compression via the Yang–Baxter equation (YBE), the resulting circuit is constant depth. This study provides a new way to simulate TDHF with the aid of a quantum device as well as provides a new direction for the application of YBE symmetry in quantum chemistry simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

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

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Strain-modulated intercalated phases of Pb monolayer with dual periodicity in SiC(0001)-graphene interface

Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.

36 MATERIALS SCIENCE↗