Search for hep solar neutrinos and the diffuse supernova neutrino background using all three phases of the Sudbury Neutrino Observatory
Not Available
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.
Not Available
This study is the first of a multi-phase effort to assess the sensitivity of light-water reactor (LWR) core-level prediction biases to changes in lattice-level calculation parameters. Prediction bias is the measured-to-predicted difference in a core-level quantity of interest (QOI) which can be estimated by comparing the simulation results with the plant-measured data for key nuclear parameters. The LWR two-step neutronics codes employed herein are the SCALE–Polaris lattice physics code (v6.3.1) and the Purdue Advanced Reactor Core Simulator (PARCS) nodal diffusion simulator (v3.4.2), both funded and used for confirmatory analysis to support licensing by the US Nuclear Regulatory Commission (NRC). Polaris–PARCS is used to model Watts Bar Unit 1 cycles 1–3 and Peach Bottom Unit 2 cycles 1–3. This study focuses on the impact of changes to Polaris input options such as scattering treatment or quadrature settings and how these input options induce changes in core-level quantities of interest (QOIs)bias. The report documents multiple bias assessments for different modeling choices and compares the bias magnitude to the QOI measurement uncertainties. Future companion reports will investigate the sensitivity of core-level LWR prediction bias to Polaris input options and Polaris-computed QOIs such as few-group assembly-homogenized cross sections to gain an understanding of the key drivers of prediction bias at lattice and core levels for application of a two-step LWR neutronics procedure in a licensing scenario.
The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.
Polymers are widely used industry, but degrade due to environmental conditions, such as elevated temperature or radiation. The degradation of polymers is dominated oxidation in oxygenated environments. Such oxidation is commonly referred to as thermal-oxidation or photo-oxidation and in many situations is heterogeneous through thickness of the polymer. While heterogeneous oxidation of polymers is a common phenomenon, the effect on the lifetime of polymers is not well understood as methods to evaluate such oxidation are time consuming and challenging to implement. Here we show a novel method for quantitative evaluation of heterogeneous oxidation in polymers using a fast and accurate color approach. The proposed color approach was observed to improve the efficiency of obtaining heterogenous oxidation profiles in half the time of common approaches, such as indention or Fourier transform infrared spectroscopy. Furthermore, the method was found to produce trends similar to that observed using common approaches and the applicability of the method was verified by mapping oxygen concentration through thickness using energy dispersive X-ray spectroscopy. Lastly, the utility of the proposed method was demonstrated through the generation of an oxidation design space, which is envisioned to assist with designing experiments towards estimation of polymer lifetime. Finally, we anticipate this work to be a starting point for further evaluation of heterogeneous oxidation in polymers using optical properties.
High thermal conductivity nuclear fuels offer important potential advantages over traditional oxide-based fuels such as higher burnup, reduction in fission gas release, and better overall safety of the system. One proposed approach to high thermal conductivity fuels utilizes high thermal conductivity nonfissile additives with UO 2 fuel to lower the fuel operating temperature and thereby take advantage of the highly favorable radiation resistance of UO 2 at lower operating temperatures. However, differential swelling in the matrix and high conductivity additive phases during high dose irradiation could lead to internal cracking and poor performance. In the current study, ceria (CeO 2 ) and zirconia (ZrO 2 ) surrogate matrices were used to model UO2 behavior. Additives of 10 vol. % Al 2 O 3 or SiC in the form of short fibers or platelets were used for the high conductivity second phase. The nuclear fuel surrogates were sintered to achieve densities greater than 93% of the ideal values. Scanning electron microscopy (SEM) imaging and X-ray diffraction confirmed the uniform distribution of the second phase and that no intermetallic second phase was formed during sintering. The thermal conductivity of the sintered samples was measured from 50 °C to 900 °C and confirmed the desirable increase compared to pure CeO 2 /ZrO 2 pellets. Samples were irradiated with 20 MeV Ni 6+ ions at midrange doses ranging from 1 to 15 displacements per atom (dpa) and temperatures from 300 °C to 700 °C. Post irradiation characterization revealed a good stability of the samples at low to medium doses with matrix lattice parameter swelling of < 0.14 % but showed a significant microstructural deterioration and decrease of the mechanical properties at 15 dpa.
The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.
C–H terminated nanometer scale diamonds ( d = 1 to 15 nm) are synthesized from 1-fluoroadamantane at high pressure (6–8 GPa) and high temperature (500–1500 °C) in a multianvil press. High resolution transmission electron microscopy, X-ray diffraction, Raman, diffuse reflectance Fourier transform infrared, and X-ray absorption spectroscopies demonstrate the excellent crystallinity and atomically flat C–H terminated surfaces of nanodiamonds with (111) and (110) facets. The importance of hydrogen to the synthesis of nanodiamond and its faceting is discussed. Following vacancy generation, annealing and oxidation of the nanodiamonds, optically detected magnetic resonance and electron spin resonance coherence times ( T 2 = 0.9 and 2.1 μs) of nitrogen vacancy (NV) centers are measured. Finally, the obtained T 2 values are equivalent to the shallow NV centers (depth <10 nm) in bulk diamond crystals and larger nanocrystals prepared by mechanical milling.
Thermal Energy Storage (TES) is a key component for solar thermal applications to bridge the gap between the demand for thermal energy and the supply of solar energy, whose availability depends on the time of day and season. Thus, cost-effective packed-bed thermal containers filled with a solid storage medium have been proposed for high-temperature sensible heat storage as materials are abundant and relatively cheap. Thus, it is necessary to investigate their performance and temperature profiles during the charge-discharge cycle. Several models are available for this purpose. Typically, the more detailed a model, the greater the computational effort required to solve it, and hence a time-efficient model is needed to prevent excessively long computation times for long-term analysis. At the more basic level, the common Hughes E-NTU model and the less realistic simplified Infinite-NTU model are very important for their less time and computational effort. In this paper, the appropriateness of employing the Infinite-NTU model was evaluated to investigate the performance of a typical and scalable rock-bed TES as a case study. The results presented provide a methodology to quickly test the validity of the model and predict the temperature profile for the case under study. Accordingly, such simple charge-discharge cycle thermal performance predictions are important to plan, design, and rapidly deploy a reliable and economical solar thermal system for the supply of valuable heat to high-temperature demanding applications of power generation and industrial processes as part of a rapid shift towards non-polluting renewable energy. Keywords: Solar Thermal, TES, Packed-bed, NTU model, Temperature profile
Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.
We present the results of Monte Carlo simulations of the microchannel plate (MCP) response to x rays in the 250 eV to 25 keV energy range as a function of both x-ray energy and impact angle and their comparisons with the experimental results from the X8A beamline at the National Synchrotron Light Source at Brookhaven National Laboratory. Incoming x rays interact with the lead glass of the microchannel plate, producing photoelectrons. Transport of the photoelectrons is neglected in this model, and it is assumed that photoelectrons deposit all their energy at the point they are created. This deposition leads to the generation of many secondary electrons, some fraction of which diffuse to the MCP pore surface where they can initiate secondary electron cascades in the pore under an external voltage bias. X-ray penetration through multiple MCP pore walls is increasingly important above 5 keV, and the effect of this penetration on MCP performance is studied. In agreement with past measurements, we find that the dependence of MCP sensitivity with angle relative to the pore bias changes from a cotangent dependence to angular independence and then proceeds to a secant dependence as the x-ray energy increases. We also find that with the increasing x-ray energy, the MCP gain sensitivity as a function of bias voltage decreases. The simulations also demonstrate that for x rays incident normal to the MCP surface, spatial resolution shows little dependence on the x-ray energy but degrades with the increasing x-ray energy as the angle of incidence relative to the surface normal increases. This agrees with experimental measurements. Simulation studies have also been completed for MCPs gated with a subnanosecond voltage pulse. In this work, we find that the optical gate profile width increases as the x-ray energy is increased above 5 keV, a consequence of increased x-ray penetration at energies >5 keV. Simulations of the pulsed dynamic range show that the dynamic range varies between ~100 and 1000 depending on x-ray energy and peak voltage.
The diffusion of defects (e.g., vacancies and interstitials) and elements used for dating (e.g., He and Pb) in a mineral structure is a thermal process: This is the primary assumption used to determine the age and thermal history of minerals. For instance, thermal history reconstruction, through the number and length distribution of tracks produced by spontaneous fission of 238 U, is obtained by assuming a thermal event to be the only energy source for shortening of fission tracks. Here, we report a new, non-thermal energy source that induces additional shortening of fission tracks by the irradiation of alpha-recoils from the alpha-decay of 238 U and 232 Th. We simulate alpha-decay induced track-shortening by combining ion accelerator irradiations with transmission electron microscopy. This allows for the first observation of track-shrinkage during in situ ion irradiation. We show that rather than alpha-particles, alpha-recoils induce a significant shortening of fission tracks by nuclear-collisions. The shortening of track-length can be quantified as a function of alpha-decay event dose. However, apatite is less sensitive than zircon to this non-thermal process. Furthermore, the findings exemplify the interactions among different types of self-irradiation from alpha-particles, alpha-recoils and fission-fragment nuclei in single mineral grains and have important implications for the use of zircon and apatite for radiometric dating and thermochronology.
This progress report (Level 3 Milestone Number M3SF-24LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Work Package Number SF-24LL01030205. The research is focused on actinide and radionuclide sequestration in hydrothermally altered repository materials. In FY24, we completed a rigorous analysis of Se sorption to iron oxide phases using our L-SCIE sorption database. This effort explicitly accounts for surface titration behavior of oxide surfaces that was the subject of a recently published manuscript (Han et al., 2023). With this in mind, our Se sorption analysis now yields a more robust workflow for developing self-consistent surface complexation modeling approaches that can be adapted to specific SCM conceptual and numerical approaches (i.e. non-electrostatic, diffuse layer, triple layer models). In late FY24, we plan to publish the results of our comprehensive surface complexation modeling of Se(IV) and Se(VI) sorption to iron oxide mineral phases and provide a path forward to developing robust radionuclide sorption models for use in performance assessment. In FY24, we also submitted a manuscript summarizing our approach to integrating radionuclide sorption and coprecipitation phenomena and evaluation of radionuclide partitioning values across a range of radionuclides relevant to performance assessment. We demonstrated our approach in detail using Se sorption and coprecipitation with iron oxide minerals as a test case. This manuscript was recently accepted for publication in Applied Geochemistry (Balboni et al., Accepted). We also continued experiments to identify radionuclide interaction with hydrothermally altered crystalline repository and backfill materials. Recent research performed at Los Alamos National Laboratory (LANL) and Sandia National Laboratory (SNL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are now examining how mineral alteration affects retardation behavior of a suite of radionuclides of interest to repository performance assessment. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database, developed using our L-SCIE approach, can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. A short update of results to date is presented below.
Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.
Canopy spectral reflectance collected measured in 2015-2018 along the NGEE-Arctic tram located on the Barrow Environmental Observatory (BEO), Barrow, Alaska. Measurements were made using a PP-Systems UniSpec-DC spectrometer mounted on an automatically operated instrument cart. Spectra were collected at 137 stop points at 0.5 m intervals along the tram transect. Transect runs were collected daily, or multiple times per day, from April or May through to August or September each year. Downwelling radiance was measured using a 2 meter fiber optic cable connected to a cosine diffuser. Upwelling (i.e. reflected) radiance was measured using a 2 meter cable connected to a 12 degree field-of-view (FOV) lens. Canopy reflectance was calculated using the ratio of upwelling to downwelling radiance measured over a 99.99% reflectance Spectralon standard measured at the start of each measurement set. Unprocessed spectral data are presented in .spu files. Metadata are provided in csv files, and the instrument manual (pdf). The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter includes the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. Two models are traditionally used to predict the defect concentration evolution, the standard rate theory (SRT) and the atomistic kinetic Monte Carlo (AKMC). We show in this work a large discrepancy in the defect concentrations between both methods when the average number of defects in the AKMC simulation box is close or lower than one. The reason is that AKMC naturally captures strong space and time correlations between vacancies and interstitials generated by the finite size of the periodic simulation box. These correlations strongly affect the recombination rate and the point defect concentrations. SRT fails to predict such correlations, and the corresponding solution deviates from the more accurate solution given by the AKMC simulation under similar conditions. These finite size correlation effects are strong when the elimination of point defects occurs by recombination only, but can still be significant in the presence of sinks. In order to account for the spacio-temporal correlation in a continuum framework, we introduce a Correlated Pair Theory (CPT). This theory fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations. Inversely, if the goal is for the AKMC to reproduce defect concentrations in bulk using small simulation boxes, the naturally occurring correlations need to be corrected. We show here that the CPT can be used to modify the elimination rates in the AKMC simulations, so as to yield point defect concentrations in agreement with SRT.
Concentrated hydrogen-bonded electrolytes (CoHBEs) are structured, electrochemically stable, less-volatile alternatives to aqueous and dilute nonaqueous electrolytes, however, with high viscosities that limit molecular diffusion. This work provides an understanding of the proton conduction mechanism in CoHBEs based on mixtures of acids and azoles and establishes a link between the structurally dictated transport properties and the proton-coupled electron transfer (PCET) reaction rates that can be leveraged for enhancing electrochemical reactions. Diffusion and relaxation NMR studies suggest a breaking of the viscosity–conductivity tradeoff, where at high azole concentrations (>45 mol%), Grotthuss transport is more likely with lowered proton transfer energy barriers between the azole and the acid according to the machine learning (ML) accelerated ab initio path integral MD (AI-PIMD) simulations. Proton conduction pathways are found to be switchable between the hydrogen bonding networks of the acid and the azole, with imidazole chain forming structures better facilitating Grotthuss hopping. Supported by small-angle neutron scattering studies, the chains are found to have six member molecules on average with maximum of 3 to 4 imidazole/imidazoliums at 50 to 60 mol%. Despite their high viscosities, the measured PCET rates for quinones and phenazines measured in the protic CoHBEs present relatively high electron transfer rate constants (k 0 ~ 10 −4 cm/s), validated by rotating disc electrode and scanning electrochemical microscopy measurements. The results demonstrate that strategic tuning of hydrogen-bond donor–acceptor interactions enables the decoupling of proton transport and viscosity, thereby impacting PCET reactions.
Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.
In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.