Search NASA⌕ Search

SEARCH · Search NASA

Results for “Theoretical Kinetics”

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

Conformer-Dependent Chemistry: Experimental Product Branching of the Vinyl Alcohol + OH + O 2 Reaction

The concentration of formic acid in Earth's tropo-sphere is underestimated by detailed chemical models compared to field observations. Phototautomerization of acetaldehyde to its less stable tautomer vinyl alcohol, followed by the OH-initiated oxidation of vinyl alcohol, has been proposed as a missing source of formic acid that improves the agreement between models and field measurements. Theoretical investigations of the OH + vinyl alcohol reaction in excess O 2 conclude that OH addition to the alpha carbon of vinyl alcohol produces formaldehyde + formic acid + OH, whereas OH addition to the β site leads to glycoaldehyde + HO 2 . Furthermore, these studies predict that the conformeric structure of vinyl alcohol controls the reaction pathway, with the anti-conformer of vinyl alcohol promoting α OH addition, whereas the syn-conformer promotes β addition. However, the two theoretical studies reach different conclusions regarding which set of products dominate. We studied this reaction using time-resolved multiplexed photoionization mass spectrometry to quantify the product branching fractions. Our results, supported by a detailed kinetic model, conclude that the glycoaldehyde product channel (arising mostly from syn-vinyl alcohol) dominates over formic acid production with a 3.6:1.0 branching ratio. Additionally, this result supports the conclusion of Lei et al. that conformer-dependent hydrogen bonding at the transition state for OH-addition controls the reaction outcome. As a result, tropospheric oxidation of vinyl alcohol creates less formic acid than recently thought, increasing again the discrepancy between models and field observations of Earth's formic acid budget.

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

OH Roaming during the Ozonolysis of α-Pinene: A New Route to Highly Oxygenated Molecules?

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. Here, in this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ~20% at room temperature up to ~70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

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

Auger electron angular distributions following excitation or ionization from the Xe 3d and F 1s levels in xenon difluoride

Linearly polarized synchrotron radiation has been used to record polarization dependent, non-resonant Auger electron spectra of XeF 2 , encompassing the bands due to the xenon M 45 N 1 N 45 , M 45 N 23 N 45 , M 45 N 45 N 45 and M 45 N 45 V and fluorine KVV transitions. Resonantly excited Auger spectra have been measured at photon energies coinciding with the Xe 3d 5/2 → σ* and the overlapped Xe 3d 3/2 /F 1s → σ* excitations in XeF 2 . The non-resonant and resonantly excited spectra have enabled the Auger electron angular distributions, as characterized by the β A parameter, to be determined for the M 45 N 45 N 45 transitions. In the photon energy range over which the Auger electron angular distributions were measured, theoretical results indicate that transitions into the εf continuum channel dominate the Xe 3d photoionization in XeF 2 . In this limit, the theoretical value of the atomic alignment parameter (A 20 ) characterizing the core ionized state becomes constant. This theoretical value has been used to obtain the Auger electron intrinsic anisotropy parameters (α 2 ) from the β A parameters extracted from our non-resonant Auger spectra. For a particular Auger transition, the electron kinetic energy measured in the resonantly excited spectrum is higher than that in the directly ionized spectrum, due to the screening provided by the electron promoted into the σ* orbital. The interpretation of the F KVV Auger band in XeF 2 has been discussed in relation to previously published one-site populations of the doubly charged ions (XeF 2 2+ ). Here, the experimental results show that the ionization energies of the doubly charged states predominantly populated in the decay of a vacancy in the F 1s orbital in XeF 2 tend to be higher than those populated in the decay of a vacancy in the Xe 4d level in XeF 2 .

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

Modeling Hydrogen Evolution Reaction Kinetics through Explicit Water–Metal Interfaces

Despite the apparent simplicity of the hydrogen evolution reaction (HER) and the decades of research into it, controversy remains in the literature regarding the identity of the active site and the competition between the Heyrovsky and Tafel steps. In this work, we use charge-extrapolated ab initio simulations with explicit water in conjunction with mean-field microkinetic modeling to explore the mechanism for HER on both close-packed (111) and stepped (211) transition metals. First, we show that atop H*, beyond a monolayer of hollow H*, is unlikely to play a role in the HER mechanism, given its very positive adsorption energies. The energetics suggests the Volmer–Heyrovsky mechanism to predominate on fcc transition metals under typical operating conditions. Here, we evaluate our theoretical results vs several experimental observations. We show that the Volmer–Heyrovsky mechanism predicts an activity volcano with its peak at a H* binding ΔG H* ≈ 0 eV, consistent with experiment. In contrast, the Volmer–Tafel volcano shows a broad rate plateau between ΔG H* ≈ 0 eV and ΔG H* ≈ – 0.4 eV. We find our theoretical Tafel slopes to be consistent with experimental ones on a range of transition metals. We show that, in line with experimental observations, the introduction of a CO(g) atmosphere shifts the strong binding metals toward the weak binding leg. Our study suggests that the simple thermodynamic approach to HER activity still holds, even when a detailed kinetic picture is considered.

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

Design of Ru-Ni diatomic sites for efficient alkaline hydrogen oxidation

Anion exchange membrane fuel cells are limited by the slow kinetics of alkaline hydrogen oxidation reaction (HOR). Here, we establish HOR catalytic activities of single-atom and diatomic sites as a function of *H and *OH binding energies to screen the optimal active sites for the HOR. As a result, the Ru-Ni diatomic one is identified as the best active center. Guided by the theoretical finding, we subsequently synthesize a catalyst with Ru-Ni diatomic sites supported on N-doped porous carbon, which exhibits excellent catalytic activity, CO tolerance, and stability for alkaline HOR and is also superior to single-site counterparts. In situ scanning electrochemical microscopy study validates the HOR activity resulting from the Ru-Ni diatomic sites. Furthermore, in situ x-ray absorption spectroscopy and computational studies unveil a synergistic interaction between Ru and Ni to promote the molecular H 2 dissociation and strengthen OH adsorption at the diatomic sites, and thus enhance the kinetics of HOR.

25 ENERGY STORAGE↗
DOE OSTI2022

In situ recording of Mars soundscape

Abstract Before the Perseverance rover landing, the acoustic environment of Mars was unknown. Models predicted that: (1) atmospheric turbulence changes at centimetre scales or smaller at the point where molecular viscosity converts kinetic energy into heat 1 , (2) the speed of sound varies at the surface with frequency 2,3 and (3) high-frequency waves are strongly attenuated with distance in CO 2 (refs. 2–4 ). However, theoretical models were uncertain because of a lack of experimental data at low pressure and the difficulty to characterize turbulence or attenuation in a closed environment. Here, using Perseverance microphone recordings, we present the first characterization of the acoustic environment on Mars and pressure fluctuations in the audible range and beyond, from 20 Hz to 50 kHz. We find that atmospheric sounds extend measurements of pressure variations down to 1,000 times smaller scales than ever observed before, showing a dissipative regime extending over five orders of magnitude in energy. Using point sources of sound (Ingenuity rotorcraft, laser-induced sparks), we highlight two distinct values for the speed of sound that are about 10 m s −1 apart below and above 240 Hz, a unique characteristic of low-pressure CO 2 -dominated atmosphere. We also provide the acoustic attenuation with distance above 2 kHz, allowing us to explain the large contribution of the CO 2 vibrational relaxation in the audible range. These results establish a ground truth for the modelling of acoustic processes, which is critical for studies in atmospheres such as those of Mars and Venus.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fundamental limitations on photoisomerization from thermodynamic resource theories

Small, out-of-equilibrium, and quantum systems defy simple thermodynamic expressions. Such systems are exemplified by molecular switches, which exchange heat with a bath. These molecules can photoisomerize, or change conformation, or switch, on absorbing light. The photoisomerization probability depends on kinetic details that couple the molecule’s energetics to its dissipation. Therefore, a simple, general, thermodynamic-style bound on the photoisomerization probability seems out of reach. We derive such a bound using a resource theory. The resource-theory framework is a set of mathematical tools, developed in quantum information theory, used to generalize thermodynamics to small and quantum settings. From this toolkit has been derived a generalization of the second law, the thermomajorization preorder. We use thermomajorization to upper-bound the photoisomerization probability. Then, we compare the bound with an equilibrium prediction and with a Lindbladian model. We identify a realistic parameter regime in which the Lindbladian evolution saturates the thermomajorization bound. We also quantify the energy coherence in the electronic degree of freedom, and we argue that this coherence cannot promote photoisomerization. Overall, this work illustrates how quantum-information-theoretic thermodynamics can elucidate complex quantum processes in nature, experiments, and synthetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2023

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

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

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

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

Regression Analysis with the Directed Infusion of Data

Integrating artificial intelligence and machine learning tools into industry necessitates large-scale collaborative efforts that ensure the robust and accurate execution of downstream analytics such as time series prediction, uncertainty quantification, grid optimization, and condition monitoring. However, concerns related to data privacy pervade the nuclear industry due to the proprietary nature of its data and the possibility of data leakage. Legacy techniques such as encryption often require the explicit transmission of data to trustworthy parties, thereby inviting data leakage concerns. The ideal collaboration scenario avoids the explicit dissemination of data/code while maintaining experimental fidelity, which is currently accomplished using various techniques such as trusted execution environments, homomorphic encryption, differential privacy, and multimatrix masking. These techniques, however, often necessitate a trade-off between trust, efficiency, and utility. This article extends a previously proposed technique called the directed infusion of data (DIOD) that ensures data privacy, allows for scalable obfuscation, and combats the risk of data leakage without compromising utility. The experiments discussed in this article examine a regression-type scenario using DIOD with the goal of preserving the inferential link between two variables. Using the point-kinetics equations, regression experiments compare the performance of a model trained using the original data to that of a model trained using the obfuscated data, which produced identical results. Our claim is further strengthened by an information theoretic proof and experiment, which showed that the inferential content between variables remains the same after obfuscation, thereby avoiding the required communication of the proprietary data.

47 - OTHER INSTRUMENTATION↗
DOE OSTI2024

Photoelectron spectra in circularly and elliptically polarized laser pulses

We present results of numerical simulations of the time-dependent Schrödinger equation and theoretical analysis concerning the interaction of a rare gas atom with circularly and elliptically polarized laser pulses. In agreement with recent observation in circularly polarized fields, the photoelectron energy spectra for counterrotating electrons are peaked at lower kinetic energy than those for corotating electrons. We show that this difference can be interpreted as being due to the additional pathways to ionization that are available for the counterrotating electrons only. Furthermore, our results show, in agreement with earlier work, that the offset angle by which the emission of electrons is rotated in an elliptically polarized field increases with each successive above-threshold ionization (ATI) order and is larger for the emission from counterrotating states as compared to that from corotating states. As a result, a simple model based on the interference and relative phase difference between just three continuum states provides remarkable agreement and once again emphasizes the importance of the additional ionization pathways available for the counterrotating electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗
DOE OSTI2024

Exploring the potential and impact of single-crystal active materials on dry-processed electrodes for high-performance lithium-ion batteries

Roll-to-roll powder-to-film dry processing (DP) and single-crystal (SC) active materials (AMs) with many advantages are two hot topics of lithium-ion batteries (LIBs). However, DP of SC AMs for LIBs is rarely reported. Consequently, the impact of SC AMs on dry-processed LIBs is not well understood. Herein, for the first time, via a set of experimental and theoretical studies of the conventional polycrystalline-AM- and SC-AM-based DPed electrodes (DPEs), this work not only reports a high-performance dry SC-AM cathode for LIB manufacturing, but also establishes some fundamental understanding of SC-based dry-processed electrodes, including their morphology, structure, mechanical strength, electronic conductivity and LIB electrochemical behavior. Further, the results suggest that DP of SC AMs is promising, which can dramatically improve the electrochemical kinetics at electrode level and particle level. Specifically, for the rate capability and long-term cyclability in full cells, SC DPEs exhibit a discharge specific capacity of 152.1 mAh g -1 at 1C and a capacity retention rate of 79.9 % at C/3 over 500 cycles, which are superior to those of PC DPEs (135.6 mAh g -1 and 68.3 %) at the same conditions and are further confirmed by the simulation data from the theoretical modelling study. Therefore, this comprehensive work marks a significant milestone for DP strategy and SC AMs, enlightening future research and development of LIB manufacturing.

25 ENERGY STORAGE↗
DOE OSTI2020

A single-component water-lean post-combustion CO 2 capture solvent with exceptionally low operational heat and total costs of capture – comprehensive experimental and theoretical evaluation

A comprehensive evaluation of a recently developed water-lean amine-based solvent, namely N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (2-EEMPA), has been performed to analyze its post-combustion CO 2 capture performance. This evaluation comprises (1) fundamental characterization of the solvent–CO 2 interaction using vapor–liquid equilibria, kinetics and viscosity measurements; (2) process characterization of the CO 2 capture performance as measured in a laboratory scale continuous flow system and via Aspen Plus® simulation using a flue gas simulant; as well as (3) a full techno economic analysis of the capture process at industrial scale with corresponding projections of critical metrics. This paper summarizes the many parts of this comprehensive evaluation and shows how the various parts come together to empower validated conclusions about its process performance. Notably, it is projected that this solvent can operate at a regeneration heat rate of 2.0 GJ per tonne CO 2 for post-combustion capture, and at a total cost of capture of $50.6/tonne CO 2 . Overall, with further process optimization significant reductions in the capture cost are predicted.

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

Realizing Two-Electron Transfer in Ni(OH) 2 Nanosheets for Energy Storage

The theoretical capacity of a given electrode material is ultimately determined by the number of electrons transferred in each redox center. The design of multi-electron transfer processes could break through the limitation of one–electron transfer and multiply the total capacity, but is difficult to achieve because multiple electron transfer processes are generally thermodynamically and kinetically more complex. In this work, we report the discovery of two-electron transfer in monolayer Ni(OH) 2 nanosheets, which contrasts with the traditional one-electron transfer found in multilayer materials. The first-principles calculations predict that the first oxidation process from Ni 2+ →Ni 3+ occurs easily, whereas the second electron transfer from Ni 3+ →Ni 4+ is strongly hindered in multilayer materials by both the interlayer hydrogen bonds and the domain H structure induced by the Jahn-Teller distortion of the Ni 3+ (t 2g 6 e g 1 ) centered octahedra. In contrast, the second electron transfer can easily occur in monolayers because all H’s are fully exposed. Experimentally, the as-prepared monolayer is found to deliver an exceptional redox capacity of ~ 576 mAh/g, nearly two times the theoretical capacity of one-electron processes. In-situ experiments demonstrate that monolayer Ni(OH) 2 can transfer two electrons and mostly transform to Ni 4+ in the charging process, while bulk Ni(OH) 2 can be partially transformed . Our work reveals a new redox reaction mechanism in atomically thin Ni(OH) 2 nanosheets, and suggests a promising path towards tuning the electron transfer numbers to multiply the capacity of the relevant energy storage materials.

25 ENERGY STORAGE↗