Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LIQUID METAL COOLED REACTORS”

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.

492 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Anomalous superconductivity in twisted MoTe 2 nanojunctions

Introducing superconductivity in topological materials can lead to innovative electronic phases and device functionalities. Here, we present a unique strategy for quantum engineering of superconducting junctions in moiré materials through direct, on-chip, and fully encapsulated 2D crystal growth. We achieve robust and designable superconductivity in Pd-metalized twisted bilayer molybdenum ditelluride (MoTe 2 ) and observe anomalous superconducting effects in high-quality junctions across ~20 moiré cells. Unexpectedly, the junction develops enhanced, instead of weakened, superconducting behaviors, exhibiting fluctuations to a higher critical magnetic field compared to its adjacent Pd 7 MoTe 2 superconductor. In addition, the critical current further exhibits a notable V-shaped minimum at zero magnetic field. These features are unexpected in conventional Josephson junctions and absent in junctions of natural bilayer MoTe 2 created using the same approach. We discuss implications of these observations, including the possible formation of mixed even- and odd-parity superconductivity at the moiré junctions. Our results also demonstrate a pathway to engineer and investigate superconductivity in fractional Chern insulators.

Science & Technology - Other Topics↗
DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗
DOE OSTI2025

Quantum graph models for transport in filamentary switching

The formation of metallic nanofilaments bridging two electrodes across an insulator is a mechanism for resistive switching. Examples of such phenomena include atomic synapses, which constitute a distinct class of memristive devices the behavior of which is closely tied to the properties of the filament. Until recently, experimental investigation of the low-temperature regime and quantum transport effects has been limited. However, with growing interest in understanding the true impacts of the filament on device conductance, comprehending quantum effects has become crucial for quantum neuromorphic hardware. Here, we discuss quantum transport resulting from filamentary switching in a narrow region where the continuous approximation of the contact is not valid, and only a few atoms are involved. In this scenario, the filament can be represented by a graph depicting the adjacency of atoms and the overlap between atomic orbitals. Using the theory of quantum graphs with locally diffusive node scattering, we calculate the scattering amplitude of charge carriers on this graph and explore the interplay between filamentary formation and quantum transport effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry↗
DOE OSTI2025

Order–Disorder Phase Stabilization by Pressure‐Induced Charge Transfer Enhances the Ferroelectric Photovoltaic Effect in Multiferroic BaFe 4 O 7

Abstract Multiferroic ferroelectric photovoltaic (FPV) materials, combining magnetic and ferroelectric properties, are of paramount importance for optoelectronic and photovoltaic applications. However, optimizing both the remanent polarization and the optical bandgap—key factors for enhanced FPV performance—presents a significant challenge due to their trade‐off. This work shows that pressure‐induced charge transfer between different metal sites can break this trade‐off. Above ≈20 GPa, charge transfer between different trivalent iron (Fe) sites in the multiferroic material BaFe 4 O 7 leads to Fe valence disproportionation, FeO 4 tetrahedra disorder, and Jahn–Teller distortion of FeO 6 octahedra. These changes reduce the bandgap, lower resistivity, and enhance ferroelectric polarization, resulting in a 2.5‐fold increase in photocurrent. Upon decompression, BaFe 4 O 7 retains an order–disorder structure, optimal ferroelectric and optical properties at ambient conditions. This work provides a novel pathway to simultaneously optimizing ferroelectricity and bandgap via pressure‐induced charge transfer, overcoming the traditional trade‐off in FPV materials, and offers a promising approach for developing high polarization performance, narrow‐bandgap FPV materials.

Chemistry↗
DOE OSTI2025

In‐Plane Anisotropy in van der Waals NiTeSe Ternary Alloy

Abstract The anisotropic properties of materials profoundly influence their electronic, magnetic, optical, and mechanical behaviors and are critical for a wide range of applications. In this study, the anisotropic characteristics of Ni‐based van der Waals materials, specifically NiTe 2 and its alloy NiTeSe, utilizing a combination of comprehensive scanning tunneling microscopy (STM), angle‐resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations, are explored. Unlike 1T‐NiTe 2 , which exhibits trigonal in‐plane symmetry, the substitution of Te with Se in NiTe 2 (resulting in the NiTeSe alloy) induces a pronounced in‐plane anisotropy. This anisotropy is clear in the STM topographs, which reveal a distinct linear order of charge distribution. Corroborating these observations, ARPES measurements and DFT calculations reveal an anisotropic Fermi surface centered at the point, which is notably elongated along the k y direction, leading to directional variations in in‐plane carrier velocities. Consequently, the Fermi velocity is highest along the k x direction where the linear charge distribution aligns in real space and is lowest along the k y direction. These findings offer valuable insights into the tunability of anisotropic properties in ternary transition metal dichalcogenide systems, highlighting their potential applications in the development of anisotropic electronic and optoelectronic devices.

Chemistry↗
DOE OSTI2025

Effect of Entropic Constraints on the Thermodynamics of Molecular Adsorption in Nano‐Porous Materials

Abstract Gas separation is a critical industrial process that consumes a significant amount of energy due to the widely used techniques that are currently employed. Adsorptive materials—such as metal–organic frameworks (MOFs)—show promise as an energy‐efficient alternative. Of particular current interest are novel, temperature‐dependent separation processes in MOFs, such as the recently reported separation of ternary isomeric hydrocarbon mixtures within one and the same material. However, the mechanisms of these highly desirable separations remain poorly understood. Herein, through a combination of ab initio simulations and statistical mechanics, it is shown that the temperature dependence is the result of a constraint on the guest molecule's entropic degrees of freedom when loaded into the MOF, caused by the fortuitous tight fitting of the guest inside the pore. While the framework applies to all molecular adsorption in porous media, it is essential for the description of large molecules in small pores, which is demonstrated here using the separation of C6 isomers in Ca(H 2 tcpb) as a test case. The developed framework and analysis not only reveal the reason why separation occurs but also predict the temperatures at which it takes place, thus opening the door to newly designed MOFs with tailor‐made precision.

Chemistry↗
DOE OSTI2026

Power and particle exhaust in the ST-E1 fusion power plant

Power exhaust challenges and potential solutions for a 5 m major radius, low-aspect ratio burning tokamak have been explored. 1D edge plasma models have been used to screen for access to detachment using short and long outer divertor legs in double and single null configurations, using Ar as the primary impurity and assuming tungsten plasma-facing components (PFCs). These show that detachment access can be accessed for all but the most conservative assumptions on scrape-off layer (SOL) width and power, but that trade-offs will be required between magnet engineering and the size of the acceptable window of as-yet uncertain plasma parameters. SOLPS-ITER was used to further model selected plasma scenarios, confirming that Ar seeding can be used to achieve dissipative divertor scenarios with peak deposited heat fluxes below 15 MWm -2 . Initial scoping of first wall loads and positioning of limiters has been carried out, showing the feasibility of protecting the breeding blanket wall during steady state without impeding tritium breeding. Initial PFC technology selection is also presented, identifying this as a critical area where further work is needed to find an attractive solution for helium-cooled PFCs that can handle high heat fluxes without excessive power requirements. Key questions and trade-offs for concept development have been identified, including: how to achieve high radiation for reduction of SOL power without core performance degradation; whether power exhaust can be well-controlled in a double null plasma; mechanical design and materials challenges of high-heat flux PFCs; and control of material erosion, redeposition and tritium retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Massive red spiral galaxies in SDSS-IV MaNGA survey

ABSTRACT Massive red spiral galaxies (MRSGs) are supposed to be the possible progenitors of lenticular galaxies (S0s). We select a large sample of MRSGs ($M_*\gt 10^{10.5}\rm {\rm M}_{\odot }$) from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) DR17 using the g − r colour versus stellar mass diagram, along with control samples of blue spirals and S0s. Our main results are as follows: (1) After comparing the Sérsic index, concentration parameter, asymmetry parameter distribution, size–mass relation, and Σ1 (stellar mass surface density within the central 1 kpc)−mass relation, we find MRSGs are similar to S0s and have more compact and symmetric structures than blue spirals. MRSGs also resemble S0s in Dn4000, metallicity, Mgb/$\rm \left\langle Fe \right\rangle$, and V/σ radial profile. (2) By using MaNGA 2D spectra data, we separate the spatial regions into inner (R < 0.8Re) and outer (0.8 < R < 1.5Re) regions, and detect residual star formation in the outer regions of MRSGs. (3) When we select a sub-sample of MRSGs with NUV − r > 5, we find that they are completely star formation quenched in both inner and outer regions. Compared to optically selected MRSGs, NUV − r selected MRSGs appear to be more concentrated and have more massive dark matter haloes. The similarities between S0s and MRSGs suggest the possible evolutionary trend between MRSGs and S0s.

Astronomy & Astrophysics↗
DOE OSTI2026

Methanol mixing−controlled compression ignition enabled via homogeneous charge compression ignition of dimethyl ether through catalytic decomposition of methanol

Methanol is a potentially attractive fuel for marine and off−road engines owing to its availability at bunkering and global distribution locations. Although methanol is well−distributed worldwide, its fuel chemistry and ignition properties make it poorly suited as a direct drop−in replacement for diesel fuel in compression−ignition engines. However, industrial processes are regularly used to convert methanol, via catalytic dehydration, to dimethyl ether (DME) over nonprecious metal catalysts. This chemical conversion can occur at relatively low pressures, temperatures, and catalyst space velocities, highlighting a potential opportunity to generate DME via onboard catalytic dehydration of methanol. DME’s fuel kinetic and ignition properties for compression ignition are much more favorable than those of methanol or even diesel fuel, but DME is more challenging than diesel fuel or methanol to pump, store, and deliver through conventional diesel fueling injection hardware. Thus, a potential opportunity exists to use the ignition and kinetic properties of DME, with the transportation and delivery advantages of methanol, in a methanol−fueled mixing−controlled compression−ignition engine. The present work explores performance, combustion behavior, and emissions reduction opportunities for methanol mixing−controlled combustion, enabled by a HCCI of DME that represents a small fraction of the total fuel energy that can be generated onboard via catalytic dehydration of methanol.

Splitter, Derek [ORNL] (ORCID:0000000174044047)↗
DOE OSTI2026

DNA Crystals as a Template for Patterned Functional Materials

DNA nanotechnology offers a wide toolkit of molecular functionalities and scales, including intricate motifs less than 10 nm and periodic structures exceeding 100 µm. At larger scales, however, there are often significant tradeoffs for DNA structures, namely stability and mechanical strength. This work describes the design, synthesis, and characterization of a functionalized DNA crystal. Using a ligated DNA crystal grants significant freedom for various functional materials to be applied, in this case, semiconducting cadmium sulfide and palladium metal. Properties investigated in this study include stability, mechanical strength, and optoelectronic properties such as photoluminescence (PL) and electric conductivity. Significant changes are observed based on the functional material applied to DNA crystals. The Young's modulus of the crystal varies by about five orders of magnitude when functionalized with palladium. PL and semiconductive behaviors were observed when cadmium sulfide was attached. These crystals represent an expansion of the capabilities of DNA structures at these length scales, and additionally a platform for future studies exchanging the materials or altering the ligation scheme.

CdS↗
DOE OSTI2026

Development of in-situ polymerized intrinsically conductive resin and low-cost carbon pigments offering high conductivity for sensing, EMI shielding and lighting protection

Electrically conductive composites are emerging across diverse industries such as electronic, automotive, aerospace, advanced air mobility, biomedical, infrastructure, defense and security offering static charge dissipation, electromagnetic interference shielding, lighting protection, sensing, dicing, corrosion monitoring, etc. Conductivity enhanced composites provide several advantages compared to conventional metals including weight reduction, corrosion resistance, energy efficient processability, tunable properties and multifunctionality. Polymers are typically insulating in nature and require conducting filler for electron transport. However, dispersion and polymer-filler interphases are critical and often disrupt conducting pathways. Besides, conductive fillers such as graphene, carbon nanotube, MXene, silver nanowire, etc. are expensive, limiting their wide adoption in composite industry. On the other hand, a limited number of intrinsically conductive polymers are available among which polyaniline (PANI) has been widely studied due to its high conductivity, thermal and chemical stability. However, PANI is difficult to process and exhibits weak mechanical properties. In brief, there is a significant demand for electrically conductive polymer formulation with cost-effective conducting fillers that offer processability in scale to expand the market of conductivity enhanced materials.

Kumar, Vipin [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

Wetter, Michael [Lawrence Berkeley National Labora↗
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2026

Single-crystal NMC622 from combustion-assisted solid-state synthesis: Tunable particle size, voltage-dependent stability, and Li-ion diffusion

Cost-effective and scalable synthesis of single-crystal lithium–metal-oxide cathodes with controlled size and morphology remains challenging for advancing high-performance lithium-ion batteries. Here, we demonstrate a solution combustion–assisted solid-state method for producing single-crystal NMC622 with tunable particle size at various lithium stoichiometries. Structural and electronic characteristics were probed using synchrotron-XANES, EXAFS, XRD and XPS, confirming well-defined coordination environments and phase purity up to Li ≈ 1.3. The electrochemical behavior of LixNMC (x = 1.0 and 1.3) single crystals was evaluated across 2.8–4.3 V and 2.8–4.7 V windows to elucidate the effects of lithium content and operating voltage. At 2.8–4.3V, Li1.3NMC exhibits higher initial capacity due to its larger lithium inventory, while Li1.0NMC shows superior long-term retention driven by its larger crystal size and reduced structural distortion. Increasing the cutoff voltage to 4.7 V enhances the initial capacity of both compositions by 17–25%, but long-term capacity retention ultimately converges to values comparable to those at 4.3 V due to voltage-induced degradation, where Li+ diffusion constant remains lying ranging from 10−12 –10−11 cm2s−1. Overall, this study establishes a scalable combustion-assisted route for synthesizing tunable single-crystal NMC622 and highlights the interplay between lithium stoichiometry, particle size, and voltage window in governing cathode stability.

Roy, Subrata C [Jackson State University]↗
DOE OSTI2024

Real-space visualization of a defect-mediated charge density wave transition

Here, we study the coupled charge density wave (CDW) and insulator-to-metal transitions in the 2D quantum material 1T-TaS 2 . By applying in situ cryogenic 4D scanning transmission electron microscopy with in situ electrical resistance measurements, we directly visualize the CDW transition and establish that the transition is mediated by basal dislocations (stacking solitons). We find that dislocations can both nucleate and pin the transition and locally alter the transition temperature T c by nearly ~75 K. This finding was enabled by the application of unsupervised machine learning to cluster five-dimensional, terabyte scale datasets, which demonstrate a one-to-one correlation between resistance—a global property—and local CDW domain-dislocation dynamics, thereby linking the material microstructure to device properties. This work represents a major step toward defect-engineering of quantum materials, which will become increasingly important as we aim to utilize such materials in real devices.

4D-STEM↗
DOE OSTI2026

Computational discovery of a tetragonal MBene phase: diamond-shaped boron lattices, viable MAB precursors, and selective CO 2 reduction

MBenes are a class of two-dimensional transition-metal borides derived from layered MAB phases. Here, we use first-principles calculations to predict a new tetragonal MB (1 : 1) MBene phase featuring a diamond-shaped B–B lattice and systematically assess its phase stability, synthetic accessibility, and catalytic performance towards the CO 2 reduction reaction (CO 2 RR). Four tetragonal configurations are benchmarked against known hexagonal and orthorhombic phases, and a multi-tier screening identifies six robust members (CrB, FeB, MoB, WB, IrB, and PtB). Exploration of the I4/mmm MAB space reveals 11 viable precursors, including the experimentally synthesized Ir 2 ZnB 2 , which validates our screening approach. CO 2 RR free-energy diagrams show that CrB favors CH 3 OH while FeB, MoB, WB, IrB, and PtB preferentially yield HCOOH. The corresponding limiting potentials are competitive with, and in some cases superior to, those reported for state-of-the-art orthorhombic and hexagonal MBenes. These findings expand the MBene landscape by establishing a tetragonal phase with structural, synthetic, and catalytic promise.

Lu, Linguo [University of Puerto Rico, San Juan, P↗
DOE OSTI2023

Challenges and Prospects of Sodium‐Ion and Potassium‐Ion Batteries for Mass Production

The exponential growth of the lithium‐ion (LIB) market is causing a significant disparity between the supply chain and demand for its resources. In this regard, sodium‐ion and potassium‐ion batteries are promising alternatives to LIBs due to their low cost. However, the larger sizes of Na + and K + ions create challenges that prevent them from achieving energy densities comparable to LIBs while maintaining an acceptable cycle life. Here, in this perspective, the aim is to evaluate the status of Na‐ion and K‐ion batteries and the challenges associated with them on both fundamental and commercial levels. The focus is on the structural instability arising from phase transitions during cycling, intricate chemical degradation processes, and potential avenues for enhancing their performance with a specific goal of improving their viability for grid‐scale energy storage. Materials production and abundance limitations for the chemistries of the state‐of‐the‐art materials and account for critical parameters from both the perspective of researchers and investors are analyzed. This analysis aims to provide insights into the strategic trade‐offs required to effectively implement the technology in real‐world applications, such as grid‐scale storage and other areas. Furthermore, the utilization of metals with low or no supply‐chain problems as an important aspect of these trade‐offs is considered.

25 ENERGY STORAGE↗