Results for “predictive state estimation” — Page 64

DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

DOE OSTI2026

Predicting the single-site and multi-site event discrimination power of dual-phase time projection chambers

Dual-phase xenon time projection chambers (TPCs) are widely used in searches for rare dark matter and neutrino interactions, in part because of their excellent position reconstruction capability in 3D. Despite their millimeter-scale resolution along the charge drift axis, xenon TPCs face challenges in resolving single-site (SS) and multi-site (MS) interactions in the transverse plane. In this paper, we build a generic TPC model with an idealized signal readout, and use Fisher Information (FI) to predict its theoretical capability of differentiating SS and MS events using the electroluminescence signal. We also demonstrate via simulation that, when only statistical photon noise is present, the theoretical limits can be approached with conventional reconstruction algorithms like maximum likelihood estimation, and with a convolutional neural network classifier. The implications of this study on future TPC experiments will be discussed.

Physics↗

DOE OSTI2025

Predictions for the Detectability of Milky Way Satellite Galaxies and Outer-Halo Star Clusters with the Vera C. Rubin Observatory

We predict the sensitivity of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) to faint, resolved Milky Way satellite galaxies and outer-halo star clusters. We characterize the expected sensitivity using simulated LSST data from the LSST Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) accessed and analyzed with the Rubin Science Platform as part of the Rubin Early Science Program. We simulate resolved stellar populations of Milky Way satellite galaxies and outer-halo star clusters over a wide range of sizes, luminosities, and heliocentric distances, which are broadly consistent with expectations for the Milky Way satellite system. We inject simulated stars into the DC2 catalog with realistic photometric uncertainties and star/galaxy separation derived from the DC2 data itself. We assess the probability that each simulated system would be detected by LSST using a conventional isochrone matched-filter technique. We find that assuming perfect star/galaxy separation enables the detection of resolved stellar systems with $M_V$ = 0 mag and $r_{1/2}$ = 10 pc with >50% efficiency out to a heliocentric distance of ~250 kpc. Similar detection efficiency is possible with a simple star/galaxy separation criterion based on measured quantities, although the false positive rate is higher due to leakage of background galaxies into the stellar sample. When assuming perfect star/galaxy classification and a model for the galaxy-halo connection fit to current data, we predict that 89 +/- 20 Milky Way satellite galaxies will be detectable with a simple matched-filter algorithm applied to the LSST wide-fast-deep data set. Different assumptions about the performance of star/galaxy classification efficiency can decrease this estimate by ~7%-25%, which emphasizes the importance of high-quality star/galaxy separation for studies of the Milky Way satellite population with LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Polymers from Levoglucosenone and Other Anhydrosugars (CRADA Final Report)

The goal of this work is to identify promising drop-in or novel polymers that can be produced from levoglucosenone (LGO), CyreneTM, or other anhydrous sugars. LGO is a potentially important feedstock for polymers produced from renewable sources and can be produced in high yields by the company Circa from the pyrolysis of cellulose and biomass. NLR has a long history with anhydrosugars, pyrolysis and biomass-based polymers and has recently developed machine learning (ML) tools for predicting polymer properties. The DOE has recently supported efforts to identify coproducts from pyrolysis of biomass and their use in the production of performance advantaged polymers. DOE is also supporting the development of ML and molecular modeling for predicting and understanding polymer properties. NLR will leverage this expertise and experience to map out the pathways to polymers from LGO, use ML to identify the properties of novel polymers and estimate the market application space. From this effort NLR will deliver a summary report that outlines the most promising pathways based on quantitative and qualitative metrics for producing polymers from LGO.

09 BIOMASS FUELS↗

NASA NTRS1980

Oxidation state of the atmosphere and crust of Venus from Pioneer Venus results

Recent spacecraft observations of Venus permit a detailed model of sulfur chemistry in the atmosphere-lithosphere system. Pioneer Venus experiments confirm that, as predicted, COS and H2S are dominant over SO2 in the lower atmosphere, and that the equilibrium concentrations of S2 and S3 are significant. Many criteria serve to bracket the oxidation state of the crust: it is nearly certain that the S2(2-)/SO4(2-) buffer regulates the oxygen fugacity, and that FeO is at least as abundant as F2O3 in crustal silicates. A highly oxidized crust (as, for example, would result from O2 absorption complementary to escape of vast amounts of H2) is incompatible with the gas-phase sulfur chemistry. If the Pioneer Venus mass spectrometer estimates of the abundance of sulfur gases are correct, earth-like models for the bulk composition of Venus are seriously in error, and a far lower FeO content is required for Venus.

Lewis, J. S.↗

DOE OSTI2025

Large Neutrino "Collider"

We propose using current and future large-volume neutrino telescopes as ``Large Neutrino Colliders" (L$ν$Cs) to explore TeV-scale physics beyond the Standard Model. Cosmic neutrinos with energies above 100 PeV colliding with nucleons in the detector reach center-of-mass energies beyond the 14 TeV limit of the Large Hadron Collider (LHC). Using recently predicted and measured high-energy and ultra-high-energy neutrino fluxes from IceCube and KM3NeT, we estimate mass-scale sensitivities for representative new physics scenarios at 1--30 km$^3$ L$ν$Cs. Our results demonstrate that L$ν$Cs provide a novel avenue to probe multi-TeV particles with sensitivities comparable to, or even surpassing, those of the LHC.

Bai, Yang [Wisconsin U., Madison; Argonne]↗

DOE OSTI2025

Constraining cosmological parameters using the pairwise kinematic Sunyaev-Zel’dovich effect with CMB-S4 and future galaxy cluster surveys

We present a forecast of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) measurement that will be achievable with the future CMB-S4 experiment. CMB-S4 is the next stage for ground-based cosmic microwave background experiments, with a planned wide-area survey that will observe approximately 50% of the sky. We construct a simulated sample of galaxy clusters that have been optically selected in a Legacy Survey of Space and Time–like survey and have spectroscopic redshifts. For this cluster sample, assuming the likelihood is Gaussian, we predict that CMB-S4 will reject the null hypothesis of zero pairwise kSZ signal at 36⁢𝜎. We estimate the effects of systematic uncertainties such as scatter in the mass-richness scaling relation and cluster miscentering. We find that these effects can reduce the signal-to-noise ratio of the CMB-S4 pairwise kSZ measurement by 20%. We explore the constraining power of the measured kSZ signal in combination with measurements of the galaxy clusters’ thermal SZ emission on two extensions to the standard cosmological model. The first extension allows the dark energy equation of state 𝑤 to vary. We find the CMB-S4 pairwise kSZ measurement yields a modest reduction in the uncertainty on 𝑤 by a factor of 1.36 over the Planck’s 2018 uncertainty. The second extension tests general relativity by varying the growth index 𝛾. In conclusion, we find that CMB-S4’s pairwise kSZ measurement will yield a 28⁢𝜎 constraint on 𝛾 and strongly constrain alternative theories of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1998

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Elishakoff, Isaac↗

NASA NTRS1995

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗

NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗

NASA NTRS2006

An Intercomparison of Microphysical Retrieval Algorithms for Upper Tropospheric Ice Clouds

The large horizontal extent, location in the cold upper troposphere, and ice composition make cirrus clouds important modulators of the earth's radiation budget and climate. Cirrus cloud microphysical properties are difficult to measure and model because they are inhomogeneous in nature and their ice crystal size distribution and habit are not well characterized. Accurate retrievals of cloud properties are crucial for improving the representation of cloud scale processes in large-scale models and for accurately predicting the earth's future climate. A number of passive and active remote sensing retrieval algorithms exist for estimating the microphysical properties of upper tropospheric clouds. We believe significant progress has been made in the evolution of these retrieval algorithms in the last decade, however, there is room for improvement. Members of the Atmospheric Radiation measurement program (ARM) Cloud properties Working Group are involved in an intercomparison of optical depth(tau), ice water path, and characteristic particle size in clouds retrieved using ground-based instruments. The goals of this intercomparison are to evaluate the accuracy of state-of-the-art algorithms, quantify the uncertainties, and make recommendations for improvement.

Comstock, Jennifer M.↗

NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology) is a developing technology enabling in-situ aerosol particle sampling and analysis in a small spacecraft mission envelope. It integrates a passive aerosol sample collection system into a probe’s thermal protection system (TPS) to remove the need for heat shield separation and active descent control (parachutes, gliders, etc.). AERACEPT is being validated against the requirements of the Nephele mission concept, which targets the middle and lower Venus cloud layers. The proposed technology employs 3D Carbon-Carbon (3D-CC) at the probe’s nose to withstand the extreme environments encountered in a Venus entry without producing pyrolysis gases that would contaminate the mission sample collection. Any shape change of this 3D-CC sample inlet will affect airflow through the sampling system, thereby impacting aerosol sample collection bias and efficiency. Accurate prediction of the nose’s material response during entry is there-fore of fundamental importance to the development and practical use of this technology. As a result, the AERACEPT project is planning a test campaign in the newly built PlasmatronX facility with a novel “Open Iso-Q” test article to validate state-of-the-art material response tools for a simulated Venus entry. The PlasmatronX is a 350 kW inductively-coupled plasma facility developed and run by the Center for Hypersonics and Entry Systems Studies at the University of Illinois at Urbana-Champaign. The facility can support ground testing for a variety of planetary destinations by simulating entry conditions in Nitrogen, Air, and Carbon Dioxide. To best simulate the Venus atmosphere and match flight-like recession, the planned AERACEPT test campaign will use Carbon Dioxide as its test gas. AERACEPT has designed two test articles for this campaign: 1) a standard test article whose curved surface approximates a constant applied heat flux (“Closed Iso-Q” model) and 2) a novel inlet test article with a through-hole at the stagnation point (“Open Iso-Q” model). Both article types consist of a 3D-CC sample bonded to a graphite fixture. All articles will be laser scanned before and after testing to estimate shape change and recession, and the “Closed Iso-Q” articles will be instrumented with thermocouples to provide temperature histories. The test article temperatures and inlet shape change will be compared against simulations run with the Porous material Analysis Toolbox based on OpenFOAM (PATO) to validate the AERACEPT 3D-CC material response model. Due to the high thermal conductivity of 3D-CC, multidimensional effects have a large influence on the test article temperature and the standard 1-D material response tools are not sufficient for test planning. The test campaign features (e.g. article design, run duration, heat fluxes) are therefore informed by a series of 3-D simulations using state-of-the-art material response tools such as PATO. The test campaign is scheduled for spring of 2024. Expected results include 1) Key drivers of selected test conditions (e.g. matching mission flight atmosphere & recession), 2) Simulation results that informed test article design, 3) Acquired test data such as recession measurements, temperature histories, and photos, and 4) Comparisons between test results and simulation predictions.

AERACEPT↗

NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology) is a developing technology enabling in-situ aerosol particle sampling and analysis in a small spacecraft mission envelope. It integrates a passive aerosol sample collection system into a probe’s thermal protection system (TPS) to remove the need for heat shield separation and active descent control (parachutes, gliders, etc.). AERACEPT is being validated against the requirements of the Nephele mission concept, which targets the middle and lower Venus cloud layers. The proposed technology employs 3D Carbon-Carbon (3D-CC) at the probe’s nose to withstand the extreme environments encountered in a Venus entry without producing pyrolysis gases that would contaminate the mission sample collection. Any shape change of this 3D-CC sample inlet will affect airflow through the sampling system, thereby impacting aerosol sample collection bias and efficiency. Accurate prediction of the nose’s material response during entry is there-fore of fundamental importance to the development and practical use of this technology. As a result, the AERACEPT project is planning a test campaign in the newly built PlasmatronX facility with a novel “Open Iso-Q” test article to validate state-of-the-art material response tools for a simulated Venus entry. The PlasmatronX is a 350 kW inductively-coupled plasma facility developed and run by the Center for Hypersonics and Entry Systems Studies at the University of Illinois at Urbana-Champaign. The facility can support ground testing for a variety of planetary destinations by simulating entry conditions in Nitrogen, Air, and Carbon Dioxide. To best simulate the Venus atmosphere and match flight-like recession, the planned AERACEPT test campaign will use Carbon Dioxide as its test gas. AERACEPT has designed two test articles for this campaign: 1) a standard test article whose curved surface approximates a constant applied heat flux (“Closed Iso-Q” model) and 2) a novel inlet test article with a through-hole at the stagnation point (“Open Iso-Q” model). Both article types consist of a 3D-CC sample bonded to a graphite fixture. All articles will be laser scanned before and after testing to estimate shape change and recession, and the “Closed Iso-Q” articles will be instrumented with thermocouples to provide temperature histories. The test article temperatures and inlet shape change will be compared against simulations run with the Porous material Analysis Toolbox based on OpenFOAM (PATO) to validate the AERACEPT 3D-CC material response model. Due to the high thermal conductivity of 3D-CC, multidimensional effects have a large influence on the test article temperature and the standard 1-D material response tools are not sufficient for test planning. The test campaign features (e.g. article design, run duration, heat fluxes) are therefore informed by a series of 3-D simulations using state-of-the-art material response tools such as PATO. The test campaign is scheduled for spring of 2024. Expected results include 1) Key drivers of selected test conditions (e.g. matching mission flight atmosphere & recession), 2) Simulation results that informed test article design, 3) Acquired test data such as recession measurements, temperature histories, and photos, and 4) Comparisons between test results and simulation predictions.

AERACEPT↗

DOE OSTI2021

Search for R-parity-violating supersymmetry in a final state containing leptons and many jets with the ATLAS experiment using $\sqrt{s}=13\hbox{TeV}$ proton–proton collision data

A search for R-parity-violating supersymmetry in final states characterized by high jet multiplicity, at least one isolated light lepton and either zero or at least three b-tagged jets is presented. The search uses 139fb –1 of $\sqrt{s}=13\hbox{TeV}$ proton–proton collision data collected by the ATLAS experiment during Run 2 of the Large Hadron Collider. The results are interpreted in the context of R-parity-violating supersymmetry models that feature gluino production, top-squark production, or electroweakino production. The dominant sources of background are estimated using a data-driven model, based on observables at medium jet multiplicity, to predict the b-tagged jet multiplicity distribution at the higher jet multiplicities used in the search. Machine-learning techniques are used to reach sensitivity to electroweakino production, extending the data-driven background estimation to the shape of the machine-learning discriminant. No significant excess over the Standard Model expectation is observed and exclusion limits at the 95% confidence level are extracted, reaching as high as 2.4 TeV in gluino mass, 1.35 TeV in top-squark mass, and 320 (365) GeV in higgsino (wino) mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Gate Set Tomography

Gate set tomography (GST) is a protocol for detailed, predictive characterization of logic operations (gates) on quantum computing processors. Early versions of GST emerged around 2012-13, and since then it has been refined, demonstrated, and used in a large number of experiments. This paper presents the foundations of GST in comprehensive detail. The most important feature of GST, compared to older state and process tomography protocols, is that it is calibration-free . GST does not rely on pre-calibrated state preparations and measurements. Instead, it characterizes all the operations in a gate set simultaneously and self-consistently, relative to each other. Long sequence GST can estimate gates with very high precision and efficiency, achieving Heisenberg scaling in regimes of practical interest. In this paper, we cover GST's intellectual history, the techniques and experiments used to achieve its intended purpose, data analysis, gauge freedom and fixing, error bars, and the interpretation of gauge-fixed estimates of gate sets. Our focus is fundamental mathematical aspects of GST, rather than implementation details, but we touch on some of the foundational algorithmic tricks used in the pyGSTi implementation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Predictive Analytics for Hydropower Fleet Intelligence

A primary challenge in hydropower industry is the ability to maintain cost-competitiveness, reliability, and security of hydropower assets through evolving power system contexts and aging of the fleet. Maintaining cost-effective and reliable operations under these conditions is expected to require new modernization and maintenance paradigms for changing contexts. Changes in existing practices for O&M will require an understanding of the current state and health of hydropower assets, and the impact of changing paradigms on asset health and reliability. The Hydropower Fleet Intelligence project is developing and evaluating standardized methodologies and analysis tools for data-driven asset reliability and management technologies for hydropower, leading to eventual predictive maintenance planning, repair/replacement decision making, and asset-reliability and cost-optimized operations. A key question is the feasibility of using existing data sets at hydropower facilities to perform assessments of asset reliability. This document uses data from hydropower facilities to assess the potential for using available analytics methods for asset reliability estimates. In addition to reliability assessments, the feasibility of using existing analytics techniques for several other potential applications is discussed. Finally, a case study that a data-driven model is trained to learn nominal operations via vibration data from an asset of a certain plant, and then utilized to identify anomalies on a similar asset from a different plant, highlighting the generic use of proposed Prognostics and Health Management (PHM) approaches.

Yucesan, Yigit↗

NASA NTRS1986

Sequential estimation and satellite data assimilation in meteorology and oceanography

The central theme of this review article is the role that dynamics plays in estimating the state of the atmosphere and of the ocean from incomplete and noisy data. Objective analysis and inverse methods represent an attempt at relying mostly on the data and minimizing the role of dynamics in the estimation. Four-dimensional data assimilation tries to balance properly the roles of dynamical and observational information. Sequential estimation is presented as the proper framework for understanding this balance, and the Kalman filter as the ideal, optimal procedure for data assimilation. The optimal filter computes forecast error covariances of a given atmospheric or oceanic model exactly, and hence data assimilation should be closely connected with predictability studies. This connection is described, and consequences drawn for currently active areas of the atmospheric and oceanic sciences, namely, mesoscale meteorology, medium and long-range forecasting, and upper-ocean dynamics.

Ghil, M.↗

DOE OSTI2021

Estimation of secondary organic aerosol viscosity from explicit modeling of gas-phase oxidation of isoprene and α -pinene

Abstract. Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we apply the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model to conduct explicit chemical modeling of isoprene photooxidation and α-pinene ozonolysis and their subsequent SOA formation. The detailed gas-phase chemical schemes from GECKO-A are implemented into a box model and coupled to our recently developed glass transition temperature parameterizations, allowing us to predict SOA viscosity. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model, including autoxidation and particle-phase reactions, leading to the formation of high-molar-mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

54 ENVIRONMENTAL SCIENCES↗