Calculations of widths in the problem of decay by proton emission
Abstract
We develop a new fully quantum method for determination of widths for nuclear decay by proton emission where multiple internal reflections of wave packet describing tunneling process inside proton–nucleus radial barrier are taken into account. Exact solutions for amplitudes of wave function, penetrability and reflection are found for step barrier (at arbitrary ) which approximates the realistic barrier. In contrast to semiclassical approach and twopotential approach, we establish by this method essential dependence of the penetrability on the starting point in the internal well where proton starts to move outside (for example, for the penetrability is changed up to 200 times; accuracy is ). We impose a new condition: in the beginning of the proton decay the proton starts to move outside from minimum of the well. Such a condition provides minimal calculated halflife and gives stable basis for predictions. However, the halflives calculated by such an approach turn out to be a little closer to experimental data in comparison with the semiclassical halflives. Estimated influence of the external barrier region is up to 1.5 times for changed penetrability.
PACS numbers: 03.65.Xp, 23.50.+z, 27.70.+q
Keywords: tunneling, multiple internal reflections, wave packet, decay by proton emission, penetrability and reflection, halflife
1 Introduction
Nuclei beyond the proton drip line are groundstate proton emitters, i. e. nuclei unstable for emission of proton from the ground state. Associated lifetimes, ranging from sec to few seconds, are sufficiently long to obtain wealth of spectroscopic information. Experimentally, a number of proton emitters has been discovered in the mass region , 150, and 160 (see [1, 2, 3, 4] and references in cited papers). A new regions of proton unstable nuclei is supposed to be explored in close future using radioactive nuclear beams.
Initially, the parent nucleus is in quasistationary state, and the proton decay may be considered as a process where the proton tunnels through potential barrier. In theoretical study one can select three prevailing approaches [5]: approach with distorted wave Born approximation (DWBA), twopotential approach (TPA), and approach for description of penetration through the barrier in terms of onedimensional semiclassical method (WKBA). In systematical study these approaches are correlated between themselves, while calculation of penetrability of the barrier is keystone in successful estimation of gamma widths. While the third approach studies such a question directly, in the first and second approaches the penetrability of the barrier is not studied and the width is based on correlation between wave functions in the initial and final states. However, the most accurate information on amplitudes and phases of these wave functions and correspondence between them can be obtained from unite picture of penetration of proton through the barrier, which is used in the WKBA approach (up to the second order approximation).
The main objective of this paper is to pass from semiclassical unite description of the process of penetration of proton through the barrier used in the WKBA approach to its fully quantum analogue, to put a fully quantum grounds for determination of the penetrability in this problem. In order to provide such a formalism, we have improved method of multiple internal reflections (MIR, see Refs. [6, 7, 8, 9, 10]) generalizing it on the radial barriers of arbitrary shapes. In order to realize this difficult improvement, we have restricted ourselves by consideration of the spherical groundstate proton emitters, while nuclear deformations are supposed to be further included by standard way. This advance of the method never studied before allows to describe dynamically a process of penetration of the proton through the barrier of arbitrary shape in fully quantum consideration, to calculate penetrability and reflection without the semiclassical restrictions, to analyze abilities of the semiclassical and other models on such a basis.
This paper is organized in the following way. In Sec. 2, formalism of the method of multiple internal reflections in description of tunneling of proton through the barrier in proton decay is presented. Here, we give solutions for amplitudes, define penetrability, width and halflife. In Sec. 3, results of calculations are confronted with experimental data and are compared with semiclassical ones. Here, using the fully quantum basis of the method, we study a role of the barrier shape in calculations of widths in details. In particular, we observe essential influence of the internal well before the barrier on the penetrability. We discuss shortly possible interconnections between the proposed approach and other fully quantum methods of calculation of widths. In Sec. 4, we summarize results. Appendixes include proof of the method MIR and alternative standard approach of quantum mechanics used as test for the method MIR and for the results presented.
2 Theoretical approach
An approach for description of onedimensional motion of a nonrelativistic particle above a barrier on the basis of multiple internal reflections of stationary waves relatively boundaries has been studied in number of papers and is known (see [11, 12, 13] and references therein). Tunneling of the particle under the barrier was described successfully on the basis of multiple internal reflections of the wave packets relatively boundaries (approach was called as method of multiple internal reflections or method MIR, see Refs. [6, 7, 8, 9]). In such approach it succeeded in connecting: 1) continuous transition of solutions for packets after each reflection, total packets between the abovebarrier motion and the underbarrier tunneling; 2) coincidence of transmitted and reflected amplitudes of stationary wave function in each spatial region obtained by approach MIR with the corresponding amplitudes obtained by standard method of quantum mechanics; 3) all nonstationary fluxes in each step, are nonzero that confirms propagation of packets under the barrier (i. e. their “tunneling”). In frameworks of such a method, nonstationary tunneling obtained own interpretation, allowing to study this process at interesting time moment or space point. In calculation of phase times this method turns out to be enough simple and convenient [10]. It has been adapted for scattering of the particle on nucleus and decay in the spherically symmetric approximation with the simplest radial barriers [6, 7, 9] and for tunneling of photons [7, 10]. However, further realization of the MIR approach meets with three questions. 1) Question on effectiveness. The multiple reflections have been proved for the motion above one rectangular barrier and for tunneling under it [13, 7, 10]. However, after addition of the second step it becomes unclear how to separate the needed reflected waves from all their variety in calculation of all needed amplitudes. After obtaining exact solutions of the stationary amplitudes for two arbitrary rectangular barriers [9, 6], it becomes unclear how to generalize such approach for barriers with arbitrary complicate shape. In Ref. [14] multiple internal reflections of the waves were studied for tunneling through a number of equal rectangular steps separated on equal distances. However, the amplitudes were presented for two such steps only, in approximation when they were separated on enough large distance, and these solutions in approach of multiple internal reflections were based of the amplitudes of total wave function obtained before by standard method (see Appendix A, eqs. (7), (18), (19) in this paper). So, we come to a serious unresolved problem of realization of the approach of multiple reflections in real quantum systems with complicated barriers, and clear algorithms of calculation of amplitudes should be constructed.
2) Question on correctness. Whether is interference between packets formed relatively different boundaries appeared? Whether does this come to principally different results of the approach of multiple internal reflections and direct methods of quantum mechanics? Note that such interference cannot be appeared in tunneling through one rectangular barrier and, therefore, it could not visible in the previous papers.
3) Question on uncertainty in radial problem. Calculations of halflives of different types of decays based on the semiclassical approach are prevailing today. For example, in Ref. [15] agreement between experimental data of decay halflives and ones calculated by theory is demonstrated in a wide region of nuclei from up to nucleus with and (see Ref. [16] for some improved approaches). In review [17] methodology of calculation of halflives for spontaneousfission is presented (see eqs. (21)–(24) in p. 321). Let us consider protondecay of nucleus where proton penetrates from the internal region outside with its tunneling through the barrier. At the same boundary condition, reflected and incident waves turn out to be defined with uncertainty. How to determine them? The semiclassical approach gives such answer: according to theory, in construction of well known formula for probability we neglect completely by the second (increasing) item of the wave function inside tunneling region (see Ref. [18], eq. (50.2), p. 221). In result, equality has no any sense (where and are coefficients of penetrability and reflection). Condition of continuity for the wave function and for total flux is broken at turning point. So, we do not find reflection . We do not suppose on possible interference between incident and reflected waves which can be non zero. The penetrability is determined by the barrier shape inside tunneling region, while internal and external parts do not take influence on it. The penetrability does not dependent on depth of the internal well (while the simplest rectangular well and barrier give another exact result). But, the semiclassical approach is so prevailing that one can suppose that it has enough well approximation of the penetrability estimated. It turns out that if in fully quantum approach to determine the penetrability through the barrier (constructed on the basis of realistic potential of interaction between proton and daughter nucleus) then one can obtain answer “no”. Fully quantum penetrability is a function of new additional independent parameters, it can achieve essential difference from semiclassical one (at the same boundary condition imposed on the wave function). This will be demonstrated below.
2.1 Decay with radial barrier composed from arbitrary number of rectangular steps
Let us assume that starting from some time moment before decay the nucleus could be considered as system composite from daughter nucleus and fragment emitted. It‘s decay is described by a particle with reduced mass which moves in radial direction inside a radial potential with a barrier. We shall be interesting in the radial barrier of arbitrary shape, which has successfully been approximated by finite number of rectangular steps:
(1) 
where are constants (). We define the first region 1 starting from point , assuming that the fragment is formed here and then it moves outside. We shall be interesting in solutions for above barrier energies while the solution for tunneling could be obtained after by change . A general solution of the wave function (up to its normalization) has the following form:
(2) 
(3) 
where and are unknown amplitudes, and are unknown amplitudes of transmission and reflection, is spherical function, are complex wave vectors. We shall be looking for solution for such problem in approach of multiple internal reflections (we restrict ourselves by a case of orbital moment while its nonzero generalization changes the barrier shape which was used as arbitrary before in development of formalism MIR and, so, is absolutely non principal).
According to the method of multiple internal reflections, scattering of the particle on the barrier is considered on the basis of wave packet consequently by steps of its propagation relatively to each boundary of the barrier (idea of such approach can be understood the most clearly in the problem of tunneling through the simplest rectangular barrier, see [7, 9, 10] and Appendix A where one can find proof of the method and analysis of its properties). Each step in such consideration of propagation of the packet will be similar to on from the first steps, independent between themselves. From analysis of these steps recurrent relations are found for calculation of unknown amplitudes , , and for arbitrary step , summation of these amplitudes are calculated. We shall be looking for the unknown amplitudes, requiring wave function and its derivative to be continuous at each boundary. We shall consider the coefficients , , and , , as additional factors to amplitudes . Here, bottom index denotes number of the region, upper (top) signs “” and “” denote directions of the wave to the right or to the left, correspondingly. At the first, we calculate , … and , …:
(4) 
Using recurrent relations:
(5) 
and selecting as starting the following values:
(6) 
we calculate successively coefficients …, … and …. At finishing, we determine coefficients :
(7) 
the amplitudes of transmission and reflection:
(8) 
and corresponding coefficients of penetrability and reflection :
(9) 
We check the property:
(10) 
which should be the test, whether the method MIR gives us proper solution for wave function. Now if energy of the particle is located below then height of one step with number , then for description of transition of this particle through such barrier with its tunneling it shall need to use the following change:
(11) 
For the potential from two rectangular steps (with different choice of their sizes) after comparison between the all amplitudes obtained by method of MIR and the corresponding amplitudes obtained by standard approach of quantum mechanics, we obtain coincidence up to first 15 digits. Increasing of number of steps up to some thousands keeps such accuracy and fulfillment of the property (10) (see Appendix B where we present shortly the standard technique of quantum mechanics applied for the potential (19) and all obtained amplitudes). This is important test which confirms reliability of the method MIR. So, we have obtained full coincidence between all amplitudes, calculated by method MIR and by standard approach of quantum mechanics, and that is way we generalize the method MIR for description of tunneling of the particle through potential, consisting from arbitrary number of rectangular barriers and wells of arbitrary shape.
2.2 Width and halflife
We define width of the decay of the studied quantum system by following the procedure of Gurvitz and Kälbermann [19]:
(12) 
where is the preformation probability and is the normalization factor. is the penetrability coefficient in propagation of the particle from the internal region outside with its tunneling through the barrier, which we shall calculate by approach MIR or by approach WKB. In approach WKB we define it so:
(13) 
where and are the second and third turning points. According to [15], the normalization factor is given by simplified way by or by improved way by so:
(14) 
The halflife of the decay is related to the width by well known expression:
(15) 
For description of interaction between proton and the daughter nucleus we shall use the spherical symmetric proton–nucleus potential (at case ) in Ref. [20] having the following form:
(16) 
where , and are Coulomb, nuclear and centrifugal components
(17) 
Here, and are the nucleon and proton numbers of the daughter nucleus, is the value for the protondecay, is the strength of the nuclear component, is radius of the daughter nucleus, is the effective radius of the nuclear component, is diffuseness. All parameters are defined in Ref. [20]. Note that in this paper we are concentrating on the principal resolution of question to provide fully quantum basis for calculation of the penetrability and halflife in the problem of the proton decay, while the proton–nucleus potential can be used in simple form that does not take influence on the reliability of the developed methodology of multiple internal reflections absolutely and could be naturally included for modern more accurate models.
3 Results
Today, there are a lot of modern methods able to calculate halflives, which have been studied experimentally well. So, we have a rich theoretical and experimental material for analysis. We shall use these nuclei: , and . Such a choice we explain by that they have small coefficient of quadruple deformation and at good approximation can be considered as spherical (we have ). We shall study protondecay on the basis of leaving of the particle with reduced mass from the internal region outside with its tunneling through the barrier. This particle is supposed to start from and move outside. Using technique of the and coefficients in eqs. (4)–(6), we calculate total amplitudes of transmission and reflection by eqs. (8), the penetrability coefficient by eqs. (9). We check the found amplitudes, coefficients and comparing them with corresponding amplitudes and coefficients calculated by standard approach of quantum mechanics presented in Appendix B. We restrict ourselves by eq. (14) for and find width by eq. (12) and halflive by eq. (15). We define the penetrability by eq. (13), calculate width and halflive by eqs. (12) and (15).
3.1 Dependence of the penetrability on the starting point
The first interesting result which we have obtained is essential dependence of penetrability on the position of the first region where we localize the wave incidenting on the barrier. In particular, we have analyzed how much the internal boundary takes influence on the penetrability. Taking into account that width of each interval is 0.01 fm, we consider point as a starting point (with error up to 0.01 fm), from here proton begins to move outside and is incident on the internal part of the barrier in the first stage of the proton decay. In the Fig. 1 [left panel] one can see that halflive of the proton decay of is changed essentially at displacement of . So, we establish essential dependence of the penetrability on the starting point , where the proton starts to move outside by approach MIR.
At fm this dependence allows us to achieve very close coincidence between the halflive calculated by the approach MIR and experimental data.
3.2 Dependence of the penetrability on the external region
The region of the barrier located between turning points and is main part of the potential used in calculation of the penetrability in the semiclassical approach (up to the second correction), while the internal and external parts of this potential do not take influence on it. Let us analyze whether convergence exists in calculations of the penetrability in the approach MIR if to increase the external boundary (). Keeping width of each interval (step) to be the same, we shall increase (through increasing number of intervals in the external region), starting from the external turning point , and calculate the corresponding penetrability . In Fig. 1 [central panel] one can see how the penetrability is changed for with increasing . Dependence of the halflife on is shown in the next figure 1 [right panel]. One can see that the method MIR gives convergent values for the penetrability and halflife at increasing of . From such figures we find that inclusion of the external region into calculations changes the halflife up to 1.5 times ( sec is the minimal halflife calculated at fm, sec is the halflife calculated at fm, or 50 percents). So, error in determination of the penetrability in the semiclassical approach (if to take the external region into account) is expected to be the same as a minimum on such a basis.
3.3 Results for the proton emitters , and
So, the fully quantum study of the penetrability of the barrier for the proton decay give us its large dependence on the starting point. In order to give power of predictions of halflives calculated by the approach MIR, we need to find recipe able to resolve such uncertainty in calculations of the halflives. We shall introduce the following hypothesis: we shall assume that in the first stage of the proton decay proton starts to move outside the most probably at the coordinate of minimum of the internal well. If such a point is located in the minimum of the well, the halflive obtains minimal value. So, as criterion we could use minimum of halflive for the given potential, which has stable basis. Let us analyze which results such approach gives. We shall compare the halflives calculated by approach MIR and by the semiclassical approach with experimental data. We should take into account that the halflives obtained before are for the proton occupied ground state while it needs to take into account probability that this state is empty in the daughter nucleus. In order to obtain proper values for the halflives we should divide them on the spectroscopic factor (which we take from [5]), and then to compare them with experimental data. Results of such calculations are presented in Table 1. One can see that the calculated halflives by MIR approach turn out to be a little closer to experimental data in comparison with halflives obtained by the semiclassical approach.
Parent nucleus  Halflivevalues, sec  Points, fm  

Nucleus  , MeV  
0.947  0.66  7.43  
1.214  0.59  3.32  7.34  
1.086  0.51  3.41  7.46 
3.4 Comparison with other approaches of calculations of widths of proton decay
Halflife of the proton decay is defined on the basis of width which can be calculated by different approaches. For determination of width we shall use systematics of different approaches proposed in Ref. [5]. The proton emitters are narrow resonances with extremely small widths. Perturbative approach based on standard reaction theory could be expected to be accurate. Let us analyze two following approaches in such a direction.
3.4.1 The distorted wave Born approximation method
The resonance width can be expressed through transition amplitude, which in the distorted wave Born approximation (DWBA) is given so [5]:
(18) 
The DWBA calculations of the decay width requires knowledge of the quasistationary initial state wave function, , the final state wave function, , and interaction potential. The initial state wave function, , is written as a product of the daughternucleus wave function, , and the proton wave function, . The radial wave function of the proton is found by numerically integrating the Schrödinger equation with onebody potential, and it should be irregular part of the Coulomb wave function, , in asymptotic region. So, such wave function is complex and it defines nonzero flux. As we use condition of continuity of total flux (i.e. absence of sources inside spatial region) we cannot obtain zero wave function in whole region of its definition, and at , in particular.
In the final state the wave function of the decaying nuclear system can be written as a product of the intrinsic wave function of the proton and the daughter nucleus (an inner core). Radial part of the proton wave function is , where is the regular Coulomb function. By other words, this wave function is real, and, therefore, it gives zero flux exactly determined on the basis of the total wave function in the initial state. The total wave functions in the initial and final states correspond to different processes (with different total fluxes). This confirms that they, complete wave functions, do not take reflection from the barrier inside the internal region into account (but they are defined by different boundary conditions in the initial and final states only). Here, question about determination of the decay width is passed on successful determination of perturbation of the potential (that has another physical basis for the definition of the decay width as definition on the basis of the penetrability of the barrier). However, the question about separation of the total wave function in the internal region before the barrier into the incident and reflected waves remains unresolved in the DWBA method.
Now, if we pass from real radial potential in optical model approach to complex one, then we shall introduce new additional independent parameter into our problem while the penetrability could be calculated for real radial barrier. Essential point in determination of the decay width in the DWBA method is accurate normalization of the wave functions in the initial and final states. It could introduce some (essential) uncertainty in calculation of width also while the penetrability is independent on such normalization absolutely.
One can calculate the decay width through timereversed capture process. However, in such calculations shape of the barrier is approximated by inverse oscillator (or other potentials with knowing exact solutions of the wave function) and the penetrability for such a barrier could be calculated. It is clear that both internal well and external region do not take influence on results absolutely (like calculations in semiclassical approach). But, this is possible to resolve this problem accurately and taking whole studied shape of the potential barrier into account that we have demonstrated above in the fully quantum approach MIR.
3.4.2 The twopotential method
In the modified twopotential approach (TPA) introduced by Gurvitz and Kalbermann in [19] (details and examples can be found in [22], see also [21, 5, 23]) the decay width is defined so (see (16) in [5], and some details):
(19) 
where , is reduced mass, is radial coordinate of the barrier height, is the radial wave function for the first radial potential including internal well up to point , is the regular radial wave function for the second radial potential including external region, starting from point and without the internal well and with asymptotic behavior
(20) 
Both wave functions are real and defined at different energy levels. So, in the TPA approach we do not consider fluxes and do not calculate penetrability. We do not study possible reflection of proton from the barrier in the internal well. We escape from a problem of separation of the total wave function in the internal well into the incident and reflected waves which takes influence on the resulting penetrability essentially (for example, for the simplest rectangular barrier with rectangular well such an uncorrect separation of the same exact wave function can give infinite penetrability that is explained by increased role of interference between incident and reflected waves). Success in obtaining the resulting width is dependent on accuracy of correspondence between internal and external wave functions and which should be calculated from different Schrödinger equations with independent normalization. The correspondence between these wave functions is determined concerning only one boundary point (or it possible shift [23]) separating two potentials and boundary conditions at or at . In contrary, the correspondence between the incident, transmitted and reflected waves in the MIR approach is determined concerning the barrier as the whole potential (with needed restrictions of the radial problem) that corresponds to fully quantum and unified consideration of penetration of the proton through the barrier shown in principle of nonlocality of quantum mechanics. In particular, the transmitted wave in the MIR approach is strongly dependent on the depth of the internal well and its shape, while the external wave function in the TPA approach is absolutely independent on these depth and shape (such a dependence can be found in the wave function , but starting from the simplest WKB approach factor F directly includes it also). By other words, we have strong correspondence between incident, reflected and transmitted waves in the MIR approach and a possible week correspondence between the internal and external wave functions in the TPA approach. This plays the essential role in calculations of the decay widths and explains so large difference between the essential dependence of penetrability on the starting point in the MIR approach and practically full absence of such a dependence in the TPA approach.
The simplest example demonstrated why this dependence really exists and it could be not small, can be found in classical tasks of quantum mechanics. Let us consider definition of the penetrability in [18] (see eq. (25.3), p. 103):
(21) 
where is the penetrability, and are wave numbers of transmitted and incident waves, i.e. concerning the left asymptotic part of the potential and its asymptotic right part (see Fig. 5 in [18], p. 103), is the transmitted amplitude of the wave function. This formula demonstrates that decreasing of the left part of potential increases the wave number (as is connected with asymptotic presentations (25.1) and (25.2) of waves) and, so, changes the total penetrability . Result on the essential dependence of the penetrability of the starting point above has the similar sense, but has been obtained concerning the realistic barrier with the internal well and takes into account change of the internal amplitudes also. This contradicts with a possible little dependence of penetrability on the shape of the internal well in the TPA approach. So, these points seem to be reduction of the TPA approach, and confirm that this approach does not determine the penetrability in the fully quantum consideration in the problem of proton decay. At the same time, comparison of results obtained by such approach and results obtained by principally other fully quantum developments sometimes leads to some confusion as the TPA approach has been called as the fully quantum. So, approaches for determination of the decay widths on the basis of penetrability are physically motivated, could be more accurate and have perspective for research.
4 Conclusions
The new fully quantum method (called as the method of multiple internal reflections, or MIR) for calculation of widths for the decay of the nucleus by emission of proton in the spherically symmetric approximation and the realistic radial barrier of arbitrary shape is presented. Note the following:

Solutions for amplitudes of wave function (described motion of the proton from the internal region outside with its tunneling through the barrier), penetrability and reflection are found by the method MIR for step radial barrier at arbitrary . These solutions are exactly solvable and have been obtained in the fully quantum approach for the first time. At limit these solutions could be considered as exact ones for the realistic proton–nucleus potential with needed arbitrary barrier and internal hole. Estimated error of the achieved results is .

In contrast to the semiclassical approach and the TPA approach, the approach MIR gives essential dependence of the penetrability on the starting point inside the internal well where proton starts to move outside in the beginning of the proton decay. For example, the penetrability of the barrier calculated by MIR approach for is changed up to 200 times in dependence on position of (see Fig. 1, the left panel). The amplitudes calculated by MIR approach we compared with the corresponding amplitudes obtained (for the same potential) by independent standard stationary approach of quantum mechanics presented in Appendix B and we obtained coincidence up to first 15 digits for all considered amplitudes. This important test confirms that presence of the essential dependence of the penetrability of the starting point is result independent on the fully quantum method applied. Such a result could be connected with a possibility to introduce initial condition which could be imposed on proton decay in the fully quantum consideration. Comparison with the WKB and TPA approaches shows that such approaches have no such a perspective (having physical sense and opening a possibility to obtain a new information about the proton decay), which fully quantum study of the penetrability gives.

In order to resolve uncertainty in calculations of the halflives caused by the dependence of the penetrability on , we have introduced the hypothesis: in the first stage of the proton decay the proton starts to move outside at the coordinate of minimum of the internal well. Such condition provides minimal value for the calculated halflife and gives stable basis for predictions in the MIR approach. However, the halflives calculated by the MIR approach turn out to be a little closer to experimental data in comparison with the halflives obtained by the semiclassical approach (see Tabl. 1).

Taking the external region of the potential after the barrier into account, halflive calculated by the MIR approach is changed up to 1.5 times (see Fig. 1, the right panel).
A main advance of the MIR method developed in this paper is not a new attempt to describe experimental data of halflives more accurately than other approaches do this, but rather this method seems to be the first tools for estimation of the penetrability of any desirable barrier of the proton decay in the fully quantum consideration.
Appendix A Tunneling of packet through onedimensional rectangular step
Main ideas and formalism of the multiple internal reflections can be the most clearly analyzed in the simplest problem of tunneling of the particle through onedimensional rectangular barrier in whole axis [6, 7, 8, 9, 10]. Let us consider a problem of tunneling of a particle in a positive direction through an onedimensional rectangular potential barrier (see Fig. 2). Let us label a region I for , a region II for and a region III for , accordingly. We shall study an evolution of its tunneling through the barrier.
In standard approach, with energy less than the barrier height the tunneling evolution of the particle is described using a nonstationary propagation of WP
(22) 
where stationary WF is:
(23) 
and , , and are the total energy and mass of the particle, accordingly. The weight amplitude can be written in the standard gaussian form and satisfies to a requirement of the normalization , value is an average energy of the particle. One can calculate coefficients , , and analytically, using a requirements of a continuity of WF and its derivative on each boundary of the barrier. Substituting in eq. (22) instead of the incident , transmitted or reflected part of WF , defined by eq. (23), we receive the incident, transmitted or reflected WP, accordingly.
We assume, that a time, for which the WP tunnels through the barrier, is enough small. So, the time necessary for a tunneling of an particle through a barrier of decay in decay of a nucleus, is about seconds. We consider, that one can neglect a spreading of the WP for this time. And a breadth of the WP appears essentially more narrow on a comparison with a barrier breadth. Considering only subbarrier processes, we exclude a component of waves for abovebarrier energies, having included the additional transformation
(24) 
where function satisfies to the requirement
The method of multiple internal reflections considers the propagation process of the WP describing a motion of the particle, sequentially on steps of its penetration in relation to each boundary of the barrier [11, 12, 13]. Using this method, we find expressions for the transmitted and reflected WP in relation to the barrier. At the first step we consider the WP in the region I, which is incident upon the first (initial) boundary of the barrier. Let’s assume, that this package transforms into the WP, transmitted through this boundary and tunneling further in the region II, and into the WP, reflected from the boundary and propagating back in the region I. Thus we consider, that the WP, tunneling in the region II, is not reached the second (final) boundary of the barrier because of a terminating velocity of its propagation, and consequently at this step we consider only two regions I and II. Because of physical reasons to construct an expression for this packet, we consider, that its amplitude should decrease in a positive direction. We use only one item in eq. (23), throwing the second increasing item (in an opposite case we break a requirement of a finiteness of the WF for an indefinitely wide barrier). In result, in the region II we obtain:
(25) 
Thus the WF in the barrier region constructed by such way, is an analytic continuation of a relevant expression for the WF, corresponding to a similar problem with abovebarrier energies, where as a stationary expression we select the wave , propagated to the right.
Let’s consider the first step further. One can write expressions for the incident and the reflected WP in relation to the first boundary as follows
(26) 
A sum of these expressions represents the complete WF in the region I, which is dependent on a time. Let’s require, that this WF and its derivative continuously transform into the WF (25) and its derivative at point (we assume, that the weight amplitude differs weakly at transmitting and reflecting of the WP in relation to the barrier boundaries). In result, we obtain two equations, in which one can pass from the timedependent WP to the corresponding stationary WF and obtain the unknown coefficients and .
At the second step we consider the WP, tunneling in the region II and incident upon the second boundary of the barrier at point . It transforms into the WP, transmitted through this boundary and propagated in the region III, and into the WP, reflected from the boundary and tunneled back in the region II. For a determination of these packets one can use eq. (22) with account eq. (24), where as the stationary WF we use:
(27) 
Here, for forming an expression for the WP reflected from the boundary, we select an increasing part of the stationary solution only. Imposing a condition of continuity on the timedependent WF and its derivative at point , we obtain 2 new equations, from which we find the unknowns coefficients and .
At the third step the WP, tunneling in the region II, is incident upon the first boundary of the barrier. Then it transforms into the WP, transmitted through this boundary and propagated further in the region I, and into the WP, reflected from boundary and tunneled back in the region II. For a determination of these packets one can use eq. (22) with account eq. (24), where as the stationary WF we use:
(28) 
Using a conditions of continuity for the timedependent WF and its derivative at point , we obtain the unknowns coefficients and .
Analyzing further possible processes of the transmission (and the reflection) of the WP through the boundaries of the barrier, we come to a deduction, that any of following steps can be reduced to one of 2 considered above. For the unknown coefficients , , and , used in expressions for the WP, forming in result of some internal reflections from the boundaries, one can obtain the recurrence relations:
(29) 
Considering the propagation of the WP by such way, we obtain expressions for the WF on each region which can be written through series of multiple WP. Using eq. (22) with account eq. (24), we determine resultant expressions for the incident, transmitted and reflected WP in relation to the barrier, where one can need to use following expressions for the stationary WF:
(30) 
Now we consider the WP formed in result of sequential reflections from the boundaries of the barrier and incident upon one of these boundaries at point () or at point (). In result, this WP transforms into the WP , transmitted through boundary with number , and into the WP , reflected from this boundary. For an independent on parts of the stationary WF one can write:
(31) 
where the sign “+” (or “”) corresponds to the WP, tunneling (or propagating) in a positive (or negative) direction and incident upon the boundary with number . Using and , one can precisely describe an arbitrary WP which has formed in result of multiple reflections, if to know a “path” of its propagation along the barrier. Using the recurrence relations eq. (29), the coefficients and can be obtained.
(32) 
Using the recurrence relations, one can find series of coefficients , , and . However, these series can be calculated easier, using coefficients and . Analyzing all possible “paths” of the WP propagations along the barrier, we receive:
(33) 
where
(34) 
All series , , and , obtained using the method of multiple internal reflections, coincide with the corresponding coefficients , , and of the eq. (23), calculated by a stationary methods [18]. Using the following substitution
(35) 
where is a wave number for a case of abovebarrier energies, expression for the coefficients , , and for each step, expressions for the WF for each step, the total eqs. (33) and (34) transform into the corresponding expressions for a problem of the particle propagation above this barrier. At the transformation of the WP and the timedependent WF one can need to change a sign of argument at function. Besides the following property is fulfilled:
(36) 
Appendix B Direct method
We shall add shortly solution for amplitudes of the wave function obtained by standard technique of quantum mechanics which could be obtained if to use only condition of continuity of the wave function and its derivative at each boundary, but on the whole region of the studied potential. At first, we find functions and (from the first boundary):
(37) 
Then, using the following recurrent relations:
(38) 
we calculate next functions , , …, and by such a formula:
(39) 
the functions , , …. From and we find amplitudes , and amplitude of transmission :
(40) 
Now using the recurrent relations:
(41) 
and such a formula:
(42) 
we consistently calculate the amplitudes , , , …, . At finishing, we find amplitude of reflection :