Some Remarks on Effective Range Formula in Potential Scattering arxiv:812.347v1 [math-ph] 16 Dec 28 Khosrow Chadan Laboratoire de Physique Théorique Université de Paris XI, Bâtiment 21, 9145 Orsay Cedex, France (Khosrow.Chadan@th.u-psud.fr) Abstract In this paper, we present different proofs of very recent results on the necessary as well as sufficient conditions on the decrease of the potential at infinity for the validity of effective range formulas in 3-D in low energy potential scattering (André Martin, private communication, to appear. See Theorem 1 below). Our proofs are based on compact formulas for the phase-shifts. The sufficiency conditions are well-known since long. But the necessity of the same conditions for potentials keeping a constant sign at large distances are new. All these conditions are established here for dimension 3 and for all angular momenta l. LPT Orsay 8-41 November 28 Unité Mixte de Recherche UMR 8627 - CNRS
I - Introduction We consider potential scattering in 3-dimensions with a real, local, and spherically symmetric potential V(r) [1-3]. The spherical symmetry allows us to consider separately each angular momentum state l, and study the corresponding radial Schrödinger equation. We begin with the S-wave (l = ) in order to see clearly the main points without algebraic complications. We therefore have the following approximation, called effective range formula, for the S-wave phase-shift δ (k) [1-3] : k cotg δ (k) = 1 + 1 k small a 2 r k 2, (1) where k, in appropriate units, is the momentum of the particle. This formula was found very useful in the past for low energy scatterings of particles, especially for nucleon-nucleon scattering. The conditions for the validity of (1) are not very stringent. In general, for having a decent scattering theory, and finite phase-shifts δ l (k), continuous and bounded for all l, and all k, one needs only [1-3] rv(r) L 1 (, ). (2) Under this condition, one has then, at most, a finite number n of bound states, and the Levinson Theorem [1-3] : δ () δ ( ) = nπ if no resonance at E = k 2 =, (3) δ () δ ( ) = nπ + π a resonance at E =. 2 Usually, one chooses δ ( ) =. In the case of a resonance at E =, a becomes infinite. We shall assume (2) throughout this paper. Whatever the number of bound states n, one has then [3, 4] : a) for having a finite scattering length a, it is sufficient to have r 2 V(r) L 1 (1, ) ; (4) 2
b) and for having a finite effective range r it is sufficient to have V(r) r r s, s > 5. (5) Remark 1. Assuming only (2), the phase-shifts δ l (k) are defined only for real values of k, and without extra conditions on the decrease of the potential at infinity, they cannot be extended to complex values of k. However, if one assumes e 2µr V(r) L 1 (1, ), then the S-matrix S l (k) = exp[2iδ l (k)] is meromorphic in Im k < µ, and one can extend (1) in this strip of the k-plane, which may contain the lowest bound states k j = iγ j, γ j < µ. This is the case for neutron-proton scattering. We refer the reader to [3] for more details. Very recently, the above conditions have been made more precise (André Martin, private communication, to appear). His essential points can be summarized as follows : Theorem 1. For the S-wave, and in the absence of a resonance at E = k 2 =, if the potential V(r) keeps a constant sign beyond some finite R, then for having a finite, it is both necessary and sufficient to have r 2 V(r) L 1 (, ). Then, in order to have also r finite, it is again necessary and sufficient to have r 4 V(r) L 1 (, ). One may ask then what happens to the second term in (1) when r 4 V(r) / L 1 (1, )? This also has been studied for V(r) r s at infinity [4], and shown by André Martin on examples. Remark 2. If the effective range formula (1) is exact for all k, the corresponding potential in decreasing exponentially at infinity [5]. It corresponds to a Jost function analytic in the whole k-plane, with one zero and one pole. The purpose of the present paper is to give a different proof of Theorem 1, and generalize also to all anguler momenta l. We end this introduction by quoting the following general theorem, which will 3
be useful for our proofs : Theorem 2. (Hille, [6]). Consider the differential equation φ (r) V(r)φ (r) =. (6) If rv(r) L 1 (, ), this equation has a unique solution χ (r) such that lim χ (r) = 1, lim χ r r (r) =. (7) If V(r) is real, and ultimately keeps a constant sign, the condition on V is both necessary and sufficient for the existence of such a solution χ. There is also a solution φ (r), non-unique (modulo the addition of αχ (r)), such that φ (r) lim = 1, lim φ r r r (r) = 1. (8) If V(r) satisfies also the condition r 2 V(r) L 1 (, ), then χ (r) satisfies lim r[χ (r) 1] =, lim r 2 χ r r (r) =, (9) and there exists a unique solution ψ (r) such that lim [ψ (r) r] =, lim r[ψ r r (r) 1] =. (1) Again, if V(r) keeps a constant sign, the condition r 2 V(r) L 1 (1, ) is both necessary as well as sufficient for the existence of such a solution ψ. Obviously, φ and χ, or ψ and χ, are two independent solutions of (6) since the Wronskians at r =, are, according to (7), (8), and (1), φ χ φ χ = ψ χ ψ χ = 1. (11) Remark 3. In general, φ (), so that, the solution ϕ (r), with ϕ () =, is a linear combination of the fundamental solutions φ and χ : ϕ (r) = αφ (r)+β χ (r). (12) 4
Because of (7), ϕ obviously satisfies (8). If we have r 2 V(r) L 1 (1, ), then ϕ = ψ + αχ, and because of (9), ϕ satisfies also (1). The solution ψ being unique, the same is true for ϕ. Note that ϕ (r) is given by the solution of the Volterra integral equation (here we normalized it to ϕ () = 1!) ϕ (r) = r + whereas χ (r) is the solution of r (r t)v(t)ϕ (t)dt, (13) χ (r) = 1+ r (t r)v(t)χ (t)dt. (14) In case V(r) is positive, it is obvious on these equations, and on the basis of (6), that : a) ϕ (r) is an increasing convex function of r. Also ϕ (r) is increasing, and therefore, ϕ (r) 1 for all r. Assuming now also r2 V(r) L 1 (1, ), one has as seen on (13) ; ϕ ( ) = A <, A = rv(r)ϕ (r)dr, (15) b) χ (r) is a positive, convex, and decreasing function, with χ ( ) = 1. II - Proof of Theorem 1 The proof is based on the single formula for the phase-shift [7] ϕ 2 (k,r) δ (k) = k ϕ 2 (k,r)+k 2 ϕ 2 (k,r) V(r)dr, (16) where one assumes V(r), ϕ being the reduced radial wave function, solution of [1-4] ϕ (k,r)+k 2 ϕ(k,r) = V(r)ϕ(k,r), r [, ), ϕ(k,) =, ϕ (k,) = 1. (17) 5
In some sense, (16) is an absolute definition of the phase-shift since, for potentials satisfying (2), one has, automatically, δ ( ) =. Note that, in (16), the normalization of ϕ is irrelevant. However, we keep the normalization ϕ (k,) = 1 for convenience. The formula (16) is valid for all k >. One can prove, in fact, the following : Theorem 3 ([7]). Under the conditions V(r) and rv(r) L 1 (, ), the formula (16) is valid for all k, δ (k) is a continuous and bounded function, and δ ( ) =. It is also differentiable for k >. In order to have also differentiability at k =, with a finite derivative δ (), it is sufficient to have r 2 V(r) L 1 (, ). In making k in (16), the integral diverges, according to Theorem 2 of Hille, if we only have rv(r) L 1 (, ). However, there is the factor k in front of it, and the net result is δ () =. Remark 4. If we make k = in (16), the denominator becomes ϕ 2 (,r). If there are bound states, ϕ(,r) has, according to the nodal theorem [8], n zeros in (, ). Between these zeros, it has maxima and minima, and so the integral is meaningless. In case of a resonance at k =, one has ϕ(, ) = constant, and ϕ (, ) =, and so, again, the breakdown of (16). If V(r), there are no bound states, and no resonance at E = k 2 = [1, 2, 3, 5]. Finiteness of a. Consider now δ (k) δ (k) δ () a = lim lim k k k k = δ (k = ), (18) assuming V(r). Therefore, according to Theorem 3 above, r 2 V L(, ) is sufficient to secure that δ () is finite, that is, a is finite. So, in essence, reference [7] contains already the proof of the old result that r 2 V L(, ) a finite. From (16), we also have a = V(r) ϕ2 (r) ϕ 2 (r) dr. (19) Again, on the basis of Theorem 2 of Hille, r 2 V L(, ) is also necessary to make a finite. This completes the proof of the first part of Theorem 1. Finally, 6
let us remark that, using ϕ = Vϕ in (19), and integrating by parts, we find ( ) a = lim R ϕ (R) R ϕ also a known result. (R) Finiteness of r. We must first compare (1) with (16). From the well-known expansion of cotg x : cotg x = 1 ( ) 1 x 2x 6 + x2 9 + x4 945 + and using (δ is an odd function of k [1-3]), (2) we find easily δ (k) = nπ ka +bk 3 +, k, (21) r = 2 3 a 2b a 2. (22) We assume now, of course, that a is finite. It follows that the finiteness of r and b are completely equivalent. We can therefore concentrate ourselves on (16). Combining (16) and (19), we find b = lim k δ (k)+ka k 3 1 = lim k k 2 V(r) ϕ2 (k,r) ϕ 2 (r) ϕ2 (r) ϕ 2 (k,r) k 2 ϕ 2 ϕ2 ϕ 2 ( ϕ 2 +k 2 ϕ 2 ) dr. (23) Writing now ϕ 2 ϕ 2 ϕ2 ϕ 2 = (ϕϕ ϕ ϕ )(ϕϕ + ϕ ϕ ), we get, from (17) and ϕ = Vϕ, that (ϕϕ ϕ ϕ ) = k 2 ϕϕ. Since, at r =, ϕ and ϕ vanish, and their derivatives are one (remember (2)!), we find r ϕ(k,r) ϕ (r) ϕ (r) ϕ (k,r) = k 2 ϕ(k,t) ϕ (t) dt. (24) Therefore, taking the limit in (23), we find b = ( 2ϕ (r) ϕ (r) r V(r) ϕ2 (t) dt ϕ4 (r) ) ϕ 4 dr. (25) (r) 7
Now, from the asymptotic behaviour of ϕ (r) r, and ϕ (r) 1, for r, one sees immediately that the fraction in (25) behaves exactly as 1 3 r4. It follows that, for V(r), b is finite if and only if r 4 V(r) L 1 (, ). The same conclusions hold therefore for r. The net conclusion is : Theorem 1. If V(r), r 2 V(r) L 1 (, ) is both necessary and sufficient for having a finite. And for having also r finite, it is both necessary and sufficient to have r 4 V(r) L 1 (, ). We must now include bound states. Bound states. We shall choose, as is usually done, δ ( ) =. If there are bound states of energies γ 2 j, γ j >, j = 1, n, one has the Levinson Theorem (3) : δ() δ() δ( ) = nπ. One can define then δ (k) = δ (k) 2 n j=1 Arctg γ j k. (26) One has now, again, δ( ) =, and δ() = δ() nπ =. Using the inverse problem theory of Gel fand and Levitan [5], one can calculate the potential V(r) corresponding to δ(k). If we write V(r) = V(r)+ V(r), (27) it can be shown that the additional potential V(r) has the following asymptotic behaviour [5] V(r) = r C j e γ j r, (28) where C j are positive constants. V(r) keeps therefore a negative sign for large values of r, and is fast decreasing. We assume, of course, that V(r) is not exponentially decreasing for, otherwise, there would be no problem, and a and r would be finite. So, for the purpose of Theorem 4, V and V are equivalent. Since V has no bound states, or a resonance at E = because of δ() =, the theorem applies to δ(k). Now, as it is easily seen by using the known expansion 8
π Arctg x = x 2 1 x + 1 +, (29) 3x3 the coefficients of the expansions of δ (k) and δ (k) for small k, (21), are related to each other by (remember that, for δ(k), n = in (21)) a = a 2 j b = b 2 j 1 γ 3 j 1 γ j,. (3) They differ by finite quantities, and this completes the proof of Theorem 1. If there is a resonance at E =, in the first term in the sum in (26) there is no factor 2, and one has δ() = π +(n 1)π. In this case, as we said before (Remark 4), 2 ϕ ( ) = C, ϕ ( ) =, and so a =. III - Higher l The validity of the usual scattering theory leading to a continuous and bounded phase-shift δ l (k) is, as we said in the introduction, secured always by rv(r) L 1 (, ). Under this condition, the generalization of (16) for l > is [7] : ϕ 2 δ l (k) = k [ l (k,r) (u l ϕ l u l ϕ l )2 +(v l ϕ V(r) dr, (31) l v l ϕ l )2] where u l (kr)andv l (kr) areappropriately normalizedspherical Bessel andneumann functions, and ϕ l (k,r) the solution of the reduced radial Schrödinger equation [1-4] ϕ l (k,r)+k2 ϕ l (k,r) = ϕ l (k,r) = r r l+1 (2l+1)!! + [ V(r)+ l(l+1) r 2 The Effective Range Formula becomes now [1-3] : ] ϕ l (k,r), (32) k 2l+1 cotg δ l (k) = 1 a l + 1 2 r lk 2 + (33) 9
In order to continue further, one needs now the equivalent of Theorem 2 in the presence of the centrifugal potential l(l + 1)/r 2. One can prove very easily, by mimiking the proofs ot Theorem 2 [6], that one has : Theorem 2. Consider the equation φ (r) V(r)φ (r) = l(l+1) r 2 φ (r). (6 ) If rv(r) L 1 (, ), this equation has a unique solution χ (r) such that ( lim r rl χ (r) = 1, lim r l χ (r) ) =. (7 ) r If V(r) is real, and ultimately keeps a constant sign, the condition on V(r) is both necessary and sufficient for the existence of χ (r). There is also a solution φ (r), non-unique (modulo the addition of αχ (r)), such that ( lim r r l 1 φ (r) = 1, r lim r l φ (r) ) = 1. (8 ) If V(r) satisfies also the condition r 2 V(r) L 1 (, ), then r l χ (r) satisfies lim r[ r l χ (r) 1 ] =, lim r 2( r l χ (r) ) =, (9 ) r r and there exists a unique solution ψ (r) such that [ lim r l ψ (r) r ] [ (r =, lim r l ψ (r) ) ] 1 =. (1 ) r r The solutions φ and χ, or ψ and χ, are two independent solutions of (6 ), and their wronskians are φ χ φ χ = ψ χ ψ χ = (2l+1). (11 ) Again,ifV(r)keepsaconstantsignbeyondsomeR, theconditionr 2 V(r) L 1 (1, ) is both necessary and sufficient for the existence of ψ. In short, one gets Theorem 2 from Theorem 2 by replacing χ by r l χ, and φ and ψ by r l φ and r l ψ, and changing the right-hand side of (11) to (11 ). From these formulas (7 ) to (1 ), one can immediately obtain the asymptotic properties of χ, φ, and ψ themselves, to be used in (31). 1
The proofs are based on the Volterra integral equation for χ and φ [1, 9] : χ (r) = r l φ (r) = r l+1 r r r l+1 r l r l r l+1 V(r )χ (r )dr, (2l+1) r l+1 r l r l r l+1 V(r )φ (r )dr, (2l+1) by iterating them, starting from the zero order solutions r l and r l+1, respectively, and mimiking exactly the proof of Theorem 2 of [6]. We leave the details to the reader. It is quite standard. One can then continue the analysis, as was done for the S-wave, and one finds : (34) Theorem 4. In the absence of a bound state at E =, if the potential keeps a constant sign for r > R, then : a) for having a finite a l, it is necessary and sufficient to have r 2l+2 V(r) L 1 (, ) ; b) the effective range r l is finite if and only if r 2l+4 V(r) L 1 (, ). In case there is a bound state at zero energy (for l 1, it is a real bound state at E =, with an L 2 wave function, instead of being a resonance. It contributes by π to the Levinson Theorem), the scattering length a l is infinite [1-3], as in the case of a resonance when l =. This is obvious on (3), when making γ 1. Acknowledgement. The author is grateful to André Martin for having communicated to him his results prior to publication, and for useful correspondence. He wishes also to thank Professor Kenro Furutani, Reido and Takao Kobayashi, and the Science University of Tokyo (Noda) where this work began, for warm hospitality and financial support. 11
References [1] NewtonRG1982ScatteringTheoryofWavesandParticles2ndedn(NewYork: Springer) [2] Galindo A and Pascual P 199 Quantum Mechanics vol 1 (Berlin : Springer) [3] Mott NF and Massey HSW 1965 The Theory of Atomic Collision 3d edn (Oxford UK : Oxford) chapter II [4] Levy B R and Keller J B 1963 J. Math. Phys. 4 54 ; Handelsman R A, Pao Y-P and Lew J S 1968 Il Nuovo Cimento 55A 453 [5] Chadan K and Sabatier PC 1989 Inverse Problems in Quantum Scattering Theory 2nd edn (Berlin : Springer) [6] Hille E 1969 Lectures on Ordinary Differential Equations (Reading MA : Addison-Wesley) chapter 9 [7] Chadan K, Kobayashi R and Kobayashi T 21 J. Math. Phys. 42 431 [8] Courant R and Hilbert D 1953 Methods of Mathematical Physics (New York : Interscience) vol 1. See also Chadan K and Montes A 1968 J. Math. Phys. 9 1898 [9] Regge T and De Alfaro 1965 Potential Scattering(Amsterdam: North Holland) [1] Chadan K and Kobayashi R 26 J. Phys. A 39 1 12