Decomposition of scattering phase shifts and reaction cross sections using the complex scaling method

Decomposition of scattering phase shifts and reaction cross sections using the complex scaling method

Myagmarjav ODSUREN Meme Media Laboratory, Hokkaido University, Sapporo 060-8628, Japan Nuclear Research Center, National University of Mongolia, Ulaanbaatar 210646, Mongolia    Kiyoshi KATŌ Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan    Masayuki AIKAWA Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan    Takayuki MYO General Education, Faculty of Engineering, Osaka Institute of Technology, Osaka 535-8585, Japan Research Centre for Nuclear Physics (RCNP), Osaka University, Ibaraki 567-0047, Japan
December 3, 2018

We apply the complex scaling method to the calculation of scattering phase shifts and extract the contributions of resonances in a phase shift and a cross section. The decomposition of the phase shift is shown to be useful to understand the roles of resonant and non-resonant continuum states. As examples, we apply this method to several two-body systems: (i) a schematic model with the Gyarmati potential which produces many resonances, (ii) the system which has a Coulomb barrier potential in addition to an attractive nuclear interaction, and (iii) the system which has no barrier potential. Using different kinds of potentials, we discuss the reliability of this method to investigate the resonance structure in the phase shifts and cross sections.

PACS numbers

21.60.Gx, 24.10.-i

Valid PACS appear here
preprint: APS/123-QED

I Introduction

Nuclear scattering is the most important phenomena from which we can obtain information and knowledge on various nuclear properties. Many theoretical approaches and experimental techniques have been developed to extract physics from the scattering phenomena. Resonances observed as some peaks in scattering cross sections provide us with the precious information for understanding nuclear interactions and structures. In particular, it is indispensable in recent developments of unstable nuclear physics to investigate the resonances involved in the unbound states locating above the many-particle decay thresholds, because the unstable nuclei barely have bound states and most of the excited states are resonances. Furthermore, to understand weakly bound states in unstable nuclei, such as a halo structure, we need to investigate continuum states together with bound states because of a strong coupling between them Ho12 ().

Study of resonances in the scattering problem of light nuclei has been carried out by using various methods, one of which is the complex scaling method (CSM) Ag71 (). The theory of the complex scaling was proposed mathematically Ag71 () and it has been extensively applied to the atomic and nuclear physics Ho83 (); Mo98 (); Ao06 (). In the CSM, resonant states of the many-body systems are described using the appropriate basis functions. The resonance wave functions are obtained as eigenstates together with bound states by carrying out the diagonalization of the complex scaled Hamiltonian. Hence, we can calculate expectation values of physical quantities for resonances Ho97 (). Recently, resonances decaying up to five-body systems have been studied in the CSM and successfully compared with experiments My11 (); My12 ().

We have also discussed that the CSM is very useful in studies of not only resonant states but also scattering states Ao06 (). The complex scaling separates the resonant states with resonance energies and widths () from the continuum states obtained on a ’-line’, where is a scaling angle in the CSM Ag71 (). We can extract scattering properties from the continuum solutions on the ’-line’ together with the resonance solutions. Calculations of the phase shifts of continuum states with complex energies along the ’-line’ provide the background phase shifts exclusive of the resonances existing in the energy region between the ’-line’ and the real energy axis Ho97 (). As examples of applications of the CSM to scattering problems, response functions for external electric fields and breakup cross sections of two-neutron halo nuclei have been discussed and shown for observed data to be well described by a two-neutron-plus-core three-body model My98 (); My07 (); My01 (). In their results, it is concluded that a sequential breakup process is dominant rather than a direct three-body breakup.

In these calculations, a complex scaled Green’s function has been introduced to obtain the response functions as observables. We can also calculate scattering quantities such as phase shifts in a form of sum of resonances and background terms. It is possible to investigate the resonance contributions and to obtain a deep understanding of resonance structure by separation of a scattering quantity. Suzuki et al. Su05 () showed that scattering phase shifts can be calculated from the continuum level density (CLD) which is expressed using the complex scaled Green’s function. The CLD is also called a time delay and has a relation with the -matrix Le69 (); Sh92 (). In the CSM, we obtain the discretized continuum states as complex energy eigenstates. Using those complex eigen-energies, we can evaluate the CLD as a smoothed real-energy function.

In this paper, we discuss the explicit relation between the scattering phase shifts and complex-energy eigenvalues in the CSM via the CLD. The results provide us with deeper understanding of the role of resonant states characterised by the widths described as an imaginary part of the eigen-energy. We show the results of typical potential scattering which has many resonances near the real energy. We also analyze the several realistic systems and compare the results with the observed data. The observed scattering cross sections have various energy distributions as a result of interference between resonance and background terms of phase shifts Fa61 (). We see this interference in the and systems.

The outline of the paper is as follows. In Section II, the method is briefly explained, and the scattering phase shifts are explicitly shown by using the complex energy solutions in the CSM. In Section III, the present method is applied to several two-body systems; a typical potential problem, so called Gyarmati potential Gy90 () is discussed. The phase shifts of and systems are presented together with experimental ones. Finally, in Section IV, discussion and summary are given.

Ii Decomposition of phase shifts

ii.1 Complex scaling method

We take up two-body systems which are described by the Schrödinger equation


where the Hamiltonian consists of kinetic energy and potential for the relative motion between two bodies. The eigenvalue problem is generally solved under a boundary condition of asymptotic out-going waves for bound states and resonances. The out-going boundary condition directly enables us to solve bound states in an functional basis set because the states have negative energies and a damping behavior in the asymptotic region. Resonant states are unbound and defined as the eigenstates belonging to the complex eigen-energy which corresponds to a complex momentum value in the lower half plane (unphysical plane Ao06 ()). The resonant states cannot be solved in the functional space due to asymptotically divergent behavior. Furthermore, continuum states of arbitrary positive energies cannot also be obtained under the out-going condition.

The complex scaling has been introduced to solve resonant states within basis functions, and is defined by the following complex-dilatation transformation for relative coordinate and momentum Ag71 ();


where is a scaling angle and . The maximum value is determined so as to keep analyticity of the potential. For example, for a Gaussian potential. This transformation makes every branch cut rotated by on the complex energy plane. In the wedge region pinched by the rotated branch cut and the positive energy axis, resonance eigenstates are obtained by solving the following eigenvalue problem


within an appropriate non-orthogonal basis set . The index is to distinguish the eigenstates of the complex scaled Hamiltonian . The bound states are obtained on the negative energy axis independently from as well as the ordinary bound states. Because of a finite number of basis states, the continuum states are discretized with complex energies distributed on the rotated branch cut (-line).

The eigenvalues and eigenstates of the complex scaled Schrödinger equation (3) are classified as


where and are the number of bound states and the number of resonant states which depends on , respectively. The complex energies of resonant states are obtained as , when . The discretized energies of continuum states are -dependent and expressed as .

These three-kinds solutions of the complex-scaled Schrödinger equation construct the extended completeness relation My98 ();

where the tilde ( ) in bra-states means the bi-orthogonal states with respect to the ket-states due to non-hermitian property of . The integration of the third term is taken along the rotated branch cut . The extended completeness relation has been proven for single- and coupled-channel systems Gi03 (); Gi04 (). In the case of eigenstates within a finite number of basis states, the integration for continuum states is approximated by the summation of discretized states as My98 ()


It has been investigated that the reliability of the approximation of the continuum states are confirmed by using a sufficient large basis number of in the CSM Ao06 (); Su05 ().

ii.2 Continuum level density and scattering phase shift

The CLD, as function of the real energy is defined as Le69 (); Sh92 ()



are the full and free Green’s functions, respectively. The CLD is also related to the scattering phase shift and the relation is expressed in a single channel case as Le69 (); Sh92 ():


Using Eqs. (7) and (8), we can obtain the phase shift in terms of the eigenvalues of and by integrating the CLD over the energy . When we apply the complex scaling and obtain the complex scaled Green’s function, the CLD can be expressed as


Furthermore, we apply the extended completeness relation given in Eq. (6) to the calculation of in Eq. (9), where we use the solutions obtained by diagonalization of matrix elements of and with a finite number of basis functions. The energy eigenvalues of are given as with . All of these values are distributed on the ’-line’. The CLD of Eq. (9) is approximated as

where . It is important to note that the approximated CLD, basically has a dependence on the scaling angle because we employ a finite number of the basis states to expand the complex scaled wave functions. In the calculation we adopt a sufficiently large number of to keep the numerical accuracy and to make the -dependence negligible in the solutions Su05 (). Thus we calculate the phase shift from ;

Figure 1: The geometrical expressions for phase shifts: , and as functions of the energy . Both and cases are displayed in the upper panel a) and lower panel b), respectively. The details are explained in the text.

Performing integration of every term, we obtain the following expression:

where we assume . When we define , and as


respectively, we express the phase shift as


The geometrical expressions of , and are given for or in Fig. 1. The resonance phase shift is the angle of the th resonant pole at an energy on the real energy axis. At the resonance energy, , the relation is confirmed for every resonant pole. In addition, at and at for each resonance. Similarly the phase shifts from continuum terms including asymptotic part, are given by the angles of the discretized continuum energies. At , the continuum terms of the phase shifts go to because of the relation . Thus, for and the Levinson theorem is confirmed as


where the number of the Pauli forbidden states is included in when they exist.

The cross section is described by using these phase shifts, and we can see the contributions from both resonant poles and continuum terms. When we concentrate our interest on the contribution from a single resonant pole, other terms are described as a background phase shift. We can have the similar discussion as was done by Fano Fa61 () because the partial cross section for the orbital angular momentum is given as


where with the reduced mass . The phase shift is expressed in the form as , where and are the single resonance and the background terms including all other terms given in Eq. (14), respectively. The shape of the cross section can be investigated by evaluating the resonance () and background () phase shifts.

Iii Applications to several systems

iii.1 Typical potential scattering

We apply the method of analyzing the phase shifts to a simple schematic potential which is introduced in Ref.  Gy90 (). The explicit form of the Hamiltonian is given as


where the so-called Gyarmati potential is


It is assumed that (MeV).

This potential has an attractive pocket in the inside and a repulsive barrier in the outside. In this system, one bound state and several resonant poles for and are obtained Ho97 (); My97 (). It is interesting to see the contributions from those resonant poles to the scattering quantities. In Ref.  Ho97 (), it has been shown that -transition strengths into the resonant states exhaust the sum rule value. In Ref. My97 (), continuum states on the ’-line’ of the complex energies have been investigated through the scattering phase shifts on the rotated branch cut. The calculated phase shifts for different values suggest that the resonant states located above the -line behave like bound states. The phase shifts approach to in the higher energy region where is the number of those resonant states as will be discussed in Section IV.

To solve the eigenvalue problem of Eq. (3), we employ the Gaussian basis functions Hi03 () given as


where the range parameters are given by a geometric progression as , and the normalization factor is given as


We take and employ the optimal values of and so as to obtain stationary resonance solutions. The results of eigenvalues for bound and resonant states obtained with are shown in Table I. The results are same as those in Refs. Ho97 (); My98 (). In addition, the consistent and stable results were obtained for the CLD at different scaling angles ( and ) in Fig. 4 of Ref. Su05 (). Using those eigenvalues including the continuum states, we calculate phase shifts by Eq. (14).

wave wave
E (MeV) state E (MeV) state
-1.928 bound -0.675 bound
0.310- resonance 1.171- resonance
1.632-0.123 resonance 2.031- resonance
2.249-1.040 resonance 2.832- resonance
2.854-2.570 resonance 3.934- resonance
Table 1: Bound and resonance energies with decay widths calculated with for the and partial waves.
Figure 2: The resonance and continuum phase shifts and the sum for .
Figure 3: The phase shifts of the Gyarmati potential for and the subtraction of the resonance terms one by one.
Figure 4: Upper panel: The decomposition of scattering phase shifts of the system for (a) , (b) and (c) , respectively. The dashed-lines and dotted-lines represent the contributions of the resonance and continuum terms, respectively. The solid-lines display total scattering phase shifts. The experimental data He10 () are shown with open circles. Lower panel: The distributions of eigenvalues are given in the complex energy plane for each partial wave.
Figure 5: Upper panel: The decomposition of the scattering phase shifts of system for (a) and (b) . The dashed-lines and dotted-lines represent the contributions of the resonance and continuum terms, respectively. The solid-lines display total scattering phase shifts. The experimental data Ho66 () are shown with open circles. Lower panel: The distributions of eigenvalues are displayed in the complex energy plane.
Figure 6: Total cross sections as functions of relative energy of the system. The dotted-line shows the present calculation and open circles are experimental data  Au62 ().
Figure 7: The partial cross sections of the (a) and (b) waves. The dashed, dotted and solid-lines denote the results of resonance, continuum contributions and partial cross sections, respectively.

In Figs. 2 and 3, we show the calculated phase shifts of the state. Here, we define resonance and continuum phase shifts, and as


respectively, where , and are given in Eq. (13). The phase shifts, and have the -dependence coming from the number of resonant poles above the ’-line’. Both of two-kinds of phase shifts and do not show zero in the negative energy region. On the other hand the sum of and , , becomes -independent and zero in the negative energy region due to the cancelation of each other. Therefore, the total phase shift, , is independent of the scaling angle and finite only in the positive energy () expect for . Using this property, we can rewrite the total phase shifts given by integration from as


From the resonance term , we can distinguish the 1st and 2nd resonances of the state in Table 1 at the corresponding resonance energies (0.31 MeV and 1.63 MeV), but the structure of higher resonances are not clearly seen. The continuum term has no structure and always negative values indicating a repulsive interaction nature. Its behavior looks like of a hard sphere scattering. In fact, we can estimate the hard sphere radius as fm from the behavior of in Fig. 2.

In order to see the resonance effect on the phase shift clearly, we calculate the phase shift subtracting the resonance term from as


In Fig. 3, the results of calculation are shown for , where . It is found that the effects of the 1st and 2nd resonances are remarkable, but the 3rd and 4th resonances do not have notable effects. It is shown that every resonance definitely changes the phase shift by at higher energies.

The phase shifts of the state are also calculated and the results show the similar behavior as the case for one bound and four resonances shown in Table I.

iii.2 Scattering phase shifts in the and systems

Figure 8: The partial cross sections in (a) and (b) waves of the system. The dashed, dotted and solid-lines denote the results of resonance, continuum contributions and partial cross sections, respectively.
Figure 9: The real part of complex phase shift on the ’-line’ for the wave at different scaling angles. Energy is given by the absolute value measured from the threshold.

We apply the method to the realistic problems of the and systems comparing with the observed data. The system including a long range Coulomb interaction is described by the following orthogonality condition model (OCM) Sa69 (). Hamiltonians are given as:


where we use MeVfm, MeV, fm and fm. The term is a Pauli-potential Ku84 () to project out the Pauli forbidden states from the relative motion in the OCM. The explicit form of is given by using the harmonic oscillator wave functions as


where we use the range parameter fm and the potential strength MeV. It is noted that the Coulomb potentials between two particles are expressed by folding potentials with charge distributions of the Gaussian and point in and , respectively.

The eigenvalue problems for these Hamiltonians are solved on the basis of the Gaussian basis function method Hi03 () mentioned above. The resonance solutions for and are obtained as in MeV comparing with the experimental data Ti04 () in MeV, respectively. Using continuum solutions in addition to these resonance ones as shown in lower panels of Fig. 4, the CLD and the scattering phase shifts are calculated. Here, to see the resonance and continuum phase shifts for positive energies, we redefined them in positive energies;

It is noticed that for . The results of and are shown in upper panels of Fig. 4.

From lower panels of Fig. 4, it is seen that one resonant pole in every partial wave of and is obtained. There are no other resonant poles which can make any structure in the resonance phase shifts as seen in upper panels of Fig. 4. In the case of , the resonance term shows a sharp resonance behavior because of the very small resonance width, and the continuum term shows a rather strong repulsive behavior. The case shows that the resonance behavior is weakened by the large continuum contribution repulsively. On the other hand, in the case, the resonance behavior remains due to small contribution of the continuum term which shows a repulsive nature like other partial waves. This repulsive nature of the continuum terms for and can be understood in association with presence of the Pauli-forbidden states. Because the forbidden states act as bound states to be orthogonal to the scattering solutions, they cause a repulsive nature of the continuum phase shifts.

The system is described by using the Hamiltonians


where for the potential we use the so-called microscopic KKNN potential Ka79 () given by

The parameters are given in Table 2. The Pauli-potential is defined with the harmonic oscillator wave function as


where the oscillator size parameter is taken to be fm and MeV.

(MeV) (fm) (MeV) (fm)
-96.3 0.36 34.0 0.20
central 77.0 0.90 -85.0 0.53
51.0 2.50
(MeV) (fm) (MeV) (fm)
spin-orbit -16.8 0.52 -20.0 0.396
20.0 2.200
Table 2: Parameters of the KKNN potential Ka79 ().

Using the Gaussian basis functions, we solve the complex scaled eigenvalue problems for the Hamiltonians of Eq. (27) with and as well. The results for the and waves are presented in Fig. 5. One resonant pole of the system is obtained: MeV for and MeV for , which are compared with the experimental data MeV for and MeV for Ti02 (). Using these results and Eq. (LABEL:eq3-2-3), we calculate the resonance, continuum and total phase shifts, which are shown in upper panels of Fig. 5 together with experimental data. We can see a good agreement between theoretical and experimental results. The resonance phase shift of increases rapidly due to the small decay width. Although has a larger width, the phase shift of shows a clear resonance behavior beyond . The continuum phase shifts of both states are very similar. This trend seems due to the same wave scattering and a small effect of the force to the background states.

The property of the scattering phase shifts is determined from a sum of resonance and continuum terms. Therefore, the observed resonances depend on not only resonant states as poles but also the contribution from the continuum state. The and systems show to keep resonance behavior in spite of existence of continuum contributions. On the other hand, the resonant poles higher than the 3rd one in the Gyarmati potential can not be distinguished as resonances in the phase shifts. They are absorbed in the continuum states.

iii.3 Scattering cross sections of the and systems

The partial cross sections are given in Eq. (16), and the total cross section is expressed as


In Fig. 6, we show the comparison of our computed total cross sections with experiments of the system. The experimental data are taken from Ref. Au62 (). It is found that the calculated total cross sections are in good agreement with the experimental data in a wide energy region. A very sharp peak is observed at the low energy around 2 MeV and has a long tail distribution in higher energies. The low energy cross section dominantly comes from partial waves as seen below. The -wave has no resonance due to no barrier. A virtual state corresponding to a pole on the negative imaginary momentum axis can appear in -waves when the interaction between and neutron has an appropriate strength. However, the KKNN potential is not so strong as to produce such a virtual state Ma00 () and has a repulsive nature for -waves.

Figure 7 shows the partial cross sections of and waves which are decomposed into contributions of resonance and continuum terms, respectively. The partial cross sections are calculated using Eq. (16). Resonance and continuum cross sections are calculated using resonance and continuum phase shifts of and , respectively, in Eq. (LABEL:eq3-2-3) and are presented in Fig. 5. We see that the peak of the total cross section corresponds to the sharp resonance peak in and the resonance cross section of gives rather broad distribution. The continuum cross sections show the similar behavior in the and waves, which are presented by doted-lines in Fig. 7.

In Fig. 8, partial cross sections and their decomposition into resonance and continuum terms are shown for and waves of the system. The partial cross section for is too sharp as like the delta-function because the small decay width of the resonance. Then, we skipped the partial cross section. For and , resonance cross sections have the shapes like the Breit-Wigner form. The continuum contribution of is rather large, while that is small in . The partial cross section of is not so different from the resonance cross section as compared to the case. It is interesting that the peak energies of the partial cross section fairly shift from the position of the resonance energies.

Iv Discussion and summary

The advantage of the CSM is to decompose the unbound states into resonant and non-resonant continuum states by rotating the branch cut on the complex energy plane with a parameter . As a result, we can separate a physical quantity into two parts associated with resonant and non-resonant continuum states. To investigate directly the properties of the continuum states on the ’-line’ in the CSM, the phase shifts on the ’-line’ were calculated for the Gyarmati potential My97 (). The phase shifts on the ’-line’ are obtained as complex values and the real parts are presented for various values in Fig. 9. For , the calculated phase shift is very similar to the in Figs. 2 and 3. In Fig. 9, the values are taken for resonant poles so as to be separated one by one from the continuum states. The phase shifts show a jump of at the asymptotic energies due to existence of a resonance between neighboring values. This behavior is completely the same as shown in Fig. 3. From this result, we can understand the physical meaning of the ’-line’ and the solutions obtained along this line in the CSM. In a practical application, we can analyse the scattering quantities decomposing them into sharp resonance term and background term like the relation discussed above.

The decomposition of a cross section is not so simple, because it is not equal to a direct sum of resonance and continuum cross sections due to their interference. Since the study by Fano Fa61 (), many discussions have been done so far. Inserting into Eq. (16), we have the so-called Fano formula Sh12 ()


where and . In many discussions, and defined in Eq. (13) have been used instead of and . When , they almost coincide with each other. Using the Fano formula, we can understand how the shape of the cross section deviates from the Breit-Wigner form.

In this work we discussed scattering phase shifts and cross sections in the framework of the CSM to investigate the unbound states of two-cluster systems. Using complex energy eigenvalues of resonant and non-resonant continuum states, the analytic form of the phase shifts is derived in a form of sum of resonance and continuum terms in addition to the constant term coming from bound states. This decomposition of phase shifts are very useful to see the resonance contributions in the observed phase shifts and cross sections.

The framework was applied to a simple schematic potential, the so-called Gyarmati potential Gy90 (), producing many resonances for and states. The results indicate that resonances embedded in continuum energies are exposed by using the CSM, and their contributions to the phase shift are made distinctive. This means that the present approach is very promising to understand the role of resonances and their structures in the scattering observables. It is confirmed that the present results of calculation of the phase shift and extracting the resonance terms show a good correspondence with the previous calculations in Ref. My97 ().

Applying the present framework to the and systems, we obtained the good reproduction of the observed phase shifts and cross sections. The decomposition into resonance an continuum terms makes clear that resonance contributions are dominant but continuum terms and their interference are not negligible. To understand the behavior of observed phase shifts and the shape of the cross sections, both resonance and continuum terms are necessary to be taken into account. If the continuum term is zero, the cross section exhibits a typical Breit-Wigner form. As was discussed by Fano Fa61 (), deviation from the Breit-Wigner form can be investigated by calculating the interference between resonance and continuum terms.

The present method of analyzing the phase shifts and cross sections is also useful in nuclear data evaluations. To study a wider range of nuclear data, it is desirable to develop the method further to include multi-channel systems and to treat many-body systems.


We thank the support by ”RD Platform Formation of Nuclear Reaction Data in Asian Countries (2010-2013)”, JSPS AA Science Platform Program. One of the authors, M.O. thanks for hospitality at the Nuclear Reaction Data Centre (JCPRG) and Theoretical Nuclear Physics Laboratory of Hokkaido University.


  • (1) H. Horiuchi, K. Ikeda and K. Katō, Prog. Theor. Phys. Suppl. 192, 1 (2012).
  • (2) J. Aguilar and J.M. Combes, Commun. Math. Phys. 22, 269 (1971); E. Balslev and J. M. Combes, Commun. Math. Phys. 22, 280 (1971).
  • (3) Y. K. Ho, Phys. Rep. 99, 1 (1983).
  • (4) N. Moiseyev, Phys. Rep. 302, 211 (1998).
  • (5) S. Aoyama, T. Myo, K. Katō and K. Ikeda, Prog. Theor. Phys. 116, 1 (2006).
  • (6) M. Homma, T. Myo, and K. Katō, Prog. Theor. Phys. 97, 561 (1997).
  • (7) T. Myo, R. Ando and K. Katō, Phys. Letters B 691, 150 (2010).
  • (8) T. Myo, Y. Kikuchi and K. Katō, Phys. Rev. C 85, 034338 (2012).
  • (9) T. Myo, A. Ohnishi and K. Katō, Prog. Theor. Phys. 99, 801 (1998).
  • (10) T. Myo, K. Katō, H. Toki, and K. Ikeda, Phys. Rev. C 76, 024305 (2007).
  • (11) T. Myo, K. Katō, S. Aoyama and K. Ikeda, Phys. Rev. C 63, 054313 (2001).
  • (12) R. Suzuki, T. Myo, and K.Katō, Prog. Theor. Phys. 113, 1273 (2005).
  • (13) S. Shlomo, Nucl. Phys. A539, 17 (1992).
  • (14) R. D. Levine, Quantum Mechanics of Molecular Rate Processes, 101 (Clarendon Press, Oxford, 1969).
  • (15) U. Fano, Phys. Rev. 124, 1866 (1961).
  • (16) A. Csótó, B. Gyarmati, A. T. Kruppa, K. F. Pál and N. Moiseyev, Phys. Rev. A 41, 3469 (1990).
  • (17) B. G. Giraud and K. Katō, Ann. of Phys. 308, 115 (2003).
  • (18) B. G. Giraud, K. Katō and A. Ohnishi, J. of Phys. A 37, 11575 (2004).
  • (19) T. Myo and K. Katō, Prog. Theor. Phys. 98, 1275 (1997).
  • (20) E. Hiyama, Y. Kino and M. Kamimura, Prog. Part. Nucl. Phys. 51, 223 (2003).
  • (21) S. Saito, Prog. Theor. Phys. 40, 893 (1968); 41, 705 (1969); Prog. Theor. Phys. Suppl. 62, 11 (1977).
  • (22) V. I. Kukulin, V. M. Krasnopol’sky, V. T. Voronchev and P. B. Sazonov, Nucl. Phys. A417, 128 (1984).
  • (23) D. R. Tilley, J. H. Kelley, J. L. Godwin, D. J. Millener, J. Purcell, C. G. Sheu, and H. R. Weller, Nucl. Phys. A 745, 155 (2004).
  • (24) N. P. Heydenburg, G. M. Temmer, Phys. Rev. 104, 123 (1956); C. W. Reich, J. L. Russell, G. C. Phillips, Phys. Rev. 104, 135 (1956); C. M. Jones, G. C. Phillips, P. D. Miller, Phys. Rev. 117, 525 (1960); T. A. Tombrello, L. S. Senhouse, Phys. Rev. 129, 2252 (1963).
  • (25) H. Kanada, T. Kaneko, S. Nagata and M. Nomoto, Prog. Theor. Phys. 61, 1327 (1979).
  • (26) D. R. Tilley, C. M. Cheves, J. L. Godwin, G. M. Hale, H. M. Hofmann, J. H. Kelley, C. G. Sheu, H. R. Weller, Nucl. Phys. A708, 3 (2002).
  • (27) B. Hoop Jr and H. H. Barschall, Nucl. Phys. 83, 65 (1966); Th. Stammbach and R. L. Walter, Nucl. Phys. A180, 225 (1972).
  • (28) S. M. Austin, H. H. Barschall and R. E. Shamu, Phys. Rev. 126, 1532 (1962); R. E. Shamu and J. G. Jenkin, Phys. Rev. 135, B99 (1964); F. J. Vaughn, W. L. Imhof, R. G. Johnson and M. Walt, Phys. Rev. 118, 683 (1960); Los Alamos Physics and Gryogenics Groups, Nucl. Phys. 12, 291 (1959).
  • (29) H. Masui, S. Aoyama, T. Myo, K. Katō and K. Ikeda, Nucl. Phys. A 673, 207 (2000).
  • (30) I. Shimamura, Adv. Quant. Chem. 63, 165 (2012).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description