Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT) calculations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

DFT simulation of the X-ray diffraction pattern of aluminum-ion-intercalated graphite used as the cathode material of the aluminum-ion battery

Many breakthroughs have been achieved in rechargeable aluminum-ion battery technologies in recent years. Most recently, operando X-ray diffraction (XRD) combined with density functional theory (DFT) calculations was reported to study the chloroaluminate anion (AlCl 4 -)-intercalated graphite cathode of the battery. However, there are quite a few discrepancies between the measured and simulated XRD patterns. This work is focused on the simulation of XRD patterns of graphite intercalation compounds (GICs) with DFT calculations. Our results reveal that both the geometry of AlCl 4 - in graphite and the gallery height of GICs are dependent on the intercalant density. At low intercalant density, the gallery height keeps constant, but at high intercalant densities, the gallery height is linearly related to the intercalant density. Our simulated XRD patterns are highly consistent with the measured operando XRD patterns. Not only do the angles of the peaks match very well, but also the relative intensities and the corresponding electrode capacities show reasonable agreement with the experimental results. The DFT simulation of the XRD pattern provides significant information on the stage index and the charge capacity of the GIC electrode.

25 ENERGY STORAGE↗
DOE OSTI2025

Fermi surface reconstruction due to orthorhombic distortion in the Dirac semimetal YbMnSb 2

A Dirac semimetal with magnetic atoms as constituents delivers an interesting platform to investigate the interplay of Fermi surface (FS) topology, electron correlation, and magnetism. One such family of semimetals is YbMn⁢ 𝑃⁢𝑛 2 (𝑃⁢𝑛 = Sb,Bi), which is being actively studied due to the intertwined spin and charge degrees of freedom. Here, in this study, we investigate the relationship between the magnetic/crystal structures and FS topology of YbMnSb 2 using single crystal x-ray diffraction, neutron scattering, magnetic susceptibility, magnetotransport measurements, and complimentary density functional theory (DFT) calculations. Contrary to previous reports, the x-ray and neutron diffraction reveal that YbMnSb 2 crystallizes in an orthorhombic 𝑃⁢𝑛⁢𝑚⁢𝑎 structure with a notable antiphase displacement of the magnetic Mn ions that increases in magnitude upon cooling. As a result of this orthorhombicity, first-principles DFT calculation reveals a reduced Brillouin zone and more anisotropic FS of YbMnSb 2 compared to YbMnBi 2 . Moreover, the hole-type carrier density drops by two orders of magnitude as YbMnSb 2 orders antiferromagnetically, indicating band folding in a magnetic ordered state. In addition, the Landau level fan diagram yields a nontrivial nature of the Shubnikov–de Haas quantum oscillation frequency arising from the Dirac-like Fermi pocket. These results imply that YbMnSb 2 is an ideal platform to explore the interplay of subtle lattice distortion, magnetic order, and topological transport arising from relativistic quasiparticles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Critical Assessment of Electronic Structure Descriptors for Predicting Perovskite Catalytic Properties

The discovery and design of materials which can efficiently catalyze the oxygen reduction and evolution reactions at reduced temperatures is important for facilitating the widespread adoption of fuel cell and electrolyzer technologies. Numerous studies have produced correlations between catalytic properties, such as oxygen surface exchange or electrode area specific resistance (ASR), and properties of the catalyst material. However, correlations have historically been limited in scope (e.g., using only a few materials or at a single temperature) and it has been difficult to provide detailed assessments of their robustness. Here, in this study, we assess the ability of the O p-band center electronic structure descriptor, obtained from density functional theory (DFT) calculations, to correlate with oxygen surface exchange rates, diffusivities, and area specific resistances for a large database of perovskite oxide catalytic properties. By data mining the literature, we obtain 747 catalytic property value data points spanning 299 unique perovskite compositions from 313 studies. We assess linear correlations of each property with the O p-band center and find generally modest correlations that are qualitatively useful (prediction mean absolute errors of about 0.5 log units are typical), where the correlations are improved at higher temperatures (e.g., 800 °C vs. 500 °C) and significantly improve when considering fits to the subset of materials which have multiple independent measurements. These findings suggest that the spread of property data is significantly influenced by experimental uncertainty, and subsequent measurements of additional materials will likely improve the O p-band center correlations.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Prediction of organic homolytic bond dissociation enthalpies at near chemical accuracy with sub-second computational cost

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol-1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

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

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Spin-polarized imaging of strongly interacting fermions in the ferrimagnetic state of the Weyl candidate CeBi

CeBi has an intricate magnetic phase diagram whose fully polarized state has recently been suggested as a Weyl semimetal, though the role of f states in promoting strong interactions has remained elusive. Here we focus on the less-studied but also time-reversal symmetry-breaking ferrimagnetic phase of CeBi, where our density functional theory (DFT) calculations predict additional Weyl nodes near the Fermi level E F . We use spin-polarized scanning tunneling microscopy and spectroscopy to image the surface ferrimagnetic order on the itinerant Bi p states, indicating their orbital hybridization with localized Ce f states. We observe suppression of this spin-polarized signature at E F , coincident with a Fano line shape in the conductance spectra, suggesting the Bi p states partially Kondo screen the f magnetic moments, and this p – f hybridization causes strong Fermi-level band renormalization. Furthermore, the p-band flattening is supported by our quasiparticle interference measurements, which also show band splitting in agreement with DFT, painting a consistent picture of a strongly interacting magnetic Weyl semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

In-situ structure and catalytic mechanism of NiFe and CoFe layered double hydroxides during oxygen evolution

NiFe and CoFe (MFe) layered double hydroxides (LDHs) are among the most active electrocatalysts for the alkaline oxygen evolution reaction (OER). Herein, we combine electrochemical measurements, operando X-ray scattering and absorption spectroscopy, and density functional theory (DFT) calculations to elucidate the catalytically active phase, reaction center and the OER mechanism. We provide the first direct atomic-scale evidence that, under applied anodic potentials, MFe LDHs oxidize from as-prepared α-phases to activated γ-phases. The OER-active γ-phases are characterized by about 8% contraction of the lattice spacing and switching of the intercalated ions. DFT calculations reveal that the OER proceeds via a Mars van Krevelen mechanism. The flexible electronic structure of the surface Fe sites, and their synergy with nearest-neighbor M sites through formation of O-bridged Fe-M reaction centers, stabilize OER intermediates that are unfavorable on pure M-M centers and single Fe sites, fundamentally accounting for the high catalytic activity of MFe LDHs.

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

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications

Here, the increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH 3 ) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

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

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2022

Epitaxial growth, magnetoresistance, and electronic band structure of GdSb magnetic semimetal films

Motivated by observations of extreme magnetoresistance (XMR) in bulk crystals of rare-earth monopnictide (RE-V) compounds and emerging applications in novel spintronic and plasmonic devices based on thin-film semimetals, we have investigated the electronic band structure and transport behavior of epitaxial GdSb thin films grown on III-V semiconductor surfaces. The Gd 3 + ion in GdSb has a high spin S=7/2 and no orbital angular momentum, serving as a model system for studying the effects of antiferromagnetic order and strong exchange coupling on the resulting Fermi surface and magnetotransport properties of RE-Vs. Here, we present a surface and structural characterization study mapping the optimal synthesis window of thin epitaxial GdSb films grown on III-V lattice-matched buffer layers via molecular-beam epitaxy. To determine the factors limiting XMR in RE-V thin films and provide a benchmark for band-structure predictions of topological phases of RE-Vs, the electronic band structure of GdSb thin films is studied, comparing carrier densities extracted from magnetotransport, angle-resolved photoemission spectroscopy (ARPES), and density-functional theory (DFT) calculations. ARPES shows a hole-carrier rich, topologically trivial, semimetallic band structure close to complete electron-hole compensation, with quantum confinement effects in the thin films observed through the presence of quantum-well states. DFT-predicted Fermi wave vectors are in excellent agreement with values obtained from quantum oscillations observed in magnetic field-dependent resistivity measurements. An electron-rich Hall coefficient is measured despite the higher hole-carrier density, attributed to the higher electron Hall mobility. The carrier mobilities are limited by surface and interface scattering, resulting in lower magnetoresistance than that measured for bulk crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Buffer concentration dramatically affects the stability of S-nitrosothiols in aqueous solutions

S-nitrosothiols (RSNOs) are an important group of nitric oxide (NO)-donating compounds with low toxicity and wide biomedical applications. In this paper, we, for the first time, demonstrate that the concentration of buffer remarkably affects the stability of RSNOs including naturally occurring S-nitrosoglutathione (GSNO) and synthetic S-nitroso-N-acetylpenicillamine (SNAP). For a solution with a high concentration of GSNO (e.g., 50 mM) and an initial near-neutral pH, the optimal buffer concentration is close to the GSNO concentration under our experimental conditions. A lower buffer concentration does not have adequate buffer capacity to resist the pH drop caused by GSNO decomposition. The decreased solution pH further accelerates GSNO decomposition because GSNO is most stable at near-neutral pH according to our density functional theory (DFT) calculations. A higher-than-optimal buffer concentration also reduces the GSNO stability because buffer ingredients including phosphate, Tris base, and HEPES consume NO/N 2 O 3 . In contrast to GSNO, the highest SNAP stability is obtained when the starting solution at a neutral pH does not contain buffer species, and the stability decreases as the buffer concentration increases. Furthermore, this is because SNAP is more stable at mildly acidic pH and the SNAP decomposition-induced pH drop stabilizes the donor. When the RSNO concentration is low (e.g., 1 mM), the buffer concentration also matters because any excess buffer accelerates the donor decomposition. Since the effect of buffer concentration was previously overlooked and suboptimal buffer concentrations were often used, this paper will aid in the formulation of RSNO solutions to obtain the maximum stability for prolonged storage and sustained NO release.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Revisiting point defect thermodynamics in group IVB and VB transition metal carbides

We present a comprehensive re-examination of point defect thermodynamics in group IVB and VB transition metal carbides (TMCs) with the rocksalt structure using a combination of density functional theory (DFT) calculations and a statistical mechanical Wagner-Schottky model within the canonical ensemble. The most stable configurations of point defects were discovered using basin-hopping global optimization, driven by either a machine learning interatomic potential (MLIP) or DFT. A key finding is the identification of previously unreported dicarbon antisites—a C–C dimer occupying a metal site—as the structural (constitutional) defects on the carbon-rich side of stoichiometry in all group IVB and VB TMCs except TaC. Furthermore, dicarbon antisite-containing thermal defect complexes, such as quadruple and interbranch defects, can dominate in TMCs under specific stoichiometric and temperature conditions. In conclusion, by incorporating dicarbon antisites into the defect landscape, this work provides a revised understanding of the thermodynamics of point defects in TMCs.

Carbides↗
DOE OSTI2021

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗
DOE OSTI2025

Anion Vacancies Coupling with Heterostructures Enable Advanced Aerogel Cathode for Ultrafast Aqueous Zinc‐Ion Storage

As a potential cathode material, manganese-based sulfide has recently attracted increasing interest due to its many advantages in aqueous zinc-ion storage. Unfortunately, some challenges such as sluggish kinetics, unstable structure, and controversial phase transition mechanism during the energy storage process hinder its practical application. Herein, inspired by density functional theory (DFT) calculations, a novel 3D sulfur vacancy-rich and heterostructured MnS/MXene aerogel is designed, and used as a cathode for aqueous Zn-ion batteries/hybrid capacitors (ZIBs/ZICs) for the first time. Thanks to the synergistic modification strategy of sulfur vacancies and heterostructures, the as-constructed MnS/MXene//Zn ZIBs exhibit significantly enhanced electrochemical properties, especially outstanding rate capability and cyclic stability. More encouragingly, the as-assembled MnS/MXene//porous carbon (PC) ZICs exhibit an ultrahigh energy density, a high power density, and a splendid cycling lifespan. Most notably, systematic kinetic analyses, ex situ characterizations, and DFT calculations illustrate that MnS/MXene first irreversibly converts into MnO x @ZnMnO 3 /MXene, and then undergoes a reversible conversion from MnO x @ZnMnO 3 /MXene to MnOOH@ZnMn 2 O 4 /MXene, accompanied by the co-insertion/extraction of H + and Zn 2+ . Further, the synergistic modification strategy of sulfur vacancies and heterostructures and the thorough mechanistic study proposed in this work offer valuable guidance for designing and exploiting high-performance cathodes in aqueous zinc-based energy storage devices.

25 ENERGY STORAGE↗