-
$ {}^9_{\Lambda} $ Be within the$ \alpha+\alpha+\Lambda $ three cluster system is described by the GCM, the configurations of$ {}^5_{\Lambda} $ He and$ {}^9_{\Lambda} $ Be are shown in Fig. 1. The Pauli principle between two α clusters is taken into account with the GCM, whereas antisymmetrization between Λ hyperon and other nucleons is not required. Therefore the wave functions of$ {}^5_{\Lambda} $ He and$ {}^9_{\Lambda} $ Be can be written as:Figure 1. (color online) The coordinate systems and angular momentum labels for (a)
$ {}^5_{\Lambda}\mathrm{He} $ and (b)$ {}^9_{\Lambda}\mathrm{Be} $ $ \begin{aligned}[b] &\Psi_J({}^5_{\Lambda}\mathrm{He})=\sum_{d}f_J(d)\phi_{\alpha}\phi_J(d,r)Y_J(\hat{{r}}) \\&\Psi_J^{}({}^9_{\Lambda}\mathrm{Be})=\sum_{lL}\sum_{dD}f_{lLJ}(d,D)[\Phi_{l,\alpha\alpha}(d)\phi_L(D,R)Y_L(\hat{{R}})]_J \end{aligned} $
(1) where
$ \phi_{\alpha} $ and$ \Phi_{l,\alpha\alpha} $ represent the wave functions of α particle and$ {}^8_{} $ Be with total angular momentum l, respectively. Some details of the calculation for the$ \alpha + \alpha $ part can be found in the Appendix or in [17, 72, 73].The size parameters b appropriate for the coordinates in Fig. 1 are chosen so that combinations with the corresponding reduced massed lead to an identical value of
$ \hbar/\omega $ :$ \begin{aligned}[b] M_Nb_N^2=\;&M_{\Lambda}b_{\Lambda}^2=\mu_{\alpha \alpha}b_{\alpha \alpha}^2=\mu_{\Lambda \alpha}b_{\Lambda \alpha}^2=\mu_{\Lambda, \alpha\alpha}b_{\Lambda,\alpha\alpha}^2\\=\;&\mu_{\Lambda\alpha, \alpha}b_{\Lambda\alpha, \alpha}^2=\hbar/\omega \end{aligned} $
(2) The total Hamiltonian of
$ {}^9_{\Lambda} $ Be can be divided as:$ \begin{array}{*{20}{l}} \begin{aligned} & H^{}=H_{{}^8_{}Be}+H_R^{} \end{aligned} \end{array} $
(3) with
$ \begin{aligned}[b] & H_{{}^8_{}Be}=T-T_G+V_N+V_C\\&H_R^{}=T_R^{}+\sum_{i=1}^8V_{\Lambda N}^i \end{aligned} $
(4) where T is the total kinetic energy of 8 nucleons in
$ {}^8 $ Be,$ T_G $ is the center-of-mass kinetic energy of$ {}^8 $ Be.$ T_R $ is the kinetic energy associated with the coordinate R.$ V_N $ is the effective two-body nuclear interaction energy and$ V_C $ is the Coulomb interaction energy.$ V_{\Lambda N}^i $ denotes the$ \Lambda N $ interaction between the hyperon Λ and the i-th nucleon.Practically, we adopt the Volkov (or Minnesota) interaction as the effective two-body nuclear potential which has the general form :
$ \begin{aligned}[b] V_N=\;&\dfrac{1}{2}\sum\limits_{i\neq{j}}^8\{(W-MP_{\sigma\tau}+BP_{\sigma}-HP_{\tau})\}_{ij}\\&\times\sum\limits_{k=1}^{k_{max}}V^0_{k} \exp\left(-\dfrac{r_{ij}^2}{a_k^2}\right), \end{aligned} $
(5) where the parameters of Volkov and Minnesota interactions are all listed in Table 1. In our calculations Volkov No.1 is utilized as the
$ NN $ interaction.Interaction k $ V^0_{k} $ $ a_k $ $ W_k $ $ M_k $ $ B_k $ $ H_k $ Volkov No.1 1 -83.34 1.60 1-M M 0 0 2 144.86 0.82 1-M M 0 0 Volkov No.2 1 -60.65 1.80 1-M M 0 0 2 61.14 1.01 1-M M 0 0 Minnesota 1 200 1/ $ \sqrt{1.487} $ u/2 1-u/2 0 0 2 -178 1/ $ \sqrt{0.639} $ u/4 1/2-u/4 u/4 1/2-u/4 3 -91.85 1/ $ \sqrt{0.465} $ u/4 1/2-u/4 -u/4 u/4-1/2 Table 1. Amplitudes
$ V^0_{k} $ (in MeV) and ranges$ a_k $ (in fm) of the Volkov No.1, Volkov No.2 and Minnesota interactions. The standard values are$ M=0.6 $ and$ u=1 $ , but these parameters can be modified in order to reproduce some important properties of the system.The Coulomb interaction can be written as
$ V_C=\dfrac{1}{2}\sum\limits_{i\neq{j}}^8\left(\dfrac{1}{2}+t_{zi}\right)\left(\dfrac{1}{2}+t_{zj}\right)\dfrac{e^2}{r_{ij}}, $
(6) where the isospin z-component equals
$ t_z=+{1}/{2} $ for the proton and$ t_z=-{1}/{2} $ for the neutron.The two-body
$ \Lambda N $ interaction is chosen as the YNG interaction, which is given by:$ \begin{aligned}[b] V_{\Lambda N}(r)=\;&\sum_i\Bigg\{(V_D^0+V_{EX}^0 P_r)\exp\left[-\left(\dfrac{r}{\beta_{i}}\right)^2\right]\\&+(V_{D}^{\sigma\sigma}\sigma_{\Lambda}\sigma_{N} +V_{EX}^{\sigma\sigma}\sigma_{\Lambda}\sigma_{N} P_r)\exp\left[-\left(\dfrac{r}{\beta_{i}}\right)^2\right]\Bigg\} \end{aligned} $
(7) where
$ P_r $ is the space exchange parameter and$ \begin{array}{*{20}{l}} \begin{aligned} \begin{cases} &V_D^{0}=\dfrac{V^0(E)+V^0(O)}{2},\ V_{EX}^{0}=\dfrac{V^0(E)-V^0(O)}{2}\\ \\\; &V_{D}^{\sigma\sigma}=\dfrac{V^0_{\sigma\sigma}(E)+V^0_{\sigma\sigma}(O)}{2},\ V_{EX}^{\sigma\sigma}=\dfrac{V^0_{\sigma\sigma}(E)-V^0_{\sigma\sigma}(O)}{2}\\ \end{cases} \end{aligned} \end{array} $
(8) The parameters of the YNG model used in this work are taken from [23] and listed in Table 2. With the parameters adopted in
$ \Lambda N $ potential, the binding energy of$ {}^{5}_{\Lambda} $ He can be reproduced as 3.10 MeV, which is in good agreement with the experimental value 3.12$ \pm $ 0.02 MeV [74].$ \beta_i $ $ V^0(E) $ $ V^0(O) $ $ V^0_D $ $ V^0_{EX} $ 1.5 −9.93 −7.66 −8.795 −1.135 0.9 −227.73 −82.55 −155.140 −72.590 0.5 1021.17 717.40 869.285 151.885 Table 2.
$ \Lambda N $ interaction depth of the YNG model. The Fermi momentum$ k_F $ is 0.9 fm−1.The following restrictions on the channels and the generator coordinates are adopted in the calculation of
$ {}^{5}_{\Lambda} $ He:$ \begin{array}{*{20}{l}} \begin{aligned} &{channel}\ \ \ (J)=(0)\\ &\alpha-\Lambda:\ \ \ d=0.01, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5 \ \rm{fm} \end{aligned} \nonumber \end{array} $
In this work the Λ particle is treated as spinless, since the
$ \Lambda N $ spin-spin interaction does not contribute to the coupling potential between Λ particle and spin-saturated α cluster. In addition, the Λ one-body spin-orbit potential has been found to be very weak in the experiment.In GCM framework the complex scaling can be introduced by letting the generator coordinate be complex scaled one. Such transformation has been used in [17] for two cluster system
$ \alpha+\alpha $ . Here, for three cluster system the transformation is extended to$ d_{\theta}=e^{i\theta}d $ and$ D_{\theta}=e^{i\theta}D $ . In addition, because of the non-hermitian property caused by CSM we need to use c-product (bi-orthogonal product) [75] instead of the normal inner product to calculate the matrix elements. More detailed information regarding the implementation of the CSM in the framework of GCM can be found in [18]. -
There are many practical methods to select the resonance states among positive energy eigenstates. One of them is the analytic continuation in the coupling constant (ACCC) method, which is based on the analytic continuation from a bound-state pole to S-matrix pole of a resonant state. As a variant of ACCC, inverse ACCC (IACCC) [76, 77] method can also complete the task of analytic continuation.
Explicitly speaking, this can be achieved by adding an artificial attractive pseudo potential V to the original Hamiltonian H:
$ \begin{array}{*{20}{l}} H(\delta)=H+\delta \times V, \end{array} $
(9) where δ is a coupling constant to vary the strength of the pseudo potential, or from another point of view, if the original Hamiltonian can be written as
$ \begin{array}{*{20}{l}} \begin{aligned} H(\lambda)=H_0+\lambda V_{att}, \nonumber \end{aligned} \end{array} $
where
$ V_{att} $ is an attractive potential and$ \lambda=1 $ represents the physical value. Then we have$ \delta=\lambda-1 $ and$ V=V_{att} $ . When there is a bound state with a value of λ, the binding energy will decrease with the shrinking of the parameter λ and this bound state will reach the threshold at$ \lambda = \lambda_0 $ , where$ E(\lambda = \lambda_0)=0 $ . Beyond this critical point the bound state will become a resonant or a virtual state.For a two-body system
$ (n=2) $ it can be shown that the square root of the binding energy$ k_l=\sqrt{-E} $ , behaves near the threshold$ \lambda = \lambda_0 $ as$ \begin{array}{*{20}{l}} k_l= \begin{cases} \sqrt{\lambda-\lambda_0}, \ l\neq 0\\ \lambda-\lambda_0,\ l=0 \end{cases} , \end{array} $
(10) Therefore, one can consider
$ k_l $ as a function of complex variable x:$ \begin{array}{*{20}{l}} x= \begin{cases} \sqrt{\lambda-\lambda_0}, \ l\neq 0\\ \lambda-\lambda_0,\ l=0 \end{cases} , \end{array} $
(11) and
$ k_l $ analytically continue from the$ \lambda>\lambda_0 $ region to the$ \lambda<\lambda_0 $ region.In practice we generally use the Padé approximant to formulate
$ k_l $ :$ \begin{array}{*{20}{l}} \begin{aligned} &k_l^{[N_1,N_2]}(x)=\dfrac{a_1x+a_2x^2+\cdots+a_{N_1}x^{N_1}}{1+b_1x+b_2x^2+\cdots+b_{N_2}x^{N_2}} . \end{aligned} \end{array} $
(12) In general, for systems with
$ n>2 $ ,$ k_l $ is the relative momenta to the nearest desintegration threshold$ k_l=\sqrt{-(E_n-E_{\rm{thr}})} $ . In the calculations, we take the threshold energy as the sum of the intrinsic energies of two free α particles, namely$ E_{thr}=2E_{\alpha} $ . Contrary to the$ n=2 $ case, for$ n>2 $ systems, l can not determine anymore if a bound state turns into a resonant state or a virtual one when$ \lambda<\lambda_0 $ .In this work, the EPD-ACCC method applied is based on the effective N-N interaction adopted in microscopic cluster model:
$ \begin{array}{*{20}{l}} \begin{aligned} &V_N(r)=F(r)(W-MP_{\sigma\tau}+BP_{\sigma}-HP_{\tau}), \end{aligned} \end{array} $
(13) where
$ F(r) $ is often the form of Gaussian function.In the EPD-ACCC method the parameter M(or u) becomes the new coupling constant. In conventional ACCC method, with the increasing of the coupling constant δ (or λ) the attraction effect of the potential must be stronger and thus the binding energy of the bound state will increase gradually. This is also true if parameter u in Minnesota interaction is adopted as the coupling constant. However, on the contrary, for Volkov interaction we can find that the binding energy of the bound state will decrease with the increasing of the coupling constant M. As M decreases from its physical value (0.573 in our calculations) the system becomes increasingly bound, and the resonant energy approaches the threshold. Further decreasing M causes the system to cross the threshold and become bound, allowing us to obtain a series of bound states and their corresponding M parameters. The complex variable that we should utilize is
$ \sqrt{M-M_0} $ , where$ E(M_0)-E_{thr}=0 $ and$ M<M_0 $ for bound state.Compared to introducing an auxiliary attractive potential, using parameters from the original interaction as the coupling constant has the advantage that the results do not depend on the choice of the attractive potential. Additionally, for single α cluster the nuclear potential energy is only dependent on
$ M+W $ and$ M+W $ is generally restricted to be a constant. Therefore the internal energy of a single free α particle is independent on the choice of parameters M and u in Volkov and Minnesota interactions, respectively. This is very convenient when we use the ACCC method to estimate the resonant states in$ \alpha+\alpha+\Lambda $ system. -
In order to test the reliability of EPD-ACCC method we choose the
$ 4_1^{+} $ state of$ {}^9_{\Lambda} $ Be as a proof-of-concept example. The single channel case and the coupled channel case are both considered. For single channel case, namely$ (l,L)=(4,0) $ , the generator coordinates are taken to be$ d=1,2,3,\cdots,25\ \rm{fm} $ and$ D=1,2,3,\cdots,20\ \rm{fm} $ . In the coupled channel calculations, the following restrictions on the channels and the generator coordinates are adopted as follows:$ \begin{aligned} &{channel}:\ J^{\pi}=4^+ \\& (l,L)=(4,0), (0,4), (2,2), (2,4), (4,2), (4,4)\\ &d=1,2,3,4,5,6,7,8\ \rm{fm},\ \ D=1,2,3,4,5,6,7,8\ \rm{fm}\\ \nonumber \end{aligned} $
We solve the coupled channel Hill–Wheeler equations and reproduce the so called
$ {}^8 $ Be analog states, where the configurations are$ {}^8 $ Be($ 0^+ $ ) + Λ,$ {}^8 $ Be($ 2^+ $ ) + Λ, and$ {}^8 $ Be($ 4^+ $ ) + Λ for the$ 0_1^{+} $ ,$ 2_1^{+} $ ,$ 4_1^{+} $ states of$ {}^9_{\Lambda} $ Be, respectively. The results are listed in Table 3. It is noticeable that the binding energy of Λ particle is defined as$ B_{\Lambda} $ =$ E(^9_\Lambda Be)-E(^8 Be(0^+_1)) $ . We can find that the Λ particle occupies the$ 0s $ orbit coupling to the$ {}^8 $ Be core with a high probability, which can be seen from the large s-wave components$ 0.9735 $ ,$ 0.9680 $ and$ 0.9274 $ for the$ 0_1^{+} $ ,$ 2_1^{+} $ ,$ 4_1^{+} $ states, respectively. In addition, we can find that the$ 0_1^{+} $ and$ 2_1^{+} $ in$ {}^8 $ Be analog states are both the bound states. However, the$ 4_1^{+} $ state calculated by the bound state approximation has a positive binding energy of Λ particle, which indicates that this state might be a resonance. More specifically, to clarify the nature of the excited states we should compare their energies with the lowest$ {}^{5}_{\Lambda} $ He + α threshold. The$ B_{\Lambda} $ values alone are insufficient to determine whether the excited states are bound or unbound. The$ 0^+_{1} $ and$ 2^+_{1} $ states are below the lowest$ {}^{5}_{\Lambda} $ He + α threshold, indicating that they are bound states. No excited states are found between the$ {}^{5}_{\Lambda} $ He + α and α + α + Λ thresholds, which suggests that potential$ ^{5}_\Lambda $ He + α resonance states are not within the scope of this work. Additionally, we refer to the resonance states found above the α + α + Λ threshold as "α + α + Λ resonances".$ J^{\pi} $ Energy (MeV) $ B_{\Lambda} $ (MeV)momenta $ w_{l,L}^2 $ $ 0_1^{+} $ −61.13(−60.45) −7.056 $ (l,L) $ (0,0)0.9735 (2,2)0.02596 (4,4)5 $ \times10^{-4} $ $ 2_1^{+} $ −58.18(−57.5) −4.106 $ (l,L) $ (2,0)0.9680 (0,2)0.01248 (2,2)0.01149 (2,4)3 $ \times10^{-4} $ (4,2)0.007596 (4,4)2 $ \times10^{-4} $ $ 4_1^{+} $ −50.62(−50.02) 3.454 $ (l,L) $ (4,0)0.9274 (0,4)7 $ \times10^{-4} $ (2,2)0.06105 (2,4)6 $ \times10^{-4} $ (4,2)0.009859 (4,4)3 $ \times10^{-4} $ Table 3. The weight
$ w_{l,L}^2 $ of a channel specified by the angular momenta of α-α part$ (l) $ and Λ-$ { }^8 $ Be part$ (L) $ , defined by Eq.(1). The total binding energies E and Λ binding energy$ B_{\Lambda} $ are also listed. It should be noted that$ B_{\Lambda} $ is defined as the binding energy measured from the calculated$ { }^8 \mathrm{Be}(0^+) $ energy$ -54.074 \;\mathrm{MeV} $ . Energies in parentheses are the results of the single channel calculation of$ (l,L)=(J, 0) $ . The$ \Lambda N $ interaction of YNG form is used.Therefore we utilize the EPD-ACCC method and the CSM to obtain the complex energy of
$ 4_1^{+} $ state. At first we use EPD-ACCC method to estimate the resonance of$ 4_1^{+} $ state with single channel calculations, the order of the Padé approximant is denoted by$ [N_1,N_2] $ .Choosing
$ N_1+N_2 $ bound state energies of$ {}^9_{\Lambda} $ Be close to$ \alpha+\alpha $ threshold, the coefficients in Padé rational function can be determined. The trajectories of resonance with coupling constant λ (here λ means the Majorana parameter M) with several Padé orders$ [N_1,N_2] $ are plotted in Fig. 2(a), where the values of coupling constant M corresponding to the points on the trajectories are taken to be 0.515, 0.520, 0.525, 0.530, 0.535, 0.540, 0.545, 0.550, 0.555, 0.560, 0.565, 0.570 and 0.573 from left to right. The rightmost point of the trajectory represents the real physical situation, namely the resonance we desire. Fig. 2(b) is similar with Fig. 2(a) but for coupled channel case. The resonant energies and decay widths obtained by CSM are also marked in the figures with red stars, which shows the good agreement between our two theoretical methods.Figure 2. (color online) Positions of the resonance poles with
$ 4_1^{+} $ state for the Padé approximant of different orders$ [N_1,N_2] $ (the approximant order is indicated in square brackets). (a) Single channel case ($ (l,L)=(4,0) $ ) (b) Coupled channel case. The values of coupling constant M corresponding to the points on the trajectories are taken as 0.515, 0.520, 0.525, 0.530, 0.535, 0.540, 0.545, 0.550, 0.555, 0.560, 0.565, 0.570 and 0.573 from left to right. The rightmost point of the trajectory is the real physical situation. The positions of the resonance poles obtained by CSM are also marked in the figures with red stars.In order to further corroborate the validity of this technique, we may artificially manipulate the M parameter to a higher value (namely broader resonance). Subsequently, employing the same Padé approximant, we can estimate the resonant state corresponding to the new parameter M. The CSM method is still applied as a benchmark during this process. For instance, we can set the value of parameter M to be larger and use the original bound states obtained in the coupled channel case, the complex energies of resonant states are listed in Table 5. Through the above calculations, we can see that for resonance states with broad decay widths, this method can also extrapolate reliable values.
M $ E_{r} $ (MeV)Γ (MeV) $ E^{CSM}_{r} $ (MeV)$ \Gamma^{CSM} $ (MeV)0.575 3.8592 0.4969 3.9311 0.5205 0.585 4.5075 0.7287 4.5345 0.7822 0.595 5.1320 1.0004 5.0928 1.0822 0.605 5.7350 1.3091 5.6057 1.4107 0.615 6.3182 1.6523 6.0734 1.7549 0.625 6.8831 2.0276 6.4966 2.0981 Table 5. Resonant energies and decay widths for the
$ 4_1^{+} $ state of$ {}^9_{\Lambda} $ Be with coupled channel calculations. Several different values of the parameter M are considered in order to test the stabilization and reliability of EPD-ACCC method. The results obtained by EPD-ACCC method are listed in the second and third columns. Here, the order of the Padé approximant is taken as$ [N_1,N_2] $ =[3,3]. In the fourth and fifth columns we also display the results extracted by CSM. -
The complex scaling is introduced by using the complex scaled generator coordinates
$ d_{\theta}=e^{i\theta}d $ and$ D_{\theta}=e^{i\theta}D $ in three cluster system. We calculate both the single channel and coupled channel cases for$ 4_1^{+} $ state of$ {}^9_{\Lambda} $ Be with CSM and extract the resonances from the complex energy spectra with different angle θ.At first, we show the results in the single channel case, which dominates the
$ 4_{1}^{+} $ state seen from Table 3. In Fig. 3(a) we display the complex energy spectrum, where the scaling angles are taken to be$ {12}^{\circ} $ ,$ {14}^{\circ} $ and$ {16}^{\circ} $ . We can find that a α + α + Λ resonance appears above the$ {}^8 $ Be$ (0^+)+\Lambda $ threshold. Using the conventional stability condition, the accurate value is extracted from the θ trajectory in Fig. 3(b) and the optimal scaling angle is also marked. With this optimal angle the location of the resonance is determined to be$ E-i\Gamma/2=4.0766-0.05096i $ MeV (energy of two free α particles already subtracted), which is a pretty narrow resonance. By comparing the real part of this complex energy$ -50.06 $ MeV with the positive energy$ -50.02 $ MeV obtain by bound approximation in Table 3, the good reliability of the bound state approximation for narrow resonance is also confirmed.Figure 3. (color online) The complex energy spectrum for
$ 4_1^+ $ state with GCM wave functions in the single channel case. The blue square, red circle and orange triangle correspond to$ \theta=12^{\circ},14^{\circ},16^{\circ} $ , respectively. Due to the small resonance energy of$ {}^8 $ Be($ 0^+ $ ) above the$ \alpha+\alpha $ threshold,$ {}^8 $ Be($ 0^+ $ )$ +\Lambda $ continuum and$ \alpha+\alpha+\Lambda $ continuum are indistinguishable, which can be seen from Fig. 3(a).$ {}^8 $ Be($ 4^+ $ )$ +\Lambda $ continuum is also marked in Fig. 3(a). In Fig. 3(b) we plot the θ trajectory of the resonant complex energy. The angle θ are taken from 4$ {}^{\circ} $ to 16$ {}^{\circ} $ with the step of 1$ {}^{\circ} $ . It can be seen from it that the optimal energy of the resonant state is about at$ \theta={9.7}^{\circ} $ .Similarly in Fig.4 (a)(b) we show the complex energy spectra of the
$ 4_1^{+} $ state calculated with coupled channel and the θ trajectory, respectively. The resonance marked by the red cycle is located at$ E-i\Gamma/2=3.8050-0.2368i $ MeV with θ stabilization condition. This broader resonance indicates that there is stronger repulsive effect caused by the coupled channel, or more explicitly mainly due to higher D wave Λ particle coupling to$ {}^8 $ Be$ (2^+) $ .Figure 4. (color online) Complex energy spectrum for
$ 4_1^+ $ state with GCM wave functions in the coupled channel case. The blue square, red circle and orange triangle correspond to$ \theta=20^{\circ},23^{\circ},26^{\circ} $ respectively. The angle θ are taken from 12.5$ {}^{\circ} $ to 15$ {}^{\circ} $ with the step of 0.25$ {}^{\circ} $ . It can be seen from Fig. 4(b) that the optimal energy of the resonant state is about at$ \theta={13.6}^{\circ} $ .The results of two theoretical method EPD-ACCC and CSM for
$ 4_1^{+} $ resonance are listed in Table 4. We can find that the resonances obtain by these two methods are very consistent with each other. In addition, through the investigation of$ 4_1^+ $ resonance in the range of$ M>0.573 $ as seen in Table 5, the good stabilization of EPD-ACCC method is further confirmed by the good consistency between two theoretical methods.$ J^{\pi}=4_1^{+} $ CSM EPD-ACCC $ E_r $ (MeV)Γ (MeV) $ E_r $ (MeV)Γ (MeV) single channel 4.0766 0.1020 4.2847 0.1068 coupled channel 3.8050 0.4736 3.7265 0.4556 Table 4. Resonant energies and decay widths for the
$ 4_1^{+} $ state of$ {}^9_{\Lambda} $ Be. The theoretical values are given by CSM and EPD-ACCC, respectively.Finally, we present the energy levels of
$ {}^8 $ Be and$ {}^9_{\Lambda} $ Be in Fig. 5, including both experimental data and theoretical results. The first column shows the experimental data for$ {}^8 $ Be. The shaded areas represent the decay widths, whose values are also indicated in parentheses. The second column contains our numerical results obtained using the microscopic R-matrix method, which are in good agreement with the experimental values. The third to seventh columns display the energy spectra for$ {}^9_{\Lambda} $ Be. The third column shows the energies of the bound states$ 0_1^{+} $ and$ 2_1^{+} $ , along with the resonance state$ 4_1^{+} $ obtained using the EPD-ACCC and CSM. The fourth column contains the$ {}^8 $ Be analog states from [43]. The fifth and sixth columns present the experimental data Exp.(1) and Exp.(2) taken from [79−81] and [82, 83], respectively. The final column shows the recalibrated experimental data Exp.(2)* [74] from Exp.(2). It is important to note that the resonant energies of$ {}^8 $ Be displayed in Fig. 5 differ from those in [17, 19], this is because of that in [17, 19] we did not use the parameter$ b=1.36 $ fm for the free α particle. Instead, we used$ b^{*}=1.3748 $ fm for free α particle to maximize its binding energy, which results in a difference of about 0.0316 MeV in the resonant energies.Figure 5. (color online) Low lying energy levels of
$ {}^8 $ Be and$ {}^9_{\Lambda} $ Be. Here, we take$ \alpha+\alpha+\Lambda $ threshold as the zero point of energy. The black solid lines represent the binding energies of bound states and resonant energies of resonances, and the shaded areas indicate the corresponding decay widths. The first column corresponds to experimental data of$ {}^8 $ Be [78], the second column presents the energy levels of$ 0^{+} $ ,$ 2^{+} $ , and$ 4^{+} $ states calculated based on the microscopic R-matrix method. The third column represents the$ 0_1^{+} $ ,$ 2_1^{+} $ , and$ 4_1^{+} $ states of$ {}^9_{\Lambda} $ Be obtained through the EPD-ACCC calculations, while the fourth column displays the theoretical results taken from [43]. The observed energies of$ {}^9_{\Lambda} $ Be in column Exp.(1) are taken from [79−81] and the observed energies in column Exp.(2) are taken from [82, 83]. Additionally, in column Exp.(2)*, we illustrate the recalibration of the observed values in column Exp.(2) according to [74].Our results for the
$ 0_1^{+} $ and$ 2_1^{+} $ bound states show consistency with those obtained using the orthogonality condition model and Gaussian expansion method in [43], and they also align well with the experimental data. Additionally, our results for the$ 4_1^{+} $ resonant state are also in good agreement with those from [43]. The consistent resonances obtained from EPD-ACCC and CSM validate the good reliability of the EPD-ACCC method. -
The matrix elements required are listed as follows:
The
$ {}^5_{\Lambda} $ He caseMatrix elements of overlap:
$ \left\langle{\phi_{\alpha}\phi_J(d_1,r)Y_J(\hat{{r}})}|{\phi_{\alpha}\phi_J(d_2,r)Y_J(\hat{{r}})}\right\rangle =4\pi\exp\left(-\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{4}\right)\mathscr{J}_J\left(\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{2}\right). $
(A1) Matrix elements of kinetic:
$ \begin{aligned}[b] \left\langle{\phi_{\alpha}\phi_J(d_1,r)Y_J(\hat{{r}})}{T}|{\phi_{\alpha}\phi_J(d_2,r)Y_J(\hat{{r}})}\right\rangle =\;&4\pi\hbar\omega\exp\left(-\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{4}\right)\bigg[\left(\dfrac{3}{4}-\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{8}\right)\mathscr{J}_J\left(\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{2}\right)\\&+\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{4}\mathscr{J}^{\prime}_J\left(\dfrac{(d_1^2+d_2^2)\nu_{\Lambda,\alpha}}{2}\right)\bigg]. \end{aligned} $
(A2) Matrix elements of
$ \Lambda N $ potential:$ \left\langle{\phi_{\alpha}\phi_J(d_1,r)Y_J(\hat{{r}})}{\sum_{i=1}^{4}V_{\Lambda N}(i)}|{\phi_{\alpha}\phi_J(d_2,r)Y_J(\hat{{r}})}\right\rangle =4\left(1+\dfrac{\mu}{\nu_{\Lambda N}}\right)^{-3/2}4\pi\exp(-\gamma(d_1^2+d_2^2))(V_{\Lambda N}^D\mathscr{J}_J(\gamma^+d_1d_2)+V_{\Lambda N}^E\mathscr{J}_J(\gamma^-d_1d_2)). $
(A3) The notations used above are defined by:
$ \gamma=\dfrac{1}{4}\left(\nu_{\Lambda \alpha}+\dfrac{\mu \nu_{\Lambda N}}{\mu+\nu_{\Lambda N}}\right),\;\ \gamma^{\pm}=\dfrac{1}{2}\left(\nu_{\Lambda \alpha}+\dfrac{(-1\pm1)\nu_{\Lambda N}^2}{\mu+\nu_{\Lambda N}}-\dfrac{\mu \nu_{\Lambda N}}{\mu+\nu_{\Lambda N}}\right) $
(A4) with
$ \mu=\dfrac{1}{\beta_{\Lambda N}^2}, \ \nu_N=\dfrac{1}{b_N^2}, \ \nu_{\Lambda}=\dfrac{1}{b_{\Lambda}^2}, \ \nu_{\Lambda N}=\dfrac{\nu_{\Lambda}\nu_{ N}}{\nu_{\Lambda}+\nu_{N}},\ \nu_{\Lambda \alpha}=\dfrac{4\nu_{\Lambda}\nu_{N}}{\nu_{\Lambda}+4\nu_{N}},\ \nu_{\Lambda ,\alpha\alpha}=\dfrac{8\nu_{\Lambda}\nu_{N}}{\nu_{\Lambda}+8\nu_{N}}. $
(A5) The wave functions of two α clusters are constructed by the Brink wave functions. The matrix elements involving two α clusters are given as follows at first:
Kernels of
$\alpha+\alpha$ Norm kernel:$ \left\langle{\Phi_{l,\alpha\alpha}(d_1)}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle =(1+(-1)^l)4\pi\exp\left(-\dfrac{(d_1^2+d_2^2)\nu_N}{2}\right)\left(\mathscr{J}_l(d_1d_2\nu_N)-4\mathscr{J}_l\left(\dfrac{d_1d_2\nu_N}{2}\right)+3\delta_{l0}\right) $
(A6) Kinetic kenrel:
$ \left\langle{\Phi_{l,\alpha\alpha}(d_1)}{T-T_G}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle =\hbar\omega\left(\dfrac{3}{4}(A-1)+\dfrac{\nu_N}{2}\dfrac{d}{d\nu_N}\right)\left\langle{\Phi_{l,\alpha\alpha}(d_1)}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle $
(A7) where A is the nucleon number of
$ {}^8 $ Be.$ NN $ interaction kernel:$ \begin{aligned}[b] \left\langle{\Phi_{l,\alpha\alpha}(d_1)}{V_{NN}}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle =\;&(1+(-1)^l)4\pi\exp(-(d_1^2+d_2^2)\dfrac{\nu_N}{2})\times\{(2V_d+2V_e)[\mathscr{J}_{l}(d_1d_2\nu_N)\\&-2\mathscr{J}_{l}(d_1d_2\dfrac{\nu_N}{2})+\delta_{l0}-2[\exp(-qd_1^2\dfrac{\nu_N}{2})+\exp(-qd_2^{2}\dfrac{\nu_N}{2})](\mathscr{J}_{l}(d_1d_2\dfrac{\nu_N}{2})-\delta_{l0})]\\&+2V_d\exp[-q(d_1^2+d_2^{ 2})\dfrac{\nu_N}{2}][\mathscr{J}_{l}((2-2q)d_1d_2\dfrac{\nu_N}{2})-2\mathscr{J}_{l}((1-2q)d_1d_2\dfrac{\nu_N}{2})+\mathscr{J}_{l}(2qd_1d_2\dfrac{\nu_N}{2})]\\&+2V_e\exp[-q(d_1^2+d_2^{ 2})\dfrac{\nu_N}{2}][\mathscr{J}_{l}((1-2q)d_1d_2\dfrac{\nu_N}{2})-2\mathscr{J}_{l}(2qd_1d_2\dfrac{\nu_N}{2})+\mathscr{J}_{l}((1+2q)d_1d_2\dfrac{\nu_N}{2})]\} \end{aligned} $
(A8) with
$ \begin{aligned}[b] q=\dfrac{\mu_N}{2(\nu_N+2\mu_N)},\qquad \mu_N=\dfrac{1}{\beta_{NN}^2}\quad V_d=v_{NN}^{0}(\dfrac{2q\nu_N}{\mu_N})^{3/2}(8W+4B-4H-2M), V_e=v_{NN}^{0}(\dfrac{2q\nu_N}{\mu_N})^{3/2}(8M+4H-4B-2W) \end{aligned} $
(A9) The
$ NN $ interaction used here is of the form$ V_N(r)=v_{NN}^{0}(W-MP_{\sigma\tau}+BP_{\sigma}-HP_{\tau})\exp(-\dfrac{r^2}{\beta_{NN}^2}) $
(A10) The
$ {}^9_{\Lambda} $ Be caseMatrix elements of overlap:
$ \begin{aligned}[b] &\left\langle{[\Phi_{l,\alpha\alpha}(d_1)\phi_{L}(D_1,R)Y_{L}(\hat{{R}})]_J}|{[\Phi_{l,\alpha\alpha}(d_2)\phi_{L}(D_2,R)Y_{L}(\hat{{R}})]_J}\right\rangle \\=\;& \left\langle{\Phi_{l,\alpha\alpha}(d_1)}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle 4\pi\exp(-\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{4})\mathscr{J}_L(\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{2}) \end{aligned} $
(A11) Matrix elements of kinetic
$ T_R $ :$ \begin{aligned}[b] &\left\langle{[\Phi_{l,\alpha\alpha}(d_1)\phi_{L}(D_1,R)Y_{L}(\hat{{R}})]_J}T_R|{[\Phi_{l,\alpha\alpha}(d_2)\phi_{L}(D_2,R)Y_{L}(\hat{{R}})]_J}\right\rangle \\=\;& \left\langle{\Phi_{l,\alpha\alpha}(d_1)}|{\Phi_{l,\alpha\alpha}(d_2)}\right\rangle 4\pi\hbar\omega\exp(-\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{4})[(\dfrac{3}{4}-\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{8})\mathscr{J}_L(\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{2})\\&+\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{4}\mathscr{J}^{\prime}_L(\dfrac{(D_1^2+D_2^2)\nu_{\Lambda,\alpha\alpha}}{2})] \end{aligned} $
(A12) Matrix elements of
$ \Lambda N $ potential:$ \begin{aligned}[b] &\left\langle{[\Phi_{l_1,\alpha\alpha}(d_1)\phi_{L_1}(D_1,R)Y_{L_1}(\hat{{R}})]_J}{\sum_{i=1}^{8}V_{\Lambda N}(i)}|{[\Phi_{l_2,\alpha\alpha}(d_2)\phi_{L_2}(D_2,R)Y_{L_2}(\hat{{R}})]_J}\right\rangle \\=\;&16(1+\dfrac{\mu}{\nu_{\Lambda N}})^{-3/2}(4\pi)^2\exp[-\alpha_1(D_1^2+D_2^2)-\alpha_2(d_1^2+d_2^2)](-1)^{l_1+L_1}\sum_{m=0}^3 \begin{pmatrix} 3\\m \end{pmatrix}(-1)^{m+1}\sum_{k_1,\cdots,k_6}(-1)^{k_1+\cdots+k_6}\hat{k}_1\cdots\hat{k}_6\\&\times[V_{\Lambda N}^D\mathscr{J}_{k_1}(\alpha_4^+(m)d_1d_2)\mathscr{J}_{k_2}(\alpha_3^+D_1D_2)\mathscr{J}_{k_3}(\alpha_6^+d_1D_2)\mathscr{J}_{k_4}(\alpha_6^+d_2D_1)\mathscr{J}_{k_5}(\alpha_5 d_1D_1)\mathscr{J}_{k_6}(\alpha_5 d_2D_2)\\&+V_{\Lambda N}^E\mathscr{J}_{k_1}(\alpha_4^-(m)d_1d_2)\mathscr{J}_{k_2}(\alpha_3^-D_1D_2)\mathscr{J}_{k_3}(\alpha_6^-d_1D_2)\mathscr{J}_{k_4}(\alpha_6^-d_2D_1)\mathscr{J}_{k_5}(\alpha_5 d_1D_1)\mathscr{J}_{k_6}(\alpha_5 d_2D_2)]\\&\times\sum_{K_1, \cdots, K_4}C_1(K_1 K_2 K_3 K_4,k_5k_6,l_1L_1l_2L_2J)C_2(K_1K_2K_3K_4,k_1k_2k_3k_4,J) \end{aligned} $
(A13) where
$ C_1(K_1 K_2 K_3 K_4,k_5k_6,l_1L_1l_2L_2J)=W(K_1K_2l_1L_1;Jk_5)(k_5K_1|l_1)(k_5K_2|L_1)W(K_3K_4l_2L_2;Jk_6)(k_6K_3|l_2)(k_6K_4|L_2) $
(A14) and
$ C_2(K_1K_2K_3K_4,k_1k_2k_3k_4,J)=\sqrt{\hat{K}_1\hat{K}_2\hat{K}_3\hat{K}_4}(k_1k_3|K_1)(k_4k_2|K_2)(k_1k_4|K_3)(k_3k_2|K_4) \begin{Bmatrix} k_1 & k_3 & K_1\\k_4 & k_2 & K_2\\K_3&K_4&J \end{Bmatrix} $
(A15) with
$ (ab|c)=(a0b0|c0) $ and$ \hat{d}=2d+1 $ .The notations used above are defined by:
$ \begin{aligned}[b] &\alpha_1=\dfrac{1}{4}(\nu_{\Lambda,\alpha\alpha}+\dfrac{\mu\nu_{\Lambda N}}{\mu+\nu_{\Lambda N}}) \qquad \alpha_2=\dfrac{1}{16}(8\nu_N+\dfrac{\mu\nu_{\Lambda N}}{\mu+\nu_{\Lambda N}}) \qquad \alpha_3^{\pm}=\dfrac{1}{2}(\nu_{\Lambda,\alpha\alpha}-\nu_{\Lambda N}\pm\dfrac{\nu_{\Lambda N}^2}{\mu+\nu_{\Lambda N}})\\&\alpha_4^{\pm}=\dfrac{1}{8}(\nu_{N}-\nu_{\Lambda N}\pm\dfrac{\nu_{\Lambda N}^2}{\mu+\nu_{\Lambda N}})\qquad \alpha_4^{\pm}(m)=\alpha_4^{\pm}+\dfrac{2m-3}{4}\nu_N\\&\alpha_5=\dfrac{1}{4}(-\nu_{\Lambda N}+\dfrac{\nu_{\Lambda N}^2}{\mu+\nu_{\Lambda N}})\qquad \alpha_6^{\pm}=\dfrac{1}{4}(-\nu_{\Lambda N}\pm\dfrac{\nu_{\Lambda N}^2}{\mu+\nu_{\Lambda N}}) \end{aligned} $
(A16)
Analytic continuation in the coupling constant for resonances in $ {}^9_{\Lambda} $ Be
- Received Date: 2024-06-01
- Available Online: 2025-01-01
Abstract: The application scope of the analytic continuation in the coupling constant (ACCC) can be extended to the exchange parameters of the effective nucleon-nucleon interaction in microscopic cluster model. Based on such exchange parameter dependent ACCC (abbreviated as EPD-ACCC), we study the