SEARCH · Engineering Papers

Results for “matrix product state”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Finding ways to reduce nuclear waste: searching for the unknown one step at a time

In my home country, Venezuela, research has been stagnant. Due to the political turmoil and the crisis, many educated people have left the country in search for a better life. This has caused a deficit in any technological and scientific advances, making Venezuela one of the first South American countries to have its rate of publications decline by 29% in 2013. Currently, Venezuela lacks the infrastructure and the means to keep up with the research progress as compared to other countries in South America, such as Brazil. Since coming to the United States (US), and currently working for a national laboratory, the active research environment endorses a wide range of careers and engineering programs that allow researchers to thrive at any given field. Researchers have access to funds and tools to succeed in developing materials for the future. There are 17 national laboratories in the US, and all of these have a different research focus/objective. As examples, Los Alamos National Laboratory and Sandia National Laboratory focus is on national homeland security, weapon science, radiation effects, among others. Argonne National Laboratory focuses on nuclear energy, energy storage, high performance computing, etc. At Idaho National Laboratory (INL) the research focuses on innovating nuclear energy and clean energy resources, critical infrastructure materials, along with fuel cycle solutions to manage, dispose and find ways to recycle current and future radiological waste. Compared to other national laboratories, INL focuses slightly more on applied processes and how nuclear energy can be innovated to next reactor design and technologies. The research being conducted at INL made me apply for a Seaborg distinguished postdoctoral position. For the position itself, the researcher must submit a proposal related to actinide chemistry on a research field area. In this position, 50% of my time will be focused on my own proposal. The proposal that I am working on is focused on the innovation of nuclear energy and fuel cycle recycling, which is why I was mainly interested on working at this national laboratory. To give a bit more context of what my proposal is about, a little bit of background is necessary: After the nuclear fuel (UO2) is used in a reactor, the fuel matrix is then characterized by various fission products (FP). Among these FP (including rare earth elements, alkali/alkaline earths, and actinides), many can potentially be recovered through nuclear reprocessing technologies. In pyroprocessing, the used nuclear fuel undergoes electrochemical dissolution into a molten chloride salt mixture in an electrorefiner. Initially, uranium is reduced onto an inert cathode by applied potentials. However, numerous remaining FPs accumulate in the melt and pose challenges for recovery by an inert electrode, particularly the rare earth elements (e.g., Nd, Gd, Pr, Sm) due to their multivalent oxidation states and tendencies toward side reactions, leading to their dissolution in the electrolyte. These recovery challenges result in inefficiencies and necessitate the continual discarding of the molten chloride salt, thereby generating additional waste. Furthermore, the presence of rare earth elements and other fission products in the molten salt electrolyte alters its physical and chemical properties, affecting both uranium recovery efficiency and the longevity of the molten chloride salt. To improve the recovery efficiency of the FP, specifically rare earth elements, I am investigating the fundamental interactions between rare earth elements in the molten chloride salt and their metallic form. The kinetic pathways and the chemical reactions of these elements will give insights on how the recovery efficiency can be improved. The interactions and speciation of these elements are being studied by spectro-electrochemistry at high temperature environments in quartz and other ceramic materials (e.g., alumina crucibles). Some of the challenges I am facing specifically relates the reactivity of some of these elements with different glass and crucible materials. Although my research focuses on fundamental science, it will benefit the applied process by generating new scientific knowledge and closing the gap for an efficient recycling of the waste: one step at a time.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

FIB-ToF-SIMS characterization of irradiated U-10Zr

Post-irradiation examination (PIE) is critical for the performance assessment and qualification of nuclear fuels. Secondary ion mass spectrometry (SIMS) is a powerful materials characterization technique that allows for elemental and isotopic mapping with a depth resolution greater than EDS and EPMA. However, it has not yet been applied to PIE of metallic nuclear fuel. Here, in this work, we characterize an fast neutron spectrum irradiated U-10Zr fuel sample using a time-of-flight SIMS (ToF-SIMS) system connected to a FIB/SEM system, which allows for flexible sample analysis compared to a dedicated ToF-SIMS instrument. Analysis of the resulting hyperspectral micrograph data was aided by the development of an unsupervised machine learning (ML) algorithm that iterates on existing methods to segment the 3D micrographic datasets based on the similarity of mass spectra. The results showed that the FIB-ToF-SIMS instrument was potentially capable of spatially resolving closed fission gas bubbles in 3D by continued ion sputtering of the analyzed volume. Additionally, the ML algorithm proved useful in revealing the chemical segregation of light fission products (those with an atomic mass between approximately 85–105 amu, such as ruthenium and rhodium) plus matrix zirconium, heavy fission products (those with an atomic mass between approximately 135–150 amu, such as the lanthanides) and uranium. Future studies are planned to conduct FIB-ToF-SIMS analysis on more irradiated U-Zr samples to study the constituent redistribution.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗

DOE OSTI2026

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab↗

DOE OSTI2025

Leveraging Hamiltonian simulation techniques to compile operations on bosonic devices

Circuit quantum electrodynamics enables the combined use of qubits and oscillator modes. Despite a variety of available gate sets, many hybrid qubit-boson (i.e. qubit-oscillator) operations are realizable only through optimal control theory, which is oftentimes intractable and uninterpretable. We introduce an analytic approach with rigorously proven error bounds for realizing specific classes of operations via two matrix product formulas commonly used in Hamiltonian simulation, the Lie–Trotter–Suzuki and Baker–Campbell–Hausdorff product formulas. We show how this technique can be used to realize a number of operations of interest, including polynomials of annihilation and creation operators, namely (a) p (a † ) q for integer p, q. We show examples of this paradigm including obtaining universal control within a subspace of the entire Fock space of an oscillator, state preparation of a fixed photon number in the cavity, simulation of the Jaynes–Cummings Hamiltonian, and simulation of the Hong-Ou-Mandel effect. This work demonstrates how techniques from Hamiltonian simulation can be applied to better control hybrid qubit-boson devices.

bosonic qubits↗

DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Commodity Thermoplastic Elastomer-Enabled Templated Synthesis of Large-Pore Ordered Mesoporous Materials

Fabrication of ordered mesoporous materials (OMMs) has predominantly relied on templating-based methods. However, these methods are constrained by several limitations, especially the limited pore sizes attainable with commercially available surfactants used as structure-directing agents. To unlock the full potential of the OMMs, it is essential to develop synthetic strategies that facilitate the production of large-pore OMMs using scalable processes and cost-effective precursors. This work demonstrates the use of thermoplastic elastomer (TPE)-derived carbon replicas for synthesizing ordered mesoporous silica (OMS) and metal oxides (OMMOs) via precursor infiltration and template removal. The nanostructural evolution of the resulting inorganic materials was systematically investigated. Specifically, using tetraethyl orthosilicate (TEOS) as a silica precursor, this method can produce an OMS with relatively large pores. To establish the generalizability of this process, the fabrication approach was extended to other commercially available TPEs with varied chemical compositions and molecular weights while consistently resulting in ordered structures. Additionally, this synthetic strategy can be successfully applied to the production of OMMOs, including tin and titanium oxide matrix chemistries, yielding pore sizes of 16.0 and 19.2 nm, respectively. By developing a general method and using low-cost precursors, this work presents a scalable approach for fabricating large-pore OMMs with tunable pore textures and matrix chemistries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Semi-inclusive deep-inelastic scattering on a polarized spin-1 target. I. Cross section and spin observables

We develop the theoretical framework for semi-inclusive deep-inelastic scattering on a polarized spin-1 target and apply it to scattering on the polarized deuteron with spectator nucleon tagging. In Part I (this article) we present the general form of the semi-inclusive cross section and polarization observables for the spin-1 target. A relativistically covariant formulation in terms of four-vectors and invariant polarization parameters is employed. The target polarization is described by a spin density matrix with vector and tensor polarization. The spin and azimuthal angle dependence of the semi-inclusive cross section is derived and parametrized in terms of invariant structure functions. To validate the result, the structure functions are expressed as photon-target helicity amplitudes with known symmetry properties. The expressions presented here are kinematic (no assumptions about particle production dynamics) and valid in all regions of the deep-inelastic final state (current and target fragmentation regions). In Part II (following article), we consider deep-inelastic scattering on the polarized deuteron with spectator nucleon tagging as a special case of target fragmentation. The semi-inclusive structure functions are computed by separating nuclear and hadronic structure, and the polarization observables are explored as functions of the tagged nucleon momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

We present lattice QCD calculations of the odd Mellin moments of pion valence-quark generalized parton distribution up to fifth order ⟨𝑥 4 ⟩ and for the skewness range [−0.33, 0] using operator product expansion of bilocal quark-bilinear operators. The calculations are performed on an ensemble with lattice spacing 𝑎 = 0.04 fm and valence pion mass 300 MeV, employing boosted pion states with momenta up to 2.428 GeV and momentum transfers reaching 2.748 GeV 2 . We employ ratio-scheme renormalization and next-to-leading logarithmic resummed perturbative matching. At zero skewness, our results are consistent with previous lattice studies. By combining matrix elements at multiple values of skewness and momentum transfer, skewness-dependent moments are obtained through simultaneous polynomiality-constrained fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

An Evaluation and Qualification of U.S.-Based Research Reactors for Irradiation Capabilities Supporting Advanced Nuclear Systems

Irradiation experiments are a prerequisite for evaluating nuclear reactor system designs, analyzing the performance of these systems, and obtaining licenses. Likewise, irradiation facilities are necessary for producing the radioisotopes used in industrial and medical applications. Recent developments in modeling and simulation capabilities and advancements in computational resources have further enabled the design of irradiation experiments for evaluating radiation-induced phenomena and determining nuclear fuel, material, and system design and safety criteria pertaining to both normal and accident scenarios. These computational tools and models require comprehensive experimental datasets acquired under prototypic radiation conditions—for exploring material and system performance under the uniquely harsh environments found in nuclear reactors—to enable verification and validation for qualification and licensing purposes. However, qualification of irradiation experimental facilities, primarily research and test reactors (RTRs), necessitates that their performance be evaluated based on the irradiation environment (e.g. flux, power, testing capabilities) using an appropriate scoring matrix. Although many university campus RTRs are available for research and development (R&D) activities and initiatives, this study focuses on evaluating and qualifying the irradiation facilities (mostly RTRs) within the United States that are suitable for advanced nuclear fuel, material, and system irradiation experiments aimed at establishing operational-performance limits and informing component and fuel designs so as to improve operational efficiencies and mitigate proliferation vulnerabilities, as well as for radioisotope production aimed at multipurpose applications. As a result, the findings of the present study support the acceleration of nuclear fuel and material qualifications, thus hastening new and advanced nuclear energy system demonstrations and radioisotope production efforts by using extended R&D.

irradiation experiment↗

DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗

DOE OSTI2025

Amplitude analysis of ψ3686→γKS0KS0

Using (2712 ± 14) × 106ψ(3686) events collected with the BESIII detector, we perform the first amplitude analysis of the radiative decay ψ3686→γKS0KS0$$ \psi (3686)\to \gamma {K}_S^0{K}_S^0 $$ within the mass region MKS0KS0<2.8$$ {M}_{K_S^0{K}_S^0}<2.8 $$ GeV/c2. Employing a one-channel K-matrix approach for the description of the dynamics of the KS0KS0$$ {K}_S^0{K}_S^0 $$ system, the data sample is well described with four poles for the f0-wave and three poles for the f2-wave. The determined pole positions are consistent with those of well-established resonance states. The observed f0 and f2 states are found to be in agreement with those produced in radiative J/ψ decays. The production behaviors of f0 and f2 poles in ψ3686→γKS0KS0$$ \psi (3686)\to \gamma {K}_S^0{K}_S^0 $$ are qualified with their residues and the converted branching fractions. By comparing with J/ψ→γKS0KS0$$ J/\psi \to \gamma {K}_S^0{K}_S^0 $$ decay, the ratios Bψ3686→γf0,2BJ/ψ→γf0,2$$ \frac{\mathcal{B}\left(\psi (3686)\to \gamma {f}_{0,2}\right)}{\mathcal{B}\left(J/\psi \to \gamma {f}_{0,2}\right)} $$ are determined, which provides crucial experimental inputs on the internal structure of the f0,2 states, especially their potential mixing with glueball components.

Ablikim, M↗

DOE OSTI2020

Reduced methane recovery at high pressure due to methane trapping in shale nanopores

By 2050, shale gas production is expected to exceed three-quarters of total US natural gas production. However, current unconventional hydrocarbon gas recovery rates are only around 20%. Maximizing production of this natural resource thus necessitates improved understanding of the fundamental mechanisms underlying hydrocarbon retention within the nanoporous shale matrix. In this study, we integrated molecular simulation with high-pressure small-angle neutron scattering (SANS), an experimental technique uniquely capable of characterizing methane behavior in situ within shale nanopores at elevated pressures. Samples were created using Marcellus shale, a gas-generative formation comprising the largest natural gas field in the United States. Our results demonstrate that, contrary to the conventional wisdom that elevated drawdown pressure increases methane recovery, a higher peak pressure led to the trapping of dense, liquid-like methane in sub-2 nm radius nanopores, which comprise more than 90% of the measured nanopore volume, due to irreversible deformation of the kerogen matrix. These findings have critical implications for pressure management strategies to maximize hydrocarbon recovery, as well as broad implications for fluid behavior under confinement.

58 GEOSCIENCES↗

DOE OSTI2026

Double copy root of Hawking thermality

The Hawking radiation spectrum from a collapsing null shell can be derived via the double copy of a simpler gauge theory calculation. Analyzing the non-abelian Yang-Mills root of this process, we demonstrate that the radiation spectrum is thermal in the color charge eigenvalue 𝜆, not energy. Considering the 𝑆⁢𝑈⁡(𝑁 𝑐 ) gauge theory in the large 𝑁 𝑐 limit, we find the differential spectrum ⅆ⁢𝑁/ⅆ⁢𝜆 is a product of the gravitationally familiar Planck-like factor and the color phase space density, modeled here as the Wigner semicircle from random matrix theory. This reveals that apparent energy thermality in gravity is the direct dual of charge thermality in its underlying non-abelian gauge theory.

Carrasco, John Joseph M. [Northwestern University,↗

DOE OSTI2021

Measurement of the Higgs boson production rate in association with top quarks in final states with electrons, muons, and hadronically decaying tau leptons at $\sqrt{s} =$ 13 TeV

The rate for Higgs (${\mathrm{H}} $) bosons production in association with either one (${\mathrm{t}} {\mathrm{H}} $) or two (${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton–proton collisions recorded at a center-of-mass energy of $13\,\text {Te}\text {V} $ by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137$\,\text {fb}^{-1}$. The analysis is aimed at events that contain ${\mathrm{H}} \rightarrow {\mathrm{W}} {\mathrm{W}} $, ${\mathrm{H}} \rightarrow {\tau } {\tau } $, or ${\mathrm{H}} \rightarrow {\mathrm{Z}} {\mathrm{Z}} $ decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among ${\mathrm{t}} {\mathrm{H}} $, ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ and ${\mathrm{t}} {\mathrm{H}} $ signals correspond to $0.92 \pm 0.19\,\text {(stat)} ^{+0.17}_{-0.13}\,\text {(syst)} $ and $5.7 \pm 2.7\,\text {(stat)} \pm 3.0\,\text {(syst)} $ of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and to 1.4 (0.3) for ${\mathrm{t}} {\mathrm{H}} $ production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $y_{{\mathrm{t}}}$ of the Higgs boson to the top quark divided by its SM expectation, $\kappa _{{\mathrm{t}}}=y_{{\mathrm{t}}}/y_{{\mathrm{t}}}^{\mathrm {SM}}$, is constrained to be within $-0.9< \kappa _{{\mathrm{t}}}< -0.7$ or $0.7< \kappa _{{\mathrm{t}}}< 1.1$, at 95% confidence level. This result is the most sensitive measurement of the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ production rate to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2001

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗

DOE OSTI2023

Analysis of Unpolarised p+p¿ Photoproduction with the GlueX Experiment

This thesis presents measurements of spin-density matrix elements in unpolarised p+p? photoproduction on a proton target. The dominant resonance contribution to the dipion system is the r(770) meson. Due to the large production cross section for this resonance, a valid comparison can be made between the obtained final results and r(770) spin-density matrix elements measured previously with other experiments. The measurement was performed over the 3:0 ? 11:6 GeV beam energy regime, which is a more extensive energy range than has ever been studied previously for the r(770). Results were obtained by analysing data from the GlueX experiment based at Jefferson Lab. Extended maximum likelihood fits were applied to extract three spin-density matrix elements using Markov chain Monte Carlo based parameter estimations. This was performed using various binning configurations to probe the energy, mass, and four-momentum transfer dependence of the determined physics observables. Spin-density matrix elements are shown to be consistent with the model of s-channel helicity conservation at low ?t. The effects of pomeron and f2 exchanges are clearly visible in the energy dependence of the measured observables. These observations provide valuable insights into the relative strengths of both processes as a function of the photon energy, and may enable theorists to disentangle the f2=P coupling ratio. Spin-density matrix elements are seen to be highly dependent on the reconstructed resonance mass. This observation is likely to be a result of non-resonant S-wave background processes, and emphasises the need for a more detailed model of the p+p? angular distribution that considers all of the competing angular momentum components that contribute to the measured final state. The statistical precision of measurements performed for this thesis surpass what was achievable in previous studies of the r(770) by several orders of magnitude. Studies of the energy and four-momentum transfer dependence, and insights into the effects of non r(770) background contributions provide valuable input for production models. This will help inform the choice of wave-sets used for partial wave analyses, supporting GlueX in its search for exotic hybrid meson states.

Fitches, James↗

DOE OSTI2023

Observations of ozone, acyl peroxy nitrates, and their precursors during summer 2019 at Carlsbad Caverns National Park, New Mexico

Carlsbad Caverns National Park (CAVE) is located in southeastern New Mexico and is adjacent to the Permian Basin, one of the most productive oil and natural gas (O&G) production regions in the United States. Since 2018, ozone (O 3 ) at CAVE has frequently exceeded the 70 ppbv 8-hour National Ambient Air Quality Standard. We examine the influence of regional emissions on O 3 formation using observations of O 3 , nitrogen oxides (NO x = NO + NO 2 ), a suite of volatile organic compounds (VOCs), peroxyacetyl nitrate (PAN), and peroxypropionyl nitrate (PPN). Elevated O 3 and its precursors are observed when the wind is from the southeast, the direction of the Permian Basin. We identify 13 days during the July 25 to September 5, 2019 study period when the maximum daily 8-hour average (MDA8) O 3 exceeded 65 ppbv; MDA8 O 3 exceeded 70 ppbv on 5 of these days. The results of a positive matrix factorization (PMF) analysis are used to identify and attribute source contributions of VOCs and NO x . On days when the winds are from the southeast, there are larger contributions from factors associated with primary O&G emissions; and, on high O 3 days, there is more contribution from factors associated with secondary photochemical processing of O&G emissions. The observed ratio of VOCs to NOx is consistently high throughout the study period, consistent with NO x -limited O 3 production. Finally, all high O 3 days coincide with elevated acyl peroxy nitrate abundances with PPN to PAN ratios > 0.15 ppbv ppbv -1 indicating that anthropogenic VOC precursors, and often alkanes specifically, dominate the photochemistry.

54 ENVIRONMENTAL SCIENCES↗