Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase kinetics”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Neural network interatomic potential-driven analysis of phase stability in Ti–V alloys at the atomistic scale

The evolution of the ω phase in titanium–vanadium (Ti–V) alloys is critical for their mechanical properties, particularly in aerospace and biomedical applications. Here, this study employs a Rapid Artificial Neural Network (RANN) potential to model the ω phase evolution at the atomistic level, demonstrating a high degree of consistency with experimental observations, unlike the Modified Embedded Atom Method (MEAM), which fails to capture this phase transformation accurately. RANN simulations replicate key phenomena such as the nucleation of α precipitates at ω/β interfaces and accurate lattice orientations, enhancing our understanding of phase stability and transformation kinetics. The findings affirm that RANN potentials can significantly improve the prediction accuracy of complex material behaviors, offering a powerful tool for designing advanced materials with tailored properties such as solute effect in various stacking fault energies. This approach not only bridges the gap between theoretical predictions and empirical data but also sets a new direction for future research in materials science, emphasizing the integration of machine learning techniques in the development and optimization of new alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2022

Multiscale Interaction of Turbulence, Temperature, and Soot Formation: Measurements for Critical Assessments of Chemical Kinetics and Mechanisms

This program targeted an improved understanding of how turbulence affects soot formation chemistry (directly and indirectly) with a particular focus on detailing the multiscale coupling between flow turbulence, mixing, thermal transport, and soot formation kinetics in gas-phase reacting systems. In this project, simultaneous, quantitative multi-dimensional velocity, soot volume fraction, and gas-phase temperature measurements in turbulent non-premixed sooting flames were developed. These new diagnostics tools can be used to elucidate the relative effects of fluid kinematics, mixing, thermal processes, and turbulence/scalar time/length scales on soot formation chemistry and topology. Simultaneous multi-parameter measurements were acquired within a series of flames with different Reynolds numbers and spatial locations to investigate the inter-dependence of soot, temperature, and turbulence within the soot inception, growth, and oxidation region. Preliminary assessment of the results has been completed at the project end and detailed analysis is planned following the completed program.

02 PETROLEUM↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗
NASA NTRS1991

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗
DOE OSTI2020

Requirements and Conceptual Design of Off-gas Systems for the Reprocessing of Metallic Fuels

An assessment has been conducted to determine how key regulations regarding volatile radionuclide emissions to the atmosphere may apply to the off-gas streams associated with electrochemical reprocessing. The scope of this assessment was based upon a generic electrochemical reprocessing scheme with a throughput rate of 200 MTIHM/y applied to metallic fuel discharged from a sodium fast reactor (SFR), but the findings are able to be translated to other advanced nuclear scenarios as merited. Air dispersion modeling was performed using the EPA CAP-88 model and evaluated the uncontrolled decontamination factors (DFs) that would be required to achieve regulatory compliance with the dose-based limits set forth by EPA regulation 40 CFR 190.10(a). These DFs were compared to those required by fuel cycle–based limits set forth by EPA regulation 40 CFR 190.10(b). Two theoretical sites with disparate climatological conditions were selected for air dispersion modeling (Idaho and Tennessee). The radionuclides modeled included 3 H, 85 Kr, 129 I, and selected alpha-emitting transuranic isotopes (referred to here as 239 Pu-TRU <1y ). It was found that the fuel cycle-based limits in 40 CFR 190.10(b) are most restrictive for 85 Kr and 239 Pu-TRU <1y , with DFs of 3 and 6.1E+09, respectively. The dose-based limit as derived from 40 CFR 190.10(a) could require mitigation of tritium in some scenarios, with an estimated DF of about 3 for the reference scenarios. The fuel cycle-based limit for 129 I resulted in a DF of about 240 for the reference scenario. The need for iodine mitigation based on dose to the public depended upon the physical form of iodine as either particulate or vapor-phase species. Emission of iodine from the facility as a vapor necessitated DFs of about 2 but emission as a particulate would require DFs >6,000 to meet thyroid dose-based limits. Effects of physical form on needed iodine mitigation are significant, but the understanding of speciation of iodine both during electrochemical reprocessing and after release to the atmosphere is limited. The electrochemical processing unit operations were evaluated to identify potential release points for the volatile radionuclides and to assess the potential for retention of the radionuclides within the process (thus decreasing the need for mitigation). Mitigation strategies for 3 H, 85 Kr, 129 I, and 239 Pu-TRU <1y were identified. In all cases, there are reasonably achievable pathways to regulatory compliance, although in some cases additional R&D is merited to verify the chemical speciation of these isotopes and to develop and demonstrate potential treatment technologies for this application. Whether or not additional off-gas controls (beyond common operations such as HEPA filtration and oxygen and moisture control) are needed for any of these regulated or volatile radionuclides depends on the (a) type of facility (NRC-regulated or DOE), (b) used fuel process rate, (c) used fuel burnup and composition, (d) speciation and retention of volatile radionuclides in the process and in the cell gas cleanup system, (e) site-specific parameters such as location, meteorology, stack height, and site boundaries, and (f) levels of conservatism and safety factors used in assessing compliance to air emissions regulations. Performance of this assessment revealed several areas where information is lacking or additional research is required in order to better determine if or what kinds of off-gas control might be needed. First, and most significantly, the understanding of the chemical speciation and physical form and partitioning of iodine during electrochemical processing operations is lacking and prevents the ability to accurately assess the potential iodine mitigation requirements. Future research in this area should be multifaceted and include thermodynamic modeling of iodine speciation in different process steps, experiments to quantify the kinetics of vapor-phase and melt-phase transitions, bench-scale experiments to determine the potential chemical and physical form of iodine emissions from the electrorefining process, and verification of iodine behavior with experiments utilizing operational facilities. Similarly, an improved understanding of iodine behavior in the environment after release from the facility stack will be required to refine dose estimations, as particulate and vapor-phase emissions result in significantly different doses to the MEI.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Energizing Robust Sulfur/Lithium Electrochemistry via Nanoscale-Asymmetric-Size Synergism

Sluggish redox kinetics and dendrite growth perplex the fulfillment of efficient electrochemistry in lithium–sulfur (Li–S) batteries. The complicated sulfur phase transformation and sulfur/lithium diversity kinetics necessitate an all-inclusive approach in catalyst design. Herein, a compatible mediator with nanoscale-asymmetric-size configuration by integrating Co single atoms and defective CoTe 2–x (Co SA -CoTe 2–x @NHCF) is elaborately developed for regulating sulfur/lithium electrochemistry synchronously. Substantial electrochemistry and theoretical analyses reveal that CoTe 2–x exhibits higher catalytic activity in long-chain polysulfide transformation and Li 2 S decomposition, while monodispersed Co sites are more effective in boosting sulfur reduction kinetics to regulate Li 2 S deposition. Such cascade catalysis endows Co SA -CoTe 2–x @NHCF with the all-around service of “trapping-conversion-recuperation” for sulfur species during the whole redox reaction. Furthermore, it is demonstrated by in situ transmission electron microscopy that initially formed electronic-conductive Co and ionic-conductive Li 2 Te provide sufficient lithiophilic sites to regulate homogeneous Li plating and stripping with markedly suppressed dendrite growth. Consequently, by coupling the Co SA -CoTe 2–x @NHCF interlayer and Li@Co SA -CoTe 2–x @NHCF anode, the constructed Li–S full batteries deliver superior cycling stability and rate performance, and the flexible pouch cell exhibits stable cycling performance at 0.3 C. In conclusion, the gained insights into the synergistic effect of asymmetric-size structures pave the way for the integrated catalyst design in advanced Li–S systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Compression-rate dependence of pressure-induced phase transitions in Bi

Abstract It is qualitatively well known that kinetics related to nucleation and growth can shift apparent phase boundaries from their equilibrium value. In this work, we have measured this effect in Bi using time-resolved X-ray diffraction with unprecedented 0.25 ms time resolution, accurately determining phase transition pressures at compression rates spanning five orders of magnitude (10 –2 –10 3 GPa/s) using the dynamic diamond anvil cell. An over-pressurization of the Bi-III/Bi-V phase boundary is observed at fast compression rates for different sample types and stress states, and the largest over-pressurization that is observed is ΔP = 2.5 GPa. The work presented here paves the way for future studies of transition kinetics at previously inaccessible compression rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Unveiling the Morphological and Physical Mechanism of Burn–in Loss Alleviation by Ternary Matrix Toward Stable and Efficient All–Polymer Solar Cells

All–polymer solar cells (All–PSCs) are considered the most promising candidate in achieving both efficient and stable organic photovoltaic devices, yet the field has rarely presented an in–depth understanding of corresponding device stability while efficiency is continuously boosted via the innovation of polymer acceptors. Herein, a ternary matrix is built for all–PSCs with optimized morphology, improved film ductility and importantly, boosted efficiency and better operational stability than its parental binary counterparts, as a platform to study the underlying mechanism. The target system PQM–Cl:PTQ10:PY–IT (0.8:0.2:1.2) exhibits an alleviated burn–in loss of morphology and efficiency under light soaking, which supports its promoted device lifetime. The comprehensive characterizations of fresh and light–soaked active layers lead to a clear illustration of opposite morphological and physical degradation direction of PQM–Cl and PTQ10, thus resulting in a delicate balance at the optimal ternary system. Specifically, the enlarging tendency of PQM–Cl and shrinking preference of PTQ10 in terms of phase separation leads to a stable morphology in their mixing phase; the hole transfer kinetics of PQM–Cl:PY–IT host is stabilized by incorporating PTQ10. Finally, this work succeeds in reaching a deep insight into all–PSC's stability promotion by a rational ternary design, which booms the prospect of gaining high–performance all–PSCs.

14 SOLAR ENERGY↗
DOE OSTI2023

In Situ Neutron Scattering Study of the Structure Dynamics of the Ru/Ca 2 N:e – Catalyst in Ammonia Synthesis

NH 3 synthesis is one of the most critical industrial processes. Compared to commercial iron catalysts, Ru catalysts show high intrinsic activity in this reaction but suffer from hydrogen poisoning. By loading Ru onto supports such as electrides and hydrides, the hydrogen poisoning problem can be significantly alleviated. However, relevant studies on the structural dynamics of the Ru/electride catalysts under reaction conditions are very scarce. Taking advantage of the high sensitivity to hydrogen species, it is possible to obtain insights into the structural changes during the reaction using in situ neutron techniques. In this study, we have investigated the structural evolution of the Ru/Ca 2 N:e – catalyst during the ammonia synthesis reaction by in situ neutron scattering (inelastic neutron scattering, INS) technique. In situ INS experiments suggest that Ca 2 N:e – is likely converted to the Ca 2 NH phase during the reaction. Unlike the previously known structure where H and N atoms are intermixed, the formed Ca 2 NH exhibits a segregated structure where the H and N atoms are located in different layers separated by the Ca layer. In conclusion, density functional theory calculations of the reaction energetics reveal that there are minor changes in the barriers and thermodynamics of the first N hydrogenation step between the two phases (Ca 2 NH phase with segregated H/N layers and intermixed Ca 2 NH phase), suggesting the impact of the phases on the reaction kinetics to be relatively minimal.

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

Surface Crust Formation Controls Evaporation Kinetics of Secondary Organic Aerosols

Gas-particle partitioning is critical for the evolution of secondary organic aerosols (SOA) in the atmosphere. SOA particles evaporate more slowly than expected at nearly size-independent rates, but the underlying mechanism remains controversial. Here, in this study, we apply kinetic multilayer modeling to simulate evaporation of α-pinene SOA, demonstrating that surface crust formation, emerging from accumulation of low-volatility compounds at the particle surface, leads to slow evaporation and reduced size dependence of the evaporation rate. While evaporation induced by decomposition of oligomers would naturally lead to size-independent evaporation rates, we observe and simulate nearly size-independent slow evaporation of polyethylene glycol mixture particles containing polymeric species that do not decompose, confirming the relevance of composition-dependent diffusivity for size-independent, slow evaporation. Slow evaporation of limonene SOA was also observed in environmental chamber experiments, and model simulations demonstrate strong surface crust formation with bulk diffusivity being depressed by up to 5 orders of magnitude compared to the inner bulk. We present experimental evidence using a surface-based mass spectrometry technique that shows that the particle surface becomes enriched in high molecular weight compounds upon evaporation of monomers. Our findings imply that viscous surface crusts may also limit the growth and chemical transformation of SOA particles, influencing their impacts on air quality and climate.

heterogeneous chemistry↗
DOE OSTI2023

Vibrations and Phase Stability in Mixed Valence Antimony Oxide

α-Sb 2 O 4 (cervantite) and β-Sb 2 O 4 (clinocervantite) are mixed valence compounds with equal proportions of Sb III and Sb V as represented in the formula Sb III Sb V O 4 . Their structure and properties can be difficult to calculate owing to the Sb III lone-pair electrons. Here, we present a study of the lattice dynamics and vibrational properties using a combination of inelastic neutron scattering, Mössbauer spectroscopy, nuclear inelastic scattering, and density functional theory (DFT) calculations. DFT calculations that account for lone-pair electrons match the experimental densities of phonon states. Mössbauer spectroscopy reveals the β phase to be significantly harder than the α phase. Calculations with O vacancies reveal the possibility for nonstoichiometric proportions of Sb III and Sb V in both phases. An open question is what drives the stability of the α phase over the β phase, as the latter shows pronounced kinetic stability and lower symmetry despite being in the high-temperature phase. Since the vibrational entropy difference is small, it is unlikely to stabilize the α phase. In conclusion, our results suggest that the α phase is more stable only because the material is not fully stoichiometric.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗