Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time-resolved”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ultrafast Raman probe of the photoinduced superconducting to normal state transition in the cuprate Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝛿

Here, we report an ultrafast Time-Resolved Raman scattering study of the out-of-equilibrium photoinduced dynamics across the superconducting to normal state phase transition of the cuprate Bi 2 ⁢Sr 2⁢ CaCu 2 ⁢O 8+𝛿 . Using the polarization-resolved momentum space selectivity of Raman scattering, we track the superconducting condensate destruction and recovery dynamics with subpicoseconds time resolution in the antinodal region of the Fermi surface where the superconducting gap is maximum. Leveraging ultrafast Raman thermometry, we find a significant dichotomy between the superconducting condensate and the quasiparticle temperature dynamics near the antinodes, which cannot be framed in terms of a single effective electron temperature. The present work demonstrates the ability of Time-Resolved Raman scattering to selectively probe out-of-equilibrium pathways of different electronic subdegrees of freedom during a photoinduced phase transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

The chemRIXS Instrument for Solution-Phase Resonant Inelastic X-ray Scattering at the Linac Coherent Light Source

The chemRIXS instrument at the Linac Coherent Light Source (LCLS) offers new opportunities for studying solution-phase systems with time-resolved soft X-ray spectroscopic methods such as resonant inelastic X-ray scattering (RIXS) using the recently commissioned high-repetition-rate LCLS-II X-ray free electron laser. The orders-of-magnitude X-ray flux improvement provided by the superconducting accelerator, combined with corresponding advances in the optical laser system and the liquid jet recirculation system, enables studies on dilute systems with high signal-to-noise compared to what was possible with the LCLS-I copper accelerator. These capabilities open up time-resolved X-ray absorption spectroscopy and RIXS to entirely new classes of samples, as well as enabling the development of new soft X-ray spectroscopic methods on liquid samples. An overview of the beamline components and the first LCLS-II commissioning results are presented.

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

Electron-Deficient Phenazines: Synthesis, Structures, Redox, and Optoelectronic Properties in the Gas Phase, in Solution, and in the Solid State

High-temperature gas-phase reactions of 1,4-C4F8I2 with phenazine (PHNZ) afforded new thermally- and air-stable derivatives with cyclo-C4F8 substituents, viz. 1,2-PHNZ(C4F8), 2,3-PHNZ(C4F8), 1,2:6,7-PHNZ(C4F8)2, 1,2:7,8-PHNZ(C4F8)2, and 1,2:8,9-PHNZ(C4F8)2. Reactions in the presence of Cu powder resulted in reductive-defluorination/aromatization of some of the cyclo-C4F8 substituents to produce 2,3-PHNZ(C4F4), 2,3-PHNZ(C4HF3), and 1,2:6,7-PHNZ(C4F8)(C4F4), which are formally tri- or tetrafluoro derivatives of benzo[a]phenazine. Single-crystal X-ray diffraction structures of these compounds demonstrate the planarity of 2,3-PHNZ(C4F4) and 2,3-PHNZ(C4HF3) and ordered p-stacking in 1,2-PHNZ(C4HF3), 1,2-PHNZ(C4F4), 1,2-PHNZ(C4F8), and 2,3-PHNZ(C4F8). Low-temperature gas-phase photoelectron spectra of the PHNZ(C4F8)1,2 compounds show that they are strong electron acceptors, with the three PHNZ(C4F8)2 isomers having electron affinities exceeding that of C60. Cyclic voltammetry in acetonitrile show that the five PHNZ(C4F8)1,2 compounds exhibit reversible one-electron reductions with large positive shifts in their reduction potentials relative to unsubstituted PHNZ. Spectroelectrochemical absorption and EPR spectra of PHNZ(C4F8)1,2- • anion radicals and extensive DFT calculations revealed the significant role that the different substitution patterns have on determining the optoelectronic properties of the PHNZ(C4F8)1,2 derivatives and PHNZ(C4F8)1,2- • anion radicals. Time-resolved microwave conductivity measurements and photoluminescence spectra of blends of the PHNZ(C4F8)1,2 derivatives with the donor poly(3-hexylthiophene) (P3HT) were recorded to assess the possible use of the new compounds in optoelectronic devices.

Balser, Sebastian↗
DOE OSTI2026

Kinetic control of phase evolution and defect-mediated recombination in AACVD-grown copper antimony sulphide thin films

Copper antimony sulphide (CAS) is a multinary chalcogenide semiconductor in which small deviations from stoichiometry can drive phase competition and strong defect-mediated modulation of optoelectronic properties. However, the systematic roles of copper precursor fraction and deposition time in governing phase evolution, off-stoichiometry, and recombination dynamics in aerosol-assisted chemical vapour deposition (AACVD)-grown undoped CAS thin films remain insufficiently understood. In this work, undoped CAS thin films were deposited by AACVD using Cu(dedtc)2 and Sb(dedtc)3 single-source precursors at 550 °C and a carrier gas flow rate of 150 sccm, while the Cu(dedtc)2 mole fraction (x = 0.15 – 0.55) and deposition time (1 – 2 hours) were systematically varied to probe how growth kinetics influence phase composition, microstructure, and defect-mediated optical properties without post-deposition annealing or extrinsic doping. Increasing copper precursor content drives phase evolution toward tetrahedrite-dominant CAS films at intermediate Cu(dedtc)2 mole fractions, with the film deposited at x = 0.35 exhibiting the strongest tetrahedrite character within the parameter space examined. The films are also copper-rich, antimony-poor, and sulphur-deficient, consistent with off-stoichiometric growth and intrinsic defect formation, plausibly including copper interstitials, Cu-on-Sb antisites, and sulphur vacancies. These growth-dependent compositional deviations are accompanied by tunable indirect optical bandgaps of approximately 1.60 – 2.18 eV and weak visible photoluminescence governed by defect-mediated recombination. Time-resolved photoluminescence reveals bi-exponential decay behaviour with lifetimes of approximately 0.1 – 3.6 ns, with emission dominated by slower donor–acceptor pair recombination and a smaller contribution from faster trap-assisted pathways. Collectively, these results establish an explicit kinetic process–structure–defect–property relationship for AACVD-grown CAS thin films and provide a growth–structure–property framework relevant to future optimization of CAS-based optoelectronic and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2027

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2026

Reconstruction and Dissolution of Copper Catalysts during Electrocatalytic Nitrate Reduction

Copper-based electrocatalysts are widely explored for electrochemical nitrate remediation, yet their stability under operating conditions remains poorly understood. While nanoscale and subnanoscale Cu motifs are known to restructure during the nitrate reduction reaction (NO 3 RR), how these transformations translate into irreversible material loss remains unclear. Here, we quantify Cu dissolution during NO 3 RR as a function of catalyst architecture and electrolyte chemistry using two model systems: single-atom Cu (Cu 1 ) and Cu nanoparticles (Cu NP ). Time-resolved leaching measurements, integrated with in situ X-ray absorption spectroscopy (XAS), reveal measurable Cu loss for both catalysts during NO 3 RR. Dissolution is concentrated at the initiation of electrolysis, coinciding with rapid restructuring, and depends strongly on morphology, with Cu NP consistently exhibiting greater Cu loss than Cu 1 . In situ XAS reveals that Cu 1 forms transient metallic clusters under reduction that largely redisperse upon returning to open-circuit voltage, whereas Cu NP undergoes reduction of an oxidized surface layer accompanied by sustained Cu loss during electrolysis. Notably, the presence of nitrate significantly intensifies restructuring and Cu loss, highlighting the critical role of electrolyte composition. Furthermore, these findings establish a direct link between electrochemical restructuring and Cu dissolution, unveiling electrolyte-dependent interactions as key determinants of catalyst durability in electrochemical nitrate conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Unimolecular Dynamics of Partially Deuterated Hydroperoxyalkyl Intermediates (•QOOD) in Cyclohexane and Cyclopentane Oxidation

The hydroperoxyalkyl radicals (•QOOH) formed in cyclohexane and cyclopentane oxidation with carbon radical center at the β site relative to the hydroperoxy group have recently been observed through their infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and cyclic ether products. Partial deuteration shifts the IR transitions associated with OH/OD motion to lower frequencies, providing a new window on the IR spectroscopy and unimolecular dynamics of these transient intermediates. Specifically, the overtone OD stretch (2ν OD ) has been identified at 5215.0 cm –1 and 5216.0 cm –1 for the β-QOOD intermediates in cyclohexane and cyclopentane oxidation, respectively. Unimolecular decay rates of the •QOOD intermediates at these energies have been obtained from the time-resolved appearance of OD products. The experimental rates are compared with statistical microcanonical rates evaluated using RRKM theory, including heavy-atom tunneling associated with simultaneous O–O bond elongation and C–O–O angle contraction along the reaction pathway. The experimental rates provide a further test of the previously determined transition state (TS) barriers, which were computed utilizing a benchmarking approach that builds on higher-level reference calculations for the smaller ethane oxidation system. Here, the experimental rate measurements for both •QOOD intermediates agree well with the computed rates after accounting for small changes in zero-point energy and a minor empirical adjustment of the TS barriers in both systems.

Ethers↗
DOE OSTI2026

Large Polaron Screening Delays Mobility Roll-Off in Halide Perovskites under High Injection

The interplay between carrier density regimes and the various scattering processes is fundamental for understanding charge transport efficiency in high-performance semiconductors. Combining ultrafast time-resolved terahertz spectroscopy (TRTS) and transient microwave photoconductivity (TRMC), we measure carrier mobilities over 7 orders of magnitude (10 13 –10 20 carriers cm –3 ), reaching a regime where carrier–carrier scattering dominates and mobility begins to roll off. We find the Caughey–Thomas model can be effectively employed for transport modeling in 3D metal halide perovskites (MHPs), and the mobility roll-off near ∼10 19 cm –3 can be tuned by MHP composition. Carrier lifetime in MHPs starts decreasing at much lower carrier densities due to interband multi-carrier recombination, while intraband carrier mobility remains unchanged as carrier–carrier scattering is screened by the presence of polarons. The validated Caughey–Thomas model provides a practical input for device simulations, prevents overestimation of diffusion length under high injection, and suggests that MHPs are excellent candidates for unipolar device applications under high carrier densities.

14 SOLAR ENERGY↗
DOE OSTI2026

Charge Carrier and Spin Diffusion in a Polycrystalline and Single Crystal Lead Halide Perovskite Semiconductor

The directed movement of spin-bearing excitations in semiconductors can enable many potential schemes for spintronic applications. Understanding the mechanism of spin motion, as opposed to charge carrier or exciton motion, may require specialized techniques and sample conditions. Here we employ the noncontact and time-resolved technique of light-induced transient grating spectroscopy (LITG) to measure both carrier and spin motion in a perovskite semiconductor with varying crystallinity. The carrier motion aligns with expectations from past studies undertaken at low fluence, revealing an ambipolar diffusion coefficient on the order of 1 cm 2 /sec and diffusion length of roughly 1.5 μm for single crystals that is strongly attenuated as crystallite size decreases. The spin diffusion is measured with cross-polarized excitation beams and uncovers an intensity-dependent coefficient rising above 100 cm 2 /sec. The fast room-temperature spin relaxation limits the spin diffusion length, but the remarkable speedup of spin over carrier diffusion suggests a mechanism involving a combination of exchange-mediated and doping effects that enhance spin transport.

14 SOLAR ENERGY↗
DOE OSTI2026

Dynamic Nanocrystal-Ligand Boundaries: Reversible Photoinduced Ligand Detachment from Quantum Dots in Solution

The porosity of ligand shells of colloidal quantum dots (QDs) can influence the overall rate and yield of charge transfer processes occurring at their surfaces. However, the density of ligand shells on QDs can also influence their colloidal and photochemical stability. We used time-resolved infrared spectroscopy to show that photoinduced ligand detachment, the tendency for certain ligands to detach from QD surfaces when the nanocrystals are promoted to their excitonic excited states, can be used to transiently enhance the porosity of oleic acidpassivated CdSe QDs in solution. Furthermore, we synthesized CdSe QDs with varying ligand shell densities to examine the corresponding influence that van der Waals interactions among ligands have on the yield of photoinduced ligand detachment and the time scale on which ligands return to QD surfaces. We observed that oleic acid ligands on CdSe QDs with lower shell densities have a higher probability of escape for longer periods of time. Despite this, oleic acid ligands on fully passivated CdSe QDs are still able to photodetach, resulting in a transient increase of their ligand shell porosity. In contrast, QDs with multilayer ligand coronas exhibit negligible photoinduced ligand detachment because the outer molecular layers introduce a type of cage effect, preventing the escape of the interior ligands. Our findings suggest the intriguing possibility that photoinduced ligand detachment can be used to transiently decrease the density of ligand shells of QDs to facilitate charge transfer processes while still allowing them to be fully passivated between excitation events for photochemical and colloidal stability.

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

Kinetics vs Thermodynamics: Engineering Photoredox Reactivity from an Upper Excited State of Fe II

Ultrafast deactivation of metal-to-ligand charge-transfer (MLCT) states into low-lying metal-centered states has long limited the utility of first-row transition-metal complexes in a broad range of applications, including photoredox catalysis. Here we bypass such limitations using a Fe II -pyridinium carbene complex to enable electron transfer reactivity directly from a higher-lying MLCT state. Time-resolved spectroscopic data revealed that the higher-energy MLCT manifold relaxes to lower-lying metal-centered states with a time constant of ca. 3 ps. The nature of the metal-centered (MC) excited state as a 3 MC (i.e., S = 1) was inferred from variable-temperature transient absorption studies. These experiments also revealed that the ground-state recovery process occurs at or very near the “Marcus barrierless regime”, thereby allowing for an estimate of its excited-state redox potential. The detailed picture of the energetics of this compound that emerged enabled fine-tuning of competing thermodynamic and kinetic pathways to effect electron transfer from the higher-energy MLCT excited state prior to relaxation to the ligand-field manifold. We believe these results open an unexplored landscape for the use of earth-abundant first-row transition-metal-based chromophores for applications in light-to-chemical energy conversion.

charge transfer↗
DOE OSTI2026

pH-Dependent Vibrational Dynamics Drives Excited-State Quenching in the Phycobiliprotein Complex PC645

Phycocyanin 645 (PC645) is a closed-form lightharvesting complex found in the lumen of the photosynthetic membrane of cryptophyte algae. These peripheral antenna complexes contain bilin chromophores that absorb sunlight and transfer excitation energy to the core antenna complexes embedded in the thylakoid membrane. The location of cryptophyte antenna complex on the luminal side of the membrane is unusual. During photosynthetic activity, the pH of the lumen drops, by up to two pH units. There is little known about how this pH-change affects the light-harvesting complexes. In this study, we report multiscale simulations using a computationally efficient density functional tight-binding framework to investigate the spectroscopy and excitation energy transfer in the PC645 complex. Complementary experiments were conducted using both steady-state and time-resolved spectroscopic measurements at low, neutral, and high pH values. Our study shows that (de)protonation of specific bilin pigments, namely, the mesobiliverdins (MBVs), modulates the excitation energies, excitonic couplings, and spectral densities. These changes cause excitation transfer rates to increase by up to a factor of two to three, leading to pH-dependent energy transfer pathways in the complex. Using this model, we calculated the pH-dependent fluorescence quantum yield of the system, obtaining quantitative agreement with the experimental results. These computational simulations, supported by experiments, identify MBVs as a more prominent excitation sink than previously realized, and that this role is tuned by pH.

Maity, Sayan [Constructor Univ., Bremen (Germany);↗
DOE OSTI2026

Hidden Spin-Valley Locking Stabilizes Nanosecond Spin Polarization in 2D Perovskites

Room-temperature spin control in semiconductors is fundamental to spin-optoelectronics. Although inversion symmetry breaking offers one path for spin control in semiconductors, strong spin dephasing at elevated temperatures remains a persistent limitation. Hidden spin polarization without global symmetry breaking provides another promising material design strategy, but robust spin stabilization from this effect has yet to be experimentally realized. Here we show spin stabilization at room temperature in two-dimensional hybrid organic-inorganic perovskites through hidden spin-valley locking. We use functional non-primary ammonium cations to induce symmetry-breaking distortions in the metal-halide layers, producing giant local spin splitting while preserving global inversion symmetry. Time-resolved circular dichroism measurements reveal optically generated spin-polarized carriers that persist for 686 ps in (AzOH)2PbI4 and 4.7 ns in the lead-free analogue (AzOH)2SnI4 at room temperature, without an external magnetic field. First-principles calculations suggest that these long lifetimes arise from hidden spin valleys, enhanced dielectric screening and reduced spin-orbit-induced scattering in the Sn-based compound. Our design paradigm unlocks inversion-symmetric semiconductors with ultralong spin lifetimes, broadening the materials options for light-driven spin control.

14 SOLAR ENERGY↗
DOE OSTI2026

Elucidating Norrish type I reactive pathways by ultrafast X-ray absorption spectroscopy

Norrish type I reactions selectively cleave carbon–carbon bonds directly adjacent to carbonyl groups. Despite their broad use in combination with aromatic carbonyls for additive manufacturing and dental UV curing applications, the nature of the photochemically active state and its population mechanism remain insufficiently understood. Detailed mechanistic insight requires mapping of the photoexcited population flow involving internal conversion and intersystem crossing. We present a time-domain study of gas phase acetophenone as a prototypical aromatic carbonyl combining soft X-ray time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy at the oxygen K-edge with ab initio multiple spawning (AIMS) simulations. Exploiting the specific sensitivity of TR-NEXAFS spectroscopy to states with nπ* character, we observe population transfer from the initially excited 1 ππ* state to the 1 nπ* state with a time constant of (0.13 ± 0.02) ps after an initial induction period of (0.12 ± 0.02) ps without population transfer, in quantitative agreement with the AIMS simulations. The population in the 1 nπ* state subsequently decays via intersystem crossing, likely mediated by a 3 ππ* state, within (3.17 ± 0.66) ps to a long-lived 3 nπ* state, which is presumed to be active towards Norrish type I chemistry.

Graßl, Martin [SLAC National Accelerator Laborator↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Anomalous lattice expansion of bcc Ta in helium-loaded diamond anvil cells

Anomalous lattice expansion of bcc Ta is observed in helium-loaded diamond anvil cell (DAC) experiments near 6 GPa and above 400 K, while co-loaded W, Mo, and Au exhibit normal thermal expansion. The Ta volume expansion increases with temperature, reaching ~14% near 650 K before plateauing, and is largely retained after recovery to ambient conditions. Time-resolved measurements at 500 K reveal progressive expansion over several hours, consistent with gradual, thermally activated He incorporation. The largely retained expansion during subsequent cooling and holding at 400 K indicates that the incorporated He is predominantly kinetically trapped. The data reveal two distinct regimes: below ~550 K, modest and spatially uniform expansion consistent with interstitial He trapping at isolated, randomly distributed vacancies; above ~550 K, He-vacancy cluster growth producing dramatically larger expansion, peak broadening, and a core–shell microstructure revealed by micron-resolution X-ray diffraction mapping. These results demonstrate that He is not always a passive pressure medium in DAC experiments and reveal a previously unrecognized regime of He-mediated lattice modification in Ta, complementing conventional ion-implantation studies relevant to nuclear and fusion structural materials.

Diamond anvil cell↗
DOE OSTI2026

AgnosticStorageDynamics.jl

SAND2026-22949O AgnosticStorageDynamics.jl is a Julia software package for technology-agnostic energy storage dynamics. It provides a differentiable storage model with charge and discharge limits, efficiencies, standing losses, spill, unmet demand, and time-resolved energy state tracking. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Michelen Strofer, Carlos [Sandia National Lab. (SN↗