Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rare event”

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 163 records · Page 9

MOOSE ProbML: Parallelized probabilistic machine learning and uncertainty quantification for computational energy applications

Here, this paper presents the development and demonstration of massively parallel probabilistic machine learning (ML) and uncertainty quantification (UQ) capabilities within the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source computational platform for parallel finite element and finite volume analyses. In addressing the computational expense and uncertainties inherent in complex multiphysics simulations, this paper integrates Gaussian process (GP) variants, active learning, Bayesian inverse UQ, adaptive forward UQ, Bayesian optimization, evolutionary optimization, and Markov chain Monte Carlo (MCMC) within MOOSE. It also elaborates on the interaction among key MOOSE systems — Sampler, MultiApp, Reporter, and Surrogate — in enabling these capabilities. The modularity offered by these systems enables development of a multitude of probabilistic ML and UQ algorithms in MOOSE. Example code demonstrations include parallel active learning and parallel Bayesian inference via active learning. The impact of these developments is illustrated through five applications relevant to computational energy applications: UQ of nuclear fuel fission product release, using parallel active learning Bayesian inference; very rare events analysis in nuclear microreactors using active learning; advanced manufacturing process modeling using multi-output GPs (MOGPs) and dimensionality reduction; fluid flow using deep GPs (DGPs); and tritium transport model parameter optimization for fusion energy, using batch Bayesian optimization. These capabilities are part of the MOOSE framework.

97 - MATHEMATICS AND COMPUTING↗

Ultra-low radioactivity Kapton and copper-Kapton laminates

Polyimide-based materials, like Kapton, are widely used in flexible cables and circuitry due to their unique electrical and mechanical characteristics. This study is aimed at investigating the radiopurity of Kapton for use in ultralow background, rare-event physics applications by measuring the 238U, 232Th, and natK levels using inductively coupled plasma mass spectrometry. Commercial-off-the-shelf Kapton varieties, mea- sured at approximately 950 and 120 pg/g 238U and 232Th (1.2×104 and 490 µBq/kg), respectively, can be a significant background source for many current and next-generation ultralow background detectors. This study has found that the dominant contamination is due to the use of dicalcium phosphate (DCP), a nonessential slip additive added during manufacturing. Alternative Kapton materials were obtained that did not contain DCP and were determined to be significantly more radiopure than the commercially- available options with 12 and 19 pg/g 238U and 232Th (150 and 77 µBq/kg), respectively. The lowest radioactivity version produced (Kapton ELJ, which contains an adhesive) was found to contain single digit pg/g levels of 238U and 232Th, two-to-three orders of magnitude cleaner than commercial-off-the- shelf options. Moreover, copper-clad polyimide laminates employing Kapton ELJ as the insulator were obtained and measured and showed to be very radiopure at 8.6 and 22 pg/g 238U and 232Th (110 and 89 µBq/kg), respectively.

Arnquist, Isaac J.↗

Simulation of energy loss of fractionally charged particles using Geant4

Over the years several rare-event search experiments have shown increasing interest to search for Fractionally Charged Particles (FCPs) within the data acquired in those experiments. This, in turn, has required the need for a framework to calculate accurately the energy loss expected from FCPs as they pass through the detector material. We have developed Monte Carlo methods within Geant4 to simulate for the first time energy loss of FCPs through different physics processes. The program can be used to calculate energy deposition distributions, a major component of any FCP search analysis, as they interact with detector material. We discuss the implementation and validation of the energy loss of FCPs using a typical detector material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Removal of 210 Pb by etch of crystalline detector sidewalls

A potential source of dominant backgrounds for many rare-event searches or screening detectors is from radon daughters, specifically 210 Pb, deposited on detector surfaces, often during detector fabrication. Performing a late-stage etch is challenging because it may damage the detector. Here, this paper describes a late-stage etching technique that reduces surface 210 Pb and 210 Po by > 99x at 90% C.L.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Temperature-dependent charge barrier height of amorphous germanium contact detector

The exploration of germanium (Ge) detectors with amorphous Ge (a-Ge) contacts has drawn attention to the searches for rare-event physics such as dark matter and neutrinoless double-beta decay. The charge barrier height (CBH) of the a-Ge contacts deposited on the detector surface is crucial to suppress the leakage current of the detector in order to achieve a low-energy detection threshold and high-energy resolution. The temperature-dependent CBH of a-Ge contacts for three Ge detectors is analyzed to study the bulk leakage current (BLC) characteristics. The detectors were fabricated at the University of South Dakota using homegrown crystals. The CBH is determined from the BLC when the detectors are operated in the reverse bias mode with a guard-ring structure, which separates the BLC from the surface leakage current (SLC). Here the results show that CBH is temperature dependent. The direct relation of the CBH variation to temperature is related to the barrier inhomogeneities created on the interface of a-Ge and crystalline Ge. The inhomogeneities that occur at the interface were analyzed using the Gaussian distribution model for three detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Understanding the enhancement of scintillation light in xenon-doped liquid argon

Measuring the scintillation light in noble gases is an important detection technique in particle physics. Numerous rare event searches like neutrino beam experiments, neutrino-less double beta-decay, and dark matter searches use argon-based detectors. In liquid argon, the light yield can be enhanced by the addition of a small quantity of xenon, where ~10–1000 ppm are added. The general enhancement mechanism and its pathway via an energy transfer between argon and xenon excimers is well known, however the importance of absorption of argon excimer emission by atomic xenon has not been fully appreciated. This absorption significantly reduces the light yield in commercially available argon (extracted from air) which contains trace amounts (~0.1 ppm) of xenon. The addition of a small xenon dopant of ~10 ppm recovers this lost light resulting in an increased light yield over un-doped argon of about a factor of two. Here, in this paper, we introduce a model for the light production in xenon doped argon, including absorption and re-emission, and compare it to the measured time dependence of light emission in xenon-doped argon.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An integrated online radioassay data storage and analytics tool for nEXO

Large-scale low-background detectors are increasingly used in rare-event searches as experimental collaborations push for enhanced sensitivity. However, building such detectors, in practice, creates an abundance of radioassay data especially during the conceptual phase of an experiment when hundreds of materials are screened for radiopurity. A tool is needed to manage and make use of the radioassay screening data to quantitatively assess detector design options. We have developed a Materials Database Application for the nEXO experiment to serve this purpose. Furthermore, this paper describes this database application, explains how it functions, and discusses how it streamlines the design of the experiment.

47 OTHER INSTRUMENTATION↗

Evaluation of SNOLAB background mitigation procedures through the use of an ICP-MS based dust monitoring methodology

Dust particulate fallout on materials in use for rare-event searches is a concerning source of radioactive backgrounds due to the presence of the naturally occurring radionuclides 40 K, 232 Th, 238 U, and their progeny in dust. Much effort is dedicated to inform radioactive backgrounds from dust and evaluate the efficacy of mitigation procedures. A great portion of such effort relies on fallout models and assumed dust composition. Here, in this work, an ICP-MS based methodology was employed for a direct determination of fallout rates of radionuclides and stable isotopes of interest from dust particulate at the SNOLAB facility. Hosted in an active mine, the SNOLAB underground laboratory strives to maintain experimental areas at class 2000 cleanroom level. This work validates the mitigation procedures in place at SNOLAB and informs dust backgrounds during laboratory activities. Fallout rates of major constituents of the local rock were measured two to three orders of magnitude lower in the clean experimental areas compared to non-clean transition areas from the mine to the laboratory. An increase of approximately two orders of magnitude in the fallout rate of stable Pb isotopes was determined in an experimental area during activities involving handling of Pb bricks. Increased 40 K, 232 Th, and 238 U fallout rates were measured in clean experimental areas during activities generating particulate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗

Classification and computation of extreme events in turbulent combustion

In the design of practical combustion systems, ensuring safety and reliability is an important requirement. For instance, reliably avoiding lean blowout, flame flashback or inlet unstart is critical for ensuring safe operation. Currently, the science of predicting such events is based on prior experience, limited modeling or diagnostic tools and purely statistical approaches. Even though computational and experimental tools for studying combustion devices have vastly advanced in the last three decades, the analysis of such failure events has not been pursued widely. While the use of data for model development and calibration is being widely accepted, the extension to failure events introduces numerous challenges. In particular, the focus here is on so-called data-poor problems, where the cost of generating data is extremely high and is not easily amenable to existing computational and experimental approaches. Data-poor problems are particularly relevant when related to extreme events (also called anomalous events) that can lead to catastrophic failure of the system. It is argued that transient events that describe such failure can have different causal mechanisms. To develop the scientific inference process, a classification of such problems is used to determine specific modeling paths as well as computational tools needed. Research opportunities in the emerging field of extreme event prediction are highlighted in order to identify critical and immediate needs.

97 MATHEMATICS AND COMPUTING↗

Quantifying the Relevance of Long-Range Forces for Crystal Nucleation in Water

Understanding nucleation from aqueous solutions is of fundamental importance in a multitude of fields, ranging from materials science to biophysics. The complex solvent-mediated interactions in aqueous solutions hamper the development of a simple physical picture, elucidating the roles of different interactions in nucleation processes. In this study, we make use of three complementary techniques to disentangle the role played by short- and long-range interactions in solvent-mediated nucleation. Specifically, the first approach we utilize is the local molecular field (LMF) theory to renormalize long-range Coulomb electrostatics. Second, we use well-tempered metadynamics to speed up rare events governed by short-range interactions. Third, the deep learning-based State Predictive Information Bottleneck approach is employed in analyzing the reaction coordinate of the nucleation processes obtained from the LMF treatment coupled with well-tempered metadynamics. We find that the two-step nucleation mechanism can largely be captured by the short-range interactions, while the long-range interactions further contribute to the stability of the primary crystal state under ambient conditions. Furthermore, by analyzing the reaction coordinate obtained from the combined LMF-metadynamics treatment, we discern the fluctuations on different time scales, highlighting the need for long-range interactions when accounting for metastability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Solvation and Acidity of Carboxylic Acids Using an Ab Initio Deep Neural Network Potential

Formic and acetic acid constitute the simplest of carboxylic acids, yet they exhibit fascinating chemistry in the condensed phase such as proton transfer and dimerization. The go-to method of choice for modeling these rare events have been accurate but expensive ab-initio molecular dynamics simulations. Here, we present a deep neural network potential trained using accurate ab-initio data that can be used in tandem with enhanced-sampling methods to perform an efficient exploration of the free-energy surface of aqueous solutions of weak carboxylic acids. In particular, we show that our model captures proton dissociation and provides a good estimate of the pK a , as well as the dimerization of formic and acetic acid. This provides a suitable starting point for applications in different research areas where computational efficiency coupled with the accuracy of ab-initio methods is required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin–Landau–Verwey–Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na + , K + , and NO 3 –) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (~3–4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal–basal surfaces is more favorable than between (001) edge–edge surfaces (~5–6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (~5–10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. Further, the results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cu–Ni Oxidation Mechanism Unveiled: A Machine Learning-Accelerated First-Principles and in Situ TEM Study

Here, the development of accurate methods for determining how alloy surfaces spontaneously restructure under reactive and corrosive environments is a key, long-standing, grand challenge in materials science. Using machine learning-accelerated density functional theory and rare-event methods, in conjunction with in situ environmental transmission electron microscopy (ETEM), we examine the interplay between surface reconstructions and preferential segregation tendencies of CuNi(100) surfaces under oxidation conditions. Our modeling approach predicts that oxygen-induced Ni segregation in CuNi alloys favors Cu(100)-O c(2 × 2) reconstruction and destabilizes the Cu(100)-O (2√2 × √2)R45° missing row reconstruction (MRR). In situ ETEM experiments validate these predictions and show Ni segregation followed by NiO nucleation and growth in regions without MRR, with secondary nucleation and growth of Cu 2 O in MRR regions. Our approach based on combining disparate computational components and in situ ETEM provides a holistic description of the oxidation mechanism in CuNi, which applies to other alloy systems.

36 MATERIALS SCIENCE↗

Antimicrobial Peptides Can Facilitate Whole Blood Safety from Bacteria: A Proof of Concept

With continuous improvements to blood donor deferrals and the availability of sensitive tests for donation screening for infectious agents, bacterial contamination of whole blood (WB) and blood components stored for transfusion is a rare event. Nonetheless, it still occurs and remains a transfusion-associated risk in terms of septic transfusion reactions (STRs) and transfusion-transmitted bacterial infections with morbidity and mortality outcomes. One of the risk mitigation strategies for bacterial contamination is to implement treatment with currently available proactive pathogen reduction technologies (PRTs) for these transfusion products. Here, as a proof of concept, we tested two recently developed unique cationic antimicrobial peptides (AMPs; D-CONGA and D-CONGA-Q7) for WB safety from bacterial contamination. In this study, WB was inoculated with Escherichia coli and Staphylococcus epidermidis and treated with the two peptides to evaluate their bactericidal efficacy. The results demonstrated that D-CONGA and D-CONGA-Q7 exhibit potent inhibitory activity against the bacteria with a minimal inhibitory concentration (MIC) range of 4–8 and 1–8 μM, respectively, depending on the bacterial species tested. Time-kill kinetics further confirmed that the peptides exhibit bactericidal efficacy at 8 μM by achieving a 5-log 10 reduction (99.999%) of the bacterial load in WB with a time-dependent killing profile. Furthermore, even at 20 μM, the AMPs did not negatively impact hemolysis or hemostatic properties. We have further demonstrated using a cationic exchange resin that the cationic AMPs can be separated and removed from WB after the peptide treatments. During 35-day WB storage at 2–8 °C, 4 μM D-CONGA-Q7 one-time treatment prevented S. epidermidis growth and preserved WB quality and integrity. Overall, the results described here provide the first proof of concept that certain AMPs, such as D-CONGA and D-CONGA-Q7, can facilitate WB safety from bacteria during storage.

60 APPLIED LIFE SCIENCES↗

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗