SEARCH · Engineering Papers
Results for “nuclear shell model”
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.
Precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections
The precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections in a quasistationary state is studied. The rates in question obtained analytically are compared with their numerical counterparts resulting from a dynamical simulation based on the Smoluchowski equation. It is shown that the integral Kramers formula is applicable even in the presence of shell corrections. As the excitation energy of the nucleus being considered decreases, the inaccuracy of this formula decreases monotonically, reaching 2% (inaccuracy of the dynamical simulation) in which case the dimensionless controlling parameter (entropy barrier) becomes greater than four. At small and intermediate values of this parameter, the accuracy of the integral Kramers formula is higher than the accuracy of the approximate Kramers formula obtained from it.
Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions
This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.
Uncertainties in ab initio nuclear structure calculations with chiral interactions
We present theoretical ground state energies and their uncertainties for p -shell nuclei obtained from chiral effective field theory internucleon interactions as a function of chiral order, fitted to two- and three-body data only. We apply a Similary Renormalization Group transformation to improve the numerical convergence of the many-body calculations, and discuss both the numerical uncertainties arising from basis truncations and those from omitted induced many-body forces, as well as chiral truncation uncertainties. With complete Next-to-Next-to-Leading (N 2 LO) order two- and three-body interactions, we find significant overbinding for the ground states in the upper p -shell, but using higher-order two-body potentials, in combination with N 2 LO three-body forces, our predictions agree with experiment throughout the p -shell to within our combined estimated uncertainties. The uncertainties due to chiral order truncation are noticeably larger than the numerical uncertainties, but they are expected to become comparable to the numerical uncertainties at complete N 3 LO.
Nuclear states projected from a pair condensate
Atomic nuclei exhibit deformation, pairing correlations, and rotational symmetries. To meet these competing demands in a computationally tractable formalism, we revisit the use of general pair condensates with good particle number as a trial wave function for even-even nuclei. After minimizing the energy of the condensate, allowing for general triaxial deformations, we project out states with good angular momentum with a fast projection technique. To show applicability, we present example calculations from pair condensates in several model spaces and compare against angular-momentum projected Hartree-Fock and full configuration-interaction shell-model calculations. This approach successfully generates spherical, vibrational, and rotational spectra, demonstrating potential for modeling medium- to heavy-mass nuclei.
Formation of Hyper Hydrogen isotopes in light p-shell nuclei
Recent experiment by A1 Collaboration at Mainz, namely identification of {sub Λ}{sup 4}H hyperfragment from primary {sub Λ}{sup 9}Li hypernucleus produced in reaction (e, e′K{sup +}), stirred renewed interest in baryonic decay of hypernuclei. The important role of s{sup −1}s{sub Λ} excitation in p-shell hypernuclei as well as the existence of selection rules connected with [f] (Young diagram) was noted earlier within Translational Invariant Shell Model (TISM). The particular conditions of the present experiment (reaction is not selective and produces highly excited states) dictate extension of the simplest TISM: expansion of the harmonic oscillator basis space and compression of multiplets (to only (λμ) for L and [f] for S, T). Such modified TISM explains abundance production of {sub Λ}{sup 4}H hyperfragment and predicts production of another Hyper Hydrogen {sub Λ}{sup 6}H.
Generation of Fragment Angular Momentum in Fission
A recent analysis of experimental data [J. Wilson et al., Nature (London) 590, 566 (2021)] found that the angular momenta of nuclear fission fragments are uncorrelated. Based on this finding, the authors concluded that the spins are therefore determined only after scission has occurred. We show here that the nucleon-exchange mechanism, as implemented in the well-established event-by-event fission model freya, while agitating collective rotational modes in which the two spins are highly correlated, nevertheless leads to fragment spins that are largely uncorrelated. We find that this counterexample invalidates the conclusion in [J. Wilson et al.] that uncorrelated spins must necessarily have been generated after scission (a potentious conclusion that would rule out all models that generate the fragment spins prior to scission). Furthermore, it was reported [J. Wilson et al.] that the mass dependence of the average fragment spin has a sawtooth structure. We demonstrate that such a behavior naturally emerges when shell and deformation effects are included in the moments of inertia of the fragments at scission.
Properties of infrared extrapolations in a harmonic oscillator basis
Here, the success and utility of effective field theory (EFT) in explaining the structure and reactions of few-nucleon systems has prompted the initiation of EFT-inspired extrapolations to larger model spaces in ab initio methods such as the no-core shell model (NCSM). In this contribution, we review and continue our studies of infrared (ir) and ultraviolet (uv) regulators of NCSM calculations in which the input is phenomenological NN and NNN interactions fitted to data. We extend our previous findings that an extrapolation in the ir cutoff with the uv cutoff above the intrinsic uv scale of the interaction is quite successful, not only for the eigenstates of the Hamiltonian but also for expectation values of operators, such as r 2 , considered long range. The latter results are obtained with Hamiltonians transformed by the similarity renormalization group (SRG) evolution. On the other hand, a possible extrapolation of ground state energies in the uv cutoff when the ir cutoff is below the intrinsic ir scale is not robust and does not agree with the ir extrapolation of the same data or with independent calculations using other methods.
Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems
Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.
Microscopic optical potentials from a Green's function approach
Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.
Fourteenth Exotic Beam Summer School EBSS 2015
The Fourteenth Annual Exotic Beam Summer School EBSS 2015 was held August 2nd - August 7th, 2015, and belongs to the series of summer programs aimed at educating future workforce in nuclear physics-related areas, mostly about the challenges of radioactive ion beam physics. Through these schools the research community will be able to exploit fully the opportunities created by the exotic beam facilities. These facilities in the US include CARIBU at ANL, the NSCL and the future FRIB laboratory as well as smaller-scale university laboratories. The skill set needed by the future workforce is very diverse and a fundamental understanding of theoretical, technical, computational and applied fields are all important. Therefore, the Exotic Beam Summer Schools follow a unique approach, in which the students not only receive lectures but also participate in hands-on activities. The lectures covered broad topics in both the experimental and theoretical physics of nuclei far from stability as well as radioactive ions production and applications. The afternoons provided opportunities for "hands-on" projects with experimental equipment and techniques useful in FRIB research. Five activities were performed in groups of eight students, rotating through the activities over the five afternoons of the school. The center of the activities was an experiment at the FSU tandem accelerator, measuring the angular distribution and cross section of the 12 C(d,p) 13 C transfer reaction, measured with a silicon telescope in a scattering chamber. The experimental data were analyzed by performing a DWBA calculation with the program DWUCK, and the resulting spectroscopic factors were compared to a shell model calculation. The other activities included target preparation, digital gamma-spectroscopy and modern neutron detection methods.
Contribution of higher multiplicity collisions to elastic p{sup 6,8}He and p{sup 8,9}Li scattering within Glauber theory
Differential cross sections for elastic p{sup 6,8}He and p{sup 8,9}Li scattering at energies between 60 and 70 MeV per nucleon and at the energy of 700 MeV per nucleon were calculated. The calculations in question were performed with the wave functions found on the basis of the α–n–n (for {sup 6}He), α–t–n (for {sup 8}Li), and α–t–2n (for {sup 9}Li) three-body models and with the density from the large-scale shell model for the {sup 8}He nucleus. The respective matrix elements were derived either upon taking fully into account the multiple-scattering operator or in the optical-limit approximation. A comparison of the results of the precise and approximate calculations made it possible to estimate reliably the contribution of higher multiplicity collisions to the differential cross sections.
Toward a Fundamental Understanding of Nuclear Reactions and Exotic Nuclei
Nuclear systems near the drip lines offer an exciting opportunity to advance our understanding of the interactions among nucleons, which has so far been mostly based on the study of stable nuclei. However, this is not a goal devoid of challenges. From a theoretical standpoint, it requires the ca- pability to address within an ab initio framework not only bound, but also resonant and scattering states, all of which can be strongly coupled. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from Quantum Chromodynamics employing Hamiltonians constructed within chiral effective field theory. In this contribution, we present a brief overview of one of such methods, the ab initio no-core shell model with continuum, and its applications to nucleon and deuterium scattering on light nuclei. The first investigation of the low-lying continuum spectrum of 6He within an ab initio framework that encompasses the 4He+n+n three-cluster dynamics characterizing its lowest particle-decay channel will also be briefly presented.
Lifetime Measurements and Triple Coexisting Band Structure in S 43
Lifetime measurements of excited states in the neutron-rich nucleus 43 S were performed by applying the recoil-distance method on fast rare-isotope beams in conjunction with the Gamma-Ray Energy Tracking In-beam Nuclear Array. The new data based on γγ coincidences and lifetime measurements resolve a doublet of (3/2 - ) and (5/2 - ) states at low excitation energies. Results were compared to the π(sd)-ν(pf) shell model and antisymmetrized molecular dynamics calculations. The consistency with the theoretical calculations identifies a possible appearance of three coexisting bands near the ground state of 43 S: the K π =1/2 - band built on a prolate-deformed ground state, a band built on an isomer with a 1f$-1\atop{7/2}$ character, and a suggested excited band built on a newly discovered doublet state. The latter further confirms the collapse of the N=28 shell closure in the neutron-rich region.
Prompt-delayed γ-ray spectroscopy of neutron-rich 119.121 In isotopes
The Z = 50 shell closure, near N = 82, is unique in the sense that it is the only shell closure with the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , enclosing the magic gap. The interaction of the proton hole/particle in the above-mentioned orbitals with neutrons in the νh 11/2 orbital is an important prerequisite to the understanding of the nuclear structure near N = 82 and the νπ interaction. To explore the structural similarity between the high-spin isomeric states in In ( Z = 49), Sn ( Z = 50), and Sb ( Z = 51) isotopes from a microscopic point of view. In addition, to understand the role of a proton hole or particle in the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , respectively, with neutron holes in the ν h 11 / 2 orbital on these aforementioned isomers. The fusion and transfer induced fission reaction Be 9 ( U 238 , f ) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++, and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich In 119 , 121 isotopes. This setup enabled the prompt-delayed γ -ray spectroscopy of isotopes in the time range of 100 ns – 200 μ s . In the odd- A In 119 , 121 isotopes, indications of a short half-life 19/2 - isomeric state, in addition to the previously known 25/2 + isomeric state, were observed from the present data. Further, new prompt transitions above the 25/2 + isomer in In 121 were identified along with reevaluation of its half-life. The experimental data were compared with the theoretical results obtained in the framework of large-scale shell-model calculations in a restricted model space. Furthermore, the ( π g 9 / 2 ν h 11 / 2 ; I | H ^ | π g 9 / 2 ν h 11 / 2 ; I ) two-body matrix elements of residual interaction were modified to explain the excitation energies and the B ( E 2) transition probabilities in the neutron-rich In isotopes. The (i) decreasing trend of E (29/2 + ) - E (25/2 + ) in odd-In (with dominant configuration π g 9 / 2 - 1 ν h 11 / 2 - 2 and maximum aligned spin of 29/2 + ) and (ii) increasing trend of E (27/2 + ) - E (23/2 + ) in odd-Sb (with dominant configuration π g 7 / 2 + 1 ν h 11 / 2 - 2 and maximum aligned spin of 27/2 + ) with increasing neutron number could be understood as a consequence of hole-hole and particle-hole interactions, respectively.
SU-E-T-667: Radiosensitization Due to Gold Nanoparticles: A Monte Carlo Cellular Dosimetry Investigation of An Expansive Parameter Space
Purpose: To investigate dose enhancement to cellular compartments following gold nanoparticle (GNP) uptake in tissue, varying cell and tissue morphology, intra and extracellular GNP distribution, and source energy using Monte Carlo (MC) simulations. Methods: Models of single and multiple cells are developed for normal and cancerous tissues; cells (outer radii 5–10 µm) are modeled as concentric spheres comprising the nucleus (radii 2.5–7.5 µm) and cytoplasm. GNP distributions modeled include homogeneous distributions throughout the cytoplasm, variable numbers of GNP-containing endosomes within the cytoplasm, or distributed in a spherical shell about the nucleus. Gold concentrations range from 1 to 30 mg/g. Dose to nucleus and to cytoplasm for simulations including GNPs are compared to simulations without GNPs to compute Nuclear and Cytoplasm Dose Enhancement Factors (NDEF, CDEF). Photon source energies are between 20 keV and 1.25 MeV. Results: DEFs are highly sensitive to GNP intracellular distribution; for a 2.5 µm radius nucleus irradiated by a 30 keV source, NDEF varies from 1.2 for a single endosome containing all GNPs to 8.2 for GNPs distributed about the nucleus (7 mg/g). DEFs vary with cell dimensions and source energy: NDEFs vary from 2.5 (90 keV) to 8.2 (30 keV) for a 2.5 µm radius nucleus and from 1.1 (90 keV) to 1.3 (30 keV) for a 7.5 µm radius nucleus, both with GNPs in a spherical shell about the nucleus (7 mg/g). NDEF and CDEF are generally different within a single cell. For multicell models, the presence of gold within intervening tissues between source and target perturbs the fluence reaching cellular targets, resulting in DEF inhomogeneities within a population of irradiated cells. Conclusion: DEFs vary by an order of magnitude for different cell models, GNP distributions, and source energies, demonstrating the importance of detailed modelling for advancing GNP development for radiotherapy. Funding provided by the Natural Sciences and Engineering Council of Canada (NSERC), and the Canada Research Chairs Program (CRC)
Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: evolution to the end of core helium burning
We present a comparison of 15M ⊙ , 20M ⊙ and 25M ⊙ stellar models from three different codes|GENEC, KEPLER and MESA|and their nucleosynthetic yields. The models are calculated from the main sequence up to the pre-supernova (pre-SN) stage and do not include rotation. The GENEC and KEPLER models hold physics assumptions that are characteristic of the two codes. The MESA code is generally more flexible; overshooting of the convective core during the hydrogen and helium burning phases in MESA is chosen such that the CO core masses are consistent with those in the GENEC models. Full nucleosynthesis calculations are performed for all models using the NuGrid post-processing tool MPPNP and the key energy-generating nuclear reaction rates are the same for all codes. We are thus able to highlight the key diferences between the models that are caused by the contrasting physics assumptions and numerical implementations of the three codes. A reasonable agreement is found between the surface abundances predicted by the models computed using the different codes, with GENEC exhibiting the strongest enrichment of H-burning products and KEPLER exhibiting the weakest. There are large variations in both the structure and composition of the models—the 15M ⊙ and 20M ⊙ in particular—at the pre-SN stage from code to code caused primarily by convective shell merging during the advanced stages. For example the C-shell abundances of O, Ne and Mg predicted by the three codes span one order of magnitude in the 15M ⊙ models. For the alpha elements between Si and Fe the differences are even larger. The s-process abundances in the C shell are modified by the merging of convective shells; the modification is strongest in the 15M ⊙ model in which the C-shell material is exposed to O-burning temperatures and the γ -process is activated. The variation in the s-process abundances across the codes is smallest in the 25M ⊙ models, where it is comparable to the impact of nuclear reaction rate uncertainties. In general the differences in the results from the three codes are due to their contrasting physics assumptions (e.g. prescriptions for mass loss and convection). The broadly similar evolution of the 25M ⊙ models gives us reassurance that different stellar evolution codes do produce similar results. For the 15M ⊙ and 20M ⊙ models, however, the different input physics and the interplay between the various convective zones lead to important differences in both the pre-supernova structure and nucleosynthesis predicted by the three codes. For the KEPLER models the core masses are different and therefore an exact match could not be expected.
γ-ray measurements in fast-neutron-induced reactions on 203 Tl
Fast-neutron-induced reactions can be used to characterize reaction mechanisms, investigate nuclear structure, and impose constraints on nuclear models. Because of the difficulty in predicting such effects, experimental data are important to constrain models. Thallium isotopes located close to the doubly magic Pb 208 nucleus are important for comparison with shell-model calculations. Study the population of excited states in Tl 203 and lighter isotopes in such reactions. Our methods are as follows. γ -ray cross sections were measured. The data were taken by using the Germanium Array for Neutron-Induced Excitations spectrometer. The pulsed neutron source of the Los Alamos Neutron Science Center's Weapons Neutron Research facility provided neutrons in the energy range from 1 to 300 MeV. The time-of-flight technique was used to determine the incident neutron energies. Cross sections for emission of several γ rays in 20 reaction channels were determined. Candidates for the intruder 9 / 2 - state in Tl 203 , from the odd proton in the h 9 / 2 orbital, and the first 5 + states in Tl 202 , 204 , from admixture of configurations, were identified. The excitation energy of the candidate 9 / 2 - , π h 9 / 2 state is in good agreement with the theoretical prediction from a semi-empirical weak-coupling model and its half-life is in the range of tens of nanoseconds. The feeding of a previously known isomer in Tl 202 exhibits similarities with the feeding of other isomers in this mass region.