SEARCH · Engineering Papers

Results for “stabilization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations

The effects of plasma shaping on the H-mode pedestal structure are investigated. High fidelity kinetic simulations of the neoclassical pedestal dynamics are combined with the magnetohydrodynamic (MHD) stability conditions for triggering edge localized mode (ELM) instabilities that limit the pedestal width and height in H-mode plasmas. We use the neoclassical kinetic XGC0 code [Chang et al., Phys. Plasmas 11, 2649 (2004)] to carry out a scan over plasma elongation and triangularity. As plasma profiles evolve, the MHD stability limits of these profiles are analyzed with the ideal MHD ELITE code [Snyder et al., Phys. Plasmas 9, 2037 (2002)]. In simulations with the XGC0 code, which includes coupled ion-electron dynamics, yield predictions for both ion and electron pedestal profiles. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities are discussed. For the discharges with higher elongation, it is found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continues to evolve to higher pedestal pressures and bootstrap currents until the peeling-ballooning stability conditions are satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes are triggered at earlier times. The plasma elongation is found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency is large, and the H-mode pedestal evolves rapidly. It is found that the temperature of neutrals in the scrape-off-layer (SOL) region can affect the dynamics of the H-mode pedestal buildup. But the final pedestal profiles are nearly independent of the neutral temperature. The elongation and triangularity affect the pedestal widths of plasma density and electron temperature profiles differently. This provides a new mechanism of controlling the pedestal bootstrap current and the pedestal stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Predicting the thermodynamic stability of double-perovskite halides from density functional theory

Recently, a series of double-perovskite halide compounds such as Cs 2 AgBiCl 6 and Cs 2 AgBiBr 6 have attracted intensive interest as promising alternatives to the solar absorber material CH 3 NH 3 PbI 3 because they are Pb-free and may exhibit enhanced stability. The thermodynamic stability of a number of double-perovskite halides has been predicted based on density functional theory (DFT) calculations of compound formation energies. In this paper, we found that the stability prediction can be dependent on the approximations used for the exchange-correlation functionals, e.g., the DFT calculations using the widely used Perdew, Burke, Ernzerhof (PBE) functional predict that Cs 2 AgBiBr 6 is thermodynamically unstable against phase-separation into the competing phases such as AgBr, Cs 2 AgBr 3 , Cs 3 Bi 2 Br 9 , etc., obviously inconsistent with the good stability observed experimentally. The incorrect prediction by the PBE calculation results from its failure to predict the correct ground-state structures of AgBr, AgCl, and CsCl. By contrast, the DFT calculations based on local density approximation, optB86b-vdW, and optB88-vdW functionals predict the ground-state structures of these binary halides correctly. Furthermore, the optB88-vdW functional is found to give the most accurate description of the lattice constants of the double-perovskite halides and their competing phases. Given these two aspects, we suggest that the optB88-vdW functional should be used for predicting thermodynamic stability in the future high-throughput computational material design or the construction of the Materials Genome database for new double-perovskite halides. As a result, using different exchange-correlation functionals has little influence on the dispersion of the conduction and the valence bands near the electronic bandgap; however, the calculated bandgap can be affected indirectly by the optimized lattice constant, which varies for different functionals.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Subtlety of TiO 2 phase stability: Reliability of the density functional theory predictions and persistence of the self-interaction error

TiO 2 is an important material with broad applications that can exist in different phases with dramatically different properties. Theoretical prediction of their polymorph energetics is therefore critical for the material design and for identifying thermodynamically accessible structures. Determining TiO 2 relative phase stabilities remains challenging for first-principles methods, and density functional theory is the only approach available for studying phase stabilities at finite temperatures with acceptable computational efficiency. Here, we show that density functional theory using the recently developed efficient strongly constrained and appropriately normed (SCAN) [Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] exchange-correlation functional for the first time predicts the phase stability in qualitative agreement with the experimental results at realistic conditions. Further analysis shows that the self-interaction error intrinsic in the density functional persists in the stability prediction. By correcting the self-interaction error through an empirical approach, SCAN predicts the relative stability as well as defect properties in excellent agreement with the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Interpreting AI for fusion: An application to plasma profile analysis for tearing mode stability

Artificial intelligence models have demonstrated strong predictive capabilities for various instabilities in fusion devices such as Tokamaks, including tearing modes (TM), edge localized modes, and disruptive events, but their opaque nature raises concerns about safety and trustworthiness when applied to fusion power plants. Here, we present a physics-based interpretation framework using a TM prediction model as a demonstration that is validated through a dedicated DIII-D TM avoidance experiment. By applying Shapley analysis, we identify how profiles such as rotation, temperature, and density contribute to the model's prediction of TM stability. Our analysis shows that in our experimental scenario, core electron temperature and rotation peaking play the primary role in TM stability, while density changes have smaller effects on stability. We show that off-axis ion temperature stabilizes TMs, suggesting that off-axis neutral beam heating can further stabilize this scenario. This work presents a generalizable ML-based event prediction methodology, from training to physics-driven interpretation, bridging the gap between physics understanding and opaque ML models.

Farre-Kaga, Hiro J. [Princeton Univ., NJ (United S↗

DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗

DOE OSTI2020

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗

DOE OSTI2021

Stability of beta-induced Alfvén eigenmodes (BAE) in DIII-D

Although the stability of ellipticity, toroidal and reversed-shear Alfvén eigenmodes (EAE, TAE, RSAE) are relatively well understood, less is known about the stability of lower-frequency modes such as the beta-induced Alfvén eigenmode (BAE) but, because they are often unstable in present devices and are implicated in fast-ion transport, understanding their stability is vital. BAE stability is studied in primarily weak or reversed shear DIII-D plasmas with sub-Alfvénic deuterium beams. Modes are classified based on electron cyclotron emission, beam emission spectroscopy, magnetics, and interferometer data. The study is limited to the initial two seconds of the discharge, where the evolving q profile provides an effective scan of the dependence of stability upon q. In a dedicated experiment, BAEs are unstable at times in the discharge when the minimum of the safety factor q min is close to a rational number. The observed mode frequencies are usually close to analytic estimates of the BAE accumulation point and the eigenfunction peaks in the vicinity of q min . Unstable BAEs usually occur in bursts that chirp rapidly in frequency. To isolate the importance of thermal and beam gradients in driving the modes, the beam and electron cyclotron heating power is altered for 50–100 ms durations in reproducible discharges. As expected from the resonance condition, BAEs depend sensitively on the beam power and injection geometry. Modes only persist for ~25 ms because the anisotropic beam population only interacts strongly with the modes over a relatively narrow range of q. Lastly, a database of over 1000 beam-heated discharges shows that BAEs are more likely to be unstable when the poloidal beta exceeds 0.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Toroidal modeling of energetic passing particle drift kinetic effects on tearing mode stability

Abstract Drift kinetic effects of the neutral beam injection induced passing energetic particles (EPs) on the linear stability of the n = 1 tearing mode (TM) (with the dominant poloidal harmonic of m = 2) are numerically investigated utilizing the MARS-K code (Liu et al 2008 Phys. Plasmas 15 112503), in a tokamak plasma with finite equilibrium pressure and anisotropic thermal transport. In the low plasma pressure regime, it is found that co- (counter-) passing EPs stabilize (destabilize) the TM, agreeing with previous studies. However, as the plasma pressure increases beyond a critical value, it is found that co-passing EPs also destabilize the mode. An in-depth analysis reveals that the net effect of co-passing EPs is a result of competition between the stabilizing contribution from the non-adiabatic drift kinetic terms and the destabilizing contribution associated with adiabatic terms, with the latter becoming more dominant at higher equilibrium pressure. Non-perturbative magnetohydrodynamic-kinetic hybrid modeling also finds that co- and counter-passing EPs modify the TM eigenfunction differently, with the counter-passing EPs enhancing the sideband harmonics. Furthermore, effects of the plasma resistivity and toroidal rotation, as well as that of the equilibrium distribution of EPs in the particle pitch angle space, are also investigated, showing asymmetric results on the TM stability between the co- and counter-passing EPs. The first order finite orbit width correction is found to be stabilizing with co-passing EPs and destabilizing with counter-passing particles. Finally, drift resonances between passing EPs and the TM induce finite frequency to the mode and generate finite net torques inside the plasma, due to the neoclassical toroidal viscosity and the Reynolds stress associated with 3D perturbations.

Physics↗

DOE OSTI2023

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Stability evaluation and mitigation strategies in advanced tokamaks using 3D MHD spectroscopy

Multi-modal, active 3D MHD spectroscopy is applied in high-performance advanced tokamak scenarios to study their stability time evolution, revealing an intriguing dependence on both and . A tailored applied 3D field provides a 3D plasma response to extract the growth rate of the least stable mode. The estimated growth rate finds a decrease in stability when the minimum in the safety factor (q) passes through 2.0 and reveals inherent risks of crossing an additional rational surface at integer q min , even above the usual q = 1 sawtooth condition. Based on this result, the potential scenario in which q min ~ 2 can be safely crossed during a more stable lower β N phase was investigated, and the improved stability of this scenario is confirmed by the estimated growth rate. This shows that 3D MHD spectroscopy can offer insights into strategies for improving stability by identifying the vulnerable aspects of such scenarios. In addition, the method highlights its potential for instability avoidance by enabling early detection of multiple modes, even before magnetic coils can measure them. The measured growth rate by the 3D MHD spectroscopy shows its reliability by exhibiting a correlation with the programmed rises in plasma beta across various high β N and high q min discharges. In addition, this method is successfully applied during rapidly evolving I p ramp-up phases, a key part of the scenario development. By achieving reasonable growth rate measurements at high-performance scenario developments, this technique contributes to the development of advanced diagnostic tools for tokamak scenario stability, which will help identify an effective pathway to stable, high-performance scenarios.

Yang, S. M. [Princeton Plasma Physics Laboratory (↗

DOE OSTI2026

Clifford Hierarchy Stabilizer Codes: Transversal Non-Clifford Gates and Magic States

A fundamental problem in fault-tolerant quantum computation is the tradeoff between universality and dimensionality, exemplified by the the Bravyi-König bound for $n$-dimensional topological stabilizer codes. In this work, we extend topological Pauli stabilizer codes to a broad class of $n$-dimensional Clifford hierarchy stabilizer codes. These codes correspond to the $(n+1)$D Dijkgraaf-Witten gauge theories with non-Abelian topological order. We construct transversal non-Clifford gates through automorphism symmetries represented by cup products. In 2D, we obtain the first transversal non-Clifford logical gates including T and CS for Clifford stabilizer codes, using the automorphism of the twisted $\mathbb{Z}_2^3$ gauge theory (equivalent to $\mathbb{D}_4$ topological order). We also combine it with the just-in-time decoder to fault-tolerantly prepare the logical T magic state in $O(d)$ rounds via code switching. In 3D, we construct a transversal logical $\sqrt{\text{T}}$ gate in a non-Clifford stabilizer code at the third level of the Clifford hierarchy, located on a tetrahedron corresponding to a twisted $\mathbb{Z}_2^4$ gauge theory. Furthermore, our constructions surpass the Bravyi-König bound by achieving the logical gates in the $(n+1)$-th level of Clifford hierarchy in $n$ spatial dimension.

Kobayashi, Ryohei [Institute for Advanced Study, P↗

DOE OSTI2011

Entanglement-assisted codeword stabilized quantum codes

Entangled qubits can increase the capacity of quantum error-correcting codes based on stabilizer codes. In addition, by using entanglement quantum stabilizer codes can be construct from classical linear codes that do not satisfy the dual-containing constraint. We show that it is possible to construct both additive and nonadditive quantum codes using the codeword stabilized quantum code framework. Nonadditive codes may offer improved performance over the more common stabilizer codes. Like other entanglement-assisted codes, the encoding procedure acts only on the qubits on Alice's side, and only these qubits are assumed to pass through the channel. However, errors in the codeword stabilized quantum code framework give rise to effective Z errors on Bob's side. We use this scheme to construct entanglement-assisted nonadditive quantum codes, in particular, ((5,16,2;1)) and ((7,4,5;4)) codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Structural and ferromagnetic properties of an orthorhombic phase of MnBi stabilized with Rh additions

The article addresses the possibility of alloy elements in MnBi which may modify the thermodynamic stability of the NiAs-type structure without significantly degrading the magnetic properties. The addition of small amounts of Rh and Mn provides an improvement in the thermal stability with some degradation of the magnetic properties. The small amounts of Rh and Mn additions in MnBi stabilize an orthorhombic phase whose structural and magnetic properties are closely related to the ones of the previously reported high-temperature phase of MnBi (HT MnBi). The properties of the HT MnBi, which is stable between 613 and 719 K, have not been studied in detail because of its transformation to the stable low-temperature MnBi (LT MnBi), making measurements near and below its Curie temperature difficult. The Rh-stabilized MnBi with chemical formula Mn 1.0625–x Rh x Bi [x=0.02(1)] adopts a new superstructure of the NiAs/Ni 2 In structure family. It is ferromagnetic below a Curie temperature of 416 K. The critical exponents of the ferromagnetic transition are not of the mean-field type but are closer to those associated with the Ising model in three dimensions. The magnetic anisotropy is uniaxial; the anisotropy energy is rather large, and it does not increase when raising the temperature, contrary to what happens in LT MnBi. The saturation magnetization is approximately 3μB/f.u. at low temperatures. Thus, while this exact composition may not be application ready, it does show that alloying is a viable route to modifying the stability of this class of rare-earth-free magnet alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Resolving phase stability in the Ti-O binary with first-principles statistical mechanics methods

The Ti-O system consists of a multitude of stable and metastable oxides that are used in wide ranging applications. In this work we investigate phase stability in the Ti-O binary from first principles. We perform a systematic search for ground state structures as a function of oxygen concentration by considering oxygen-vacancy and/or titanium-vacancy orderings over four parent crystal structures: (i) hcp Ti, (ii) ω-Ti, (iii) rocksalt, and (iv) hcp oxygen containing interstitial titanium. We explore phase stability at finite temperature using cluster expansion Hamiltonians and Monte Carlo simulations. The calculations predict a high oxygen solubility in hcp Ti and the stability of suboxide phases that undergo order-disorder transitions upon heating. Vacancy ordered rocksalt phases are also predicted at low temperature that disorder to form an extended solid solution at high temperatures. Predicted stable and metastable phase diagrams are qualitatively consistent with experimental observations, however, important discrepancies are revealed between first-principles density functional theory predictions of phase stability and the current understanding of phase stability in this system.

36 MATERIALS SCIENCE↗