We derive the equations to calculate the reduced width amplitudes (RWA) of the different size clusters and deformed clusters without any approximation. These equations named Laplace expansion method are applicable to the nuclear models which uses the Gaussian wave packets. The advantage of the method is demonstrated by the numerical calculations of the ${}^{16}{rm O}+alpha$ and ${}^{24}{rm Mg}+alpha$ RWAs in $^{20}{rm Ne}$ and $^{28}{rm Si}$.
General expressions for the breakup cross sections in the lab frame for $1+2$ reactions are given in terms of the hyperspherical adiabatic basis. The three-body wave function is expanded in this basis and the corresponding hyperradial functions are obtained by solving a set of second order differential equations. The ${cal S}$-matrix is computed by using two recently derived integral relations. Even though the method is shown to be well suited to describe $1+2$ processes, there are nevertheless particular configurations in the breakup channel (for example those in which two particles move away close to each other in a relative zero-energy state) that need a huge number of basis states. This pathology manifests itself in the extremely slow convergence of the breakup amplitude in terms of the hyperspherical harmonic basis used to construct the adiabatic channels. To overcome this difficulty the breakup amplitude is extracted from an integral relation as well. For the sake of illustration, we consider neutron-deuteron scattering. The results are compared to the available benchmark calculations.
The possibility to resolve narrow structures in reaction cross sections in calculations with the Lorentz integral transform (LIT) method is studied. To this end we consider a fictitious two-nucleon problem with a low-lying and narrow resonance in the $^3P_1$ nucleon-nucleon partial wave and calculate the corresponding ``deuteron photoabsorption cross section. In the LIT method the use of continuum wave functions is avoided and one works instead with a localized function tildePsi. In this case study it is investigated how far into the asymptotic region tildePsi has to be determined in order to obtain a precise resolution of the artificially introduced E1 resonance. Comparing with the results of a conventional calculation with explicit neutron-proton continuum wave functions it is shown that the LIT approach leads to an excellent reproduction of the cross section in the resonance region and of further finer cross section details at higher energies. To this end, however, for tildePsi one has to take into account two-nucleon distances up to at least 30 fm.
We perform a detailed comparison of results of the Gamow Shell Model (GSM) and the Gaussian Expansion Method (GEM) supplemented by the complex scaling (CS) method for the same translationally-invariant cluster-orbital shell model (COSM) Hamiltonian. As a benchmark test, we calculate the ground state $0^{+}$ and the first excited state $2^{+}$ of mirror nuclei $^{6}$He and $^{6}$Be in the model space consisting of two valence nucleons in $p$-shell outside of a $^{4}$He core. We find a good overall agreement of results obtained in these two different approaches, also for many-body resonances.
The inversion of nabla Laplace transform, corresponding to a causal sequence, is considered. Two classical methods, i.e., residual calculation method and partial fraction method are developed to perform the inverse nabla Laplace transform. For the first method, two alternative formulae are proposed when adopting the poles inside or outside of the contour, respectively. For the second method, a table on the transform pairs of those popular functions is carefully established. Besides illustrating the effectiveness of the developed methods with two illustrative examples, the applicability are further discussed in the fractional order case.
In order to study structure of proto-neutron stars and those in subsequent cooling stages, it is of great interest to calculate inhomogeneous hot and cold nuclear matter in a variety of phases. The finite-temperature Hartree-Fock-Bogoliubov (FT-HFB) theory is a primary choice for this purpose, however, its numerical calculation for superfluid (superconducting) many-fermion systems in three dimensions requires enormous computational costs. To study a variety of phases in the crust of hot and cold neutron stars, we propose an efficient method to perform the FT-HFB calculation with the three-dimensional (3D) coordinate-space representation. Recently, an efficient method based on the contour integral of Greens function with the shifted conjugate-orthogonal conjugate-gradient method has been proposed [Phys. Rev. C 95, 044302 (2017)]. We extend the method to the finite temperature, using the shifted conjugate-orthogonal conjugate-residual method. We benchmark the 3D coordinate-space solver of the FT-HFB calculation for hot isolated nuclei and fcc phase in the inner crust of neutron stars at finite temperature. The computational performance of the present method is demonstrated. Different critical temperatures of the quadrupole and the octupole deformations are confirmed for $^{146}$Ba. The robustness of the shape coexistence feature in $^{184}$Hg is examined. For the neutron-star crust, the deformed neutron-rich Se nuclei embedded in the sea of superfluid low-density neutrons appear in the fcc phase at the nucleon density of 0.045 fm$^{-3}$ and the temperature of $k_B T=200$ keV. The efficiency of the developed solver is demonstrated for nuclei and inhomogeneous nuclear matter at finite temperature. It may provide a standard tool for nuclear physics, especially for the structure of the hot and cold neutron-star matters.