Search NASA⌕ Search

SEARCH · Search NASA

Results for “Renormalization Group”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ab initio multishell valence-space Hamiltonians and the island of inversion

In the shell-model framework, valence-space Hamiltonians connecting multiple major-oscillator shells are of key interest for investigating the physics of neutron-rich nuclei, which have been the subject of intense experimental activity for decades. Here we present an extension of the ab initio valence-space in-medium similarity renormalization group, which allows the derivation of such Hamiltonians nonperturbatively. Starting from initial two- and three-nucleon forces from chiral effective field theory, we then calculate properties of nuclei in the important island-of-inversion region above oxygen, so far unexplored with ab initio methods. In this work, our results in the neon and magnesium isotopes indicate the importance of neutron excitation from the $\textit{sd}$ to $\textit{pf}$ shells and ground states dominated by intruder configurations around $\textit{N}$ = 20, consistent with the conclusions from phenomenological studies. We also benchmark the excitation spectrum of 16 O with coupled-cluster theory, finding generally good agreement, and discuss implications for ground-state energies and charge radii in oxygen and calcium isotopes. Finally we outline the proper procedure for treating the longstanding issue of center-of-mass contamination, and show that with a particular choice of valence space, these spurious states can be removed successfully.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Testing isospin symmetry breaking in ab initio nuclear theory

In this work we present the first steps towards benchmarking isospin symmetry breaking in ab initio nuclear theory for calculations of superallowed Fermi β decay. Using the valence-space in-medium similarity renormalization group, we calculate b and c coefficients of the isobaric multiplet mass equation, starting from two different Hamiltonians constructed from chiral effective field theory. We compare results to experimental measurements for all T = 1 isobaric analog triplets of relevance to superallowed β decay for masses A = 10 to A = 74 and find an overall agreement within approximately 250 keV of experimental data for both b and c coefficients. A greater level of accuracy, however, is obtained by a phenomenological Skyrme interaction or a classical charged-sphere estimate. Lastly, we show that evolution of the valence-space operator does not meaningfully improve the quality of the coefficients with respect to experimental data, which indicates that higher-order many-body effects are likely not responsible for the observed discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Two-pion exchange as a leading-order contribution in chiral effective field theory

Pion exchange is the central ingredient to nucleon-nucleon interactions used in nuclear structure calculations, and one-pion exchange (OPE) enters at leading order in chiral effective field theory. In the 2S+1 L J = 1 S 0 partial wave, however, OPE and a contact term needed for proper renormalization fail to produce the qualitative, and quantitative, features of the scattering phase shifts. Cutoff variation also revealed a surprisingly low breakdown momentum Λ b ≈ 330 MeV in this partial wave. Here we show that potentials consisting of OPE, two-pion exchange (TPE), and a single contact address these problems and yield accurate and renormalization group (RG) invariant phase shifts in the 1 S 0 partial wave. Here, we demonstrate that a leading-order potential with TPE can be systematically improved by adding a contact quadratic in momenta. For momentum cutoffs Λ ≲ 500 MeV, the removal of relevant physics from TPE loops needs to be compensated by additional contacts to keep RG invariance. Inclusion of the Δ isobar degree of freedom in the potential does not change the strong contributions of TPE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Entanglement and seniority

In quantum mechanics entanglement is the most striking phenomenon which has no counterpart in classical systems. Although, different approaches have already been developed to study correlations in the case of indistinguishable particles, the exploration of the so-called mode entanglement is still in its initial stage in nuclear physics. Study of mode-entanglement in the seniority model, derivation of analytic formulas for the one-body reduced density matrix of states with seniority $ν$ = 0, 1, 2, and $ν$ = 3, and also determination of the particle number dependence of the one-body reduced density matrix for arbitrary seniority. In addition, comparison of the predictions of the seniority model with the results of the shell model to gain insight into the structure and correlations of the ground and lowest yrast states. In the seniority model the analytic results are calculated using the quasi-spin formalism. The numerical shell model calculations are carried out in the standard shell model framework using an inert core. The applied realistic effective interactions are derived earlier using the $G$-matrix formalism. The density matrix renormalization group method is also applied in order to directly calculate the mode entropies. In the seniority model simple analytical expressions are given for the mode entropies. The peculiar behavior of the half-filled shells and the seniority-zero states are revealed. Numerical results are presented for the lightest stable calcium isotopes and for 94 Ru nucleus. For 94 Ru, the seniority model accounts for the 0$g$ 9/2 mode entropies, but seniority mixing is important for certain yrast states. Furthermore, interaction induced quantum fluctuations decrease the occupation of the 0$f$ 5/2 , 1$p$ 3/2 and 1$p$ 1/2 shells, and amount in non-negligible mode entropies on these shells, too, clearly outside the scope of the simple (0$g$ 9/2 ) 4 seniority model. The 0$f$ 7/2 shell-based seniority model is more accurate for the Ca isotopes, but seniority mixing is substantial for some 44 Ca yrast states, too. Mode and one-body entanglement entropies are useful tools to investigate the structure of quantum correlations in nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Ab initio calculation of muon capture on 24 Mg

In this work we study ordinary muon capture (OMC) on 24 Mg from a first-principles perspective. Starting from a particular two- and three-nucleon interaction derived from chiral effective field theory, we use the valence-space in-medium similarity renormalization group (VS-IMSRG) framework to construct effective Hamiltonians and muon-capture operators, which nonperturbatively account for many-body physics outside the valence space. Here, the obtained nuclear matrix elements are compared against those from the phenomenological shell model. The impact of including the correlations from the nuclear shell model (NSM) as well as including the induced twobody part is studied in detail. Furthermore, the effects of realistic bound-muon wave function on the operators is studied. Finally, predictions for capture rates to the lowest excited states in 24 Na are given and compared with available data. It is found that the spectroscopic properties of 24 Mg and its OMC daughter 24 Na are fairly well described by both the NSM and VS-IMSRG, and that the effect of the hadronic two-body currents significantly reduces the OMC rates. Both models have some difficulties in matching the measured OMC rates, especially for the 2 + final states. This calls for further studies in other light nuclei with available OMC data.

20 ≤ A ≤ 38↗
DOE OSTI2023

Cold neutron-deuteron capture and Wigner-SU(4) symmetry

We calculate the cold neutron-deuteron (nd) capture cross section,σ nd to next-to-next-to leading order (NNLO) using the model-independent approach of pionless effective-field theory [EFT(π)]. At leading order we find σ nd = 0.314 ± 0.217 mb, while the experimental result is 0.508(15) mb for a laboratory neutron velocity of 2200 m/s. At next-to-leading-order (NLO), we show that σnd is sensitive to the low-energy constant (LEC) $L$$^{(0)}_{1}$ of the two-nucleon isovector current appearing at NLO. A fit of $L$$^{(0)}_{1}$ at NLO to the triton magnetic moment yields a NLO prediction of σ nd = 0.393 ± 0.164 mb, where the error comes from propagating the error from the $L$$^{(0)}_{1}$ fit. At NNLO, we find that a new three-nucleon magnetic moment counterterm is required for renormalization-group invariance of both σnd and the triton magnetic moment. Fitting the NNLO correction to $L$$^{(0)}_{1}$ (denoted $L$$^{(1)}_{1}$) to cold neutron-proton capture (σnp) yields a NNLO prediction of σ nd = 0.447 ± 0.130 mb, where the error comes from propagating the error from the $L$$^{(1)}_{1}$ fit. We also study different fittings of $L$$^{(0)}_{1}$ and $L$$^{(1)}_{1}$ to σ np , σ nd , and/or the triton magnetic moment. For example, fitting $L$$^{(0)}_{1}$ simultaneously to σ np , σ nd , and the triton magnetic moment at NLO, and fitting $L$$^{(1)}_{1}$ simultaneously to σ np and σnd at NNLO, yields σ nd = 0.480 ± 0.114 mb and 0.511 ± 0.042 mb, respectively, where errors are naively estimated from EFT(π) power counting. Additionally, we discuss how Wigner SU(4) symmetry may alter the naive EFT(π) expansion of σ nd .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Cross-shell excited configurations in the structure of 34 Si

The cross-shell excited states of 34 Si have been investigated via β decays of the 4 − ground state and the 1 + isomeric state of 34 Al. Since the valence protons and valence neutrons occupy different major shells in the ground state as well as the intruder 1 + isomeric state of 34 Al, intruder levels of 34 Si are populated via allowed β decays. Spin assignments to such intruder levels of 34 Si were established through γ-γ angular correlation analysis for the negative-parity states. The configurations of such intruder states play crucial roles in our understanding of the N = 20 shell gap evolution. A configuration interaction model derived from the FSU Hamiltonian was utilized in order to interpret the intruder states in 34 Si. Furthermore, shell model interaction derived from a more fundamental theory with the valence space in medium similarity renormalization group method was also employed to interpret the structure of 34 Si.

20 ≤ A ≤ 38↗
DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Improved structure of calcium isotopes from ab initio calculations

The in-medium similarity renormalization group (IMSRG) is a powerful and flexible many-body method to compute the structure of nuclei starting from nuclear forces. Recent developments have extended the IMSRG from its standard truncation at the normal-ordered two-body level, the IMSRG(2), to a precision approximation including normal-ordered three-body operators, the IMSRG(3)-N 7 . This improvement provides a more precise solution to the many-body problem and makes it possible to quantify many-body uncertainties in IMSRG calculations. We explore the structure of 44,48,52 Ca using the IMSRG(3)-N 7 , focusing on understanding existing discrepancies of the IMSRG(2) to experimental results. We find a significantly better description of the first 2 + excitation energy of 48 Ca, improving the description of the shell closure at N=28. At the same time, we find that the IMSRG(3)-N 7 corrections to charge radii do not resolve the systematic underprediction of the puzzling large charge radius difference between 52 Ca and 48 Ca. We present estimates of many-body uncertainties of IMSRG(2) calculations applicable also to other systems based on the size extensivity of the method.

39 ≤ A ≤ 58↗
DOE OSTI2025

Variations in the charge radii of indium isotopes between 𝑁=52 and 82

Measurements of the 5⁢𝑝 2 𝑃 3/2 → 9⁢𝑠 2 𝑆 1/2 and 5⁢𝑝 2 𝑃 1/2 → 8⁢𝑠 2 𝑆 1/2 transitions in the indium atom, combined with new atomic physics calculations, were used to extract the changes in mean-square nuclear charge radii, 𝛿⁡⟨𝑟 2 ⟩, of the indium (𝑍 = 49) isotopes 101–111 In , 113–123 In , and 125–131 In . With a proton hole in the closed nuclear shell of 𝑍 = 50, indium provides a detailed study of the effect of unpaired nucleons adjacent to the proton-shell closure, allowing investigation into the charge radii for isotopes between the two major neutron-shell closures at 𝑁 = 50 and 𝑁 = 82. A study of the variations in charge radii between neighboring isotopes with neutron number (the ‘odd-even staggering') is presented and provides further insight and challenges for the theoretical description of the size of proton-hole nuclei. Two nuclear theories, density functional theory and the valence-space in-medium similarity renormalization group method, were employed to interpret the data. The new information obtained in this work provides valuable insights into the successes and shortcomings of the theoretical approaches employed.

ab initio calculations↗
DOE OSTI2020

Renormalization scale setting for heavy quark pair production in $e^+e^–$ annihilation near the threshold region

Heavy fermion pair production in e + e – annihilation is a fundamental process in hadron physics and is of considerable interest for various phenomena. In this paper, we will apply the principle of maximum conformality (PMC) to provide a comprehensive analysis of these processes. The PMC provides a systematic, unambiguous method for determining the renormalization scales of the QCD coupling constant for single-scale and multiple-scale applications. The resulting predictions eliminate any renormalization scheme-and-scale ambiguities, eliminate the factorial renormalon divergences, and are consistent with the requirements of the renormalization group. It is remarkable that two distinctly different scales are determined by using the PMC for heavy fermion pair production near the threshold region. One scale is the order of the fermion mass m f , which enters the hard virtual corrections, and the other scale is of order v m f , where v is the quark velocity, which enters the Coulomb rescattering amplitude. The PMC scales yield the correct physical behavior and reflect the virtuality of the propagating gluons (photons) for the QCD (QED) processes. Moreover, we demonstrate the consistency of PMC scale setting from QCD to QED. Perfect agreement between the Abelian unambiguous Gell-Mann-Low and the PMC scale-setting methods in the limit of zero number of colors is demonstrated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Clock model interpolation and symmetry breaking in O(2) models

The $q$-state clock model is a classical spin model that corresponds to the Ising model when $q=2$ and to the $XY$ model when $q\to\infty$. The integer-$q$ clock model has been studied extensively and has been shown to have a single phase transition when $q=2$,$3$,$4$ and two phase transitions when $q>4$.We define an extended $q$-state clock model that reduces to the ordinary $q$-state clock model when $q$ is an integer and otherwise is a continuous interpolation of the clock model to noninteger $q$. We investigate this class of clock models in 2D using Monte Carlo (MC) and tensor renormalization group (TRG) methods, and we find that the model with noninteger $q$ has a crossover and a second-order phase transition. We also define an extended-$O(2)$ model (with a parameter $\gamma$) that reduces to the $XY$ model when $\gamma=0$ and to the extended $q$-state clock model when $\gamma\to\infty$, and we begin to outline the phase diagram of this model. These models with noninteger $q$ serve as a testbed to study symmetry breaking in situations corresponding to quantum simulators where experimental parameters can be tuned continuously.

Hostetler, Leon↗
DOE OSTI2022

Preparation of the SU(3) lattice Yang-Mills vacuum with variational quantum methods

Studying QCD and other gauge theories on quantum hardware requires the preparation of physically interesting states. The variational quantum eigensolver provides a way of performing vacuum state preparation on quantum hardware. Here in this work, variational quantum eigensolver is applied to pure SU(3) lattice Yang-Mills on a single plaquette and one dimensional plaquette chains. Bayesian optimization and gradient descent were investigated for performing the classical optimization. Ansatz states for plaquette chains are constructed in a scalable manner from smaller systems using domain decomposition and a stitching procedure analogous to the density matrix renormalization group. Small examples are performed on IBM’s superconducting Manila processor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Standard model Higgs inflation supplemented by minimal dark matter

Renormalization group analysis with the present measurements of the top quark mass m t = 172.69 ± 0.30 GeV [1] indicates that the Standard Model (SM) Higgs potential becomes unstable at energy scales ∼ 10 10 GeV . This may be interpreted as hinting at new particles at high energy. The minimal extension of the SM that can avoid this instability while leaving the SM Higgs as the sole scalar particle of the theory is obtained by adding suitable fermions to the SM. These fermions are good dark matter candidates and the model is known as the minimal dark matter model. We revisit the inflationary scenario based on the minimal dark matter model, taking into account updated parameter constraints and recent understanding of reheating dynamics. We explore the model with different values of the right-handed neutrino mass and find that the cosmological prediction is insensitive to such details. We obtained a spectral index of the cosmic microwave background n s = 0.9672 and a tensor-to-scalar ratio r = 0.0031 as a robust prediction of this scenario. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Universal features of 2 → N scattering in QCD and gravity from shockwave collisions

A remarkable double-copy relation of Einstein gravity to QCD in Regge asymptotics is Γ μ ν = 1 2 C μ C ν − 1 2 N μ N ν , where Γ μ ν is the gravitational Lipatov vertex in the 2 → 3 graviton scattering amplitude, C μ its Yang-Mills counterpart, and N μ the QED bremsstrahlung vertex. In QCD, the Lipatov vertex is a fundamental building block of the BFKL equation describing the 2 → N scattering of gluons at high energies. Likewise, the gravitational Lipatov vertex is a key ingredient in a 2D effective field theory framework describing trans-Planckian 2 → N graviton scattering. We construct a quantitative correspondence between a semiclassical Yang-Mills framework for radiation in gluon shockwave collisions and its counterpart in general relativity. In particular, we demonstrate the Lipatov double copy in a dilute-dilute approximation corresponding to R S , L , R S , H ≪ b , where R S , L , R S , H are the respective emergent Schwarzchild radii generated in shockwave collisions and b is the impact parameter. We outline extensions of the correspondence developed here to the dilute-dense computation of gravitational wave radiation in the close vicinity of one of the black holes, the construction of graviton propagators in the shockwave background, and a renormalization group approach to compute 2 → N amplitudes that incorporates graviton Reggeization and coherent graviton multiple scattering. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗