Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stability and 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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fluorine-Substituted Lithium Chloride Solid Electrolytes for High-Voltage All-Solid-State Lithium-Ion Batteries

Lithium ternary halides are promising solid electrolytes, owing to their high ionic conductivity and reasonably high oxidative and chemical stability. Recently, fluorine substitution in Li 3 MCl 6 has been suggested as a promising approach for further enhancing oxidation stability. Accordingly, this study outlines a material design strategy for F-substituted Li 3 MCl 6 through systematic theoretical analyses. Calculations reveal that the mixing limit of F in Li 3 MCl 6–x F x is in the range of 0.5–1.5, and the resulting Li 3 MCl 6–x F x phases can retain ionic conductivity above 1 mS/cm up to x = 1.0. Additionally, the calculations also predict that the formation of F-containing passivating phases could increase the oxidation potential for Li3MCl5F to ~6.3 V. The proposed material design strategy is validated through the synthesis of Li 3 YCl 5 F, which is confirmed to show both high ionic conductivity and enhanced oxidation stability. The design guidelines presented herein can accelerate the potential use of halide-based electrolyte chemistries in high-voltage all-solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Turbulence Structures in the Very Stable Boundary Layer Under the Influence of Wind Profile Distortion

Abstract In very stable boundary layers (VSBL), a “cocktail” of submeso motions routinely result in elevated mean wind speed maxima above the ground, acting as a new source of turbulence generation. This new source of turbulent kinetic energy enhances turbulent mixing and causes mean wind profile distortion (WPD). As a results, this transient distortion in the wind profile adjusts the classical log‐law. Addressing how WPD‐induced turbulence regulates flow structures, turbulent fluxes, and transitions in stability regimes across layers remains a challenge. Eddy covariance data measured at four levels on a 62‐m tower are employed to address these questions. It is shown that the WPD initiates large turbulent eddies that penetrate downward, leading to enhanced vertical mixing and comparable turbulent transport efficiencies across layers. As a consequence, turbulence intensity and fluxes are increased. As the WPD is intensified, turbulent fluxes and turbulent flux transport caused by large eddies are also enhanced, leading to a transition from very stable to weakly stable regimes. Due to the influence of WPD‐induced large eddies, the large‐eddy turbulent Prandtl number does not deviate appreciably from unity and the partitioning between turbulent kinetic and potential energies is linearly related to the gradient Richardson number.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Deciphering Chemical/Electrochemical Compatibility of Li 3 InCl 6 in 5.2 V High-Voltage LiCoO 2 All-Solid-State Batteries

Large interfacial resistance is a widely recognized impediment to the advancement of high-voltage, all-solid-state batteries. However, a comprehensive understanding of the fundamental cause behind the interfacial resistance between solid electrolytes and typical layered oxide cathodes has not yet been achieved. Here, we investigated the high-voltage stability of Li 3 InCl 6 and elucidated the underlying interfacial electrochemical reactions between LiCoO 2 and Li 3 InCl 6 . Further, the pairing of Li 3 InCl 6 with LiCoO 2 exhibited a superior capacity retention of 73.6% even at 5.2 V, much higher than 28.2% charged at 4.6 V in lithium-ion batteries after 70 cycles. The enhanced high-voltage stability of ASSBs is attributed to the stable interface formed between LiCoO 2 and Li 3 InCl 6 and the reinforced surface and bulk structure stability. On the other hand, the ultrahigh voltage still causes the partial decomposition of Li 3 InCl 6 and generates interfacial compounds such as InClO and cobalt and indium chlorides/oxides.

25 ENERGY STORAGE↗
DOE OSTI2022

Consequences of Intrapore Liquids on Reactivity, Selectivity, and Stability for Aldol Condensation Reactions on Anatase TiO 2 Catalysts

Abstract This study provides evidence and mechanistic interpretations for the significant consequences of intrapore non‐polar liquids on acetone aldol condensation turnover rates, selectivity to primary dimer products, and catalyst stability for reactions at Lewis acid‐base site pairs on TiO 2 surfaces. These non‐polar liquids confer such benefits through the preferential stabilization of transition states (TS) for adsorption (“entry”) and desorption (“exit”) steps, which place their respective reactants or products within a solvating outer sphere environment. The extent to which non‐polar fluids ( n ‐heptane) form an intrapore liquid phase within TiO 2 voids was obtained from N 2 uptakes using established formalisms that consider the different molal volume, surface tension, and volatility between N 2 and n ‐heptane. Acetone condensation rates are limited by C−H activation, an “entry” step that forms bound prop‐1‐en‐2‐olates via a TS stabilized by intrapore liquids, leading to higher aldol condensation turnover rates as n ‐heptane pressure increases and active TiO 2 surfaces become increasingly immersed within a non‐polar liquid phase. These liquids solvate the late TS structures that mediate the desorption of primary C 6 condensation products even more effectively than those involved in prop‐1‐en‐2‐olate formation or in nucleophilic attack events that later form C−C bonds. Such preferential solvation favors desorption over C−C coupling events, thus disfavoring the formation of larger oligomers that become stranded at active sites, thus leading to much slower deactivation. Moreover, solvation by non‐polar liquids also leads to C 6 alkanones as the sole products formed in a single surface sojourn. These effects of a non‐polar dense phase circumvent the inherent stability, reactivity, and selectivity hurdles that have precluded practical aldol condensation catalysis on Lewis acid‐base pairs at oxide surfaces; these consequences are demonstrated here for TiO 2 catalysts, acetone aldol condensation reactions, and n ‐heptane as the non‐polar liquid but through strategies, concepts, and mechanistic features that extend to other systems. More generally, these observations and their mechanistic origins demonstrate how a contacting liquid preferentially solvates TS structures for elementary steps that involve either reactants arriving from or products entering into an outer sphere environment that contains a dense non‐polar phase.

Kadam, Shashikant A.↗
DOE OSTI2022

Influence of Long-Range Forces on the Transition States and Dynamics of NaCl Ion-Pair Dissociation in Water

We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long-range Coulomb interactions and for a well-chosen reference system with short-range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool “State predictive information bottleneck” (SPIB), we show that the system with short-range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long-range interactions can be largely neglected for these processes because long-range forces from the direct interionic Coulomb interactions are almost completely canceled (>90%) by those from solvent interactions over the length scale where the transition takes place. Thus, for this important monovalent ion-pair system, short-range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short-range models for the dissociation event. Here, we believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long-range contributions may be more important.

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

Understanding Instability in Formamidinium Lead Halide Perovskites: Kinetics of Transformative Reactions at Grain and Subgrain Boundaries

Transformative and reconstructive reactions impart significant structural changes at particle boundaries of hybrid perovskites, which influence environmental stability and optoelectronic properties of these materials. Here, we investigate the moisture-induced transformative reactions in formamidinium based perovskites FAPbX 3 (X=I, Br) and show that the ambient stability of these materials can be adjusted from a few hours to several months. For FAPbI 3 , roles of water vapor, particle size, and light illumination on the kinetic pathways of the cubic (a) transformation to the hexagonal (d) phase are analyzed by X-ray diffraction, optical microscopy, photoluminescence and solid-state NMR spectroscopy techniques. The grain and sub-grain boundaries exhibit different α→δ-FAPbI 3 phase transformation kinetics. Our study suggests that the dynamic transformation involves the local water-induced dissolution of the cubic phase occurring at the crystal surfaces followed by precipitation of the hexagonal phase. Insights into structures and dynamics of a kinetically trapped α-|δ-FAPbI 3 are obtained by 1 H, 2 H, and 207 Pb ssNMR spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improving Electronic Conductivity of Layered Oxides through the Formation of Two-Dimensional Heterointerface for Intercalation Batteries

Synthetic strategies for the improvement in electronic conductivities and electrochemical stabilities of transition metal oxide cathodes, which are limiting factors in the performance of commercial intercalation batteries, are required for next-generation, high-performance battery systems. Furthermore, the chemical preintercalation approach, consisting of a combined sequence of a sol–gel process, extended aging, and a hydrothermal treatment, is a versatile, wet synthesis technique that allows for the incorporation of a polar species between the layers of transition metal oxides. Here, formation of a layered 2D δ-C x V 2 O 5 · n H 2 O heterostructure occurs via chemical preintercalation of dopamine molecules between bilayers of vanadium oxide followed by the hydrothermal treatment of the precipitate, leading to carbonization of the organic molecules. The presence of carbon layers within the structure has been confirmed via a combined analysis of scanning electron microscopy, X-ray diffraction, thermogravimetric analysis, Raman spectroscopy, X-ray photoelectron spectroscopy, electrochemical impedance spectroscopy, four-probe conductivity measurements, and scanning transmission electron microscopy characterization. 2D δ-C x V 2 O 5 · n H 2 O heterostructure electrodes demonstrated significantly improved electrochemical performance, particularly at higher current densities, in Li-ion cells. The heterostructure electrodes exhibited 75% of the capacity retention when the current was changed from 20 mA g –1 (206 mAh g –1 ) to 300 mA g –1 (155 mAh g –1 ), while the reference δ-V 2 O 5 · n H 2 O electrodes exhibited only 10% capacity retention in the same experiment. Remarkably, 2D δ-C x V 2 O 5 · n H 2 O heterostructure electrodes demonstrated significantly improved capacity retention (94% after 30 cycles) for bilayered vanadium oxide electrodes in Li-ion cells during galvanostatic cycling at 20 mA g –1 . The improved electrochemical performance, in both extended cycling and rate capability studies, of the 2D δ-C x V 2 O 5 · n H 2 O heterostructure electrodes in the Li-ion system is ascribed to the intermittent formation of carbon layers within the bilayered structure, which leads to increased electronic conductivity and improved structural stability of the heterostructure compared to the reference δ-V 2 O 5 · n H 2 O electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Influence of Ti‐incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H 2 O 2

Abstract Vapor‐phase propylene (C 3 H 6 ) epoxidation kinetics with hydrogen peroxide (H 2 O 2 ) strongly reflects the physical properties of Ti‐incorporated zeolite catalysts and the presence of spectating molecules (“solvent”) near active sites even without a bulk liquid phase. Steady‐state turnover rates of C 3 H 6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH) x ) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH 3 CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H 2 O 2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring‐opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C 3 H 6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH 3 CN. Thus, higher PO selectivities achieved in hydrophilic Ti‐MFI at steady‐state reflect the preferential stabilization of transition states for C 3 H 6 epoxidation (a primary reaction) relative to PO ring‐opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH 3 CN within pores and the extended pore structure.

Kwon, Ohsung↗
DOE OSTI2024

Influence of Ti‐incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H 2 O 2

Vapor-phase propylene (C 3 H 6 ) epoxidation kinetics with hydrogen peroxide (H 2 O 2 ) strongly reflects the physical properties of Ti-incorporated zeolite catalysts and the presence of spectating molecules (“solvent”) near active sites even without a bulk liquid phase. Steady-state turnover rates of C 3 H 6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH) x ) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH 3 CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H 2 O 2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring-opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C 3 H 6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH 3 CN. Thus, higher PO selectivities achieved in hydrophilic Ti-MFI at steady-state reflect the preferential stabilization of transition states for C 3 H 6 epoxidation (a primary reaction) relative to PO ring-opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH 3 CN within pores and the extended pore structure.

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

Photodetachment Dynamics and Structural Flexibility of Undercoordinated Iridium Halides IrCl n − (n = 3−5): An Experimental and Theoretical Investigation

Three undercoordinated iridium chloride anions, IrCl n − (n = 3–5), and their neutral counterparts were investigated by cryogenic anion photoelectron spectroscopy and theoretical calculations. Photodetachment of IrCl n − leads to the formation of the corresponding neutral complex, i.e., a triplet ground state for n = 3, a quartet for n = 4, and close-lying singlet and triplet for n = 5. The vertical detachment energies are determined to be 3.89, 4.98, and 5.14 eV for n = 3, 4, and 5, respectively, revealing superhalogen anion properties with increasing electron detachment energies as chloride ligands added. The IrCl 3 − spectrum features an extremely broad, lowest electron binding energy band, attributed to resonant autodetachment with prominent non-Franck–Condon profiles. In IrCl 5 − , detachment prompts a d-orbital rearrangement that drives a structural transformation from a twisted square-based pyramidal to a trigonal–bipyramidal geometry in the singlet state. In conclusion, these findings provide insights into the electronic and structural adaptability of iridium halides, advancing our understanding of ligand exchange reactions and dissociative stability in transition metal complexes.

Tang, Peng [East China Normal Univ. (ECNU), Shangh↗
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2023

Developing a roadmap for carbon capture, and storage in Oklahoma by assessing the viability of stacked storage

Abstract The Intergovernmental Panel on Climate Change concludes that CO 2 capture and storage (CCS) is critical for climate‐stabilizing energy transitions. In CCS, captured CO 2 is sequestered in saline aquifers within sedimentary basins. The CO 2 storage capacity and the rate of injection are functions of the geology of the saline aquifer, which is uncertain. To minimize impacts of this uncertainty, CCS projects could include backup plans, such as co‐locating geologic CO 2 storage (GCS) sites with or near existing CO 2 ‐enhanced oil recovery (CO 2 ‐EOR) operations. These “stacked storage” projects could hedge against uncertainty in the saline formation performance because captured CO 2 could be injected into either location in the event of unexpected events (e.g., the injectivity decreases). Here, we investigate the possibility and ramifications of developing CCS networks in Oklahoma that are amendable to stacked storage. We find that stacked storage is possible in Oklahoma but the counties with the lowest‐cost saline storage resources do not have existing CO 2 ‐EOR operations. At the systems level, we find it is slightly more expensive (e.g., $1/tCO 2 to $5/tCO 2 ) to site GCS in counties with CO 2 ‐EOR projects. This increased expense is largely due to increased CO 2 transportation costs because hundreds of km of additional pipeline is required to capture CO 2 from the lowest‐cost sources. Overall, our results suggest that it is optimal to build more pipelines and avoid injecting CO 2 in some of the lowest‐cost saline storage resources, to enable capturing CO 2 from the least‐cost sources. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

42 ENGINEERING↗
DOE OSTI2021

Flame aerosol synthesis of hollow alumina nanoshells for application in thermal insulation

Commercial implementation of hollow nanoshell materials is often limited by difficult scale-up and high cost. Here we present a scalable, continuous, and potentially low-cost route to hollow nanostructured materials using a High Temperature Reducing Jet flame aerosol reactor. The key advantage of this process over traditional flame aerosol synthesis methods is separation of flame chemistry and particle formation, which allows production of hollow nanoshells from an aqueous salt solution without any template or surfactant. We produced amorphous hollow alumina with an average diameter of 181 nm and shell thickness below 10 nm, and demonstrated the generality of the method by also producing hollow SiO 2 , ZrO 2 , Ga 2 O 3 , and Co 2 O 3 . We characterized the morphology, phase transition, thermal stability, and thermal conductivity of the flame-made hollow alumina by electron microscopy, x-ray diffraction, and thermal analysis methods. As a prototypical application, we incorporated hollow alumina into thermal insulation, combining it with common fiberglass insulation to reduce the thermal conductivity by 30% relative to pure fiberglass, to 0.028 W m -1 K -1 while maintaining robust fire-resistance and elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2022

Acetone to isobutene conversion on Zn x Ti y O z : Effects of TiO 2 facet

In this study, liquid-phase chemical grafting method was used to graft Zn onto TiO 2 with preferentially exposed (1 0 1) or (0 0 1) facet. The obtained Zn x Ti y O z materials were characterized using various techniques (e.g., XRD, Raman, DRIFTS etc.) and evaluated for the acetone-to-isobutene reaction. It was found that over TiO 2 (0 0 1), both terminal and bridging hydroxyls were readily titrated by Zn deposition, whereas a substantial amount of bridging hydroxyls on TiO 2 (1 0 1) remained. Although dominant Zn-O-terminated surface was obtained on two Zn x Ti y O z samples, bridging hydroxyls with high H-D exchange reactivity were observed on Zn x Ti y O z (1 0 1) compared with Zn x Ti y O z (0 0 1). The bridging hydroxyls showing rapid proton transfer efficiently stabilizes a transition state of diacetone alcohol intramolecular rearrangement for isobutene production as opposed to the diacetone alcohol dehydration.

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

Linear Free Energy Relationships Associated with Hydride Transfer From [(6,6'-R 2 -bpy)Re(CO) 3 H]: A Cautionary Tale in Identifying Hydrogen Bonding Effects in the Secondary Coordination Sphere

Here, six rhenium hydride complexes, [(6,6'-R 2 -bpy)Re(CO) 3 H] (bpy = 2,2'-bipyridine, R = OEt, OMe, NHMe, Me, F, Br), were synthesized. These complexes insert CO 2 to form rhenium formate complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 {OC(O)H}]. All the rhenium formate species were characterized using X-ray crystallography, which revealed that the bpy ligand is not coplanar with the metal coordination plane containing the two nitrogen donors of the bpy ligand but tilted. A solid-state structure of [(6,6'-Me 2 -bpy)Re(CO) 3 H] determined using MicroED also featured a tilted bpy ligand. The kinetics of CO 2 insertion into complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 H] were measured experimentally and the thermodynamic hydricities of [(6,6'-R 2 -bpy)Re(CO) 3 H] species were determined using theoretical calculations. A Brønsted plot constructed using the experimentally determined rate constants for CO 2 insertion and the calculated thermodynamic hydricities for [(6,6'-R 2 -bpy)Re(CO) 3 H] revealed a linear free energy relationship (LFER) between thermodynamic and kinetic hydricity. This LFER is different to the previously determined relationship for CO 2 insertion into complexes of the type [(4,4'-R 2 -bpy)Re(CO) 3 H]. At a given thermodynamic hydricity, CO 2 insertion is faster for complexes containing a 6,6'-substituted bpy ligand. This is likely in part due to the tilting observed for systems with 6,6'-substituted bpy ligands. Notably, the 6,6'-(NHMe) 2 -bpy ligand could in principle stabilize the transition state for CO 2 insertion via hydrogen bonding. This work shows that if only the rate of CO 2 insertion into [(6,6'-(NHMe) 2 -bpy)Re(CO) 3 H] is compared to [(4,4'-R 2 -bpy)Re(CO) 3 H] systems, the increase in rate could be easily attributed to hydrogen bonding, but in fact all 6,6'-substituted systems lead to faster than expected rates.

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

Improving Estimation of the Koopman Operator with Kolmogorov–Smirnov Indicator Functions

It has become common to perform kinetic analysis using approximate Koopman operators that transform high-dimensional timeseries of observables into ranked dynamical modes. The key to the practical success of the approach is the identification of a set of observables that form a good basis on which to expand the slow relaxation modes. Good observables are, however, difficult to identify a priori and suboptimal choices can lead to significant underestimations of characteristic time scales. Leveraging the representation of slow dynamics in terms of Hidden Markov Models (HMM), we propose a simple and computationally efficient clustering procedure to infer surrogate observables that form a good basis for slow modes. Here, we apply the approach to an analytically solvable model system as well as on three protein systems of different complexities. We consistently demonstrate that the inferred indicator functions can significantly improve the estimation of the leading eigenvalues of Koopman operators and correctly identify key states and transition time scales of stochastic systems, even when good observables are not known a priori.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗