Search NASA⌕ Search

SEARCH · Search NASA

Results for “state variable”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Studying phonon coherence with a quantum sensor

Nanomechanical oscillators offer numerous advantages for quantum technologies. Their integration with superconducting qubits shows promise for hardware-efficient quantum error-correction protocols involving superpositions of mechanical coherent states. Limitations of this approach include mechanical decoherence processes, particularly two-level system (TLS) defects, which have been widely studied using classical fields and detectors. In this manuscript, we use a superconducting qubit as a quantum sensor to perform phonon number-resolved measurements on a piezoelectrically coupled phononic crystal cavity. This enables a high-resolution study of mechanical dissipation and dephasing in coherent states of variable size ($\overline{n}$ ≃ 1 – 10 phonons). We observe nonexponential relaxation and state size-dependent reduction of the dephasing rate, which we attribute to TLS. Using a numerical model, we reproduce the dissipation signatures (and to a lesser extent, the dephasing signatures) via emission into a small ensemble (N = 5) of rapidly dephasing TLS. Our findings comprise a detailed examination of TLS-induced phonon decoherence in the quantum regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

First-Principles Evaluation of Proton Hopping in Tetrahedral Oxide Motifs

Proton-conducting oxides (PCOs) are important materials used as ionic conductors for energy conversion technologies. Existing research efforts on PCO optimization and discovery generally focus on complex perovskite-based oxides that require doping and alloying to engineer oxygen deficiency and high proton conductivity. However, the variety of chemical compositions and coordination environments in oxides poses challenges for efficient materials design. In this computational study, we construct a database of simplified motifs to elucidate the relationship between fundamental materials chemistry and proton kinetics. Specifically, we focus on the zincblende crystal structure as a proxy for tetrahedral metal–oxide (M–O) coordination environments. We systematically quantified the effects of cation type, oxidation states, and M–O bond lengths on the proton hopping barrier, and found that strong M–O bonds and metal cations with large and variable oxidation states (e.g., Mo 6+ , V 5+ ) lead to smaller proton hopping barriers. By mapping the candidate cations and their preferred bond geometries onto materials databases such as the Inorganic Crystal Structure Database (ICSD) and Materials Project, we identified real materials containing the corresponding metal–oxide units. In general, we observed good agreement between the calculated proton hopping barriers obtained in real crystal structures and those predicted by our motif database. We also discuss the limitations of our model and possible future extensions to improve its predictive capabilities. Overall, our model provides a first step for the rational design and quick screening of energy-efficient PCOs.

organic↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2026

Generation of scalable genuine multipartite Gaussian entanglement with a parametric amplifier network

Genuine multipartite entanglement is a valuable resource in quantum information science, as it exhibits inseparability across all possible bipartitions of a multimode quantum state. The large degree of quantum correlations can be exploited in various quantum information protocols, such as teleportation, dense coding, and quantum interferometry. Here, in this study, we propose a scheme to generate scalable genuine multipartite continuous-variable entangled states of light using a parametric amplifier network. We verify the presence of genuine quadripartite, hexapartite, and octapartite entanglement through a violation of the positive partial transpose criteria. Additionally, we use 𝛼-entanglement of formation to demonstrate the scalability of our approach to an arbitrary number of 2⁢𝑁 genuinely entangled parties by taking advantage of the symmetries present in our scheme.

Kim, Saesun [Univ. of Oklahoma, Norman, OK (United↗
DOE OSTI2026

ARM Data-Oriented Metrics and Diagnostics Package for Climate Model Evaluation

A Python-based metrics and diagnostics package is currently being developed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Infrastructure Team at Lawrence Livermore National Laboratory (LLNL) to facilitate the use of long-term, high-frequency measurements from the ARM Facility in evaluating the regional climate simulation of clouds, radiation, and precipitation. This metrics and diagnostics package computes climatological means of targeted climate model simulation and generates tables and plots for comparing the model simulation with ARM observational data. The Coupled Model Intercomparison Project (CMIP) model data sets are also included in the package to enable model intercomparison as demonstrated in Zhang et al. (2017). The mean of the CMIP model can serve as a reference for individual models. Basic performance metrics are computed to measure the accuracy of mean state and variability of climate models. The evaluated physical quantities include cloud fraction, temperature, relative humidity, cloud liquid water path, total column water vapor, precipitation, sensible and latent heat fluxes, and radiative fluxes, with plan to extend to more fields, such as aerosol and microphysics properties. Process-oriented diagnostics focusing on individual cloud- and precipitation-related phenomena are also being developed for the evaluation and development of specific model physical parameterizations. The version 1.0 package is designed based on data collected at ARM’s Southern Great Plains (SGP) Research Facility, with the plan to extend to other ARM sites. The metrics and diagnostics package is currently built upon standard Python libraries and additional Python packages developed by DOE (such as CDMS and CDAT). The ARM metrics and diagnostic package is available publicly with the hope that it can serve as an easy entry point for climate modelers to compare their models with ARM data. In this report, we first present the input data, which constitutes the core content of the metrics and diagnostics package in section 2, and a user's guide documenting the workflow/structure of the version 1.0 codes, and including step-by-step instruction for running the package in section 3.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil Carbon and Nitrogen Stocks Across Hillslopes Underlain by Continuous Permafrost in the Northern Arctic Foothills, Alaska, United States

Constraining the variability of soil organic carbon (SOC) and total nitrogen (TN) stocks across hillslopes in Low Arctic permafrost-affected landscapes remains a significant challenge for improving global estimates of permafrost SOC stocks. We investigated SOC and TN stocks across hillslopes at two sites in the Arctic Foothills of Alaska, United States (Happy Valley and Sagwon Hills). Average SOC and TN stocks for the 0–1-m depth interval were high (52.0 ± 15.1 kg C m −2 and 2.74 ± 0.82 kg N m −2 ) and linearly related ( R 2 = 0.74, p < 0.0001). Unlike soils of other permafrost and nonpermafrost landscapes, variability was greatest within rather than between hillslope positions. Furthermore, SOC and TN stocks in the surface 1 m did not exhibit strong patterns by hillslope position and were only weakly associated with major geomorphic parameters that typically predict SOC and TN stocks well in other landscapes. Although sampling at upper hillslope positions was largely limited to depths of less than 1.5 m due to the presence of coarse fragments in reworked glacial till, deeper observations at lower hillslope positions (footslopes, toeslopes, and basins) revealed significantly larger SOC stocks (92.0 ± 18.0 kg C m −2 at 2 m; 117.1 ± 10.4 kg C m −2 at 3 m). The unique small-scale variability in ice content, cryoturbation, patterned ground, and organic layer thickness on these broad, Low Arctic sites contributes to the relatively homogeneous distribution of SOC and TN stocks across hillslope positions in the top 1 m, but a future focus on deeper sampling may reveal greater differences in SOC and TN stocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Long duration battery sizing, siting, and operation under wildfire risk using progressive hedging

Battery sizing and siting problems are computationally challenging due to the need to make long-term planning decisions that are cognizant of short-term operational decisions. This paper considers sizing, siting, and operating batteries in a power grid to maximize their benefits, including price arbitrage and load shed mitigation, during both normal operations and periods with high wildfire ignition risk. Here we formulate a multi-scenario optimization problem for long duration battery storage while considering the possibility of load shedding during Public Safety Power Shutoff (PSPS) events that de-energize lines to mitigate severe wildfire ignition risk. To enable a computationally scalable solution of this problem with many scenarios of wildfire risk and power injection variability, we develop a customized temporal decomposition method based on a progressive hedging framework. Extending traditional progressive hedging techniques, we consider coupling in both placement variables across all scenarios and state-of-charge variables at temporal boundaries. This enforces consistency across scenarios while enabling parallel computations despite both spatial and temporal coupling. The proposed decomposition facilitates efficient and scalable modeling of a full year of hourly operational decisions to inform the sizing and siting of batteries. With this decomposition, we model a year of hourly operational decisions to inform optimal battery placement for a 240-bus WECC model in under 70 min of wall-clock time.

25 ENERGY STORAGE↗
DOE OSTI2025

Criegee Intermediates Compete Well with OH as a Cleaning Agent for Atmospheric Amides

Elucidating the chemical kinetics of stabilized Criegee intermediates (sCIs) in the atmosphere is critically important for climate modeling. Here, we report a class of very rapid bimolecular reactions of two sCIs, namely, CH 2 OO and syn-CH 3 CHOO, with amides. We used electronic structure calculations and kinetics calculations to elucidate a universal mechanism by which the oxygen atom of the carbonyl group in an amide is added to the carbon atom of the COO group in the sCI with simultaneous transfer of the amide hydrogen to the terminal oxygen atom of the carbonyl oxide of the sCI. The barriers for the mechanism are submerged below reactants by ∼9 kcal/mol, which means that the tight transition state is not the rate-determining step, and the rate constants are determined by loose free energy bottlenecks between the reactants and the precursor complexes. We calculate the rate constants due to these loose and barrierless transition states by variable-reaction-coordinate variational transition-state theory. We find that bimolecular reactions of sCIs with amides are very rapid, with rate constants 1 to 5 × 10 –10 cm 3 molecule –1 s –1 , which is a factor of 2 faster than the rate constants of carboxylic acids with sCIs and 10 2 faster than the rate constants of OH reactions with amides. This has the consequence that, under a wide range of conditions, sCIs are the major sink for atmospheric amides, in contrast to the usual assumption that OH radicals are the dominant sinks.

Amides↗
DOE OSTI2024

Rapid Optimization of Total Variation with Applications in Imaging, Additive Manufacturing, and Qualification

Total Variation optimization penalizes the gradient of a control variable or state. While this work focuses on image processing in particular, it has also found applications in inverse problems and topology optimization. In image processing, the goal is to maintain faithfulness to the original image while denoising and/or deblurring. Additionally, bilevel optimization over the spatially varying regularization weights can illuminate interfaces such as damage regions and other anomalies. We will address two fundamental challenges with TV-optimization: (i) the typical slow convergence of existing TV-optimization methods, and (ii) the selection of spatially varying TV parameters to promote interface detection. Additionally, we will apply such techniques to image data collected in additive manufacturing. In said context, stochasticity in build events induces flaws in the manufactured piece, compromising the integrity of said part. There is a critical need for in-situ monitoring to spot anomalies once they form, and in this setting we apply our total variation and hyperparameter solvers. We will develop a customized algorithm based on for extreme-scale TV-optimization that achieves super-linear or quadratic-convergence, a critical property for real-time, image-by-image analysis. A worst-case outcome is a preprocessing step that enhances image quality in-situ, specifically for out-of-focus and noisy images.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Aircraft in-situ measurements from SOCRATES constrain the anthropogenic perturbations of cloud droplet number

Aerosol-cloud interactions (ACI) in warm clouds alter reflected shortwave radiation by influencing cloud microphysical and macrophysical properties. The variable of state controlling ACI is the cloud droplet number concentration (N d ). Here, we examine the perturbations in N d due to anthropogenic aerosols (ΔN d, PD-PI ) using a perturbed parameter ensemble (PPE) hosted in the sixth Community Atmosphere Model (CAM6). Surrogate models are created for the CAM6 PPE outputs and are used to generate 1 million model variants of CAM6 by sampling 45 sources of parameter uncertainty. The range of uncertain physical parameters related to ACI are constrained with observations of aerosol and cloud properties from SOCRATES. The likely range of uncertain parameters and the associated range of ΔN d, PD-PI are more strongly constrained with observations of N d relative to observations of cloud condensation nuclei. We conduct sensitivity tests of how constraints on ΔN d, PD-PI are affected by systematic uncertainties in observations and our limitations in our surrogate models created for CAM6 PPE outputs. Based on this, we provide guidance on the impact of reducing systematic uncertainty in airborne microphysical observations and in surrogate models.

Song, Ci [Univ. of Wyoming, Laramie, WY (United St↗
DOE OSTI2026

Timeseries Photos of a Variably Inundated Stream: Umtanum Creek, Washington, United States

This dataset is associated with a broader study using game camera timeseries photos collected to evaluate stream variable inundation via changes in width (i.e. wet fraction). Four game cameras were deployed along Umtanum Creek (Washington, United States) to track changes in stream inundation over time. Drone imagery was collected at the same location on October 18, 2024 which was used to construct a digital elevation model (DEM) of the streambed topography. The associated paper and data can be found at https://doi.org/10.1016/j.envsoft.2025.106715 (Bao et al., 2025a)) and https://doi.org/10.15485/2589885 (Bao et al., 2025b), respectively. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to this readme, this data package also includes a file-level metadata (FLMD) files that describes each file and a data dictionaries (DD) that describe all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) readme; (4) field metadata; (5) field protocol; and (5) folders containing game camera photos. Game camera photos are organized into folders for each camera (CDL, CUL, CDR, CUR; see readme for information on camera naming) by the month photos were collected. All files are .csv, .jpg, or .pdf.

AI image segmentation↗
DOE OSTI2025

Investigating Particle Phase State Dependencies during a Historic Pacific Northwest Wildfire Event

The function and lifetime of an atmospheric particle are greatly dependent upon its phase state. Studying particle phase states from specific ambient source types is often challenging due to uncontrolled environmental conditions and multiple source contributions. However, a historic Pacific Northwest wildfire event occurred in September 2020, providing a rare opportunity to study a pseudo-single particle type ambient sample (i.e., wildfire) over time. Sampling was performed in Richland, Washington, capturing aged particles from shifting wildfire sources. In addition to source variability, temporal variation in wildfire particle phase state was observed from single particle analyses. Despite this, single particle elemental compositions were homogenous over time (=87% carbonaceous) with no significant temporal variation in organic functional group abundances within individual particles either (e.g., 4% relative standard deviation in alkene contributions). While particle composition is indeed known to affect phase state, water uptake appeared to be the most significant contributor to phase state variability amongst single particles here, as particle aspect ratios (related to phase state) linearly correlated with ambient relative humidity during sampling (R2 = 0.65). Therefore, while limited to a case study, meteorological considerations are found to convey greater comparative importance than particle composition for aged particle phase state predictions.

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

Weather and climate extremes in a changing Arctic

Weather and climate extremes are increasingly occurring in the Arctic. Here, in this Review, we evaluate historical and projected changes in rare Arctic extremes across the atmosphere, cryosphere and ocean and elucidate their driving mechanisms. Clear shifts occur in mean and extreme distributions after ~2000. For instance, pre-2000 to post-2000 observational probabilities of 1.5 standard deviation events increase by 20% for atmospheric heat waves, 76.7% for Atlantic layer warm events, 83.5% for Arctic sea ice loss and 62.9% for Greenland Ice Sheet melt extent — in many cases, low probability, rare extreme events in the early period become the norm in the latter period. These observed changes can be explained using a ‘pushing and triggering’ concept, representing interplay between external forcing and internal variability: long-term warming destabilizes the climate system and ‘pushes’ it to a new state, allowing subsequent variability associated with large-scale atmosphere–ocean–ice interactions and synoptic systems to ‘trigger’ extreme events over different timescales. Ongoing anthropogenic warming is expected to further increase the frequency and magnitude of extremes, such that simulated probabilities of 1.5 standard deviation events increase by 72.6% for atmospheric heat waves, 68.7% for Atlantic layer warm events and 93.3% for Greenland Ice Sheet melt rate between historic (1984–2014) and future (2069–2099) periods under a very high emission scenario. Future research should prioritize the development of physically based metrics, enhance high-resolution observation and modelling capabilities and improve understanding of multiscale Arctic climate drivers.

atmospheric dynamics↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

COMPASS-FME Synoptic Sites Level 1 Sensor Data v2-1

This is the version 2-1 Level 1 (L1) data release for COMPASS-FME environmental sensors located at our synoptic field sites. COMPASS-FME is studying sites in two distinct regions, the Chesapeake Bay and the Western Lake Erie Basin. We established the network at seven "synoptic" (observational) sites along the Chesapeake Bay and Lake Erie coastlines, collectively generating over three million observations per month, to track and comprehend environmental changes where land and water intersect. Additionally, the two regions provide an interesting contrast of saltwater and freshwater coasts that allow us to differentiate the impacts of inundation and coastal water chemistries in two nationally important coastal systems. L1 data are close to raw, but are units-transformed and have out-of-instrument-bounds, out-of-service, and outlier flags added. Duplicates and missing data are removed but otherwise these data are not filtered, and have not been subject to any additional algorithmic or human QA/QC. Any scientific analyses of L1 data should be performed with care. **This dataset will be updated quarterly with new data for the duration of the project** This dataset includes: - An overall dataset README file that describes the current version, gives citation and contact information, etc. - Site- and year-specific folders, each holding up to 12 CSV (comma separated value) data files for each site and plot in that year. - Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, as well as a general description of the site. - Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are normally logged every 15 minutes. Please see v2-0 Synoptic L1 Sensor Package Quick Start.pdf for detailed information on data package structure, temporal coverage, and versioning. This dataset was updated 2026-03-12: (i) data now go through 2025-12-31 (previous end was 2025-06-30) and (ii) dataset and file names updated to “…v2-1” (previously was “v2-0”).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Haze processing of atmospheric particles during wintertime in the Indo-Gangetic Plains

The impacts of haze on visibility, air quality, and climate are not well quantified due to a lack of understanding of the evolution of the mixing state and phase state of atmospheric particles during haze processing. The variability of the mixing state of atmospheric particles contributes significantly to uncertainties associated with the estimated aerosol radiative forcing. We collected particle samples in January 2018 from the highly polluted Indo-Gangetic plain during hazy conditions to study haze-processed particles. Single particle analysis using multi-modal micro-spectroscopy techniques revealed an abundance (40–70% by number) of potassium-rich sulfate particles from biomass-burning influenced smoke. Tilted view imaging showed that most of the organic particles that had inorganic potassium and sulfate inclusions were liquid-like while those without inclusions were more semi-solid. High-resolution mass spectrometry analysis revealed a significant presence of organosulfates and nitroxy-organosulfates in the morning samples (24%) compared to the afternoon samples (9%), despite higher relative humidity in the afternoon. Overall, our results highlight the significant contribution of both organic and inorganic sulfate to the total particulate sulfur budget during haze processing in winter, when anthropogenic emissions such as household burning, agricultural burning, and vehicular emissions are major contributors to particle mass.

Mathai, Susan [Pacific Northwest National Laborato↗
DOE OSTI2026

Oil price states and drivers: An analysis of the second-month spot-futures price differential

Oil remains a dominant component of global energy use, and its price, characterized by frequent changes and an ever-present potential for large swings, continues to be a focus of industry participants, policymakers and analysts attention. Here, this study examines the behavior of future spot oil prices using a continuous-time hidden Markov model (HMM) and daily price data spanning years 2007 to 2024. We identify six states in the second-month WTI spot-futures price differential and assess the roles of eleven futures price, quantity, financial market, and geopolitical/volatility variables in each state. The model yields several insights into the workings of the oil market and the relative roles of these drivers. We find support for several theoretical and empirical findings in the oil market literature, including the role of inventory, volatility/risk, and market responses to contango/backwardation in futures markets. A novel finding is that “normal contango” conditions represent a significant portion of next-day states in our in-sample data. Under the most volatile normal contango state, many of the oil market drivers differ markedly in coefficient signs and magnitudes from those in other states. The resulting model also performed well out-of-sample and would, in addition to understanding the impact of market drivers, be useful for short-term forecasting. Overall, the findings highlight the highly non-linear, regime-dependent interactions of the oil price and its drivers, emphasizing the importance of detailed information to market stakeholders.

Oladosu, Gbadebo A. [Oak Ridge National Laboratory↗