Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An ab initio molecular dynamics study of varied compositions of the $\text{LiF-NaF-KF}$ molten salt

With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

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

A combined experimental and machine learning exploration of Ti 2-x Zr x MnCrFeNi high entropy Laves hydrides

A series of high entropy AB 2 -type Ti 2-x Zr x MnCrFeNi alloys (x = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2) were synthesized to investigate their potential for hydrogen storage and chemical compression. The influence of the Ti/Zr ratio was explored in terms of structural, microstructural and thermodynamic properties. The storage capacity together with the reaction enthalpy and entropy changes of the synthesized high entropy alloys were compared to predictions from Machine Learning (ML) to investigate changes in these properties across the explored composition space. The results revealed that a decreasing Zr content consistently lowered the hydride formation enthalpy and increased the plateau pressure from 8 to >90 bar H 2 at 25 °C, in good agreement with ML predictions. Selected compositions (x = 1.0 and 1.2) demonstrated reversible hydrogen storage capability over 150 cycles, with capacities of 1.34–1.40 wt % H 2 and remarkable reaction kinetics (<4 min) at ambient temperature. These experimental and computational findings highlight the potential of this Laves-HEA system as tuneable, stable, and cost-effective materials suitable for long-term operations in stationary hydrogen storage and compression applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

Utilizing metal cation reactions with carbonyl sulfide to remove isobaric interferences in tandem inductively coupled plasma mass spectrometry analyses

Carbonyl sulfide (OCS) was used as a reaction gas to investigate gas phase metal (M + ) ion-molecule reactions using the Agilent 8900 inductively coupled plasma tandem mass spectrometer (ICP-MS/MS) to yield insight on how this gas may be used to remove isobaric interferences in analytical measurements. The experimental work was paired with density functional theory (DFT) calculations of the reaction enthalpy to predict whether M + will react with OCS. A multi-element standard containing 46 elements ranging from 9 to 208 u was analyzed in the presence and absence of OCS. When a reaction was observed, the dominant product was the sulfide (MS + ). Oxide products were also observed for many M + but formation was less efficient with OCS than previously observed with other reaction gases. This is likely due to the weaker OC-S bond that makes MS + formation more favorable. Increasing the flow rate from 0.1 to 0.2 mL/min (corresponding to a change in reaction gas pressure from 0.35 to 0.53 Pa (2.6 to 4.0 mTorr)) generally resulted in greater MS + production, including the secondary product MS 2 + for a few cations. The early lanthanide series ions (La + , Ce + , Pr + and Nd + ) produced greater quantities of MO + at the higher pressure, although MS + products were still the dominant product. The DFT-predicted reaction enthalpies were consistent with the observed sulfide formation, with an accuracy >90%; however, model predictions were less accurate for the minor and higher order products (< 77% for MO + ). Finally, the work presented here continues a systematic study of ion-molecule reactions in ICP-MS/MS to understand and develop new and novel ways to analyze complex mixtures with minimal pre-analysis treatment.

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

Determination of the regeneration energy of direct air capture solvents/sorbents using calorimetric methods

Regeneration energy and recyclability are important parameters that dictate the viability of a CO 2 direct air capture (DAC) process. In this study, differential scanning calorimetry (DSC) is used in conjunction with thermogravimetric analysis (TGA) and Fourier Transform Infrared Spectroscopy (FTIR) to evaluate the regeneration energy required for CO 2 -loaded 3 M aqueous potassium sarcosinate (K-SAR) solvent and crystalline methyl-glyoxal-bisiminoguanidine (MGBIG) sorbent. Based on calorimetric measurements, for aqueous K-SAR, the sensible heat amounts to 1.5 GJ/tCO 2 and the enthalpy of desorption is estimated at 3.68 GJ/tCO 2 . The total energy of regeneration will also include the enthalpy of solvent vaporization, which will depend on the thermodynamic conditions and the level of regeneration. For MGBIG, the total regeneration energy required was measured at approximately 7.0 GJ/tCO 2 . Further, FTIR measurements were used to observe CO 2 and H 2 O release during the thermal regeneration process. For K-SAR, CO 2 release is seen to take place at temperatures greater than 80 °C and, for MGBIG, the CO 2 release takes place at temperatures greater than 100 °C. The moderate temperatures required for K-SAR and MGBIG regenerations compared to other DAC solvents/sorbents suggest that low-cost sources of heat, such as geothermal heat, could potentially be used for regeneration to drive down the cost of direct air capture.

42 ENGINEERING↗
DOE OSTI2021

Importance of Engineered and Learned Molecular Representations in Predicting Organic Reactivity, Selectivity, and Chemical Properties

Machine-readable chemical structure representations are foundational in all attempts to harness machine learning for the prediction of reactivities, selectivities, and chemical properties directly from molecular structure. The featurization of discrete chemical structures into a continuous vector space is a critical phase undertaken before model selection, and the development of new ways to quantitatively encode molecules is an active area of research. Here, we highlight the application and suitability of different representations, from expert-guided “engineered” descriptors to automatically “learned” features, in different prediction tasks relevant to organic and organometallic chemistry, where differing amounts of training data are available. These tasks include statistical models of stereo- and enantioselectivity, thermochemistry, and kinetics developed using experimental and quantum chemical data. The use of expert-guided molecular descriptors provides an opportunity to incorporate chemical knowledge, domain expertise, and physical constraints into statistical modeling. In applications to stereoselective organic and organometallic catalysis, where data sets may be relatively small and 3D-geometries and conformations play an important role, mechanistically informed features can be used successfully to obtain predictive statistical models that are also chemically interpretable. We provide an overview of several recent applications of this approach to obtain quantitative models for reactivity and selectivity, where topological descriptors, quantum mechanical calculations of electronic and steric properties, along with conformational ensembles, all feature as essential ingredients of the molecular representations used. Alternatively, more flexible, general-purpose molecular representations such as attributed molecular graphs can be used with machine learning approaches to learn the complex relationship between a structure and prediction target. This approach has the potential to out-perform more traditional representation methods such as “hand-crafted” molecular descriptors, particularly as data set sizes grow. One area where this is particularly relevant is in the use of large sets of quantum mechanical data to train quantitative structure–property relationships. A general approach toward curating useful data sets and training highly accurate graph neural network models is discussed in the context of organic bond dissociation enthalpies, where this strategy outperforms regression using precomputed descriptors. Finally, we describe how graph neural network predictions can be incorporated into mechanistically informed statistical models of chemical reactivity and selectivity. Once trained, this approach avoids the expensive computational overhead associated with quantum mechanical calculations, while maintaining chemical interpretability. We illustrate examples for which fast predictions of bond dissociation enthalpy and of the identities of radicals formed through cleavage of a molecule’s weakest bond are used in simple physical models of site-selectivity and reactivity.

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

Exploring Ca–Ce–M–O (M = 3d Transition Metal) Oxide Perovskites for Solar Thermochemical Applications

Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.

14 SOLAR ENERGY↗
DOE OSTI2020

Favorable Redox Thermodynamics of SrTi 0.5 Mn 0.5 O 3–δ in Solar Thermochemical Water Splitting

Two-step, solar thermochemical splitting of water using nonstoichiometric redox-active metal oxides has emerged as an intriguing approach for large-scale hydrogen production. Perovskites have been proposed as alternatives to state-of-the-art fluorite CeO 2–δ because of their potential for lowering reduction temperature while maintaining high fuel productivity. Guided by computational insights, we explore the thermodynamic properties and water splitting efficacy of the cubic perovskite SrTi 0.5 Mn 0.5 O 3–δ (STM55). Thermogravimetric analysis is performed under controlled oxygen partial pressures (pO 2 ) and temperatures up to 1500 °C, from which both the enthalpy and entropy of reduction as a function of oxygen nonstoichiometry are determined. STM55 provides an attractive combination of moderate enthalpy, 200–250 kJ (mol-O) –1 , and high entropy, with unusual δ dependence. Using a water splitting cycle in which the material is thermally reduced at 1350 °C (pO 2 , ~10 –5 atm) and subsequently exposed to steam at 1100 °C (steam partial pressure of pH 2 O = 0.4 atm), we demonstrate a hydrogen yield of 7.4 mL g –1 . Through both half-cycles, the material remains largely in quasi-equilibrium with the gas phase, as reflected in the agreement of the measured data with predicted profiles based on the thermodynamic data. This behavior indicates rapid surface and bulk diffusion kinetics. Cyclic operation showed the material to be free of degradation and always resulted in a 2:1 yield of H 2 /O 2 . Altogether, STM55 provides outstanding performance characteristics for thermochemical hydrogen production.

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

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Vibrational Signatures of Unrealized Phosphorus Suboxide Intermediates in White Phosphorus Oxidation Reactions

White phosphorus ignition notoriously produces the phosphoric acid anhydride P 4 O 10 , yet the intermediate oxidation steps remain undetermined. Here, we report the first geometric and vibrational characterization of two P 4 O 2 isomers, P 3 OPO and P 3 PO 2 , and substantiate a previously proposed cyclic P 4 O 2 isomer. We formally assign the infrared bands observed by Andrews and Mielke at 898 and 891 cm −1 to the antisymmetric P− O−P vibrations of P 3 OPO species. Additional bands corresponding to terminal P=O and −PO 2 stretches of P 3 OPO and P 3 PO 2 discussed herein also went unrecognized due to the peculiar bonding of P 4 O 2 species compared to oxo-bridged P 4 O x (x = 3−6) species. Sequential addition of oxygen atoms to the P 4 tetrahedron appears to form P 3 OPO and P 3 PO 2 , while cyclic P 4 O 2 is formed from P 2 O dimerization. CCSD(T) geometries, CCSD(T) + MP2[δVPT2] fundamental frequencies, and enthalpies of formation extrapolated using focal-point analysis are reported. The predicted enthalpies of formation relative to tetrahedral P 4 plus 3 Σ g − O 2 for bent-P 3 OPO, extended-P 3 OPO, P 3 PO 2 , and cyclic-P 4 O 2 are −93, −88, −85, and −63 kcal mol −1 , respectively.

Molecular structure↗
DOE OSTI2023

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

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

Energetics of Ag Adsorption on and Adhesion to Rutile TiO 2 (100) Studied by Microcalorimetry

The adsorption and adhesion energies of vapor-deposited Ag on rutile TiO 2 (100) films have been measured using single crystal adsorption calorimetry (SCAC) and He + low energy ion scattering spectroscopy (LEIS). Ag grows as three-dimensional (3D) nanoparticles at 300 and 100 K. The saturation particle density is 8 × 10 16 particles/m 2 at 300 K and 2.5 × 10 17 particles/m 2 at 100 K. The differential heat of adsorption starts low, increases with Ag coverage, and finally approaches the sublimation enthalpy of Ag at both 300 and 100 K. At 300 K, the differential heat of adsorption starts from 208 kJ/mol and rises rapidly to 265 kJ/mol by 1 ML. At 100 K, Ag grows nanoparticles with smaller particle size on the terraces, while at 300 K, the Ag nanoparticles have bigger particle size and grow on step edges with stronger binding strength to the nanoparticles. Thus, the heat of Ag adsorption at 100 K starts at 141 kJ/mol (near the monomer size limit on terrace sites) and remains lower than that at 300 K until 1.5 ML. Here, the adhesion energy of Ag(solid) to rutile TiO 2 (100) is found to be large (~2.44 J/m 2 ), supporting a trend of decreasing adhesion energy with the enthalpy of oxide reduction.

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

High-Pressure Investigation of 2,4,6-Trinitro-3-bromoanisole (TNBA): Structural Determination and Piezochromism

Understanding phase transitions in energetic materials is crucial for developing predictive models of detonation. 2,4,6-Trinitro-3-bromoanisole (TNBA), an energetic material, was studied in its single-crystal form up to pressures of 45 GPa in a diamond anvil cell. The material was characterized by using X-ray, Raman, and optical transmission measurements. From single-crystal X-ray diffraction, the ambient structure of TNBA was determined which crystallizes in the P2 1 /c space group having four molecular units per unit cell. The X-ray data up to 9.2 GPa were fitted to a third-order Birch–Murnaghan equation of state by using the parameters K o = 13.2(2.4) GPa and K p = 5.1(1.4). Between 6.8 and 7.3 GPa, a phase transition was inferred in TNBA from concurrent fading of X-ray diffraction, disappearance of Raman peaks, increase in sample fluorescence, and discontinuous color change. The new phase was consistent with an amorphous state of at least partially intact molecules judging from the presence of higher-order Raman modes and irreversibility of the Raman spectra upon release. Piezochromism was observed with the translucent yellow TNBA gradually darkening and becoming opaque black at ~25 GPa. This correlated to the absorption edge gradually shifting to the red in the visible spectrum. Signs of two possible additional structural transitions were detected in the 32.4–41.0 GPa range as suggested by a jump in the absorption edge, the irreversible changes in the absorption spectrum upon release to ambient pressure, and by the lack of Raman modes in recovered samples. Here, the crystal and electronic structures of TNBA were also investigated up to 10 GPa by using DFT calculations and crystal structure prediction (CSP) simulations. In agreement with the experimentally observed transition at 7 GPa, the simulations at 10 GPa found a bevy of polymorphs lower in enthalpy and higher in density than P2 1 /c. The lowest calculated enthalpy structure was determined to be P2 1 2 1 2 1 , being in a different space group than the ambient experimental result.

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

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface Equilibration Mechanism Controls the Stability of a Model Codeposited Glass Mixture of Organic Semiconductors

While previous work has identified the conditions for preparing ultrastable single-component organic glasses by physical vapor deposition (PVD), little is known about the stability of codeposited mixtures. Here, we prepared binary PVD glasses of organic semiconductors, TPD (N,N'-Bis(3-methylphenyl)-N,N'-diphenylbenzidine) and m-MTDATA (4,4',4"-Tris[phenyl(m-tolyl)amino]triphenylamine), with a 50:50 mass concentration over a wide range of substrate temperatures (T sub ). The enthalpy and kinetic stability are evaluated with differential scanning calorimetry and spectroscopic ellipsometry. Binary organic semiconductor glasses with exceptional thermodynamic and kinetic stability comparable to the most stable single-component organic glasses are obtained when deposited at T sub = 0.78–0.90T g (where T g is the conventional glass transition temperature). When deposited at 0.94T g , the enthalpy of the m-MTDATA/TPD glass equals that expected for the equilibrium liquid at that temperature. Thus, the surface equilibration mechanism previously advanced for single-component PVD glasses is also applicable for these codeposited glasses. Furthermore, these results provide an avenue for designing high-performance organic electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Charge Separation at Organic Interfaces with Near-Zero Energy Offset: A Step toward Designing Inorganic-like Organic Semiconductors

For organic semiconductors, it is often presumed that an energy level offset at the donor–acceptor (D-A) interface is required to provide the driving force for charge separation (CS). This energy level offset unavoidably leads to a voltage loss in organic photovoltaics. In this work, by using zinc phthalocyanine (ZnPc) and fluorinated zinc phthalocyanine (F 4 ZnPc) as a model D/A interface, we found that spontaneous CS, with an enthalpy increase of ∼0.3–0.4 eV, can occur even with an interfacial energy offset as small as ∼0.1 eV. This enthalpy-increase CS process is driven by entropy. The entropic driving force can be enhanced by two factors: (1) a point-like spatial contact between the delocalized electron and hole wave function in the charge transfer exciton; (2) a small band bending near the interface originating from long-range electrostatic interaction. Our work demonstrates that effective CS can occur at interfaces with a near-zero energy level offset, which means that the energy loss at the D/A interface can be avoided.

36 MATERIALS SCIENCE↗