Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiative decay kinetics”

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

Microphysically modified magnetosonic modes in collisionless, high-β plasmas

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-β collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma β is increased. For mode amplitudes |δB/B 0 |≳0.3, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies |δB/B 0 |≳2β -1 . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-β plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Analysis of the Impact of Particle Growth on Transport Losses in Polymer-Electrolyte Fuel Cells

Voltage cycling causes catalyst nanoparticles in polymer-electrolyte fuel cells to grow. The concomitant loss of interfacial area results in larger kinetic, transport, and possibly ohmic overpotentials. This paper uses recently published experimental data and mathematical models to investigate the evolutions of transport and ohmic resistances to platinum nanoparticles located on the surface and inside the micropores of carbon black supports. Resistance to oxygen transport rises as surface area declines primarily because the flux to each remaining larger platinum particle increases. The path lengths governing oxygen diffusion to surface and buried platinum sites also increase as nanoparticles grow. Platinum nanoparticles on the surface become relatively less favorable as voltage cycling proceeds because they grow faster than platinum in micropores. Because voltage cycling causes total interfacial area to decline, and the fraction located inside micropores to increase, nanoscale ohmic losses increase as a catalyst layer decays. The practical importance of this effect is difficult to discern because proton conductivity in carbon micropores is not well characterized.

07 ISOTOPE AND RADIATION SOURCES↗

Persistent radicals in irradiated imidazolium ionic liquids probed by EPR spectroscopy

Long-lived radicals were observed in irradiated room temperature ionic liquids (RTILs) composed of the bis(trifluoromethylsulfonyl)imide anion (Tf 2 N - ) with 1-hexyl-3-methylimidazolium (hmim + ) and1-butyl-3-methylimidazolium (bmim+) cations. The EPR signal intensity increases within hours and the spectral pattern changes during the time after irradiation. The kinetics data obtained indicate the existence of at least two distinct radical cations that have different formation and decay rates. We demonstrate that oxygen does not react rapidly with the radical species formed. The bmim + Tf 2 N - was selectively deuterated at various positions which allowed us to suggest possible structures of the stable radicals formed. The first structure is a radical cation containing 2 nitrogens and 6 protons all with ca. 8.2 Gauss hyperfine splittings, and 2 protons with 2.8 Gauss splitting. The structure of the second radical cation is similar, but characterized by an odd number of hydrogens with ca. 8.2 Gauss splitting. Extensive quantum chemistry calculations were performed to attempt to identify the structure of these persistent radicals. The most likely candidates are imidazole radical cations having methyl groups at the nitrogen atoms and substituents in the second position of the imidazole ring.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Calculated fission-fragment mass yields and average total kinetic energies of heavy and superheavy nuclei

Fission-fragment mass and total-kinetic-energy (TKE) distributions following fission of even-even nuclides in the region \(74 \le Z \le 126\) and \(92 \le N \le 230\) , comprising 896 nuclides have been calculated using the Brownian shape-motion method. The emphasis is the region of superheavy nuclei. To show compatibility with earlier results the calculations are extended to include earlier studied regions. An island of asymmetric fission is obtained in the superheavy region, \(106\le Z\le 114\) and \(162\le N\le 176\) , where the heavy fragment is found to be close to \(^{208}\) Pb and the light fragment adjusts accordingly. Most experimentally observed \(\alpha \) -decay chains of superheavy nuclei with \(Z > 113 \) terminate by spontaneous fission in our predicted region of asymmetric fission. In these cases, the pronounced large asymmetry is accompanied by a low TKE value compatible with measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

The Extension of the Hauser-Feshbach Fission Fragment Decay Model to Multi-chance Fission and its Application to 239 Pu

The Hauser-Feshbach fission fragment decay model, HF3D, calculates the statistical decay of fission fragments through both prompt and delayed neutron and γ-ray emissions in a deterministic manner. While previously limited to the calculation of only first-chance fission, the model has recently been extended to include multi-chance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities–fission probabilities, pre-fission neutron energies, and the average energy causing fission– and post-scission quantities–yields in mass, charge, total kinetic energy, spin, and parity. From those fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions. The evaporation of the prompt neutrons and γ rays is calculated through the Hauser-Feshbach statistical theory, taking into account the competition between neutron and γ-ray emission, conserving energy, spin, and parity. The delayed emission is taken into account using time-independent calculation using decay data. This whole formulation allows for the calculation of prompt neutron and γ-ray properties, such as multiplicities and energy distributions, both independent and cumulative fission yields, and delayed neutron observables, in a consistent framework. Here, we describe the implementation of multi-chance fission into the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 239 Pu. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Subcritical Californium Source Drive Noise Analysis Measurements With Unreflected Uranium (93.15) Hydride

On July 6, 1989, subcritical californium source-driven noise analysis (CSDNA) measurements were performed with bare uranium hydride cylindrical assemblies at the Los Alamos National Laboratory Critical Experiments Facility in KIVA 1. Three configurations of ~7.5 cm radius-enriched (93.15 wt. % 235 U) uranium hydride with a density of ~10 g/cm 3 cylinders were assembled with uranium hydride heights of ~11, ~14, and ~16 cm. The uranium hydride (~2 and ~3 cm high) was in thin-welded stainless steel cans. The neutron multiplication factors obtained on-line from the measured ratios of spectral densities were 0.942 ± 0.002, 0.917 ± 0.003, and 0.867 ± 0.004 for uranium hydride heights of ~16, ~14, and ~11 cm., respectively. Neutron multiplication factors for these three configurations calculated by Los Alamos National Laboratory Monte Carlo methods were 0.952 ± 0.005, 0.922 ± 0.005, and 0.860 ± 0.005, all of which are in good agreement with the measurements. The break frequency noise analysis data could be fitted to obtain the prompt neutron decay constant at all subcritical states. Some data presented in this report are from notes that are not in the logbooks. A discussion of the two different point kinetics theories of the measurements that illustrated the deficiencies of the rigorous theory at low neutron multiplication factors (k eff < 0.80) is given in an Appendix. These measurements may be the only configuration of unreflected, highly enriched uranium hydride configuration assembled. The presence of hydrogen diluted the uranium density, decreasing the reactivity, but this is offset by slowing down the neutron energy to the energy ranges of higher-fission cross sections. Extrapolations to delayed critical indicate that adding one more 3 cm can and one more 2 cm can would result in the unreflected cylindrical system being slightly above critical with a critical mass of ~37 kg, which is considerably less than the critical mass of an unreflected and unmoderated highly enriched uranium metal sphere. An unreflected, highly enriched uranium hydride sphere would have a critical mass less than 37 kg.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization

While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Actinide Complexation on Ligand Radiolysis

Although the actinides boast many unique physical and chemical properties, their inherent susceptibility to radioactive decay are what make them truly interesting elements to study. The absorption of ionizing radiation from actinide decay leads to the formation of a variety of transient and steady-state radicals, ions, and molecular radiolysis products that can lead to significant changes in the surrounding environment, and ultimately dictate steady-state actinide redox distributions and the longevity of molecules designed for actinide complexation. Radiolysis of the latter leads to complexant destruction and the concomitant formation of degradation products that can complicate actinide studies and processes. However, the radiation chemistry of most actinide complexants have been studied in the absence of the actinides they were designed to complex, which can lead to inaccurate conclusions on longevity and degradation product distributions, as metal ion complexation has been historically shown to influence a ligand’s radiolytic behavior. Consequently, bridging this knowledge gap is important for actinide science. Here, I will discuss the impact of actinide complexation on the steady-state and time-resolved radiation-induced reactivity of a variety of complexants, including, tributyl phosphate (TBP) and N,N,N',N'-tetraoctyl diglycolamide (TODGA) .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Auger electron angular distributions following excitation or ionization from the Xe 3d and F 1s levels in xenon difluoride

Linearly polarized synchrotron radiation has been used to record polarization dependent, non-resonant Auger electron spectra of XeF 2 , encompassing the bands due to the xenon M 45 N 1 N 45 , M 45 N 23 N 45 , M 45 N 45 N 45 and M 45 N 45 V and fluorine KVV transitions. Resonantly excited Auger spectra have been measured at photon energies coinciding with the Xe 3d 5/2 → σ* and the overlapped Xe 3d 3/2 /F 1s → σ* excitations in XeF 2 . The non-resonant and resonantly excited spectra have enabled the Auger electron angular distributions, as characterized by the β A parameter, to be determined for the M 45 N 45 N 45 transitions. In the photon energy range over which the Auger electron angular distributions were measured, theoretical results indicate that transitions into the εf continuum channel dominate the Xe 3d photoionization in XeF 2 . In this limit, the theoretical value of the atomic alignment parameter (A 20 ) characterizing the core ionized state becomes constant. This theoretical value has been used to obtain the Auger electron intrinsic anisotropy parameters (α 2 ) from the β A parameters extracted from our non-resonant Auger spectra. For a particular Auger transition, the electron kinetic energy measured in the resonantly excited spectrum is higher than that in the directly ionized spectrum, due to the screening provided by the electron promoted into the σ* orbital. The interpretation of the F KVV Auger band in XeF 2 has been discussed in relation to previously published one-site populations of the doubly charged ions (XeF 2 2+ ). Here, the experimental results show that the ionization energies of the doubly charged states predominantly populated in the decay of a vacancy in the F 1s orbital in XeF 2 tend to be higher than those populated in the decay of a vacancy in the Xe 4d level in XeF 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy migration and scintillation kinetics in compositionally complex (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce single crystal scintillator

It is well-established that compositional tuning through binary admixture can improve scintillation performance in several materials systems, including Ce-activated garnets. Although recent work on ternary or quaternary cation admixture shows promise, the impact of this increased compositional complexity on thermal stability and carrier-defect dynamics has not been addressed. Here, we investigate a compositionally complex garnet, (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), grown by the Czochralski method using temperature-dependent photoluminescence (PL), PL decay, and thermoluminescence (TL). PL and PL decay measurements support a thermally activated Tb 3+ -Ce 3+ energy transfer, where Tb 3+ emission dominates below 60 K, but Ce 3+ emission increases from 20-300 K. Thermal quenching of Ce 3+ emission occurs around T 50 = 508 K, with an activation energy of 0.6 eV. TL and wavelength-resolved TL spectra from 20–500 K show that GYTLAG contains similar trap groups to LuAG but with a broader distribution of glow peaks below room temperature, possibly caused by quaternary cation mixing. A combination of dose dependence, partial cleaning and initial rise, and glow curve fitting to a first order continuous Gaussian distribution model are used to understand the contribution of electronic point defects to scintillation decay and afterglow at room temperature. Furthermore, these results inform how increased compositional complexity influences recombination dynamics in garnet scintillators.

Compositionally complex↗

Recombination of 𝐵 𝑐 mesons in ultrarelativistic heavy-ion collisions

High-energy heavy-ion collisions have been suggested as a favorable environment for the production of 𝐵 𝑐 mesons, due to a much larger abundance of charm and bottom quarks compared to elementary reactions. Motivated by recent CMS data for 𝐵$^{+}_{𝑐}$ production in Pb-Pb(5.02 TeV) collisions at the Large Hadron Collider (LHC), we deploy a previously developed transport approach for charmonia and bottomonia to evaluate the kinetics of 𝐵 𝑐 mesons throughout the fireball formed in these reactions. The main inputs to our approach are two transport parameters: the 𝐵 𝑐 's reaction rate and equilibrium limit. Both quantities are determined by previous calculations via a combination of charm and bottom sectors. In-medium binding energies of 𝐵 𝑐 mesons are calculated from a thermodynamic 𝑇 matrix with a lattice-QCD constrained potential, and figure in their inelastic reaction rates. Temperature-dependent equilibrium limits include charm- and bottom-quark fugacities based on their initial production. We compute the centrality dependence of inclusive 𝐵 𝑐 production and transverse-momentum (𝑝 𝑇 ) spectra using two different recombination models: instantaneous coalescence and resonance recombination. The main uncertainty in the resulting nuclear modification factors, 𝑅 AA , is currently associated with the 𝐵 𝑐 cross section in elementary 𝑝⁢𝑝 collisions, caused by the uncertainty in the branching ratio for the 𝐵$^{−}_{𝑐}$ → 𝐽/𝜓⁢𝜇 − ⁢$\overline{𝜈}$ decay. Our results indicate a large enhancement of the 𝑅 AA at low 𝑝 𝑇 , with significant regeneration contributions up to 𝑝 𝑇 ≃ 20GeV. Comparisons to CMS data are carried out but firm conclusions will require a more accurate value of the branching ratio, or alternative channels to measure the 𝐵 𝑐 production in 𝑝⁢ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Delayed Critical and Subcritical Experiments with Polyethylene Moderated Unreflected Thin 15 in. Diameter HEU Metal Plates

The thin ~15 in. diameter highly enriched uranium (HEU) metal plates were assembled to delayed criticality at the Oak Ridge Critical Experiments Facility (ORCEF) in 1969 with various thicknesses of polyethylene (varying from 1/16 to 2$\frac{3}{8}$ inches) between uranium metal plates. The average 235 U enrichment was 93.27 wt. %. These unreflected critical configurations contained 4$\frac{2}{3}$ to 20$\frac{5}{6}$ thin 15 in. diameter HEU metal plates (on loan from Los Alamos National Laboratory [LANL] and shipped to Oak Ridge National Laboratory [ORNL] on June 3, 1969). Depending on the thickness of polyethylene, the enriched uranium masses varying from 28,053 to 135,148 grams. Fractional plate sections consisted of the appropriate number of 60° pie sections. In addition to the measurement at delayed criticality, subcritical measurements were also performed by the inverse kinetic rod drop method. Prompt neutron decay constant measurements were also performed by the Rossi alpha and randomly pulsed neutron method using a time-tagged spontaneous fission californium neutron source; these are briefly reported here. At the time of these measurements in 1969, the thin HEU metal plates were in near-pristine condition with extremely little oxidation, allowing better descriptions of the uranium plates than the use of these plates in a heavily oxidized and deteriorated condition in recent reflected benchmark experiments at the LANL facility at the Nevada Test Site with these same thin highly enriched uranium metal plates. This report documents the experimental information for the measurements performed so that later researchers can perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Evaluation Program (ICSBEP) or a Nuclear Energy Agency (NEA) benchmark. Data from the experiments described should be acceptable for use as criticality safety benchmark experiments for the ICSBEP and the NEA nuclear criticality safety benchmark program once the uncertainty analysis on the measured neutron multiplication factors is completed. Additional data—such as the dimensional inspection reports, uranium isotopic information, and other relevant particulars—should be retrieved from the Y-12 Plant or LANL and incorporated in the final ICSBEP benchmark. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, the uncertainties in $k_{eff}$ could be as low as ± 0.0002 for some configurations. Other experiments with smaller-diameter than 15 in. diameter HEU metal plates have been benchmarked in HEU-METFAST-001. The prompt neutron time decay measurements could be the basis for an International Reactor Physics Benchmark Program. Preparation of the present report is part of an effort at ORNL to document more than 15 undocumented critical and subcritical experiments enumerated in ORNL/TM-2019/18 and performed by ORNL at ORCEF and other US Department of Energy critical experiments facilities using more than 500 operational days of critical facility time. This work for this report publication was supported by the Nuclear Criticality, Radiation Transport and Safety NCSP Program at ORNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alpha Decay Chains as Thermal Power Sources: Analysis and Applications for RTGs

Radioactive sources can provide power in remote and environmentally harsh locations such as the arctic or space. The generators powered by such sources are rugged and can withstand extreme temperatures, lack of sunlight, and require no human intervention for multiple years. Radioisotopes are used in thermoelectric generators to provide power at remote sites and deep in space. Isotopes like Pu-238, Cm-244, and Am-241 are used in these generators by NASA for power in space probes and spacecrafts. These power sources deliver a steady supply of energy over extended periods of time. Alpha particles created during decay do not travel far in a material. Their kinetic energy is transferred to heat that we can then convert into energy. Unlike beta and gamma decay, the slower-moving alpha particles stop in the material, making their energy available for use. Energy from these natural decay processes provides a reliable source of power. Spontaneous fission is rare and unreliable, and unlike induced fission processes, alpha decay occurs naturally and does not require external management or ignition. The ideal properties of an isotope for use as a power source depend upon the intended use. For use in an Arctic research base over a period of several years, but less than a decade, an isotope that provides high power output over a shorter lifespan may be the most suitable option. Whereas, for deep space missions where a consistent power source for decades or perhaps more than 100 years is needed that would require a very different isotope. One with a much longer half-life that would provide consistent power throughout that time and survive in that state in for these extended periods of time. These examples represent two extreme sides in terms of time frames. By analyzing the power produced by different radioactive decay processes over time, we can evaluate the suitability of various isotope decay chains for specific uses. Some unstable isotopes undergo a series of radioactive decays, transforming into different isotopes at each step and resulting in a stable isotope. The lists of isotopes in these decay processes are known as decay chains. Some of these chains, illustrated in the figures below, are currently being investigated for use in radioisotope thermoelectric generators (RTGs) designed for a range of operational durations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗