Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic analyses”

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 55 records · Page 3

Uncertainty analysis of correlated parameters in automated reaction mechanism generation

Abstract Uncertainty analysis is a useful tool for inspecting and improving detailed kinetic mechanisms because it can identify the greatest sources of model output error. Owing to the very nonlinear relationship between kinetic and thermodynamic parameters and computed concentrations, model predictions can be extremely sensitive to uncertainties in some parameters while uncertainties in other parameters can be irrelevant. Error propagation becomes even more convoluted in automatically generated kinetic models, where input uncertainties are correlated through kinetic rate rules and thermodynamic group values. Local and global uncertainty analyses were implemented and used to analyze error propagation in Reaction Mechanism Generator (RMG), an open‐source software for generating kinetic models. A framework for automatically assigning parameter uncertainties to estimated thermodynamics and kinetics was created, enabling tracking of correlated uncertainties. Local first‐order uncertainty propagation was implemented using sensitivities computed natively within RMG. Global uncertainty analysis was implemented using adaptive Smolyak pseudospectral approximations as implemented in the MIT Uncertainty Quantification Library to efficiently compute and construct polynomial chaos expansions to approximate the dependence of outputs on a subset of uncertain inputs. Cantera was used as a backend for simulating the reactor system in the global analysis. Analyses were performed for a phenyldodecane pyrolysis model. Local and global methods demonstrated similar trends; however, many uncertainties were significantly overestimated by the local analysis. Both local and global analyses show that correlated uncertainties based on kinetic rate rules and thermochemical groups drastically reduce a model's degrees of freedom and have a large impact on the determination of the most influential input parameters. These results highlight the necessity of incorporating uncertainty analysis in the mechanism generation workflow.

Gao, Connie W.↗

Developing new alkaline ceramics as possible CO 2 chemisorbents at high temperatures: The lithium and sodium yttriates (LiYO 2 and NaYO 2 ) cases

In this work, lithium and sodium yttriates (LiYO 2 and NaYO 2 ) were synthesized, characterized and tested as possible carbon dioxide (CO 2 ) captors. All the experimental syntheses and analyses were supported by the theoretical thermodynamic calculations, showing that both ceramics would be able to chemisorb CO 2 in a wide temperature range. Therefore, both ceramics were prepared by solid-state reaction and structurally characterized. In fact, the structural characterization evidenced that sodium atoms are located in a octahedral close-packed structure, while lithium ions are not so packed. The lithium ions are, indeed, in a highly distorted tetrahedron, tending to square-planar coordination. Then, based on the theoretical and structural analysis, LiYO 2 and NaYO 2 samples were investigated for the CO 2 capture process through dynamic and isothermal thermogravimetric experiments. All these experiments showed that LiYO 2 has better CO 2 capture properties than NaYO 2 , which was further correlated with the crystal structures of each ceramic. Furthermore, both ceramics presented similar CO 2 capture efficiencies regardless of the CO 2 concentration. Finally, LiYO 2 sample was tested for cyclic CO 2 chemisorption-desorption processes. The obtained results showed that the CO 2 capture efficiency is maintained high at least through 10 cycles. Therefore, these alkaline yttriates seem to present interesting alternatives as high and moderate temperature CO 2 captors, applicable for post-combustion industrial capture processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deterioration mechanism of alkali-activated materials in sulfuric acid and the influence of Cu: A micro-to-nano structural, elemental and stable isotopic multi-proxy study

In this study, a multi-proxy approach combining {sup 29}Si, {sup 27}Al and {sup 1}H MAS-NMR, FEG-EPMA, XANES at the Cu K-edge and XRD analytics with hydrochemical tools such as ICP-OES analyses, oxygen-isotope signatures, and thermodynamic modelling was applied to K-silicate-activated metakaolin specimens - with and without CuSO{sub 4}·5H{sub 2}O addition - exposed to sulfuric acid at pH = 2 for 35 days. The results revealed a multistage deterioration mechanism governed by (i) acid diffusion, (ii) leaching of K-A-S-H, (iii) microstructural damage related to precipitation of expansive (K,Ca,Al)-sulfate-hydrate phases (iv) complete dissolution of the K-A-S-H framework, (v) and formation of silica gel in the outermost corroded regions. Copper ions were mainly located in layered spertiniite-chrysocolla-like phases in the as-cured materials. The results demonstrate an overall negative effect of Cu addition on chemical material durability, implying that the reported higher durability of Cu-doped AAM in biocorrosion environments can be best explained by bacteriostatic effects.

36 MATERIALS SCIENCE↗

Semimetal–Mott insulator quantum phase transition of the Hubbard model on the honeycomb lattice

We take advantage of recent improvements in the grand canonical hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of U c /κ = 3.834(14) and the critical exponent zν = 1.185(43). Under the assumption that this corresponds to the expected antiferromagnetic Mott transition, we are also able to provide a preliminary estimate β = 1.095(37) for the critical exponent of the order parameter. We consider our findings in view of the SU(2) Gross-Neveu, or chiral Heisenberg, universality class. Here, we also discuss the computational scaling of the hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Radionuclide Vaporization from Sodium Pools for SFR Mechanistic Source Term Analysis

To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the Multiphysics Applications technical area. As part of this effort, an MST mod/sim development pathway was developed in FY21 under the NEAMS program, which outlines the high-level objectives and near-term tasks necessary to achieve the project objectives. One such recommended task was the assessment of a previously developed thermodynamic database to model radionuclide vaporization from liquid sodium pools (specifically in the context of MST analyses), which is the focus of this topic report. As such, a thermodynamic database was created in FactSage for liquid sodium pools containing oxygen and relevant elements which represent radionuclides pertinent to MST analyses. The database was demonstrated in calculating vapor fraction curves for each element in the system at concentrations relevant to SFR MST.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗

Ecological theory applied to environmental metabolomes reveals compositional divergence despite conserved molecular properties

Stream and river systems transport and process substantial amounts of dissolved organic matter (DOM) from terrestrial and aquatic sources to the ocean, with global biogeochemical implications. However, the underlying mechanisms affecting the spatiotemporal organization of DOM composition are under-investigated. To understand the principles governing DOM composition, we leverage the recently proposed synthesis of metacommunity ecology and metabolomics, termed ‘meta-metabolome ecology.’ Applying this novel approach to a freshwater ecosystem, we demonstrated that despite similar molecular properties across metabolomes, metabolite identity significantly diverged due to environmental filtering and variations in putative biochemical transformations. We refer to this phenomenon as ‘thermodynamic redundancy,’ which is analogous to the ecological concept of functional redundancy. We suggest that under thermodynamic redundancy, divergent metabolomes can support equivalent biogeochemical function just as divergent ecological communities can support equivalent ecosystem function. As these analyses are performed in additional ecosystems, potentially generalizable principles, like thermodynamic redundancy, can be revealed and provide insight into DOM dynamics.

Danczak, Robert E.↗

Evolution of intermetallic phases in an Al–Si–Ti alloy during solution treatment

A cast Al–Si–Ti alloy was solution treated at 540 °C for different periods between 0 and 72 h to understand the evolution of intermetallic phases. Only an (Al,Si) 3 Ti intermetallic phase with a low Si content was found in the as-cast alloy. The (Al,Si) 3 Ti particles were converted partly into a lamellar structure, a eutectoid phase, consisting of a Si-rich phase and an Al phase during solution treatment. The amount of the lamellar structure increased with the solution treatment time, but the composition of either constitute was kept almost unchanged regardless of the solution treatment times. The lamellar Si-rich phase is (Al,Si) 2 Ti (τ 2 ) with a TiSi 2 (C49-type) structure based on thermodynamic calculations and high-resolution TEM analyses. FCC Al phase is the product residing between the τ 2 lamellae after the completion of transformation from (Al,Si) 3 Ti to τ 2 phase. A near-rational orientation relationship (OR) between the Al and τ 2 phases is determined as Al [110]//τ 2 [100], Al ()//τ 2 [060]. The phase transformation from (Al,Si) 3 Ti to τ 2 being a result of the diffusion of Ti and Si within the original (Al,Si) 3 Ti particulates as well as the Si diffusion from the Al matrix during solution treatment is proposed. The formation of the lamellar structure in the microstructure is attributed mainly to the limited diffusivity of Ti element.

36 MATERIALS SCIENCE↗

Clumped isotope constraints on changes in latest Pleistocene hydroclimate in the northwestern Great Basin: Lake Surprise, California

Abstract During the Last Glacial Maximum (LGM) and subsequent deglaciation, the Great Basin in the southwestern United States was covered by numerous extensive closed-basin lakes, in stark contrast with the predominately arid climate observed today. This transition from lakes in the Late Pleistocene to modern aridity implies large changes in the regional water balance. Whether these changes were driven by increased precipitation rates due to changes in atmospheric dynamics, decreased evaporation rates resulting from temperature depression and summer insolation changes, or some combination of the two remains uncertain. The factors contributing to these large-scale changes in hydroclimate are critical to resolve, given that this region is poised to undergo future anthropogenic-forced climate changes with large uncertainties in model simulations for the 21st century. Furthermore, there are ambiguous constraints on the magnitude and even the sign of changes in key hydroclimate variables between the Last Glacial Maximum and the present day in both proxy reconstructions and climate model analyses of the region. Here we report thermodynamically derived estimates of changes in temperature, precipitation, and evaporation rates, as well as the isotopic composition of lake water, using clumped isotope data from an ancient lake in the northwestern Great Basin, Lake Surprise (California). Compared to modern climate, mean annual air temperature at Lake Surprise was 4.7 °C lower during the Last Glacial Maximum, with decreased evaporation rates and similar precipitation rates to modern. During the mid-deglacial period, the growth of Lake Surprise implied that the lake hydrologic budget briefly departed from steady state. Our reconstructions indicate that this growth took place rapidly, while the subsequent lake regression took place over several thousand years. Using models for precipitation and evaporation constrained from clumped isotope results, we determine that the disappearance of Lake Surprise coincided with a moderate increase in lake temperature, along with increasing evaporation rates outpacing increasing precipitation rates. Concomitant analysis of proxy data and climate model simulations for the Last Glacial Maximum are used to provide a robust means to understand past climate change, and by extension, predict how current hydroclimates may respond to expected future climate forcings. We suggest that an expansion of this analysis to more basins across a larger spatial scale could provide valuable insight into proposed climate forcings, and aid in climate model process depiction. Ultimately, our analysis highlights the importance of temperature-driven evaporation as a mechanism for lake growth and retreat in this region.

Geology↗

Studies on the Encapsulation of Brine Generated from a Process Using Selective Electrodialysis Reversal

This topical report summarizes the laboratory testing, geochemical modeling, and mineralogical testing conducted to examine the efficacy of wastewater encapsulation for the brine generated from the Saltworks Flex electrodialysis reversal technology. Within the laboratory testing, the short- and long-term physical properties of seven separate mix designs were examined along with the leaching of each mix design. Additionally, geochemical thermodynamic modeling was used along with mineralogical analyses to help gain a fundamental understanding of the chemistry involved in the encapsulation process. Encapsulation is a process that results in the immobilization of constituents of interest by creating an engineered material that contains the waste within its matrix, both physically and chemically, in such a manner that the receiving environment is protected from adverse effects. Within the context of this study, encapsulation mix designs involved mixing the Saltworks brine with fly ash and various additives, principally lime (calcium oxide) and cement, to form a paste that eventually sets and hardens into a monolith upon placement in a landfill. Moisture conditioning fly ash with wastewater sources for on-site landfill disposal is a process that is in use at some sites and has been presented as a viable option to dispose of liquid wastes. However, because chemical reactions that can lead to sequestration of constituents of interest are not optimized, leaching can occur. Therefore, in this study, mixes were created to mimic moisture conditioning, as well as engineered pastes of varying mix designs. The study found that engineered encapsulation mixes provide significant resistance to leaching compared to moisture conditioning alone. The best performing encapsulation mix used lime as the additive, which proved superior to the use of cement. After 90 days of curing, this encapsulated mix was 5.5 times stronger than the moisture conditioned mix. Additionally, this mix exhibited an approximate 2.5 order of magnitude drop in hydraulic conductivity compared to the moisture conditioned mix. When comparing leaching, nearly all the chloride was found to evolve from the moisture conditioned mix, whereas only 6.7% of the chlorine contained in the best performing encapsulated mix leached out. Mineralogical analysis via X-ray diffraction (XRD) showed a significant concentration of the chloride sequestering mineral, Friedel’s salt, in the best performing encapsulation mix, whereas none was observed in the moisture conditioned mix, pointing to it being the primary chemical mechanism for chloride retention, along with the physical impermeability of the material. With respect to the modeling study, the model was successful in predicting the presence of most mineral phases observed in the XRD analysis.

01 COAL, LIGNITE, AND PEAT↗

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanistic Insights into Acetate Selectivity on Intermetallic CuPd(110) in CO Reduction

Experimental studies demonstrate that CuPd(110) uniquely favors acetate formation during CO reduction (CORR), contrasting with the preference for ethylene on Cu surfaces. To elucidate this selectivity, we employed explicit solvation density functional theory (DFT) calculations to investigate the reaction mechanism from both thermodynamic and kinetic angles. Here, our findings reveal that on CuPd(110), the acetate-pathway intermediate H 2 CCO is thermodynamically favored at experimental conditions, while CHCHO─a precursor to ethylene─is preferred on Cu(111). Beyond thermodynamics, we find that H 2 CCO is kinetically accessible under the experimental conditions on CuPd(110), aiding acetate formation. Electron density difference analyses further corroborate distinct protonation preferences supporting this mechanism. We propose a thermodynamic screening parameter based on the Gibbs free energy, G H 2 CCO < G CHCHO , as a guide for designing Cu-based catalysts with enhanced acetate selectivity. These results offer critical mechanistic insights into the CORR product distribution and a rational framework for future catalyst design.

Acetate↗

Contrasting Effects of Urban Land Cover Change and Anthropogenic Heat on Summer Precipitation Over the Yangtze River Delta of China: Analyses From an Atmospheric Moisture Budget Perspective

Both urban land cover (ULC) change and anthropogenic heat (AH) emission are important causes of urban heat island, but their relative contributions to the changes in urban precipitation and the related mechanism remain unclear. Based on numerical simulations utilizing the latest realistic urban fraction and AH data over the Yangtze River Delta urban agglomeration, we found that ULC and AH resulted in nearly opposite effects on precipitation. Various dynamical and thermodynamic processes were involved according to the atmospheric moisture budget analyses. AH increased precipitation particularly during afternoon, and the increases were stronger during heavy precipitation events because of the enhanced moisture convergence effect together with the release of moisture storage previously accumulated in the atmosphere. Differently, ULC reduced mean precipitation mainly due to suppressed evaporation. During weak precipitation events, the suppressed evaporation was largely balanced by the intensified moisture convergence, but during heavy events, ULC caused more pronounced precipitation reduction because the moisture convergence response disappeared and failed to offset the evaporation effect. Here, the relative contributions of different dynamical and thermodynamic processes such as those related to circulation, moisture gradient, and background moisture availability to the temporal variation in the total moisture convergence were further quantified. Overall, our results help better understand the relative roles of different aspects of urbanization on precipitation, and suggest that compared to ULC, reduction in AH emission that is tightly related to the energy consumption structure could be more efficient for mitigating the risk of extreme precipitation.

54 ENVIRONMENTAL SCIENCES↗

Experimental and Theoretical Thermokinetic Analysis of the Na 3 VO 4 –CO 2 Reaction System under Different Physicochemical Conditions

Sodium vanadate (Na 3 VO 4 ) was synthesized, structural and microstructurally characterized as well as tested as possible carbon dioxide (CO 2 ) captor through thermodynamic calculations in addition to thermogravimetric dynamic and isothermal analyses. Ceramic structural characterization evidenced the formation of the Na 3 VO 4 crystal phase, while microstructural features evidenced dense agglomerates with a poor specific surface area, although synthesis temperature (600 °C) was not as high as that for other sodium ceramics. Na 3 VO 4 was investigated for the CO 2 capture process through dynamic and isothermal experiments, in the presence or absence of oxygen. All these experiments showed that Na 3 VO 4 possesses interesting CO 2 capture properties in a specific temperature range (520 and 600 °C). Isothermal products characterization allowed to elucidate the reaction pathway, implying the Na 4 V 2 O 7 formation, as possible reactive intermediate. In fact, reaction path evolution was supported with theoretical thermodynamic data. Moreover, the kinetic analysis, performed using the Avrami-Erofeev mathematical model and subsequent thermodynamic data obtention, probed the positive influence of the temperature, in both CO 2 superficial and bulk chemical sorption processes, while oxygen addition did not enhance the superficial process. Additionally, as it would be expected, kinetics were enhanced as a function of the CO 2 concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗

Alkyl-Substituted Polycaprolactone Poly(urethane-urea)s as Mechanically Competitive and Chemically Recyclable Materials

We report the mechanical performance and chemical recycling advantages of implementing alkyl-substituted poly(ε-caprolactones) (PCLs) as soft segments in thermoplastic poly(urethane-urea) (TPUU) materials. Poly(4-methylcaprolactone) (P4MCL) and poly(4-propylcaprolactone) (P4PrCL) were prepared, reacted with isophorone diisocyanate, and chain-extended with water to form TPUUs. The resulting materials’ tensile properties were similar or superior to a commercially available polyester thermoplastic poly(urethane) and had superior elastic recovery properties compared to a PCL analogue due to the noncrystalline nature of P4MCL and P4PrCL. Additionally, monomers were recovered from the TPUU materials in high yields via ring-closing depolymerization using a reactive distillation approach at an elevated temperature and a reduced pressure (240–260 °C, 25–140 mTorr) with zinc chloride (ZnCl 2 ) as the catalyst. The thermodynamics of polymerization were estimated using Van’t Hoff analyses for 4MCL and 4PrCL; these results indicated that the propyl group in 4PrCL results in a lower practical ceiling temperature (T c ) for P4PrCL.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Medium-independent hydrogen atom binding isotherms of nickel oxide electrodes

Adsorption and transfer of hydrogen atoms at solid/solution interfaces are fundamental to heterogeneous catalysis for chemical energy transformations and other processes. Reported here are electrochemical and spectroelectrochemical measurements of the thermodynamics of H-atom binding to nickel oxide electrodes, both the average and the distribution of NiO–H bond dissociation free energies (BDFEs). These are perhaps the first measurements of binding isotherms at non-metal electrodes. Remarkably, both the BDFEs and the non-Langmuirian isotherms are the same in water, acetonitrile, and dimethylformamide, and with different buffers and proton activities. Such medium independence of the BDFEs and isotherms has not been previously reported for any binary material. The medium independence supports the common use of computed hydrogen binding energies as intrinsic descriptors of surface reactivity, while the broadened isotherms add a level of complexity to such analyses. Furthermore, this work demonstrates the capability to derive key thermodynamic parameters at chemically reactive solid-liquid interfaces.

13 HYDRO ENERGY↗

Procedures for Assessing the Stability of Proposed Topological Materials

Here, we investigate the stability of MnPb 2 Bi 2 Te 6 (MPBT), which is predicted to be a magnetic topological insulator (TI), using density functional theory calculations. Our analysis includes various measures such as enthalpies of formation, Helmholtz free energies, defect formation energies, and dynamical stability. Our thermodynamic analysis shows that the phonon contribution to the energy gain from finite temperature is estimated to be less than 10 meV/atom, which may not be sufficient to stabilize MPBT at high temperatures, even with the most favorable reactions starting from binaries. While MPBT is generally robust against the formation of various defects, we find that the anti-site defect formation of MnPb is the most likely to occur, with a corresponding energy less than 60 meV. This can be attributed to the significant energy cost from compressive strain at the PbTe layer. Our findings suggest that MPBT is on the brink of stability in terms of thermodynamics and defect formation, underscoring the importance of conducting systematic analyses of the stability of proposed TIs, including MPBT, for their practical utilization. This study offers valuable insights into the design and synthesis of desirable magnetic TI materials with robust stabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗