Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface instabilities”

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 OSTI2021

Mode spectrum characteristics and onset of the low-shear MHD stability regime

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Growth of concomitant laser-driven collisionless and resistive electron filamentation instabilities over large spatiotemporal scales

Collective processes in plasmas often induce microinstabilities that play an important role in many space or laboratory plasma environments. Particularly of note is the Weibel-type current filamentation instability, which is believed to drive the creation of collisionless shocks in weakly magnetized astrophysical plasmas. In this paper, this instability class is studied through interactions of ultraintense and short laser pulses with solid foils, leading to localized generation of megaelectronvolt electrons. Proton radiographic measurements of both low- and high-resistivity targets show two distinct, superimposed electromagnetic field patterns arising from the interpenetration of the megaelectronvolt electrons and the background plasma. Particle-in-cell simulations and theoretical estimates suggest that the collisionless Weibel instability building up in the dilute expanding plasmas formed at the target surfaces causes the observed azimuthally symmetric electromagnetic filaments. For a sufficiently high resistivity of the target foil, an additional resistive instability is triggered in the bulk target, giving rise to radially elongated filaments. The data reveal the growth of both filamentation instabilities over large temporal (tens of picoseconds) and spatial (hundreds of micrometres) scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tailoring Interfaces in Solid-State Batteries Using Interfacial Thermochemistry and Band Alignment

Solid-state lithium ion batteries have enhanced thermal stability compared to batteries using liquid electrolytes. Although their practical implementation is limited by undesired side reactions between the electrode and electrolyte, different modification strategies have been proposed to stabilize these reactive interfaces. These approaches have been primarily based upon bulk materials properties, however, which may not necessarily be representative of the chemistry at the solid/solid interface. In this work, first principles calculations and X-ray scattering experiments are used to elucidate molecular-level reactivity and develop design principles for tailored interfaces based on surface and interfacial properties. Because of its well-known instability, the interface between the Li metal anode and the lithium lanthanum titanate (LLTO) solid electrolyte is applied as a model system. Density functional theory calculations are used to describe bulk, surface, and interfacial thermochemistry of the Li/LLTO system, and interfacial reconstruction is probed using ab initio molecular dynamics. These simulations of the Li/LLTO interface demonstrate facile decomposition concomitant with reduction of Ti 4+ ions. Based on further insights from computational analysis of the surface band edge positions, La 2 O 3 is proposed as an interlayer coating and is shown to provide an energetic barrier for interfacial charge transfer and decomposition reactions from X-ray reflectivity measurements and theoretical calculations. The findings in this work suggest that design strategies based on surface and interfacial properties can be used to kinetically stabilize electrode/electrolyte interfaces in solid-state batteries.

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

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Stemflow Hydrodynamics

Stemflow hydrodynamics is the study of water movement along the exterior surface area of plants. Its primary goal is to describe water velocity and water depth along the stem surface area. Its significance in enriching the rhizosphere with water and nutrients is not in dispute. Yet, the hydrodynamics of stemflow have been entirely overlooked. This review seeks to fill this knowledge gap by drawing from thin film theories to seek outcomes at the tree scale. The depth‐averaged conservation equations of water and solute mass are derived at a point. These equations are then supplemented with the conservation of momentum that is required to describe water velocities or relations between water velocities and water depth. Relevant forces pertinent to momentum conservation are covered and include body forces (gravitational effects), surface forces (wall friction), line forces (surface tension), and inertial effects. The inclusion of surface tension opens new vistas into the richness and complexity of stemflow hydrodynamics. Flow instabilities such as fingering, pinching of water columns into droplets, accumulation of water within fissures due to surface tension and their sudden release are prime examples that link observed spatial patterns of stemflow fronts and morphological characteristics of the bark. Aggregating these effects at the tree‐ and storm‐ scales are featured using published experiments. The review discusses outstanding challenges pertaining to stemflow hydrodynamics, the use of dynamic similarity and 3D printing to enable the interplay between field studies and controlled laboratory experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

On the Relationship Between the Marine Cold Air Outbreak M Parameter and Low-Level Cloud Heights in the Midlatitudes

Focusing on conditions of subsidence when low clouds are present, ground-based observations in both the North Atlantic and the Southern Ocean reveal strong relationships between cloud boundary (base and top heights) and different measures of lower tropospheric instability. The difference in potential temperature between the surface and 800 hPa (a metric called M) provides a stronger relationship than measures of inversion strength such as the lower tropospheric stability and estimated inversion strength. This is because (1) inversion strength itself does not correlate well with cloud boundaries, and (2) M contains information that appears important for cloud boundaries. These include the surface forcing through the use of sea surface rather than near-surface air temperature and an upper level close to the real cloud top. These results expand upon previous work on the importance of M as a predictor of cloud morphology. However, important differences are found in low-cloud conditions for the North Atlantic as compared to the Southern Ocean (for a given value of M): stronger inversions, deeper boundary layers, and much larger sea level pressures. Therefore, the relationship between cloud boundaries and M differs between the two regions. A general circulation model provides similar relationships as observed between M and both cloud top height and temperature but tends to place clouds higher and at colder temperatures than observed for a given M. This might cause issues with the representation of precipitation, cloud cover and radiation in the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
DOE OSTI2021

A Constitutive Structural Parameter for the Work Hardening Behavior of Additively Manufactured Ti-6Al-4V

The mechanical behavior of Ti-6Al-4V produced by additive manufacturing processes is assessed based on a formulation developed from the Kocks-Mecking relationship. A constitutive parameter derived for the microstructure is characteristic of the work hardening behavior determined by the plastic strain between the proportional limit and the strength at the instability point. The varied plasticity behavior associated with surface and build direction effects can be evaluated with this approach as presented for the case of Ti-6Al-4V under uniaxial tension.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Post-Shot Report for OMEGA Double Cylinders (CylDRT 22B)

The direct-drive double cylinder experimental platform is a high-energy-density (HED) science platform designed to image an imploding cylindrical target. The target consists of a directly-driven outer cylinder and a shock-and-collision driven inner cylinder. The purpose of this platform is to study hydrodynamic instability growth on the inner cylinder, the outer surface of which is classically Rayleigh-Taylor unstable during the acceleration phase. We present results from recent experiments at the OMEGA laser facility. These experiments were designed as a proof-of-principle for the platform, using the same cylinder exterior dimensions and direct-drive beam configuration as previous single cylinder experiments. In these experiments, three sets of targets were fielded: no machined perturbations (smooth), a sinusoidal mode-10 perturbation on the outer surface of the inner cylinder, and a mode-20 perturbation on the outer surface of the inner cylinder. The primary diagnostic was a gated x-ray framing camera which imaged the backlit inner cylinder on-axis for sixteen frames over a time window of 1 ns for each shot. A second side-lit radiograph captured one image per shot, diagnosing axial uniformity. In this report we include an overview of the results from both the backlighter and the sidelighter diagnostics. We discuss at length the experimental analysis process. Finally, we present the results of the smooth target implosion trajectory and compare them to post-shot simulations. We see favorable agreement between simulation and experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of surface roughness on phase transition timing in megaampere pulsed-power–driven exploding conductors

An understanding of material phase transitions in megaampere pulsed-power–driven exploding conductors is important for predicting the growth of hydrodynamic instabilities in magneto-inertial fusion concepts. Here, this study analyzes phase transitions in electrical conductor explosions using 1D Lagrangian and 2D arbitrary Lagrangian–Eulerian resistive magnetohydrodynamic simulations to show that micrometer-scale surface roughness can lead to the electrothermal instability (ETI), a feedback effect that concentrates resistive heating and leads to early melting and ablation. Simulations of the Mykonos electrothermal instability II (METI-II) experiment show melting begins 19% sooner for machined rods with micrometer-scale surface roughness than for rods without these features. The surface magnetic field is 41 T around the initial region of melt, representing a lower magnitude than both the 86 T from 1D simulations and the 85 T threshold reported elsewhere. In 2D simulations with micrometer-scale surface roughness, temperature measurements indicate the critical point temperature of aluminum is reached 17% faster in comparison with 1D simulations. Values from 2D simulations with surface roughness align with predictions from ETI theory, and the observed temperature redistribution further supports the ETI as an underlying mechanism. Simulation results are validated against experimental photonic Doppler velocimetry data. This study shows 1D simulations are adequate to model conductors with sub-micrometer-scale surface roughness in this high-energy-density regime; however, 2D or 3D simulations are required to capture the full range of physics for accurately describing phase transitions in conductors with micrometer-scale or larger surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Defect-induced electronic smectic state at the surface of nematic materials

Due to the intertwining between electronic nematic and elastic degrees of freedom, lattice defects and structural inhomogeneities commonly found in crystals can have a significant impact on the electronic properties of nematic materials. Here, we show that defects commonly present at the surface of crystals generally shift the wave-vector of the nematic instability to a non-zero value, resulting in an incommensurate electronic smectic phase. Such a smectic state onsets above the bulk nematic transition temperature and is localized near the surface of the sample. We argue that this effect may explain not only recent observations of a modulated nematic phase in iron-based superconductors, but also several previous puzzling experiments that reported signatures consistent with nematic order before the onset of a bulk structural distortion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Laboratory Study of the Torus Instability Threshold in Solar-relevant, Line-tied Magnetic Flux Ropes

Coronal mass ejections (CMEs) occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. This destabilization is often described in terms of an ideal magnetohydrodynamic instability called the torus instability. It occurs when the external magnetic field decreases sufficiently fast such that its decay index, ${n}_{}=-z\,\partial (\mathrm{ln}{B}_{})/\partial z$, is larger than a critical value, $n\gt {n}_{\mathrm{cr}}^{}$, where ${n}_{\mathrm{cr}}^{}=1.5$ for a full, large aspect ratio torus. However, when this is applied to solar MFRs, a range of conflicting values for ${n}_{\mathrm{cr}}^{}$ is found in the literature. To investigate this discrepancy, we have conducted laboratory experiments on arched, line-tied flux ropes and applied a theoretical model of the torus instability. Our model describes an MFR as a partial torus with foot points anchored in a conducting surface and numerically calculates various magnetic forces on it. This calculation yields better predictions of ${n}_{\mathrm{cr}}^{}$ that take into account the specific parameters of the MFR. Here, we describe a systematic methodology to properly translate laboratory results to their solar counterparts, provided that the MFRs have a sufficiently small edge safety factor or, equivalently, a large enough twist. After this translation, our model predicts that ${n}_{\mathrm{cr}}^{}$ in solar conditions falls near ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim 0.9$ and within a larger range of ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim (0.7,1.2)$, depending on the parameters. The methodology of translating laboratory MFRs to their solar counterparts enables quantitative investigations of CME initiation through laboratory experiments. These experiments allow for new physics insights that are required for better predictions of space weather events but are difficult to obtain otherwise.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bidentate Lewis Base Ligand-Mediated Surface Stabilization and Modulation of the Electronic Structure of CsPbBr 3 Perovskite Nanocrystals

The desorption of conventional ligands from the surface of halide perovskite nanocrystals (NCs) often causes their structural instability and the deterioration of the optoelectronic properties. To address this challenge, we present an approach of using a bidentate Lewis base ligand, namely, 1,4-bis(diphenylphosphino)butane (DBPP), for the synthesis of CsPbBr 3 NCs. The phosphine group of DBPP has a strong interaction with PbBr 2 precursor, forming a highly crystalline intermediate complex during the reaction. In the presence of oleic acid, the uncoordinated phosphine group of DBPP is converted into the phosphonium cation, which strongly binds with the surface bromide of the formed CsPbBr 3 NCs through hydrogen bonding. Density functional theory calculations suggest that DBPP can strongly bind to the undercoordinated lead and surface bromide ions of CsPbBr 3 NCs through its unprotonated and protonated phosphine groups, respectively. The robust binding of DBPP with the surface of perovskite NCs helps to preserve their structural integrity under various environmental stresses. Moreover, the electron density and energy levels are regulated in DBPP-capped CsPbBr 3 NCs via the efficient transfer of electrons from the ligands to the NCs, resulting in their improved photocatalytic CO 2 reduction performance. Furthermore, our study highlights the potential of using bidentate ligands to stabilize the surface of perovskite NCs and modulate their optical and electronic properties.

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

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Creation of Multi-Principal Element Alloy NiCoCr Nanostructures via Nanosecond Laser-Induced Dewetting

The multi-principal element alloy nanoparticles (MPEA NPs), a new class of nanomaterials, present a highly rewarding opportunity to explore new or vastly different functional properties than the traditional mono/bi/multimetallic nanostructures due to their unique characteristics of atomic-level homogeneous mixing of constituent elements in the nanoconfinements. Here, in this paper, the successful creation of NiCoCr nanoparticles, a well-known MPEA system is reported, using ultrafast nanosecond laser-induced dewetting of alloy thin films. Nanoparticle formation occurs by spontaneously breaking the energetically unstable thin films in a melt state under laser-induced hydrodynamic instability and subsequently accumulating in a droplet shape via surface energy minimization. While NiCoCr alloy shows a stark contrast in physical properties compared to individual metallic constituents, i.e., Ni, Co, and Cr, yet the transient nature of the laser-driven process facilitates a homogeneous distribution of the constituents (Ni, Co, and Cr) in the nanoparticles. Using high-resolution chemical analysis and scanning nanodiffraction, the environmental stability and grain arrangement in the nanoparticles are further investigated. Thermal transport simulations reveal that the ultrashort (≈100 ns) melt-state lifetime of NiCoCr during the dewetting event helps retain the constituent elements in a single-phase solid solution with homogenous distribution and opens the pathway to create the unique MPEA nanoparticles with laser-induced dewetting process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The passivity of lithium electrodes in liquid electrolytes for secondary batteries

Rechargeable Li metal batteries are currently limited by electrolyte decomposition and rapid Li consumption. Li plating and stripping greatly depend on the solid electrolyte interphase formed at the Li metal-liquid electrolyte interface. This Review discusses the reactions occurring at this interface from a corrosion science perspective, highlighting the requirements for an ideal passivation layer. Rechargeable Li metal batteries are currently limited by safety concerns, continuous electrolyte decomposition and rapid consumption of Li. These issues are mainly related to reactions occurring at the Li metal-liquid electrolyte interface. Additionally, the formation of a passivation film (that is, a solid electrolyte interphase) determines ionic diffusion and the structural and morphological evolution of the Li metal electrode upon cycling. In this Review, we discuss spontaneous and operation-induced reactions at the Li metal-electrolyte interface from a corrosion science perspective. We highlight that the instantaneous formation of a thin protective film of corrosion products at the Li surface, which acts as a barrier to further chemical reactions with the electrolyte, precedes film reformation, which occurs during subsequent electrochemical stripping and plating of Li during battery operation. Finally, we discuss solutions to overcoming remaining challenges of Li metal batteries related to Li surface science, electrolyte chemistry, cell engineering and the intrinsic instability of the Li metal-electrolyte interface.

25 ENERGY STORAGE↗