SEARCH · Engineering Papers
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.
Where does stuff come from? Observations of the sun one major way to measure abundances (but there are many others). Stars with not too many heavy elements show an interesting, repeatable pattern, which shows up in our sun too!
The main scientific contributions of this project are the following novel concepts for multidimensional arrays: shape-based similarity join (SIGMOD 2016), incremental view maintenance (SIGMOD 2017), user-defined stencil functions (HPDC 2017), and distributed caching for in-situ processing (SSDBM 2018). Building on our collaboration with the astrophysics group at LBNL, we applied these techniques to the data generated in the Palomar Transient Factory (PTF) astronomical survey. They played a pivotal role in the first-ever observation of a neutron star merger, which produces gravitational waves and turns out to be the origin of heavy elements, including gold. This has lead to a Science magazine article that has received extensive media coverage on ACM TechNews, Slashdot, FiveThirtyEight, and Quanta Magazine, among others. Additionally, two other articles detailing related aspects of the same discovery have been published in the Astrophysical Journal Letters journal. These publications have more than 3,000 citations according to Google Scholar (as of February 2022). This cross-disciplinary collaboration provided very good opportunities to apply database techniques to real-life scientific problems. The fact that they facilitated major discoveries in astrophysics proves the importance of our research. In addition to the work on multidimensional array databases, this project has also developed stochastic gradient descent (SGD) optimization algorithms for training large scale machine learning models, methods for querying in-situ data, and a database query optimizer based on sketch synopses.
The goal of this project was to further our theoretical understanding of the extreme physics and heavy element nucleosynthesis occurring in stellar explosions, in particular supernovae and neutron star mergers. To study these systems, we further developed and applied numerical simulations (radiation hydrodynamical models and nuclear reaction networks) run on high performance computers. The numerical calculations and companion theoretical analysis helped to better define the sites of astrophysical nucleosynthesis, and make predictions of the associated observables (nucleosynthetic yields and electromagnetic signals) that can be used to anticipate and interpret experimental data. This project enhances the scientific value of the nuclear data obtained from current and future experimental facilities (e.g., rare isotope beams) by improving our understanding of the astrophysical context and specific physical conditions under which we expect such reactions to take place in Nature.
Density and electron fraction of neutron star stuff leftover after a black hole eats the neutron star. Events like this one are believed to power gamma ray bursts, some of the most energetic events in the universe, and be the source of heavy elements like gold and platinum.
Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.
Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.
Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.
The Radiochemical Engineering Development Center (REDC) generates various liquid organic wastes from processing irradiated targets to recover heavy elements. In the past, these organics were discharged to the Oak Ridge National Laboratory (ORNL) liquid low-level waste (LLLW) system, along with the aqueous waste. Because of a reduction in aqueous LLLW from other ORNL generators and an increase in the radionuclide concentration, particularly 238 Pu, in the REDC organic waste stream, the organics can no longer be discharged to the LLLW system. The plan is to solidify the liquid organic waste for disposal as solid waste. Solidifying the organic liquids using PM-199 Organoclay ® would produce solid waste forms with no free liquid, which should qualify for disposal at the Waste Isolation Pilot Plant. Granular Organoclay can be added to the organic waste solutions until there is a dry layer on top of the liquid and then allowed to cure for a few hours to produce a solid wasteform.
This document is the final report for the project with the title, “Uranyl Capture with Lewis Acids and Macrocyclic Hosts.” Nuclear power is attractive for meeting the current and future energy needs of society, in that it does not release carbon dioxide or other pollutants into the atmosphere during routine use. This motivated the project, because preparation of nuclear fuel, recovery of useful components from used fuel, and handling of waste materials remain significant impediments to further deployment of important nuclear technologies. In part, these problems arise from limited availability of chemical reactions that can reliably interconvert forms of uranium and other heavy elements during preparation and processing. For example, harsh and expensive chemicals are often required for making and breaking chemical bonds to uranium, and the reactions involved are inefficient. The overall objective of this research was to harness knowledge of chemical structure and bonding to develop a useful and predictive und
Metal halide perovskite (MHP) multiple quantum wells which consist of multilayers of alternate organic and inorganic layers exhibit large exciton binding energies due to the dielectric confinement between the inorganic and organic layers. These naturally formed multiple quantum wells have strong spin-orbit coupling (SOC) due to the presence of heavy elements in their crystal structures. Although the fundamental properties of 2D MHPs are far from being entirely understood, it is widely accepted that their band edge absorption coefficient results from strong exciton interactions. However, studies demonstrating how different exciton interactions and doping effects influence electronic traps and disorder on the band edge absorption coefficient of 2D MHPs have not been demonstrated. Understanding these interactions in MHPs will allow us to access low energy optical transitions for the fabrication of solution processable short-to-mid-wavelength IR photodetectors (1 – 8 μm). Moreover, upon doping, it is possible to move the Fermi energy into the conduction band (CB) to favorably promote the transport of charges in a working device. Herein, we study the development of 2D MHPs having strong SOC, high carrier mobility, and tunable quantum well structures. Our studies shed light on the design and modulation of fundamental physical phenomena by carefully elucidating the role of dopants (n-type and p-type), exciton heterogeneity, orientation, structure, and bias stress effects on the performance of MHPs as potential IR photodetectors.
The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.
For applications related to nuclear security, safeguards, and nonproliferation, it is often critical to know the molecular compositions of lanthanide- and actinide-containing samples. Spectroscopy is a widely used tool that looks at the interaction between light and matter: Different species absorb or emit light at unique wavelengths which act as signatures. However, there is a limited number of tools that can achieve high-sensitivity, accurate measurements of lanthanide and actinide molecular compositions. Candidate methods include mass spectrometry, which usually destroys at least part of the sample and requires complicated stoichiometry to guess the original sample’s molecular compositions; optical spectroscopies, which have great atomic but limited molecular sensitivities or other drawbacks which make sensing molecules difficult like limited light sources or strong absorption in the atmosphere; and nuclear spectroscopies (gamma, neutron) which also have limited sources and long (>minute) collection times. As such, the purpose of our research is to develop a new tool to better distinguish between subtle differences in molecules containing lanthanides and actinides. Soft x-ray spectroscopy is sensitive to molecular form and is minimally intrusive/nondestructive to the sample. However, soft x-ray light with sufficient brightness for spectroscopy is typically limited to user-facilities like synchrotrons or free electron lasers, where beamtimes are competitive, and work with radiological materials may be difficult or entirely prohibited. To overcome this issue, our Team has developed a custom tabletop laser-driven, soft x-ray light source which employs high harmonic generation (HHG). Soft x-ray spectroscopy can distinguish between subtly different molecules, in the spectral range which we need to study these heavy elements. A tabletop system provides an effective and affordable tool to find both the elemental and chemical specificity of lanthanide and samples. Creating a light source in the soft x-ray spectrum is difficult because these wavelengths in the range 5-20 nm (20-350 eV photon energies) only penetrate several 100s of nm in most solid materials and only reflect well in shallow, grazing incident angles. The results are applicable to nuclear forensics, because molecular fingerprinting of lanthanide and actinide samples can be used to back out the origin and processing method of nuclear materials (Skrodzki, et al.).
In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.
From the lightest Hydrogen isotopes up to the recently synthesized Oganesson (Z = 118), it is estimated that as many as about 8,000 atomic nuclei could exist in nature. Most of these nuclei are too short-lived to be occurring on Earth, but they play an essential role in astrophysical events such as supernova explosions or neutron star mergers that are presumed to be at the origin of most heavy elements in the Universe. Understanding the structure, reactions, and decays of nuclei across the entire chart of nuclides is an enormous challenge because of the experimental difficulties in measuring properties of interest in such fleeting objects and the theoretical and computational issues of simulating strongly-interacting quantum many-body systems. Nuclear density functional theory (DFT) is a fully microscopic theoretical framework which has the potential of providing such a quantitatively accurate description of nuclear properties for every nucleus in the chart of nuclides. Thanks to high-performance computing facilities, it has already been successfully applied to predict nuclear masses, global patterns of radioactive decay like β or γ decay, and several aspects of the nuclear fission process such as, e.g., spontaneous fission half-lives. Yet, predictive simulations of nuclear spectroscopy—the low-lying excited states and transitions between them—or of nuclear fission, or the quantification of theoretical uncertainties and their propagation to basic or applied nuclear science applications, would require several orders of magnitude more calculations than currently possible. However, most of this computational effort would be spent into generating a suitable basis of DFT wavefunctions. Such a task could potentially be considerably accelerated by borrowing tools from the field of machine learning and artificial intelligence. In this paper, we review different approaches to applying supervised and unsupervised learning techniques to nuclear DFT.
Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks, producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down star formation and ultimately curtailing the growth of structure after the peak at redshift 2–3. To understand the complex interplay between gravity and feedback, we must resolve both the key physics within galaxies and map the impact of these processes over large scales, out into the cosmic web. The Advanced X-ray Imaging Satellite (AXIS) is a proposed X-ray probe mission for the 2030s with arcsecond spatial resolution, large effective area, and low background. AXIS will untangle the interactions of winds, radiation, jets, and supernovae with the surrounding interstellar medium across the wide range of mass scales and large volumes driving galaxy evolution and trace the establishment of feedback back to the main event at cosmic noon. This white paper is part of a series commissioned for the AXIS Probe mission concept; additional AXIS white papers can be found at the AXIS website.
Radial velocity (RV) surveys have discovered giant exoplanets on au-scale orbits with a broad distribution of eccentricities. Those with the most eccentric orbits are valuable laboratories for testing theories of high-eccentricity migration. However, few such exoplanets transit their host stars, thus removing the ability to apply constraints on formation from their bulk internal compositions. We report the discovery of Kepler-1704 b, a transiting 4.15 M {sub J} giant planet on a 988.88 day orbit with an extreme eccentricity of 0.921{sub -0.015}{sup +0.010}. Our decade-long RV baseline from the Keck I telescope allows us to measure the orbit and bulk heavy-element composition of Kepler-1704 b and place limits on the existence of undiscovered companions. A failed hot Jupiter, Kepler-1704 b was likely excited to high eccentricity by scattering events that possibly began during its gas accretion phase. Its final periastron distance was too large to allow for tidal circularization, so now it orbits its host from distances spanning 0.16–3.9 au. The maximum difference in planetary equilibrium temperature resulting from this elongated orbit is over 700 K. A simulation of the thermal phase curve of Kepler-1704 b during periastron passage demonstrates that it is a remarkable target for atmospheric characterization from the James Webb Space Telescope, which could potentially also measure the planet’s rotational period as the hot spot from periastron rotates in and out of view. Continued characterization of the Kepler-1704 system promises to refine theories explaining the formation of hot Jupiters and cool giant planets like those in the solar system.
We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.
Abstract Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars or black holes) and the synthesis and injection of many heavy elements into their host galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of one-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infrared light curves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine. Join us in improving and evolving this document through active community engagement. Instructions are provided at this link: https://github.com/clfryer/MM-SNe .