Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetically constrained models”

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 73 records · Page 4

Searching for Dark Matter with Superconducting Qubits

Dark matter is an enduring mystery in our quest to understand the fundamental constituents of our universe. Low mass bosons, such the axion or hidden photon, are compelling dark matter candidates. We leverage their potential interactions with electromagnetic fields, whereby the dark matter (of unknown mass) on rare occasion converts into a single photon, to devise a method of detecting these candidates. Current dark matter searches operating at microwave frequencies use a resonant cavity to coherently accumulate the field sourced by the dark matter and a near standard quantum limited (SQL) linear amplifier to read out the cavity signal. To further increase sensitivity to the dark matter signal and enable future searches, sub-SQL detection techniques are required.In this thesis, I report the development of a novel microwave photon counting technique and a new exclusion limit on hidden photon dark matter. We operate a superconducting qubit to make repeated quantum non-demolition measurements of cavity photons and apply a hidden Markov model analysis to reduce the noise to 15.7 dB below the quantum limit, with overall detector performance limited by a residual background of real photons. With the present device, we perform a hidden photon search and constrain the kinetic mixing angle to less than 1.68 x 10-15 in a band around 6.011 GHz (24.86 micro-eV) with an integration time of 8.33 s. This demonstrated noise reduction technique enables future dark matter searches to be sped up by a factor of 1300. By coupling a qubit to an arbitrary quantum sensor, more general sub-SQL metrology is possible with the techniques presented in this work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measuring the oscillator strength of intercombination lines of helium-like V ions in a laser-produced-plasma

In this work, we present results of measurements of the oscillator strength of intercombination lines of He-like Vanadium ions in high energy density (HED) laser-produced-plasmas and compare them with the simulations from commonly used codes and data from the NIST database. Whilst not yet sufficiently accurate to constrain different trusted atomic-physics models for the particular system studied, our results are in agreement with the available data within experimental error bars, yet differ from cruder approximations of the oscillator strength used in certain atomic-kinetics packages, suggesting that this general method could be further extended to be used as a measurement of the oscillator strength of additional atomic transitions under the extreme conditions that are achieved in HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement of the muon neutrino charged-current mesonless cross section in the NOvA near detector

NOvA is a long-baseline accelerator neutrino experiment at Fermilab. Its physics goals include precision neutrino oscillation measurements, neutrino interaction cross-section measurements and beyond Standard Model explorations. We present a measurement of muon neutrino charged-current cross section with zero mesons in the final state at the NOvA near detector. This measurement is performed as a function of the kinematics of the final state muon. Our chosen interaction channel is especially sensitive to quasielastic and meson exchange current interactions and it provides handles for constraining the cross section systematic uncertainties in oscillation analyses in present and future experiments. For particle identification, we use a convolutional neural network (CNN) trained on individual particles simulated in the NOvA Near detector. This allows us to select the desired signal while reducing the potential bias from neutrino interaction modeling. We study strategies for constraining the remaining charged-pion background via Michel electron information in a template fitting approach. The main experimental result is a two-dimensional differential cross section as a function of final-state muon kinetic energy and polar angle. The parameters of this measurement, including binning and unfolding, were optimized to reduce the expected systematic uncertainty in the total cross section. The final result shows good agreement with the main GENIE-based simulation framework that was independently fine-tuned in NOvA. We finally propose improvements and subsequent steps that build on this analysis and further dissect the final states of neutrino interactions. This work has been supported by US DOE grant DE-SC0015684.

Sánchez Falero, Sebastián Jesús [Iowa State U.]↗

Institutional Computing: Annual Progress Report, Precision early universe simulations to constrain nuclear reactions and beyond standard model physics

The current one-year allocation (w20 qburst) is the first year in a two year allocation of the same name. The project has supported two publications in peer-review journals, another manuscript in preparation for publication, a white paper for the Astro2020 Decadal Survey, and a number of talks (see below). Code development continues on implementing a neutrino quantum kinetic equation solver into BURST, now called QBURST. While implementing coherent and forward-scattering processes into the dynamical neutrino evolution equations, multiple numerical issues have appeared. We have remedied a number of these issues through extensive debugging and testing, which include streamlining interpolation and implementing new ODE integration techniques. Further testing is required to solve the problem at a precision of better than one part in 10 6 . We expect that we will understand the numerical difficulties and implement solutions by end of year 2021.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Kinetic-Energy–Momentum–Mass 5-Flux of a Baryon Fluid in Bargmann Spacetimes

A Bargmann spacetime is a constrained five-dimensional setting that, while introducing no new physical degrees of freedom beyond those of ordinary four-dimensional spacetime, permits Galilei physics to be expressed with a tensor formalism that respects the distinction between mass and energy while affording the conceptual and technical advantages of a spacetime metric. This framework offers a route to a strong-field ‘Galilei general relativity’ approximating the usual Poincaré general relativity introduced by Einstein. In preparation for modeling core-collapse supernovae, where such an approximation would be useful, this work generalizes the kinetic-energy–momentum–mass 5-flux 𝒯 and its associated spacetime tensor law from a simple fluid of constant particle mass to a baryon fluid whose multiple nuclear species can interconvert rest mass and internal energy. The spacetime tensor law on Bargmann–Galilei spacetime 𝐵𝒢 and its decompositions relative to comoving (‘Lagrangian’) and fiducial (‘Eulerian’) observers are derived in detail. The formalism is rendered more suitable for core-collapse supernova modeling by an extension from strict 𝐵𝒢 to a regime that might be denoted as 𝐵𝒢+: microscopically Poincaré yet macroscopically Galilei. This extension accommodates energy generation by nuclear composition changes and allows comoving energy density and pressure to contribute relative to mass density, while preserving the simplifications of Galilei bulk fluid flow and the streamlined geometry governed by the Bargmann–Galilei spacetime metric.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗

Particle acceleration during classical phase transitions on a spherical lattice

Abstract When compressed, certain lattices undergo phase transitions that may allow nuclei to gain significant kinetic energy. To explore the dynamics of this phenomenon, we develop a methodology to study Coulomb coupled N -body systems constrained to a sphere, as in the Thomson problem. We initialize N total Boron nuclei as point particles on the surface of the sphere, allowing them to equilibrate via Coulomb scattering with a viscous damping term. To simulate a phase transition, we remove N rm particles, forcing the system to rearrange into a new equilibrium. With this model, we consider the Thomson problem as a dynamical system, providing a framework to explore how non-zero temperature affects structural imperfections in Thomson minima. We develop a scaling relation for the average peak kinetic energy attained by a single particle as a function of N and N rm . For certain values of N , we find an order of magnitude energy gain when increasing N rm from 1 to 6. The model may help to design a lattice that maximizes the energy output.

Bachmann, Aidan M. (ORCID:0000000294706404)↗

Microkinetic Modeling of the Transient CO 2 Methanation with DFT–Based Uncertainties in a Berty Reactor

The transient operation of methanation reactors can become desirable when coupled with fluctuating renewable energies in Power-to-Gas scenarios, which requires suitable kinetic approaches that can describe the transient catalytic phenomena. A combined experimental and theoretical investigation of the transient CO 2 methanation is conducted using concentration forcing to derive a suitable microkinetic model. Methanation experiments are performed with a Ni/SiO 2 catalyst in a Berty-type reactor at industrially relevant conditions. The microkinetics are based on previous work and were automatically constructed for the Ni(111) facet using the Reaction Mechanism Generator. A feasible set of energetic parameters of the microkinetic models was identified in a theory-constrained optimization procedure within the DFT uncertainty space that can accurately reproduce the experimental results on a first-principles basis. The microkinetic model unravels that the formation of H 2 O* and CH 3 * control the activity and selectivity of Ni(111) under the investigated conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum

Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.

36 MATERIALS SCIENCE↗

A dynamic kinetic model captures cell-free metabolism for improved butanol production

Cell-free systems are useful tools for prototyping metabolic pathways and optimizing the production of various bioproducts. Mechanistically-based kinetic models are uniquely suited to analyze dynamic experimental data collected from cell-free systems and provide vital qualitative insight. However, to date, dynamic kinetic models have not been applied with rigorous biological constraints or trained on adequate experimental data to the degree that they would give high confidence in predictions and broadly demonstrate the potential for widespread use of such kinetic models. In this work, we construct a large-scale dynamic model of cell-free metabolism with the goal of understanding and optimizing butanol production in a cell-free system. Using a combination of parameterization methods, the resultant model captures experimental metabolite measurements across two experimental conditions for nine metabolites at timepoints between 0 and 24 hours. Here, we present analysis of the model predictions, provide recommendations for butanol optimization, and identify the aldehyde/alcohol dehydrogenase as the primary bottleneck in butanol production. Sensitivity analysis further reveals the extent to which various parameters are constrained, and our approach for probing valid parameter ranges can be applied to other modeling efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First Measurement of $π^+$-Ar and $p$-Ar Total Inelastic Cross Sections in the Sub-GeV Energy Regime with ProtoDUNE-SP Data

The ProtoDUNE-SP detector, a kiloton-scale prototype for the Deep Underground Neutrino Experiment (DUNE), is the largest liquid argon time projection chamber built to date. Operated at CERN from 2018 to 2020, it collected both cosmic-ray data and a beam consisting of positively-charged particles with discrete momentum settings across a range of 0.3 GeV/$c$ to 7 GeV/$c$. In this letter, we report the total inelastic cross section measurements for $π^+$-Ar and $p$-Ar interactions using selected $π^+$ and proton samples from the 1 GeV/$c$ beam data. These results provide the first measurement of the total inelastic cross sections for $π^+$-Ar in the 500-900 MeV kinetic energy range and for $p$-Ar below 450 MeV, both of which are directly relevant to the DUNE energy range. The measured cross sections are consistent with predictions and provide a dataset that was previously unavailable for argon targets. These measurements are essential for constraining neutrino-argon interaction models, which are crucial for the precision physics goals of the upcoming DUNE experiment.

Abbaslu, Saeed [IPM, Tehran] (ORCID:00000003356771↗

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 μs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Control of Generic Second Order Traffic Models by Driver-Assist Vehicles

We study the derivation of generic high order macroscopic traffic models from a follow-the-leader particle description via a kinetic approach. First, we recover a third order traffic model as the hydrodynamic limit of an Enskog-type kinetic equation. Next, we introduce in the vehicle interactions a binary control modeling the automatic feedback provided by driver-assist vehicles and we upscale such a new particle description by means of another Enskog-based hydrodynamic limit. The resulting macroscopic model is now a generic second order model (GSOM), which contains in turn a control term inherited from the microscopic interactions. We show that such a control may be chosen so as to optimize global traffic trends, such as the vehicle flux or the road congestion, constrained by the GSOM dynamics. By means of numerical simulations, we investigate the effect of this control hierarchy in some specific case studies, which exemplify the multiscale path from the vehiclewise implementation of a driver-assist control to its optimal hydrodynamic design.

GSOM↗

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An Assessment of Kinetic Models for Ammonia Flame Extinction

The growing interest in using ammonia as a carbon-neutral replacement fuel has prompted numerous recent research efforts towards understanding fundamental combustion characteristics and chemical kinetics of ammonia oxidation. Despite these recent efforts, there remains a large deficiency in extinction limit data for ammonia flames, particularly under elevated pressure and temperature conditions relevant to gas turbines. Even the literature reaction mechanisms developed specifically for ammonia have been untested under gas turbine relevant conditions, and existing kinetic models have shown large variations in laminar flame speed, extinction strain rate, and speciation predictions for NH3–H2 mixtures at ambient conditions. Recent work by Thomas and co-workers, for example, compared different ammonia kinetic models in the literature for predicting their extinction strain rate measurements for non-premixed NH3-H2 counterflow flames, and again have highlighted the discrepancies in model predictions. Here, our work has undertaken to further highlight important reaction pathways for ammonia oxidation that may be responsible for these discrepancies for extinction strain rate predictions, as well as elucidate sensitive reaction steps with large uncertainties that require further attention. This work points to a need for an improved understanding of ammonia combustion chemistry as well as new datasets at elevated pressure and temperature conditions to better constrain rate expressions for key reaction steps.

ammonia, flame, combustion, kinetic, kinetics↗

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗