Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation software”

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.

346 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Self‐Propelling Macroscale Sheets Powered by Enzyme Pumps

Nanoscale enzymes anchored to surfaces act as chemical pumps by converting chemical energy released from enzymatic reactions into spontaneous fluid flow that propels entrained nano‐ and microparticles. Enzymatic pumps are biocompatible, highly selective, and display unique substrate specificity. Utilizing these pumps to trigger self‐propelled motion on the macroscale has, however, constituted a significant challenge and thus prevented their adaptation in macroscopic fluidic devices and soft robotics. Using experiments and simulations, we herein show that enzymatic pumps can drive centimeter‐scale polymer sheets along directed linear paths and rotational trajectories. In these studies, the sheets are confined to the air/water interface. With the addition of appropriate substrate, the asymmetric enzymatic coating on the sheets induces chemically driven, buoyancy flows that controllably propel the sheet's motion on the air/water interface. The directionality and speed of the motion can be tailored by changing the pattern of the enzymatic coating, type of enzyme, and nature and concentration of the substrate. This work highlights the utility of biocompatible enzymes for generating motion in macroscale fluidic devices and robotics and indicates their potential utility for in vivo applications.

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

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2025

Investigating the Combustion Performance of Dual Fuel Combustion with Diesel and Port Injected Hydrogen in a Large Bore Locomotive Engine

The heavy-duty transportation sector has primarily relied on conventional diesel combustion engines given their reliability and high thermal efficiency relative to spark ignition engines, but increased focus on reducing greenhouse gas emissions has led to investigation into alternative fuels. Gaseous hydrogen fuel has garnered a great deal of recent interest in the engine community given it has zero carbon, but hydrogen is not available at the scale and cost that petroleum fuels are currently available, and this is a barrier to adoption for industries that are looking to decarbonize their operations. Because of the fuel flexibility provided, dual fuel technology offers a pathway for some industries to adopt hydrogen as a fuel source while maintaining sufficient flexibility in times and locations where the new fuel is not yet available. This computational study investigates dual fuel combustion in a large bore locomotive engine architecture using direct injected diesel and port injected gaseous hydrogen fuel. With an optimal port fuel injection configuration from previous work, simulations of varying substitution ratio, compression ratio, manifold air temperature, diesel injection timing, and diesel injection pressure were performed to understand their effect on combustion performance. Results indicated that both increased substitution ratio and higher intake air temperature accelerates hydrogen flame propagation and can result in high peak cylinder pressures. Additionally, diesel injection timing and injection pressure were demonstrated as effective methods for controlling dual fuel combustion heat release rates.

ODonnell, Patrick Christopher↗
DOE OSTI2025

Novel Zwitterionic Polyurethane-in-Salt Electrolytes with High Ion Conductivity, Elasticity, and Adhesion for High-Performance Solid-State Lithium Metal Batteries

This study presents a novel polymer-in-salt (PIS) zwitterionic polyurethane-based solid polymer electrolyte (zPU-SPE) that offers high ionic conductivity, strong interaction with electrodes, and excellent mechanical and electrochemical stabilities, making it promising for high-performance all solid-state lithium batteries (ASSLBs). The zPU-SPE exhibits remarkable lithium-ion (Li+) conductivity (3.7 × 10⁻⁴ S cm−1 at 25 °C), enabled by exceptionally high salt loading of up to 90 wt.% (12.6 molar ratio of Li salt to polymer unit) without phase separation. It addresses the limitations of conventional SPEs by combining high ionic conductivity with a Li+ transference number of 0.44, achieved through the incorporation of zwitterionic groups that enhance ion dissociation and transport. The high surface energy (338.4 J m−2) and elasticity ensure excellent adhesion to Li anodes, reducing interfacial resistance and ensuring uniform Li+ flux. When tested in Li||zPU||LiFePO₄ and Li||zPU||S/C cells, the zPU-SPE demonstrated remarkable cycling stability, retaining 76% capacity after 2000 cycles with the LiFePO4 cathode, and achieving 84% capacity retention after 300 cycles with the S/C cathode. Molecular simulations and a range of experimental characterizations confirm the superior structural organization of the zPU matrix, contributing to its outstanding electrochemical performance. The findings strongly suggest that zPU-SPE is a promising candidate for next-generation ASSLBs.

Wang, Kun↗
DOE OSTI2024

Trust-Based Detection and Mitigation of Cyber Attacks in Distributed Cooperative Control of Islanded AC Microgrids

In this study, we address the challenge of detecting and mitigating cyber attacks in the distributed cooperative control of islanded AC microgrids, with a particular focus on detecting False Data Injection Attacks (FDIAs), a significant threat to the Smart Grid (SG). The SG integrates traditional power systems with communication networks, creating a complex system with numerous vulnerable links, making it a prime target for cyber attacks. These attacks can lead to the disclosure of private data, control network failures, and even blackouts. Unlike machine learning-based approaches that require extensive datasets and mathematical models dependent on accurate system modeling, our method is free from such dependencies. To enhance the microgrid’s resilience against these threats, we propose a resilient control algorithm by introducing a novel trustworthiness parameter into the traditional cooperative control algorithm. Our method evaluates the trustworthiness of distributed energy resources (DERs) based on their voltage measurements and exchanged information, using Kullback-Leibler (KL) divergence to dynamically adjust control actions. We validated our approach through simulations on both the IEEE-34 bus feeder system with eight DERs and a larger microgrid with twenty-two DERs. The results demonstrated a detection accuracy of around 100%, with millisecond range mitigation time, ensuring rapid system recovery. Additionally, our method improved system stability by up to almost 100% under attack scenarios, showcasing its effectiveness in promptly detecting attacks and maintaining system resilience. These findings highlight the potential of our approach to enhance the security and stability of microgrid systems in the face of cyber threats.

Computer Science↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

Plasma flow generation and particle acceleration from expanding magnetic bubbles

Impulsive plasma dynamics in the laboratory are often driven by rising electric currents, yet their quantitative plasma response has not been well established. By means of fully kinetic particle-in-cell simulations and laser-driven capacitor-coil experiments, we show that a rising current expels plasma, forming an expanding magnetic bubble and accelerating particles. The expansion front velocity scales with the Alfvén speed determined by the magnetic field at its inner edge and the plasma density at its outer edge. This mechanism establishes impulsive current drive as a fundamental way that generates plasma flows and accelerates particles in laboratory plasmas, with potential relevance to astrophysics.

Zhang, Yang [Princeton University, NJ (United Stat↗
DOE OSTI2025

Detection of visible-wavelength aurora on Mars

Mars hosts various auroral processes despite the planet’s tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near–real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Science & Technology - Other Topics↗
DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗
DOE OSTI2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Abstract Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integratesin-situmeasurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697µmolCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Environmental Sciences & Ecology↗
DOE OSTI2026

Update on Radiochemical Assessment of High Burnup Commercially Irradiated Fuel

This work documents an effort to collect burnup measurements on a high burnup rod, designated 6XV, and first cycle accident tolerant fuel (ATF) rod, designated 47I, to enable benchmarking of fuel performance codes and neutronics codes. In addition to measurements, Virtual Environment for Reactor Applications (VERA) full-core-depletion analysis was also performed for the rods that were experimentally analyzed to provide an opportunity for code validation. This effort focuses on collecting data from rods irradiated at Byron Generating Station and shipped to the Oak Ridge National Laboratory (ORNL) hot-cells. This data will also anchor non-destructive examination evaluations of burnup of the various fuel rods undergoing postirradiation examination (PIE) at ORNL. Previous PIE of these fuel rods provides some guidance on the burnup trend across the fuel. Axial gamma spectroscopy scans provide a measure of relative changes in burnup across a fuel pin. Mass spectrometry based burnup measurements performed for this work at specific axial locations in the fuel are fully quantitative. By combining the mass spectrometry data with the gamma scans it is possible to more quantitatively evaluate axial variations in burnup across the entire fuel pin [1]. The combined set of burnup evaluations will be made available to other organizations that have an interest in high burnup radiochemistry data for validation of neutronic simulations and source term evaluation such as the Nuclear Regulatory Commission (NRC).

Harp, Jason [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2026

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2024

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry↗