Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Effective Theory”

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 487 records · Page 27

Measurement and microscopic description of odd–even staggering of charge radii of exotic copper isotopes

Nuclear charge radii globally scale with atomic mass number A as A 1/3 , and isotopes with an odd number of neutrons are usually slightly smaller in size than their even-neutron neighbours. This odd–even staggering, ubiquitous throughout the nuclear landscape1, varies with the number of protons and neutrons, and poses a substantial challenge for nuclear theory. Here, we report measurements of the charge radii of short-lived copper isotopes up to the very exotic 78 Cu (with proton number Z = 29 and neutron number N = 49), produced at only 20 ions s –1 , using the collinear resonance ionization spectroscopy method at the Isotope Mass Separator On-Line Device facility (ISOLDE) at CERN. We observe an unexpected reduction in the odd–even staggering for isotopes approaching the N = 50 shell gap. To describe the data, we applied models based on nuclear density functional theory and A-body valence-space in-medium similarity renormalization group theory. Through these comparisons, we demonstrate a relation between the global behaviour of charge radii and the saturation density of nuclear matter, and show that the local charge radii variations, which reflect the many-body polarization effects, naturally emerge from A-body calculations fitted to properties of A ≤ 4 nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleonic shells and nuclear masses

The binding energy of an isotope is a sensitive indicator of the underlying shell structure as it reflects the net energy content of a nucleus. Since magic nuclei are significantly lighter, or more bound, compared to their neighbors, the presence of nucleonic shell structure makes an imprint on nuclear masses. In this study, using a carefully designed binding-energy indicator, we catalog the appearance of spherical and deformed shell and subshell closures throughout the nuclear landscape. After presenting experimental evidence for shell and subshell closures as seen through the lens of nuclear masses, we study the ability of global nuclear mass models to predict local binding-energy variations related to shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Many-particle theory of nuclear systems with application to neutron star matter

The energy-density relation was calculated for pure neutron matter in the density range relevant for neutron stars, using four different hard-core potentials. Calculations are also presented of the properties of the superfluid state of the neutron component, along with the superconducting state of the proton component and the effects of polarization in neutron star matter.

Chakkalakal, D. A.↗

Chromium isotope fractionation during reduction of Chromium(VI) by Iron(II/III)-bearing clay minerals

Chromium stable isotope ratios are used to trace the reduction of Cr(VI) to Cr(III) in both ancient and modern systems. However, quantitative interpretation of Cr isotopic signatures has been stymied by the large variability in isotopic fractionation factors for Cr(VI) reduction by different reductants. Here we determine Cr isotope fractionation factors during Cr(VI) reduction by Fe(II/III)-bearing clay minerals, which are abundant in subsurface environments. Several variables were tested: pH, total Fe content of the clay, and the fraction of reduced Fe within the clay (Fe(II)/Fe(total)). The latter controls the standard reduction potential of the clay. Our results demonstrate that neither pH nor total Fe content of the clay have major effects on isotopic fractionation. In contrast, as the effective standard reduction potential of the clay and thus the standard free energy of Cr(VI) reduction become more negative, Cr isotope fractionation factors decrease in magnitude from - 4.9 to -1.3‰ according to a linear free energy relationship. This linear free energy relationship can be predicted from Marcus electron transfer theory and allows first-order predictions of Cr isotope fractionation factors to be made from the standard reduction potential or Fe(II)/Fe(total) of a clay, potentially improving our ability to model Cr isotope signatures in geochemical systems. Chromium is the first isotope system to show such a linear free energy relationship over a diverse range of reductants, including both aqueous and solid-phase reductants, and may provide a model for determining other redox-driven kinetic isotope effects in environmentally important isotope systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Application of the bRapid fission matrix burnup methodology to the JSI TRIGA Mark II research reactor

The novel burnup methodology, bRAPID, of the RAPID code system has been applied to the TRIGA Mark II research reactor at the 'Jozef Stefan' Institute. This paper presents the first step in our goal for on-line determination of isotope inventory and burnup, that the RAPID code system with its Fission Matrix (FM) method enables. Burnup analysis of changes in FM coefficients due to TRIGA fuel burnup at different irradiation times and specific power distribution in the reactor core is presented. Our results demonstrate that primary contribution for changes of FM coefficients is the burnup of the destination s, while the spectrum differences, that arise due to fuel's core location, have negligible effect on FM coefficients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

QCD resummation of dijet azimuthal decorrelations in pp and pA collisions

We study the azimuthal angular decorrelations of dijet production in both proton-proton (pp) and proton-nucleus (pA) collisions. By utilizing soft-collinear effective theory, we establish the factorization and resummation formalism at the next-to-leading logarithmic accuracy for the azimuthal angular decorrelations in the back-to-back limit in pp collisions. We propose an approach where the nuclear modifications to dijet production in pA collisions are accounted for in the nuclear modified transverse momentum dependent parton distribution functions (nTMDPDFs), which contain both collinear and transverse dynamics. This approach naturally generalizes the well-established formalism related to the nuclear modified collinear parton distribution functions (nPDFs). We demonstrate strong consistency between our methodology and the CMS measurements in both pp and pA collisions, and make predictions for dijet production in the forward rapidity region in pA collisions at LHC kinematics and for mid-rapidity kinematics at sPHENIX. Throughout this paper, we focus on the application of this formalism to a simultaneous fit to both collinear and transverse momentum dependent contributions to the transverse momentum dependent distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A new metric to control nucleation and grain size distribution in hybrid organic–inorganic perovskites by tuning the dielectric constant of the antisolvent

In perovskite research, there is a widely exploited but poorly explained phenomenon in which the addition of “antisolvents (ATS)” to precursor solutions results in higher-quality films. We explain the mechanism and driving force underlying an antisolvent-driven solvent extraction process. Density functional theory calculations uncover the defining effects of antisolvent choice on the extent of complexation between a lead salt and a methylammonium cation in solution. We experimentally validate the computational results using ultraviolet-visible spectroscopy and 207 Pb nuclear magnetic spectroscopy of methylammonium lead iodide solutions, containing both a processing solvent and an antisolvent. Furthermore, we uncover, and subsequently identify, the appearance of new species in solution as a result of the addition of the antisolvent. Furthermore, we observe that the choice of antisolvent has a substantial effect on the nature of the complexation of the methylammonium lead iodide (MAPbI 3 ) precursor species, whose origin we explain at an atomic level; specifically, the lower the dielectric of the antisolvent, the stronger the intermolecular binding energy between methylammonium cation (MA + ) cation and PbI 3 - plumbate, independent of the solvent or antisolvent interaction with the lead salt. Thin films were characterized using scanning electron microscopy; images of the films show how the addition of an antisolvent influences and, importantly, can be used to alter thin-film grain size. Grain size and distribution in thin films is reflected by the choice of antisolvent, promoting slower nucleation rates, a lower nucleation density, and hence larger final grain size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The smallest fluid on earth

High energy heavy ion collisions create quark gluon plasmas that behave like almost perfect fluids. Very similar features to those that led to this insight have also been observed in experimental data from collisions of small systems, involving protons or other light nuclei. Here we describe recent developments aimed at understanding whether, and if so how, systems that produce relatively few particles (orders of magnitude less than in typical heavy ion collisions) and are only one to a few times the size of a proton, can behave like fluids. This involves a deeper understanding of fluid dynamics and its applicability, improvements of our understanding of the initial geometry of the collisions by considering fluctuations of the proton shape, as well as advancements in the calculation of initial state effects within an effective theory of quantum chromodynamics, which can affect the observables that are used to study fluid behavior. We further address open questions and discuss future directions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrodynamic generators in relativistic kinetic theory

Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detailed examination of astrophysical constraints on the symmetry energy and the neutron skin of 208 Pb with minimal modeling assumptions

We report the symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, $R^{^{208}Pb}_{skin}$, implying a symmetry-energy slope parameter L of 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory (χEFT) and PREX-II constraints yields S 0 = 33.0$^{+2.0}_{-1.8}$ MeV, L = 53$^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin}$ = 0.17 $^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual χEFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between χ EFT, astrophysical data, and PREX-II's $R^{^{208}Pb}_{skin}$ measurement (p value = 12.3%) and that there is excellent agreement between χEFT, astrophysical data, and other nuclear experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗

Solving nuclear structure problems with the adaptive variational quantum algorithm

Here we use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the sd and pf shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interference effects for 0 ν β β decay in the left-right symmetric model

Various mechanisms may contribute to neutrinoless double beta decay in the left-right symmetric model. The interference between these mechanisms also contribute to the overall decay rate. The analysis of the contributions of these interference terms is important for disentangling different mechanisms. Here, we study interference effects contributing to the decay rate for neutrinoless double-β decay in the left-right symmetric model. The numerical values for maximum interference for several nuclides are calculated. It is observed that, for most of the interference terms, the contribution is smaller than 20% for all the nuclei considered in the study. However, the interference between the mass mechanisms (light and heavy) and η mechanism is observed to be in the range 30%–50%. The variation of the interference effect with the Q values is also studied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗