Search

The complex scaling method is commonly used to describe decaying states, but its applications are limited because the Hamiltonian operator must contain only relative coordinates. This has hindered the use of complex scaling in shell model, defined with laboratory single-particle coordinates. In particular, one of the most important model in low-energy nuclear physics, the no-core shell model, had never been considered in the frame of complex scaling. We will then present a straightforward procedure for introducing complex scaling in the no-core shell model in order to calculate nuclear resonance states. For that matter, the complex scaled two-body matrix elements must firstly be determined, and the resulting many-body Hamiltonian complex symmetric matrix must be diagonalized afterwards. Applications pertain to the bound ground states of the lightest nuclei $^2{\rm H}$, $^3{\rm H}$, $^3{\rm He}$, and $^4{\rm He}$, as well as the resonance ground states of $^5{\rm He}$ and $^5{\rm Li}$, whereby the realistic interaction Daejeon16 is utilized., Comment: 6 pages, 2 figures, submitted to Physical Review Letters

In this study, we explore the entanglement and correlation in two-nucleon systems using isospin formalism. With the help of Slater decomposition, we derive analytical expressions for various entanglement measures. Specifically, we analyse the one- and two-mode entropies, mutual informations, and a basis-independent characteristic known as the one-body entanglement entropy. To understand the impact of pairing, we consider interactions involving isovector and isoscalar L=0 pairing terms. Our findings show that certain pairing interactions can maximize one-body entanglement entropy of ground states when both total angular momentum and total isospin have zero projections. We provide numerical examples for the sd shell and explore the mutual informations in LS coupled and jj coupled single-particle bases. We find that the shell structure and angular momentum coupling significantly impact the measures of entanglement. We outline the implications of conserving angular momentum and isospin on one-mode entropies, irrespective of particle number., Comment: 32 pages, 6 figures

We propose a novel many-body framework combining the density matrix renormalization group (DMRG) with the valence-space (VS) formulation of the in-medium similarity renormalization group. This hybrid scheme admits for favorable computational scaling in large-space calculations compared to direct diagonalization. The capacity of the VS-DMRG approach is highlighted in ab initio calculations of neutron-rich nickel isotopes based on chiral two- and three-nucleon interactions, and allows us to perform converged ab initio computations of ground and excited state energies. We also study orbital entanglement in the VS-DMRG, and investigate nuclear correlation effects in oxygen, neon, and magnesium isotopes. The explored entanglement measures reveal nuclear shell closures as well as pairing correlations., Comment: 5 pages, 4 Figures, version published at Phys. Lett. B

We study mode-entanglement in the seniority model, derive analytic formulas for the one-body reduced density matrix of states with seniority $\nu = 0,\;1,\;2$, and $\nu=3$, and also determine the particle number dependence of the one-body reduced density matrix for arbitrary seniority. We carry out numerical calculations for the lightest calcium isotopes and for $^{94}{\rm Ru}$ nucleus, and investigate the structure of their ground and low energy yrast states. We investigate the fulfillment of several predictions of the seniority model regarding the behavior of one-mode entropies, which we compare with the results of configuration interaction (CI) and density matrix renormalization group (DMRG) computations. For $^{94}{\rm Ru}$, the seniority model accounts for the $0g_{9/2}$ mode entropies, but seniority mixing is important for certain yrast states. Interaction induced quantum fluctuations decrease the occupation of the $0f_{5/2}$, $1p_{3/2}$ and $1p_{1/2}$ shells, and amount in finite mode entropies on these shells, too, clearly outside the scope of the simple $(0g_{9/2})^4$ seniority model. The $0f_{7/2}$ shell based seniority model is more accurate for the light ${\rm Ca}$ isotopes, but seniority mixing is substantial for some $^{44}{\rm Ca}$ yrast states, too., Comment: 13 pages, 7 figures

We examine the mode entanglement and correlation of two fermionic particles. We study the one- and two-mode entropy and a global characteristic, the one-body entanglement entropy. We consider not only angular momentum coupled states with single configuration but use the configuration interaction method. With the help of the Slater decomposition, we derive analytical expressions for the entanglement measures. We show that when the total angular momentum is zero specific single configurations describe maximally entangled states. It turns out that for a finite number of associated modes the one- and two-mode entropies have identical values. In the shell model framework, we numerically study two valence neutrons in the $sd$ shell. The one-body entanglement entropy of the ground state is close to the maximal value and the associated modes have the largest mutual information., Comment: 20 pages, 1 figure

New sets of functions with arbitrary large finite cardinality are constructed for two-electron atoms. Functions from these sets exactly satisfy the Kato's cusp conditions. The new functions are special linear combinations of Hylleraas- and/or Kinoshita-type terms. Standard variational calculation, leading to matrix eigenvalue problem, can be carried out to calculate the energies of the system. There is no need for optimization with constraints to satisfy the cusp conditions. In the numerical examples the ground state energy of the He atom is considered.

Experimental data on Coulomb breakup and neutron removal indicate that 31Ne is one of the heaviest halo nuclei discovered so far. The possible ground state of 31Ne is either 3/2- coming from p-wave halo or 1/2+ from s-wave halo. In this work, we develop a treatable model to include deformed wave functions and a dynamical knockout formalism which includes the dependence on the nuclear orientation to study the neutron removal from 31Ne projectiles at energies around E=200 MeV/nucleon. A detailed account of the effects of deformation on cross sections and longitudinal momentum distributions is made. Our numerical analysis indicates a preference for the 31Ne ground state with spin parity 3/2-., Comment: 22 pages, 13 figures, accepted for publication in the Physical Review C

We study radiative capture reactions using quasi-separable potentials. This procedure allows an easier treatment of non-local effects that can be extended to three-body problems. Using this technique, we calculate the neutron and proton radiative capture cross sections on $^{12}$C. The results obtained are shown to be in good agreement with the available experimental data., Comment: 12 pages, 4 figures, accepted for publication in Journal of Physics G: Nuclear and Particle Physics

In coupled-channel models the poles of the scattering S-matrix are located on different Riemann sheets. Physical observables are affected mainly by poles closest to the physical region but sometimes shadow poles have considerable effect, too. The purpose of this paper is to show that in coupled-channel problem all poles of the S-matrix can be calculated with properly constructed complex-energy basis. The Berggren basis is used for expanding the coupled-channel solutions. The location of the poles of the S-matrix were calculated and compared with an exactly solvable coupled-channel problem: the one with the Cox potential. We show that with appropriately chosen Berggren basis poles of the S-matrix including the shadow ones can be determined., Comment: 11 pages, 4 figures, 59 references

The physics of open quantum systems is an interdisciplinary area of research. The nuclear "openness" manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape towards the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a Complex Scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow Shell Model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6He viewed as an $\alpha$+ n + n cluster system. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly-bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method, if it can be applied, can provide an additional information about partial energy widths associated with individual thresholds., Comment: 15 pages, 16 figures

Partial wave expansion of the Coulomb-distorted plane wave is determined and studied. Dominant and sub-dominant asymptotic expansion terms are given and leading order three-dimensional asymptotic form is derived. The generalized hypergeometric function $_2F_2(a,a;a+l+1,a-l;z)$ is expressed with the help of confluent hypergeometric functions and the asymptotic expansion of $_2F_2(a,a;a+l+1,a-l;z)$ is simplified.

The two-body Coulomb scattering problem is solved using the standard complex scaling method. The explicit enforcement of the scattering boundary condition is avoided. Splitting of the scattering wave function based on the Coulomb modified plane wave is considered. This decomposition leads a three-dimensional Schr\"odinger equation with source term. Partial wave expansion is carried out and the asymptotic form of the solution is determined. This splitting does not lead to simplification of the scattering boundary condition if complex scaling is invoked. A new splitting carried out only on partial wave level is introduced and this method is proved to be very useful. The scattered part of the wave function tends to zero at large inter-particle distance. This property permits of easy numerical solution: the scattered part of the wave function can be expanded on bound-state type basis. The new method can be applied not only for pure Coulomb potential butin the presence of short range interaction too., Comment: submitted to Phys. Rev. A

Within the complex-energy configuration interaction framework, we study correlations of valence neutrons to explain the behavior of charge radii in the neutron halo nuclei $^{6,8}$He. We find that the experimentally observed decrease of the charge radius between $^6$He and $^8$He is caused by a subtle interplay between three effects: dineutron correlations, a spin-orbit contribution to the charge radius, and a core swelling effect. We demonstrate that two-neutron angular correlations in the $2^+_1$ resonance of $^6$He differ markedly from the ground-state correlations in $^{6,8}$He., Comment: 5 pages, 5 figures

The quasi-particle energy spectrum of the Hartree-Fock-Bogoliubov (HFB) equations contains discrete bound states, resonances, and non-resonant continuum states. We study the structure of the unbound quasi-particle spectrum of weakly bound nuclei within several methods that do not rely on imposing scattering or outgoing boundary conditions. Various approximations are examined to estimate resonance widths. It is shown that the stabilization method works well for all HFB resonances except for very narrow ones. The Thomas-Fermi approximation to the non-resonant continuum has been shown to be very effective, especially for coordinate-space HFB calculations in large boxes that involve huge amounts of discretized quasi-particle continuum states., Comment: 12 pages,11 figures,submitted to PRC

Mass calculations carried out by Strutinsky's shell correction method are based on the notion of smooth single particle level density. The smoothing procedure is always performed using curvature correction. In the presence of curvature correction a smooth function remains unchanged if smoothing is applied. Two new curvature correction methods are introduced. The performance of the standard and new methods are investigated using harmonic oscillator and realistic potentials., Comment: 4 figures, submitted to Journal of Physics G: Nuclear and Particle Physics

A new method is presented for calculation of the shell correction with the inclusion of the continuum part of the spectrum. The smoothing function used has a finite energy range in contrast to the Gaussian shape of the Strutinski method. The new method is specially useful for light nuclei where the generalized Strutinski procedure can not be applied.

Two methods the complex energy shell model (CXSM) and the complex scaling (CS) approach were used for calculating isobaric analog resonances (IAR) in the Lane model. The IAR parameters calculated by the CXSM and the CS methods were checked against the parameters extracted from the direct numerical solution of the coupled channel Lane equations (CC). The agreement with the CC results was generally better than 1 keV for both methods and for each partial waves concerned. Similarities and differences of the CXSM and the CS methods are discussed. CXSM offers a direct way to study the configurations of the IAR wave function in contrast to the CS method., Comment: 3 figures

Proton emission from deformed nuclei is described within the non-adiabatic weak coupling model which takes into account the coupling to $\gamma$-vibrations around the axially-symmetric shape. The coupled equations are derived within the Gamow state formalism. A new method, based on the combination of the R-matrix theory and the oscillator expansion technique, is introduced that allows for a substantial increase of the number of coupled channels. As an example, we study the deformed proton emitter $^{141}$Ho., Comment: submitted to Phys. Rev. C

For nuclei with N = Z, the isospin degree of freedom is important and, for deformed systems, rotational bands of different isospin may be expected at low excitation energies. We have investigated, in a simple model space, the influence of the isospin-breaking Coulomb interaction on the degree of chaoticity of these rotational bands. The statistical measures used rely on an analysis of level-spacing distributions, which are extremely difficult to measure experimentally. We show, however, that the overlap intergrals between states of similar frequency reflect well the degree of chaoticity. This quantity is closely related to the experimentally more accessible gamma-decay ``spreading width''., Comment: 13 pages, 9 figures, Elsevier

The newly developed nonadiabatic method based on the coupled-channel Schroedinger equation with Gamow states is used to study the phenomenon of proton radioactivity. The new method, adopting the weak coupling regime of the particle-plus-rotor model, allows for the inclusion of excitations in the daughter nucleus. This can lead to rather different predictions for lifetimes and branching ratios as compared to the standard adiabatic approximation corresponding to the strong coupling scheme. Calculations are performed for several experimentally seen, non-spherical nuclei beyond the proton dripline. By comparing theory and experiment, we are able to characterize the angular momentum content of the observed narrow resonance., Comment: 12 pages including 10 figures

Shell corrections of the finite deformed Woods-Saxon potential are calculated using the Green's function method and the generalized Strutinsky smoothing procedure. They are compared with the results of the standard prescription which are affected by the spurious contribution from the unphysical particle gas. In the new method, the shell correction approaches the exact limit provided that the dimension of the single-particle (harmonic oscillator) basis is sufficiently large. For spherical potentials, the present method is faster than the exact one in which the contribution from the particle continuum states is explicitly calculated. For deformed potentials, the Green's function method offers a practical and reliable way of calculating shell corrections for weakly bound nuclei., Comment: submitted to Phys. Rev. C, 12 pages, 7 figures

The coupled-channel Schrodinger equation with outgoing wave boundary conditions is employed to study the fine structure seen in the proton decay of deformed even-N, odd-Z rare earth nuclei 131-Eu and 141-Ho. Experimental lifetimes and proton-decay branching ratios are reproduced. The comparison with the standard adiabatic theory is made., Comment: 4 pages, 2 figures submitted to Phys. Rev. Lett

Mathematically rigorous theory of the two-body contact interaction in three dimension is reviewed. Local potential realizations of this proper contact interaction are given in terms of Poschl-Teller, exponential and square-well potentials. Three body calculation is carried out for the halo nucleus 11Li using adequately represented contact interaction., Comment: submitted to Phys. Rev. C

Shell corrections to the nuclear binding energy as a measure of shell effects in superheavy nuclei are studied within the self-consistent Skyrme-Hartree-Fock and Relativistic Mean-Field theories. Due to the presence of low-lying proton continuum resulting in a free particle gas, special attention is paid to the treatment of single-particle level density. To cure the pathological behavior of shell correction around the particle threshold, the method based on the Green's function approach has been adopted. It is demonstrated that for the vast majority of Skyrme interactions commonly employed in nuclear structure calculations, the strongest shell stabilization appears for Z=124, and 126, and for N=184. On the other hand, in the relativistic approaches the strongest spherical shell effect appears systematically for Z=120 and N=172. This difference has probably its roots in the spin-orbit potential. We have also shown that, in contrast to shell corrections which are fairly independent on the force, macroscopic energies extracted from self-consistent calculations strongly depend on the actual force parametrisation used. That is, the A and Z dependence of mass surface when extrapolating to unknown superheavy nuclei is prone to significant theoretical uncertainties., Comment: 14 pages REVTeX, 8 eps figures, submitted to Phys. Rev. C

A FORTRAN77 program is presented that calculates fusion cross sections and mean angular momenta of the compound nucleus under the influence of couplings between the relative motion and several nuclear collective motions. The no-Coriolis approximation is employed to reduce the dimension of coupled-channels equations. The program takes into account the effects of non-linear couplings to all orders, which have been shown to play an important role in heavy-ion fusion reactions at subbarrier energies., Comment: 14 pages, submitted to CPC, the program is available upon request from hagino@phys.washington.edu

Shell corrections of finite, spherical, one-body potentials are analyzed using a smoothing procedure which properly accounts for the contribution from the particle continuum, i.e., unbound states. Since the plateau condition for the smoothed single-particle energy seldom holds, a new recipe is suggested for the definition of the shell correction. The generalized Strutinsky smoothing procedure is compared with the results of the semi-classical Wigner-Kirkwood expansion. A good agreement has been found for weakly bound nuclei in the vicinity of the proton drip line. However, some deviations remain for extremely neutron-rich systems due to the pathological behavior of the semi-classical level density around the particle threshold., Comment: 18 pages, 8 figures

Ground state energies and decay widths of particle unstable nuclei are calculated within the Hartree-Fock approximation by performing a complex scaling of the many-body Hamiltonian. Through this transformation, the wave functions of the resonant states become square integrable. The method is implemented with Skyrme effective interactions. Several Skyrme parametrizations are tested on four unstable nuclei: 10He, 12O, 26O and 28O., Comment: 5 pages, LaTeX, submitted to Phys. Rev. Lett

$\alpha-$decay through a deformed potential barrier produces significant mixing of angular momenta when mapped from the nuclear interior to the outside. Using experimental branching ratios and either semi-classical or coupled-channels transmission matrices, we have found that there is a set of internal amplitudes which are essentially constant for all even--even actinide nuclei. These same amplitudes also give good results for the known anisotropic $\alpha-$particle emission of the favored decays of odd nuclei in the same mass region. PACS numbers: 23.60.+e, 24.10.Eq, 27.90.+b, Comment: 5 pages, latex (revtex style), 2 embedded postscript figures uuencoded gz-compressed .tar file To appear in Physical Review Letters