SEARCH · Engineering Papers
Results for “Kinetic phase transition”
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.
High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions
Abstract Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12// ω B97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol −1 ) relative to those calculated at ω B97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients $${k}_{\infty }(T)$$ k ∞ ( T ) between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.
Enabling Structure/Interface Regulation for High Performance Ni-Rich Cathodes
Here, further commercialization of Ni-rich layered cathodes is hindered by severe structure/interface degradation and kinetic hindrance that occur during electrochemical operation, which leads to safety risks and reduced range in electric vehicles (EVs). Herein, by selecting elements with different solubility properties, a multifunctional strategy that synchronously fabricates perovskite-type SrZrO 3 coating and Sr/Zr co-doping is employed to strengthen the structure/interface stability and the Li + transport mobility of LiNi 0.85 Co 0.10 Mn 0.05 O 2 (NCM). Perovskite-type SrZrO 3 protective layers formed on the particle surface can substantially mitigate the unexpected interfacial side reactions and surface phase transitions. In addition, a robust crystal framework is constructed by optimizing local O coordination through the introduction of strong Zr-O bonds. Notably, Li + diffusion kinetics is effectively improved due to expanded cell parameters and O-Li-O slab spacing with the incorporation of large-diameter Sr pillar ions, as revealed by X-ray diffraction. As a result, the Sr/Zr-modified NCM achieves a remarkable capacity retention of 99.4% after 200 cycles at 1 C, and a high rate capacity of 168.9 mAh g -1 at 10 C. This work opens new avenues to develop high-performance NCM cathodes with high energy and high power for EVs with long calendar life.
Statistically averaged molecular dynamics simulations of hydrogen diffusion in magnesium and magnesium hydrides
Magnesium has a different crystal structure from its dihydride with hydrogenation leading to a phase transition from the hexagonal closely packed Mg into a tetragonal α-MgH 2 rutile type structure. Such materials exhibit complex hydrogen uptake and release kinetics because hydrogen diffusivities significantly change when the crystal structure changes. To provide a foundational understanding of (de)hydrogenation kinetics that is applicable to all stages of the reaction, we performed statistically averaged molecular dynamics simulations to derive hydrogen diffusivities as a function of temperature and hydrogen content for both magnesium and magnesium hydride. Furthermore, our studies confirm that hydrogen diffusivities in magnesium hydride are much lower than in magnesium, in agreement with experimental data. Additionally, we observe that in either magnesium or magnesium hydride, higher hydrogen compositions result in reduced diffusivities. The latter was not revealed by prior experiments, which were conducted at fixed hydrogen composition. Finally, we discover a non-Arrhenius behavior in magnesium hydride. The physical origin of this behavior is also discussed.
Solubility-limited depolymerization kinetics in the glycolysis of carbonyl-containing polymers
Chemical recycling of condensation polymers is often rationalized on the basis of the intrinsic reactivity of ester and carbonate functional groups. However, under heterogeneous conditions relevant to plastic waste processing and environmental degradation, bulk depolymerization rates often diverge from trends predicted by homogeneous chemistry. Here, we investigate how polymer–solvent compatibility, catalyst strength, and phase behavior govern the heterogeneous glycolysis of carbonyl-containing polymers. Using poly(ethylene terephthalate) (PET), glycol-modified PET (PETG), and bisphenol-A polycarbonate (PC) as model systems, we examine depolymerization kinetics at 180 °C with ethylene glycol and bisphenol A as diols under both amphoteric organosalt (TBD : MSA) and strong base (TBD) catalysis. Despite substantial differences in crystallinity and glycol uptake, PET and PETG depolymerize at comparable rates under organosalt catalysis, while PC depolymerizes significantly more slowly under identical conditions. Time-resolved molecular weight analysis and thermal characterization demonstrate that these rate differences do not arise from crystallinity, swelling, or inherent carbonyl reactivity, but instead reflect solubility-limited kinetics that constrain the transition from heterogeneous to homogeneous reaction regimes. When polymer solubility is low, depolymerization remains heterogeneous and slow; when solubility is enhanced—either through increased polymer–diol compatibility or stronger base catalysis—rapid homogeneous depolymerization is observed, reversing apparent reactivity trends. These results establish solubility and phase behavior as primary determinants of depolymerization kinetics in heterogeneous polymer recycling systems. By demonstrating how catalyst selection and solvent compatibility can expose or overcome solubility limitations, this work provides mechanistic insight to design more energy-efficient and selective chemical recycling processes. More broadly, these findings suggest that polymers with limited solvent or water compatibility may resist chemical degradation in the environment, favoring fragmentation and persistence as micro- and nanoplastics. Understanding solubility-controlled depolymerization offers a pathway toward more sustainable polymer design and end-of-life chemical recovery.
Direct Study of Changes in Catalyst Structure-Kinetic Properties During Redox Transitions
The Temporal Analysis of Products (TAP) pulse response methodology is a transient technique that provides the time resolution needed to deconvolve reaction steps from the complex networks typical in industrial catalytic processes. Traditionally, TAP measurements observe gas phase dynamics at the reactor exit but lack direct measurements of changes in the catalyst itself. Recently, a new operando technique was developed that couples gas phase transients to dynamic changes in metal centers with the precise TAP methodology for nanomole titration. Using an industrial CrOx/Al2O3 catalyst used for propane dehydrogenation, we demonstrate the capabilities of this unique device to reveal key catalytic processes: 1) total propane oxidation not associated with chromia centers, 2) reduction of Cr6+ to Cr3+ correlated with selective product formation, and 3) subsequent carbon accumulation. By utilizing incremental pulsing in a diffusion-only transport regime, the spectrokinetic device allows us to resolve detailed changes in catalyst structure, composition, and kinetic function that are otherwise indistinguishable in conventional operando devices. The unification of the TAP methodology with time-resolved spectroscopic measurements offers new and unique insights into the complex kinetic phenomena regulated by solid catalyst surfaces.
Phase-field modeling of alloy oxidation at high temperatures
Oxide growth is a complex process involving transport of reactive species, heterogeneous reactions, and microstructure evolution. Predicting oxidation kinetics and especially the oxide morphological change has been a long-standing challenge. Here we develop a phase-field model for predicting the oxide growth kinetics of a multicomponent alloy during high temperature oxidation, focusing on internal oxidation (non-protective) and its transition to external oxidation (protective). In this work, the predicted kinetics and oxide morphology are analyzed and compared to the classical Wagner’s theory and an existing analytical model by Zhao and Gleeson. Some assumptions used in the analytical models and the limitation are discussed. In addition, it is demonstrated that the morphology and distribution of the initial oxide nuclei play an important role in the later stage oxide connectivity and thus the transition to external oxidation.
A limit on the total lepton number in the Universe from BBN and the CMB
At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.
The mechanism driving a solid–solid phase transition in a biomacromolecular crystal
A solid-solid phase transition (SSPT) occurs between distinguishable crystalline forms. Because of its importance in application and theory in material science and condensed matter physics, SSPT has been studied most extensively in metallic alloys, inorganic salt or small organic molecular crystals, but much less so in biomacromolecular crystals. In general, the mechanism of SSPT at atomic and molecular levels is not well understood. Here, we describe the ordered molecular rearrangements in biomacromolecular crystals of the adenine riboswitch (riboA) aptamer using real-time serial crystallography and solution atomic force microscopy (AFM). The large, ligand-induced conformational changes drive the initial phase transition from the apo unit cell (AUC) to the trans unit cell 1 (TUC1). During this transition, coaxial stacking of P1 duplexes becomes the dominant packing interface, whereas P2-P2 interactions are almost completely disrupted, resulting in “floating” layers of molecules. The coupling points in TUC1 and their local conformational flexibility allow the molecules to reorganize to achieve the more densely packed and energetically favorable bound unit cell (BUC). Our study thus reveals the interplay between the conformational changes and the crystal phases—the underlying mechanism that drives the phase transition. Using polarized video microscopy (PVM) to monitor the SSPT in small crystals at high ligand concentration, we have identified the time window during which the major conformational changes take place, and simulated the in crystallo kinetics. Together, these results provide the spatiotemporal information necessary for informing time-resolved crystallography (TRX) experiments. Moreover, this study illustrates a practical approach to characterize SSPT in transparent crystals.
Beyond Volume Variation: Anisotropic and Protrusive Lithiation in Bismuth Nanowire
Materials storing energy via an alloying reaction are promising anode candidates in rechargeable lithium-ion batteries (LIBs) due to their much higher energy density than the current graphite anode. Until now, the volumetric expansion of such electrode particles during lithiation has been considered as solely responsible for cycling-induced structural failure. Here, we report different structural failure mechanisms using single-crystalline bismuth nanowires as the alloying-based anode. The Li–Bi alloying process exhibits a two-step transition, that is, Bi–Li 1 Bi and Li 1 Bi–Li 3 Bi. Interestingly, the Bi–Li 1 Bi phase transition occurs not only in the bulk Bi nanowire but also on the particle surface showing its characteristic behavior. The bulk alloying kinetics favors a Bi-(012)-facilitated anisotropic lithiation, whose mechanism and energetics are further studied using the density functional theory calculations. More importantly, the protrusion of Li 1 Bi nanograins as a result of anisotropic Li–Bi alloying is found to dominate the surface morphology of Bi particles. The growth kinetics of Li 1 Bi protrusions is understood atomically with the identification of two different controlling mechanisms, that is, the dislocation-assisted strain relaxation at the Bi/Li 1 Bi interface and the short-range migration of Bi supporting the off-Bi growth of Li 1 Bi. As loosely rooted to the bulk substrate and easily peeled off and detached into the electrolyte, these nanoscale protrusions developed during battery cycling are believed to be an important factor responsible for the capacity decay of such alloying-based anodes at the electrode level.
Is the Orbital-Selective Mott Phase Stable against Interorbital Hopping?
The localization-delocalization transition is at the heart of strong correlation physics. Recently, there is great interest in multiorbital systems where this transition can be restricted to certain orbitals, leading to an orbital-selective Mott phase (OSMP). Theoretically, the OSMP is widely studied for kinetically decoupled orbitals, but the effect of interorbital hopping remains unclear. Here, we show how nonlocal interorbital hopping leads to local hybridization in single-site dynamical mean-field theory (DMFT). Under fairly general circumstances, this implies that, at zero temperature, the OSMP, involving the Mott-insulating state of one orbital, is unstable against interorbital hopping to a different, metallic orbital. We further show that the coherence scale below which all electrons are itinerant is very small and gets exponentially suppressed even if the interorbital hopping is not overly small. Within this framework, the OSMP with interorbital hopping may thus reach down to extremely low temperatures T, but not to T = 0. Accordingly, it is part of a coherence-incoherence crossover and not a quantum critical point. Finally, we present analytical arguments supported by numerical results using the numerical renormalization group as a DMFT impurity solver. We also compare our findings with previous slave-spin studies.
Phase transformation mechanism in irradiation-induced superlattice formation
Atomic kinetic Monte Carlo simulations were used to model void superlattice formation under irradiation in molybdenum, driven by anisotropic diffusion of self-interstitial atoms. A change in the phase transformation mechanism from nucleation and growth to spinodal decomposition occurred with increasing dose rate, with both mechanisms leading to superlattice formation. Analysis of a rate-theory based analytical model showed that an observed change in the kinetics of vacancy accumulation, the appearance of a region of positive second derivative in the plot of average vacancy concentration versus time, was caused by the onset of spinodal instability. Further, the analytical model showed that for molybdenum and several other metals where void superlattice formation is commonly observed, the phase transformation likely occurs by nucleation and growth. However, nickel may offer the possibility of experimental observation of the transition between phase transformation mechanisms.
Modulating Operational Conditions to Mitigate Deactivation in Formate Dehydrogenation on Pd Phases
In heterogeneous catalysis, poisoning by surface-bound intermediates poses a major barrier to sustained catalyst performance in (de)hydrogenation reactions. Formate/bicarbonate systems, as liquid organic hydrogen carriers (LOHCs), offer a CO 2 -integrated, low-temperature pathway for hydrogen storage and release, making them attractive for circular energy applications. However, their lower hydrogen density and susceptibility to catalyst deactivation limit their competitiveness compared to conventional LOHCs like methylcyclohexane. Here, this study investigates the mechanistic origins of formate (HCOO – ) dehydrogenation and associated deactivation on Pd interfaces. Using density functional theory (DFT) simulations, we show that under thermocatalytic conditions, strongly bound formate accumulates on the catalyst surface (Pd(111)), blocking active sites, raising activation barriers, and leading to progressive performance loss. Because formate adsorption involves charge transfer, we exploit its sensitivity to electronic structure by modulating the electrochemical potential of the catalyst. Our results reveal that hydrogen transfer from water and formate exhibits opposing potential dependencies, providing insights into the opposing driving forces behind both catalytic activity and poisoning. To further probe this phenomenon, we examine electrochemically induced phase transitions in Pd, focusing on PdO(100) and PdH(110), which are stable under oxidizing and reducing potentials, respectively, and demonstrate enhanced dehydrogenation activity between −0.4 and 0.2 V vs standard hydrogen electrode (SHE). Complementary thermal treatments help decouple kinetic and thermodynamic contributions to intermediate binding. These findings underscore the critical role of the catalyst phase and external stimuli in dictating poison-active site interactions and highlight phase engineering as a promising strategy to mitigate deactivation. This work offers mechanistic insights and design principles for developing more resilient and efficient catalysts for LOHC applications under realistic operating conditions.
Dual phase patterning during a congruent grain boundary phase transition in elemental copper
The phase behavior of grain boundaries can have a strong influence on interfacial properties. Little is known about the emergence of grain boundary phases in elemental metal systems and how they transform. Here, we observe the nanoscale patterning of a grain boundary by two alternating grain boundary phases with distinct atomic structures in elemental copper by atomic resolution imaging. The same grain boundary phases are found by computational grain boundary structure search indicating a first-order transformation. Finite temperature atomistic simulations reveal a congruent, diffusionless transition between these phases under ambient pressure. The patterning of the grain boundary at room temperature is dominated by the grain boundary phase junctions separating the phase segments. Our analysis suggests that the reduced mobility of the phase junctions at low temperatures kinetically limits the transformation, but repulsive elastic interactions between them and disconnections could additionally stabilize the pattern formation.
Phase transformation kinetics model for metals
We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.
Enhancing Biopolyester Backbone Rigidity with an Asymmetric Furanic Monomer
Biobased furanic polyesters can exhibit performance advantages over petroleum-derived polyesters, primarily due to their rigid furan-containing backbones. Herein, we develop two strategies to polymerize methyl 5-hydroxymethyl furanoate to poly(5-hydroxymethyl furanoate) (PHMF), a furan-based polyester with even greater backbone rigidity than poly(ethylene furanoate). Thermal, spectroscopic, and computational investigations of PHMF alongside analogous furan-based and phenyl-based polyesters suggest that the high furan content of PHMF leads to its high glass transition temperature, slow crystallization kinetics, and low amorphous mobility. Molecular dynamics simulations suggest that while the backbone of PHMF is exceptionally rigid, its amorphous phase is denser than its phenyl analog due to noncovalent interchain interactions. Together, these results highlight how asymmetric furan-based monomers can modulate key properties in biobased polyesters.
In situ formed partially disordered phases as earth-abundant Mn-rich cathode materials
Abstract Earth-abundant cathode materials are urgently needed to enable scaling of the Li-ion industry to multiply terawatt hours of annual production, necessitating reconsideration of how good cathode materials can be obtained. Irreversible transition metal migration and phase transformations in Li-ion cathodes are typically believed to be detrimental because they may trigger voltage hysteresis, poor kinetics and capacity degradation. Here we challenge this conventional consensus by reporting an unusual phase transformation from disordered Li- and Mn-rich rock salts to a new phase (named δ ), which displays partial spinel-like ordering with short coherence length and exhibits high energy density and rate capability. Unlike other Mn-based cathodes, the δ phase exhibits almost no voltage fade upon cycling. We identify the driving force and kinetics of this in situ cathode formation and establish design guidelines for Li- and Mn-rich compositions that combine high energy density, high rate capability and good cyclability, thereby enabling Mn-based energy storage.
Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon
The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.