Search NASA⌕ Search

Engineering topics

Stetcu, Ionel

Publications and source records attributed to Stetcu, Ionel.

At least 19 records

Consistent $\overline{ν}$ evaluation for minor Pu isotopes with CGMF

This report follows up on the evaluation work done in 2023 where we developed a procedure and performed the first consistent evaluation of the average prompt neutron multiplicity for minor Pu isotopes, including 238,240-242 Pu, using CGMF. Details on the necessary updates to the release version of CGMF, the experimental data and experimental uncertainty quantification that went into the evaluation, and this first evaluation effort are documented in and will not be repeated in this report. Instead, we discuss here the further investigations into the CGMF model space to improve the evaluation results from FY23. Additionally, we note that these evaluation results have been transformed into the ENDF format for mean values and covariances, and validation with critical assemblies has been performed; that work is documented elsewhere. In this short report, we discuss the updated evaluation efforts performed during FY24 in Sec. 2, show results from CGMF for other prompt fission observables in Sec. 3, and then briefly conclude in Sec. 4.

07 ISOTOPE AND RADIATION SOURCES↗

Neck Rupture and Scission Neutrons in Nuclear Fission

Just before a nucleus undergoes fission, a neck is formed between the emerging fission fragments. It is widely accepted that this neck undergoes a rather violent rupture, despite the absence of unambiguous experimental evidence. The main difficulty in addressing the neck rupture and saddle-to-scission stages of fission is that both are highly nonequilibrium processes. Here, we present the first fully microscopic characterization of the scission mechanism, along with the spectrum and the spatial distribution of scission neutrons (SNs), and some upper limit estimates for the emission of charged particles. The spectrum of SNs has a distinct angular distribution, with neutrons emitted in roughly equal numbers in the equatorial plane and along the fission axis. They carry an average energy around 3±0.5 MeV for the fission of 236 U, 240 Pu, and 252 Cf, and a maximum of 16–18 MeV. We estimate a conservative lower bound of 9%–14% of the total emitted neutrons are produced at scission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep quantum circuit simulations of low-energy nuclear states

Numerical simulation is an important method for verifying the quantum circuits used to simulate low-energy nuclear states. However, real-world applications of quantum computing for nuclear theory often generate deep quantum circuits that place demanding memory and processing requirements on conventional simulation methods. Here, we present advances in high-performance numerical simulations of deep quantum circuits to efficiently verify the accuracy of low-energy nuclear physics applications. Our approach employs novel methods for accelerating the numerical simulation including management of simulated mid-circuit measurements to verify projection based state preparation circuits. In this study, we test these methods across a variety of high-performance computing systems and our results show that circuits up to 21 qubits and more than 115,000,000 gates can be efficiently simulated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microscopic evidence for scission neutrons

We discuss the neck rupture stage of fission and the emission of particles afterwards, two extremely rapid and highly non-equilibrium processes. Currently, the neck rupture cannot be directly probed by experiment, highlighting the importance of reliable theoretical predictions for this stage of fission. Time-dependent density functional theory (TDDFT) is used to simulate the spontaneous fission of 252 Cf. In conjunction with statistical models, inputs from microscopic calculations can be used to make predictions for prompt neutron and gamma emission spectra, quantities which can be experimentally measured. Here we characterize the scission mechanism within TDDFT and estimate the number of scission neutrons and their kinetic energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Noniterative finite amplitude methods for giant resonances and the application to the neutron radiative capture cross sections

We calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances with noniterative finite amplitude methods and demonstrate how the fully microscopic density functional theory predicts the giant resonances without any phenomenological parameters. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the result of E1 and M1 transitions and find that the capture cross sections for deformed nuclei are enhanced due to the contribution from the low energy M1 scissors mode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ENDF/B-VIII.1 Evaluation Updates for 234,236 U [Slides]

This presentation goes over the general evaluation procedure. Additionally covered are the fission and capture updates, reaction rates, correlation studies and finally a summary that provides closing remarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microscopic Evidence for Scission Neutrons

We discuss the neck rupture stage of fission and the emission of particles afterwards, two extremely rapid and highly non-equilibrium processes. Currently, the neck rupture cannot be directly probed by experiment, highlighting the importance of reliable theoretical predictions for this stage of fission. Time-dependent density functional theory (TDDFT) is used to simulate the spontaneous fission of 252 Cf. In conjunction with statistical models, inputs from microscopic calculations can be used to make predictions for prompt neutron and gamma emission spectra, quantities which can be experimentally measured. Here we characterize the scission mechanism within TDDFT and estimate the number of scission neutrons and their kinetic energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatial orientation of the fission fragment intrinsic spins and their correlations

New experimental and theoretical results obtained in 2021 made it acutely clear that more than 80 years after the discovery of nuclear fission we do not understand the generation and dynamics of fission fragment (FF) intrinsic spins well, in particular their magnitudes, their spatial orientation, and their correlations. The magnitude and orientation of the primary FFs have a crucial role in defining the angular distribution and correlation between the emitted prompt neutrons, and subsequent emission of statistical (predominantly $E1$) and stretched $E2$ $γ$ rays, and their correlations with the final fission fragments. Here, we present detailed microscopic evaluations of the FF intrinsic spins, for both even- and odd-mass FFs, and of their spatial correlations. In conclusion, these point to a well-defined three-dimensional FF intrinsic spin dynamics, characteristics absent in semiphenomenological studies, due to the presence of the twisting spin modes, which artificially were suppressed in semiphenomenological studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Baseline CGMF Calculations for Pu Isotopes

Updates were made to CGMF to run neutron-induced fission calculations on Pu isotopes of A =238, 239, 240, 241, and 242 (the release version of CGMF can calculate 239,241 Pu(n,f) and 238,240,242 Pu(sf) only). Here, we show the baseline calculations from CGMF for these neutron-induced fission reactions, including updates made for total kinetic energy systematics and fission barriers fit to be more consistent with the neutron-induced fission cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Consistent Nuclear Data Evaluations for Criticality Safety

Evaluations of nuclear data are based on statistical analysis of available experimental data and their uncertainties plus model calculations and their uncertainties. As the models are currently rather limited, the evaluations are heavily biassed toward experimental data, with the caveat that a thorough analysis is also required to understand possible discrepancies between data sets. Hence, as new experimental data become available, they are incorporated into the evaluation procedure. Recent measurements of the 233 U capture to fission cross section ratio at the Los Alamos Neutron Science Center have prompted a re-evaluation of the capture cross section in the resonance and fast regions up to 250 keV. We will discuss the challenges of including the new fast neutron experimental data in an evaluation that is consistent with the resonance region. We will also discuss our consistent evaluation procedure based on the Hauser-Feshbach statistical model for nuclear reactions and its application to the evaluations of 239 Pu and 139 La neutron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Consideration of memory of spin and parity in the fissioning compound nucleus by applying the Hauser-Feshbach fission fragment decay model to photonuclear reactions

Prompt and β-delayed fission observables, such as the average number of prompt and delayed neutrons, the independent and cumulative fission product yields, and the prompt γ-ray energy spectra for the photonuclear reactions on 235,238 U and 239 Pu, are calculated with the Hauser-Feshbach fission fragment decay (HF 3 ⁢D) model and compared with available experimental data. Further, in the analysis of neutron-induced fission reactions to the case of photo-induced fission, an excellent reproduction of the delayed neutron yields supports a traditional assumption that the photo fission might be similar to the neutron-induced fission at the same excitation energies regardless of the spin and parity of the fissioning systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collective enhancement in the exciton model

The preequilibrium reaction mechanism is considered in the context of the exciton model. A modification to the one-particle–one-hole state density is studied which can be interpreted as a collective enhancement. The magnitude of the collective enhancement is estimated by simulating the Lawrence Livermore National Laboratory (LLNL) pulsed-spheres neutron-leakage spectra. Here the impact of the collective enhancement is explored in the context of the highly deformed actinide, 239 Pu. A consequence of this enhancement is the removal of fictitious levels in the distorted-wave Born approximation often used in modern nuclear reaction codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Correlations between energy and $\gamma$-ray emission in 239 Pu($n, f$)

Here, we study $\gamma$-ray emission following $^{239}\mathrm{Pu}(n,\mathrm{f})$ over an incident neutron energy range of $2 < E_i < 40$ MeV. We present the first experimental evidence for positive correlations between the total angular momentum generated in fission and the excitation energy of the compound nucleus prior to fission. The $\gamma$-ray multiplicity increases linearly with incident energy below the 2\textsuperscript{nd}-chance fission threshold with a slope of $0.085 \pm 0.010$~MeV$^{-1}$. This linear trend appears to hold for the average excitation energy of the compound nucleus between $9 < \langle E_x \rangle < 19$ MeV. Most of the multiplicity increase comes from an enhancement around a $\gamma$-ray energy of 0.7 MeV, which we interpret as stretched quadrupole $\gamma$ rays that indicate an increase in total fission-fragment angular momentum with excitation energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alternative approach to quantum imaginary time evolution

There is increasing interest in quantum algorithms (QAs) that are based on the imaginary time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy postprocessing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept on the order of n , the number of qubits, by controlling the precision (number of significant digits) and the imaginary time increment. The number of quantum gates per imaginary time step is estimated to be polynomial in n . Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. Here, we illustrate our algorithm through numerical implementation on an IBM quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Open-source release of CGMF 1.1 and Integration into the MCNP6.3 ® Code [Slides]

As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗