W. Horiuchi

Erosion of N = 28 shell closure: Shape coexistence and monopole transition

Background: In neutron-rich nuclei neighboring 42Si, the quenching of the N = 28 shell gap occurs and is expected to induce the shape coexistence in their excitation spectra. Purpose: We show that different nuclear shapes coexist in N = 28 isotones 40Mg, 42Si, and 44S, and investigate observables to probe it. Method: Antisymmetrized molecular dynamics with Gogny D1S density functional is applied to describe the shape coexistence phenomena without ad hoc assumption of the nuclear shape. Results: We find that rigid shapes with different deformations coexist in the ground and the first excited 0+ states of 40Mg and 42Si, while in 44S the states exhibit large-amplitude collective motion, which does not have any particular shape. These characteristics are reflected well in the monopole transition strengths. Conclusion: The quenching of the N = 28 shell gap leads to the unique shape coexistence in the N = 28 isotones, which can be probed by the monopole transition strengths.

Large enhancement of total reaction cross sections at the edge of the island of inversion in Ti, Cr, and Fe isotopes

A systematic analysis of nuclear deformation is made for neutron-rich Ti, Cr, and Fe isotopes to explore the nuclear structure in the island of inversion near N = 40, where strong nuclear deformation is predicted. The nuclear ground states are obtained by the Skyrme Hartree-Fock method in three-dimensional coordinate space, which properly describes any nuclear shape. Three types of Skyrme interactions are employed to generate various deformed states in its isotopic chain. We find that in the island of inversion the occupation of highly elongated intruder orbits induces not only large quadrupole deformation but also large hexadecapole deformation. This appears as a sizable enhancement of the nuclear matter radius, showing the characteristic shell effect of the density profile near the nuclear surface. We show that the edge of the island of inversion, where the intruder orbit starts being occupied, can be determined by measuring the enhancement of the total reaction cross section at high incident energy. The possibility of constraining the hexadecapole deformation by a measurement of the total reaction cross sections is discussed.

Pairing core swelling effect in Pb isotopes at $N>126$

We revisit a sudden increase of the isotope shift of the charge radius of Pb isotope at $N>126$ based on a Skyrme Hartree-Fock-Bogoliubov theory. New parametrizations of the pairing interaction optimized for selected four Skyrme interactions greatly improve a description of this phenomenon. The density-dependent spin-orbit interaction is also investigated and further increases the charge radius. The pairing correlations significantly change the properties of the neutron orbits near the Fermi level and play a vital role in pulling out the well-bound protons in the Pb isotopes at $N>126$ . Regarding $^{208}$ Pb as a "core" nucleus, a novel pairing core swelling effect is proposed: The pairing interaction reduces the radius of "valence" neutron orbits by the shrinkage of diffused $1g_{9/2}$ orbit and the mixing of sharp $0i_{11/2}$ orbit. Simultaneously, the core nucleus swells, leading to the sudden enhancement of the charge radius at $N > 126$ . This characteristic behavior appears in the density profile near the nuclear surface and its measurement is highly desirable.

Possible halo structure of $^{62,72}$ Ca by forbidden-state-free locally peaked Gaussians

In order to efficiently describe nucleon orbits around a heavy core nucleus, we propose locally peaked Gaussians orthogonalized to the occupied bound states in the core. We show the advantage of those functions in both numerical stability and fast convergence by taking examples of touchstone calcium isotopes $^{62,72}$ Ca in $^{60,70}{\rm Ca}+n+n$ three-body models. Both weakly bound configurations and continuum coupling effect are taken into account. We evaluate the neutron radii and the occupation probabilities of two-neutron configurations not only for the ground state but also for some particle-bound excited states by varying the strength of the core-neutron interaction. The emergence of the halo structure in the ground state depends on the energy difference between $2s_{1/2}$ and $0g_{9/2}$ orbits. Two-neutron [consisting of $(s_{1/2})^2$ configuration] and one-neutron [consisting of $(g_{9/2}s_{1/2})$ configuration] halo structure of $^{62}$ Ca can coexist in narrow energy spacing provided that both of $2s_{1/2}$ and $0g_{9/2}$ orbits are almost degenerate and barely bound. The ground-state structure of $^{72}$ Ca is likely to be a two-neutron halo, although its emergence depends on the position of the $2s_{1/2}$ level.

preprint2021arXiv

Structure of two- and three-alpha systems in cold neutron matter

We present stability and structure of two- and three-alpha systems embedded in dilute cold neutron matter. By solving a few-alpha Schrödinger equation with quasiparticle properties, i.e., effective mass and induced two- and three-alpha interactions, which are evaluated in terms of Fermi polarons, it is shown that $^8 \mathrm{Be}$ and the Hoyle state become bound at densities of about $10^{-4}$ and $10^{-3}$ of the saturation density, respectively. It is also seen that, under cold neutron matter environment, both systems become smaller than the corresponding systems in vacuum. Our results would affect astrophysical models for stellar collapse and neutron star mergers, as well as relevant reaction rates for nucleosynthesis.

A two-neutron halo is unveiled in $^{29}$ F

We report the measurement of reaction cross sections ( $σ_R^{\rm ex}$ ) of $^{27,29}$ F with a carbon target at RIKEN. The unexpectedly large $σ_R^{\rm ex}$ and derived matter radius identify $^{29}$ F as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the $2p_{3/2}$ orbital, thereby vanishing the shell closure associated with the neutron number $N = 20$ . The results are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of $^{27}$ F but are challenged for $^{29}$ F.

Core swelling in spherical nuclei: An indication of the saturation of nuclear density

Background: Nuclear radius is one of the most important and basic properties of atomic nuclei and its evolution is closely related to the saturation of the nuclear density in the internal region but the systematics of the nuclear radii for the neutron-rich unstable nuclei is not well known. Purpose: Motivated by the recent interaction cross section measurement which indicates the 48Ca core swelling in the neutron-rich Ca isotopes, we explore the mechanism of the enhancement of the neutron and proton radii for spherical nuclei. Methods: Microscopic Hartree-Fock calculations with three sets of Skyrme-type effective interactions are performed for the neutron-rich Ca, Ni and Sn isotopes. The total reaction cross sections for the Ca isotopes are evaluated with the Glauber model to compare them with the recent cross section data. Results: We obtain good agreement with the measured cross sections and charge radii. The neutron and proton radii of the various "core" configurations are extracted from the full Hartree-Fock calculation and discuss the core swelling mechanism. Conclusions: The core swelling phenomena occur depending on the properties of the occupying valence single-neutron states to minimize the energy loss that comes from the saturation of the densities in the internal region, which appears to be prominent in light nuclei such as Ca isotopes.

Coulomb screening correction to the $Q$ value of the triple alpha process in thermal plasmas

The triple alpha reaction is a key to $^{12}$ C production and is expected to occur in weakly-coupled, thermal plasmas as encountered in normal stars. We investigate how Coulomb screening affects the structure of a system of three alpha particles in such a plasma environment by precise three-body calculations within the Debye-Hückel approximation. A three-alpha model that has the Coulomb interaction modified in the Yukawa form is employed. Precise three-body wave functions are obtained by a superposition of correlated Gaussian bases with the aid of the stochastic variational method. The energy shifts of the Hoyle state due to the Coulomb screening are obtained as a function of the Debye screening length. The results, which automatically incorporate the finite size effect of the Hoyle state, are consistent with the conventional result based on the Coulomb correction to the chemical potentials of ions that are regarded as point charges in a weakly-coupled, thermal plasma. We have given a theoretical basis to the conventional point-charge approach to the Coulomb screening problem relevant for nuclear reactions in normal stars by providing the first evaluation of the Coulomb corrections to the $Q$ value of the triple alpha process that produces a finite size Hoyle state.

Electric dipole response of low-lying excitations in the two-neutron halo nucleus $\boldsymbol{^{29}}$ F

The neutron-rich $^{28,29}$ F isotopes have been recently studied via knockout and interaction cross-section measurements. The $2n$ halo in $^{29}$ F has been linked to the occupancy of $pf$ intruder configurations. We investigate bound and continuum states in $^{29}$ F, focusing on the $E1$ response of low-lying excitations and the effect of dipole couplings on nuclear reactions. $^{29}\text{F}$ ( $^{27}\text{F}+n+n$ ) wave functions are built within the hyperspherical harmonics formalism, and reaction cross sections are calculated using the Glauber theory. Continuum states and $B(E1)$ transition probabilities are described in a pseudostate approach using the analytical THO basis. The corresponding structure form factors are used in CDCC calculations to describe low-energy scattering. Parity inversion in $^{28}$ F leads to a $^{29}$ F ground state characterized by 57.5% of $(p_{3/2})^2$ intruder components, a strong dineutron configuration, and an increase of the matter radius with respect to the core radius of $ΔR=0.20$ fm. Glauber-model calculations for a carbon target at 240 MeV/nucleon provide a total reaction cross section of 1370 mb, in agreement with recent data. The model produces also a barely bound excited state corresponding to a quadrupole excitation. $B(E1)$ calculations into the continuum yield a total strength of 1.59 e $^2$ fm $^2$ up to 6 MeV, and the $E1$ distribution exhibits a resonance at $\approx$ 0.85 MeV. Results using a standard shell-model order for $^{28}$ F lead to a considerable reduction of the $B(E1)$ distribution. The four-body CDCC calculations for $^{29}\text{F}+^{120}\text{Sn}$ around the Coulomb barrier are dominated by dipole couplings, which totally cancel the Fresnel peak in the elastic cross section. These results are consistent with a two-neutron halo and may guide future experimental campaigns.

Exploring two-neutron halo formation in the ground-state of $^{29}$ F within a three-body model

Background $\colon$ The $^{29}$ F system is located at the lower-N boundary of the "island of inversion" and is an exotic, weakly bound system. Little is known about this system beyond its two-neutron separation energy ( $S_{2n}$ ) with large uncertainties. A similar situation is found for the low-lying spectrum of its unbound binary subsystem $^{28}$ F. Purpose $\colon$ To investigate the configuration mixing, matter radius and neutron-neutron correlations in the ground state of $^{29}$ F within a three-body model, exploring the possibility of $^{29}$ F to be a two-neutron halo nucleus. Method $\colon$ The $^{29}$ F ground-state wave function is built within the hyperspherical formalism by using an analytical transformed harmonic oscillator basis. The Gogny-Pires-Tourreil (GPT) nn interaction with central, spin-orbit and tensor terms is employed in the present calculations, together with different core $+n$ potentials constrained by the available experimental information on $^{28}$ F. Results $\colon$ The $^{29}$ F ground-state configuration mixing and its matter radius are computed for different choices of the $^{28}$ F structure and $S_{2n}$ value. The admixture of d-waves with pf components are found to play an important role, favoring the dominance of dineutron configurations in the wave function. Our computed radii show a mild sensitivity to the $^{27}$ F $+n$ potential and $S_{2n}$ values. The relative increase of the matter radius with respect to the $^{27}$ F core lies in the range 0.1-0.4 fm depending upon these choices. Conclusions $\colon$ Our three-body results for $^{29}$ F indicate the presence of a moderate halo structure in its ground state, which is enhanced by larger intruder components. This finding calls for an experimental confirmation.

The $\boldsymbol{^{29}}$ F nucleus as a lighthouse on the coast of the island of inversion

The exotic, neutron-rich and weakly-bound isotope $^{29}$ F stands out as a waymarker on the southern shore of the island of inversion, a portion of the nuclear chart where the effects of nuclear forces lead to a reshuffling of the single particle levels and to a reorganization of the nuclear structure far from stability. This nucleus has become very popular, as new measurements allow to refine theoretical models. We review the latest developments and suggest how to further assess the structure by proposing predictions on electromagnetic transitions that new experiments of Relativistic Coulomb Excitation should soon become able to measure.

Two-neutron halo structure of $^{31}$ F and a novel pairing anti-halo effect

Background: A newly identified dripline nucleus $^{31}$ F offers a unique opportunity to study the two-neutron ( $2n$ ) correlation at the east shore of the island of inversion where the $N = 28$ shell closure is lost. Purpose: We aim to present the first three-body theoretical results for the radius and total reaction cross sections of $^{31}$ F. This will further help to investigate how the pairing and breakdown of the $N = 28$ shell closure influence the formation of the $2n$ -halo structure and the anti-halo effect in this mass region. Methods: A $^{29}$ F $+n+n$ three-body system is described by the cluster orbital shell model, and its total reaction cross section is calculated by the Glauber theory. Results: Our three-body calculations predict 3.48-3.70 fm for the root-mean-square radius of $^{31}$ F, which corresponds to the total reaction cross section of 1530 (1410)-1640 (1500) mb for a carbon target at 240 (900) MeV/nucleon. The binding mechanism and halo formation in $^{31}$ F are discussed. Conclusions: The present study suggests a novel anti-halo effect in this mass region: When the pairing overcome the energy gap between the $p_{3/2}$ and $f_{7/2}$ orbits, the inversion of the occupation number of these orbits takes place, and it diminishes the $2n$ -halo structure.

preprint2019arXiv

Impact of uncertainties of unbound $^{10}$ Li on the ground state of two-neutron halo $^{11}$ Li

Recently, the energy spectrum of $^{10}$ Li was measured upto 4.6 MeV, via one-neutron transfer reaction d( $^{9}$ Li, p) $^{10}$ Li. Considering the ambiguities on the $^{10}$ Li continuum spectrum with reference to new data, we report the configuration mixing in the ground state of the two-neutron halo nucleus $^{11}$ Li for two different choices of the $^{9}$ Li $+$ n potential. For the present study, we employ a three-body (core $+$ n $+$ n) structure model developed for describing the two-neutron halo system by explicit coupling of unbound continuum states of the subsystem (core $+$ n), and discuss the two-neutron correlations in the ground state of $^{11}$ Li.

preprint2019arXiv

Swelling of doubly magic $^{48}$ Ca core in Ca isotopes beyond $N=28$

Interaction cross sections for $^{42\textrm{-}51}$ Ca on a carbon target at 280 MeV/nucleon have been measured for the first time. The neutron number dependence of derived root-mean-square matter radii shows a significant increase beyond the neutron magic number $N=28$ . Furthermore, this enhancement of matter radii is much larger than that of the previously measured charge radii, indicating a novel growth in neutron skin thickness. A simple examination based on the Fermi-type distribution, and the Mean-Field calculations point out that this anomalous enhancement of the nuclear size beyond $N=28$ results from an enlargement of the core by a sudden increase in the surface diffuseness of the neutron density distribution, which implies the swelling of the bare $^{48}$ Ca core in Ca isotopes beyond $N=28$ .

preprint2019arXiv

Three-body Description of $2n$ -Halo and Unbound $2n$ -Systems: $^{22}$ C and $^{26}$ O

We study the two-neutron correlations in the ground state of the weakly-bound two-neutron halo nucleus $^{22}$ C sitting at the edge of the neutron-drip line and also in the unbound nucleus $^{26}$ O sitting beyond the neutron dripline. For the present study, we employ a three-body (core $+n+n$ ) structure model designed for describing the two-neutron halo system by explicit inclusion of unbound continuum states of the subsystem (core $+n$ ). We use either a density-independent or a density-dependent contact-delta interaction to describe the neutron-neutron interaction and its strength is varied to fix the binding energy. We report the configuration mixing in the ground state of these systems for different choices of pairing interactions.

preprint2016arXiv

Neutrino Induced 4He Break-up Reaction -- Application of the Maximum Entropy Method in Calculating Nuclear Strength Function

The maximum entropy method is examined as a new tool for solving the ill-posed inversion problem involved in the Lorentz integral transformation (LIT) method. As an example, we apply the method to the spin-dipole strength function of 4He. We show that the method can be successfully used for inversion of LIT, provided the LIT function is available with a sufficient accuracy.

preprint2016arXiv

Proton distribution radii of $^{12-19}$ C illuminate features of neutron halos

Proton radii of $^{12-19}$ C densities derived from first accurate charge changing cross section measurements at 900 $A$ MeV with a carbon target are reported. A thick neutron surface evolves from $\sim$ 0.5 fm in $^{15}$ C to $\sim$ 1 fm in $^{19}$ C. The halo radius in $^{19}$ C is found to be 6.4 $\pm$ 0.7 fm as large as $^{11}$ Li. Ab initio calculations based on chiral nucleon-nucleon and three-nucleon forces reproduce well the radii.

Correlated-basis description of $α$ -cluster and delocalized $0^+$ states in $^{16}$ O

A five-body calculation of $^{12}$ C+ $n$ + $n$ + $p$ + $p$ is performed to take a step towards solving an outstanding problem in nuclear theory: The simultaneous and accurate description of the ground and first excited 0 $^+$ states of $^{16}$ O. The interactions between the constituent particles are chosen consistently with the energies of bound subsystems, especially $^{12}$ C $+n$ , $^{12}$ C $+p$ , and $α$ -particle. The five-body dynamics is solved with the stochastic variational method on correlated Gaussian basis functions. No restriction is imposed on the four-nucleon configurations except the Pauli principle excluding the occupied orbits in $^{12}$ C. The energies of both the ground and first excited states of $^{16}$ O are obtained in excellent agreement with experiment. Analysis of the wave functions indicates spatially localized $α$ -particle-like cluster structure for the excited state and shell-model-like delocalized structure for the ground state.

Electric dipole response of $^6$ He: Halo-neutron and core excitations

Electric dipole ( $E1$ ) response of $^{6}$ He is studied with a fully microscopic six-body calculation. The wave functions for the ground and excited states are expressed as a superposition of explicitly correlated Gaussians (CG). Final state interactions of three-body decay channels are explicitly taken into account. The ground state properties and the low-energy $E1$ strength are obtained consistently with observations. Two main peaks as well as several small peaks are found in the $E1$ strength function. The peak at the high-energy region indicates a typical macroscopic picture of the giant dipole resonance, the out-of-phase proton-neutron motion. The transition densities of the lower-lying peaks exhibit in-phase proton-neutron motion in the internal region, out-of-phase motion near the surface region, and spatially extended neutron oscillation, indicating a soft-dipole mode (SDM) and its vibrationally excited mode.

Harmonic-oscillator excitations of precise few-body wave functions

A method for calculating the occupation probability of the number of harmonic oscillator (HO) quanta is developed for a precise few-body wave function obtained in a correlated Gaussian basis. The probability distributions of two- to four-nucleon wave functions obtained using different nucleon- nucleon (NN) interactions are analyzed to gain insight into the characteristic behavior of the various interactions. Tensor correlations as well as short-range correlations play a crucial role in enhancing the probability of high HO excitations. For the excited states of 4He, the interaction dependence is much less because high HO quanta are mainly responsible for describing the relative motion function between the 3N+N (3H+p and 3He+n) clusters.

Mean-field analysis of ground state and low-lying electric dipole strength in $^{22}$ C

Properties of neutron-rich $^{22}$ C are studied using the mean-field approach with Skyrme energy density functionals. Its weak binding and large total reaction cross section, which are suggested by recent experiments, are simulated by modifying the central part of Skyrme potential. Calculating $E1$ strength distribution by using the random-phase approximation, we investigate developments of low-lying electric dipole ( $E1$ ) strength and a contribution of core excitations of $^{20}$ C. As the neutron Fermi level approaches the zero energy threshold ( $\varepsilon_F >\sim -1$ MeV), we find that the low-lying $E1$ strength exceeds the energy-weighted cluster sum rule, which indicates an importance of the core excitations with the $1d_{5/2}$ orbit.

Probing neutron-skin thickness with total reaction cross sections

We analyze total reaction cross sections, $σ_R$ , for exploring their sensitivity to the neutron-skin thickness of nuclei. We cover 91 nuclei of O, Ne, Mg, Si, S, Ca, and Ni isotopes. The cross sections are calculated in the Glauber theory using the density distributions obtained with the Skyrme-Hartree-Fock method in 3-dimensional coordinate space. Defining a reaction radius, $a_R=\sqrt{σ_R/π}$ , to characterize the nuclear size and target (proton or $^{12}$ C) dependence, we find an empirical formula for expressing $a_R$ with the point matter radius and the skin thickness, and assess two practical ways of determining the skin thickness from proton-nucleus $σ_R$ values measured at different energies or from $σ_R$ values measured for different targets.

preprint2013arXiv

Spin-dipole strength functions of $^4$ He with realistic nuclear forces

Both isoscalar and isovector spin-dipole excitations of $^4$ He are studied using realistic nuclear forces in the complex scaling method. The ground state of $^4$ He and discretized continuum states with $J^π=0^-, 1^-, 2^-$ for A=4 nuclei are described in explicitly correlated Gaussians reinforced with global vectors for angular motion. Two- and three-body decay channels are specifically treated to take into account final state interactions. The observed resonance energies and widths of the negative-parity levels are all in fair agreement with those calculated from both the spin-dipole and electric-dipole strength functions as well as the energy eigenvalues of the complex scaled Hamiltonian. Spin-dipole sum rules, both non energy-weighted and energy-weighted, are discussed in relation to tensor correlations in the ground state of $^4$ He.

Ab initio study of the photoabsorption of $^4$ He

There are some discrepancies in the low energy data on the photoabsorption cross section of $^4$ He. We calculate the cross section with realistic nuclear forces and explicitly correlated Gaussian functions. Final state interactions and two- and three-body decay channels are taken into account. The cross section is evaluated in two methods: With the complex scaling method the total absorption cross section is obtained up to the rest energy of a pion, and with the microscopic $R$ -matrix method both cross sections $^4$ He( $γ, p$ ) $^3$ H and $^4$ He( $γ, n$ ) $^3$ He are calculated below 40\,MeV. Both methods give virtually the same result. The cross section rises sharply from the $^3$ H+ $p$ threshold, reaching a giant resonance peak at 26--27\,MeV. Our calculation reproduces almost all the data above 30\,MeV. We stress the importance of $^3$ H+ $p$ and $^3$ He+ $n$ cluster configurations on the cross section as well as the effect of the one-pion exchange potential on the photonuclear sum rule.

From nucleons to nuclei to fusion reactions

Nuclei are prototypes of many-body open quantum systems. Complex aggregates of protons and neutrons that interact through forces arising from quantum chromo-dynamics, nuclei exhibit both bound and unbound states, which can be strongly coupled. In this respect, one of the major challenges for computational nuclear physics, is to provide a unified description of structural and reaction properties of nuclei that is based on the fundamental underlying physics: the constituent nucleons and the realistic interactions among them. This requires a combination of innovative theoretical approaches and high-performance computing. In this contribution, we present one of such promising techniques, the ab initio no-core shell model/resonating-group method, and discuss applications to light nuclei scattering and fusion reactions that power stars and Earth-base fusion facilities.

Glauber-model analysis of total reaction cross sections for Ne, Mg, Si, and S isotopes with Skyrme-Hartree-Fock densities

A systematic analysis is made on the total reaction cross sections for Ne, Mg, Si, and S isotopes. The high-energy nucleus-nucleus collision is described based on the Glauber model. Using the Skyrme-Hartree-Fock method in the three-dimensional grid-space representation, we determine the nuclear density distribution for a wide range of nuclei self-consistently without assuming any spatial symmetry. The calculated total reaction cross sections consistently agree with the recent cross section data on Ne $+^{12}$ C collision at 240 $A$ \,MeV, which makes it possible to discuss the radius and deformation of the isotopes. The total reaction cross sections for Mg $+^{12}$ C, Si $+^{12}$ C and S $+^{12}$ C cases are predicted for future measurements. We also find that the high-energy cross section data for O, Ne, and Mg isotopes on a $^{12}$ C target at around 1000\, $A$ MeV can not be reproduced consistently with the corresponding data at 240\, $A$ MeV.