Holomorphic Extensions Associated with Fourier–Legendre Series and the Inverse Scattering Problem

In this paper, we consider the inverse scattering problem and, in particular, the problem of reconstructing the spectral density associated with the Yukawian potentials from the sequence of the partial-waves f` of the Fourier–Legendre expansion of the scattering amplitude. We prove that if the partial-waves f` satisfy a suitable Hausdorff-type condition, then they can be uniquely interpolated by a function f̃ (λ) ∈ C, analytic in a half-plane. Assuming also the Martin condition to hold, we can prove that the Fourier–Legendre expansion of the scattering amplitude converges uniformly to a function f (θ) ∈ C (θ being the complexified scattering angle), which is analytic in a strip contained in the θ-plane. This result is obtained mainly through geometrical methods by replacing the analysis on the complex cos θ-plane with the analysis on a suitable complex hyperboloid. The double analytic symmetry of the scattering amplitude is therefore made manifest by its analyticity properties in the λand θ-planes. The function f (θ) is shown to have a holomorphic extension to a cut-domain, and from the discontinuity across the cuts we can iteratively reconstruct the spectral density σ(μ) associated with the class of Yukawian potentials. A reconstruction algorithm which makes use of Pollaczeck and Laguerre polynomials is finally given.


Introduction
In the mathematical theory of the potential scattering, the inverse problem (at fixed energy) may be formulated as follows (see [1]): be the normalized Fourier-Legendre expansion of the scattering amplitude, P (·) being the Legendre polynomials.
Then, given a class V of potentials, we have the following questions: (i) Does there exist a potential V ∈ V that generates the sequence of partial-waves { f } ∞ =0 ? (ii) Is this potential unique? (iii) Does this potential depend continuously upon the given data { f } ∞ 0 ? (iv) Can we set up an algorithm for recovering this potential from the observations? Can we show its convergence accounting also for that the actual data are either affected by noise and finite in number? (v) Can we provide stability estimates for this algorithm?
Here we approach this problem assuming that V is the class of Yukawian potentials. A Yukawian potential is represented by a spherically symmetric function on R 3 of the following form: In the past sixties, various authors, in particular Martin and Targonski [2], were able to give a positive answer to the second question (i.e., the question of the uniqueness) for this class of potentials. The same authors gave also a necessary, but in general not sufficient, condition on the coefficients f to guarantee the existence of the potential; hereafter, this condition will be called the M.T. condition.
The question of the existence of the potential, i.e., point (i) above, can be split into two parts: (a) The existence of the spectral density σ(µ) in (2), associated with the partial-waves f and the Legendre expansion (1). (b) The convergence of integral (2), which gives V µ 0 (r) when σ(µ) has been determined from data.
Part (a) of this question can actually be solved by seeking for conditions on the partialwaves { f } ∞ = 0 which are sufficient to guarantee a holomorphic extension of series (1) to a suitable slit-domain in C. This connection becomes clear if we recall that the convergence properties of expansion (1), as described, for instance, by Walsh [3], are strictly analogous to the well-known convergence properties of Taylor series. In the Fourier-Legendre expansions, the region of convergence, instead of being the disk of radius ρ of the Taylor series, is an ellipse E ρ with foci in ±1 and radius ρ . =(semi-minor+semi-major)-axis. The Legendre expansion converges inside the largest ellipse within which the scattering amplitude f (θ) is holomorphic (θ ∈ C being the complexified scattering angle). The main aim of this paper is to prove that if the partial-waves { f } ∞ =0 (or a subset of them) satisfy a suitable Hausdorff-type condition, then the scattering amplitude f (θ) = f (cos θ) admits an analytic continuation to the whole complex θ-plane cut from µ 0 to +∞. Then, the spectral density σ(µ) can be iteratively reconstructed by starting from the jump function which describes the discontinuity of f (θ) across the cut. Furthermore, it is worth observing that, adopting this setting, it is possible to give a suggestive probabilistic interpretation of the Hausdorff conditions and of the partial-waves f . More precisely, we will see in Section 2.2 that the partial-waves f can be regarded as the values of a harmonic function associated with a suitable Markov process. Part (b) of question (i) remains still unanswered.
Concerning question (iii) the answer is negative: the inverse scattering problem (at fixed energy) is known to be ill-posed in the sense of Hadamard, with the data f being necessarily affected by noise and, furthermore, the experiments can provide only a finite number of them. Nevertheless, in [4] (see also [5]) it has been proposed an algorithm, based on the Pollaczek polynomials, which is able to reconstruct iteratively the spectral density σ(µ) starting from a finite set of noisy data , where ε characterizes some kind of bound on the noise. It was furthermore possible to evaluate stability estimates and to prove that the approximation error, associated with the numerical algorithm, tends to zero in the topology of the L 2 -norm, as the number of data tends to infinity and the bound ε on the noise tends to zero. We can thus say that the questions (iii) and (v) listed above, have received a quite satisfactory answer in [4]. However, that paper presents some limits and essentially two main defects: (i) The condition imposed on the partial-waves was unnecessarily restrictive. (ii) The holomorphic extension associated with the Fourier-Legendre series (1) was only conjectured but not rigorously proved.
One of the purposes of the present paper is indeed to remedy these defects. First (see Section 2.1), we introduce an appropriate Hausdorff-type condition on the partial-waves which allows us to prove that the set of the partial-waves admit a unique Carlsonian interpolant. Then, we study the Hardy space to which this interpolant belongs. In Section 2.2 we give the probabilistic interpretation of the Hausdorff condition in terms of Markov processes. Next, we obtain the holomorphic extension of the scattering amplitude f (θ) to the cut cos θ-plane, whose proof is performed in various steps. The first step consists in replacing the complex cos θ-plane (θ ∈ C) with a complex hyperboloid X (c) . The hyperboloid X (c) contains the following submanifolds: the Euclidean sphere S 2 = (iR × R 2 ) ∩ X (c) where the SO(3, R) harmonic analysis is done, and the real onesheteed hyperboloid X = R 3 ∩ X (c) that contains the cut support. The second step is a fibration on a meridian hyperbola of X (c) , implemented by means of a Radon-type transformation where the role played by the planes in the ordinary Radon transformation is now played by the horocycles. By this fibration, we can thus to reduce the harmonic analysis on the complex hyperboloid to that on a complex one-dimensional hyperbola, which contains the Euclidean circle and the real hyperbola. In this way, series (1) can be regarded as a trigonometrical series on the Euclidean circle. This last step is proved in Appendix A, while the geometrical method we used is illustrated in Appendix B. Once the trigonometrical series is obtained, we can prove that it admits a holomorphic extension to a complex cut-plane. This is proved in Section 3. Finally, in Section 4 we return back to the complex cos θ-plane (i.e., the complex hyperboloid X (c) ) by inverting the Radon transformation, and finally prove the analytic extension of the Fourier-Legendre series (1). The inversion of the Radon transform is analyzed in Appendix C. The unique interpolation of the coefficients { f } ∞ =0 turns out to be a Laplace-type transform (precisely, the composition of Laplace and Radon transformations) of the jump across the cut. Therefore, by inverting this transformation, we can reconstruct the discontinuity function which leads to the spectral density σ(µ). This latter step is described in Section 5. In the same section we briefly show how the jump function across the cuts can be reconstructed by means of the Pollaczek polynomials. Finally, Section 6 is devoted to the analysis of the connection between the first order Born approximation of the scattering amplitude and the Laplace-composed-with-Radon transform, introduced in Appendix C.

Hausdorff Moments and Hardy Spaces
Let N 0 denote the non-negative integers: N 0 = N ∪ {0} and C (+) −1/2 the complex halfplane: C Consider a sequence {a } ∞ =0 of (real) numbers a , and denote by ∆ the difference operator Then we have: (for every k 0); ∆ 0 is the identity operator, by definition. Now, suppose that there exists a positive constant M such that: We shall refer to (5) as the Hausdorff condition in view of its relevance for the solution of the Hausdorff moment problem. In fact, it can be proved [6] that condition (5) is necessary and sufficient to represent the numbers a , ∈ N 0 , as moments of a suitable function, that is: where u(x) ∈ L 2+ε (0, 1). Next, the tool we use to guarantee the uniqueness of the interpolation of a sequence of numbers is Carlson's theorem [7].
We shall refer to conditions (i) and (ii) as Carlson's bound. Moreover, an analytic function which interpolates a sequence of numbers { f } ∞ =0 and satisfy conditions (i) and (ii) above will be called a Carlsonian interpolant. Now, we can prove the following proposition (see also [8] for a preliminary version of this proposition).

Hausdorff Moments and Markov Processes
In this subsection we follow the paper of Watanabe [9]. Let (Ω, F , P) be an abstract probability field. If {y (ω); 1} is a sequence of mutually independent random variables on (Ω, F , P) and each random variable satisfies then this sequence is called a Bernoulli sequence and will be denoted by B(p). In particular, in what follows we shall consider B(1/2). Let E be the set of all couples ( , i) such that i; , i ∈ N 0 . Next, we consider the discrete Markov process x , conservative over E, which corresponds to the stochastic matrix of the transition probabilities: Π (( , i), ( + 1, i)) = Π (( , i), ( + 1, i + 1) = 1 2 ). It is called the discrete Markov process x attached to B(1/2) since this process can be derived from the Bernoulli sequence B(1/2). Let us note that where k 0, ( , i), (m, j) ∈ E. The kernel K(( , i), (m, j)) is given by: σ(·) denoting the Markov time. Now, consider an infinite sequence (m k , j k ) having no limit points in E such that for a suitable 0 b 1. Then, by the Stirling formula and taking into account equality (15), from (14) we obtain: Thus, the Martin boundary M (see [10,11]) which is induced by the process x coincides with the interval [0, 1] as a set. Hence, the generalized Poisson kernel K(( , i), b) is given by 2 b ( −i) (1 − b) i . We finally observe that for a function u defined on E the expectation E (0,0) is given by: Then, the following propositions proved by Watanabe [9] can be stated.
if and only if u is x -harmonic and E (0,0) (|u(x )|) is bounded with respect to .
then the function u( , i) is a x -harmonic function and can be represented by Formula (19).

Holomorphic Extension Associated with the Trigonometrical Series
In Appendix A we show that a series expansion in terms of Legendre polynomials P (cos u) can be regarded as a trigonometrical series, i.e., as a series of the type 1 2π ∑ ∞ =0 f e −i t (see Formulae (A3), (A4) and (A19)). Then, it is quite natural to introduce the complex plane C of the variable τ = t + iw, (t, w ∈ R), and consider in this plane the following domains (see also [12]). For ξ 0 0 we define: I = {τ ∈ C | Im τ < ±ξ 0 }. Correspondingly, we introduce the cut-domains: Figure 1a) and Figure 1b). Moreover, we will use the notationȦ = A \ 2πZ to denote every subset A of C which is invariant under the translation group 2πZ.
We can then prove the following theorem (see also [8]). Theorem 2. Consider the following series: and suppose that the set of numbers {a } ∞ =0 , a = ( (3) The jump functionF(w) describing the discontinuity of if (τ) across the cutsΞ (0) + ), is a function of class C (p−1) and satisfies the bound: wheref (σ + iν), (ν ∈ R), is the unique Carlsonian interpolant of the coefficients f , and (4)f (σ + iν) is the Laplace transform of the jump functionF(w): The following Plancherel equality holds true: Proof. Since the set {a } ∞ =0 satisfies condition (5), then, given an arbitrary non-negative number C, there exists an integer m such that for > m, | f | ≤ C. Therefore, we can write: the series on the r.h.s. of (32) being uniformly convergent for w w 0 < 0. On the other hand, series (27) can be rewritten as where T m (τ) is a trigonometric polynomial. Then, for the Weierstrass theorem on the uniformly convergent series of analytic functions, we can conclude that the series (27) converges uniformly on any compact subdomains of I Then, given an arbitrary constant C , there exists an integer 1 such that for w = 0 we have: Then, using once again the Weierstrass theorem on the uniformly convergent series of continuous functions, the series 1 2π ∑ ∞ =0 f e −i t converges to a continuous functionf (t). Hence, the first statement is proved. To prove the other statements, we introduce the following integral:f wheref (λ), (Re λ − 1 2 ), is the unique Carlsonian interpolant of the sequence { f } ∞ 0 , which exists in view of the fact that the set {a } ∞ 0 satisfies condition (5). The contour C is contained in the half-plane C (+) −1/2 , encircles the real semi-axis with Re λ > − 1 2 (or a part of it), and cross this axis at a point σ > − 1 2 , σ ∈ N 0 , as illustrated in Figure 2. Let us now consider the following inequalities, which hold true for any σ, ν ∈ R: Let us recall that λ pf (λ) (p 1) tends uniformly to zero as λ → ∞ inside any fixed half-plane Re λ δ > −1/2 and, moreover, λ (p−1)f (λ) belongs to L 1 (−∞, +∞) for Re λ −1/2. For these properties and by the use of bound (38), the integralf + (t), (t ∈ [0, 2π]), exists and the contour C can be deformed and replaced by the line Figure 2). Hence, by applying the Watson resummation method [13], we obtain for t ∈ [0, 2π]: Proceeding analogously for the integralf − (t) (Formulae (35) with η = −), we distort similarly the contour of integration and finally obtain for t ∈ [−2π, 0]: (39) we substitute for t the complex variable τ = t + iw, and see that the resulting integral provides an analytic continuation off + (t) in the strip {τ = t + iw, 0 < t < 2π, w ∈ R + }, continuous in the closure of the latter. We have: then, using bound (38), and since λ (p−1)f (λ) ∈ L 1 (−∞, +∞) for any fixed value of Re λ = σ −1/2 (see Proposition 1), the statement is proved. Similarly, the analytic continuation of the functionf − is defined in the strip {τ = t + iw, −2π < t < 0}. We thus have proved that the functionf (τ), (τ = t + iw), admits a holomorphic extension to the cut- across the cuts can be computed by replacing t with (iw) in integrals (39) and (40). We obtain From (43) we derive the following bound for the jump functionF(w): where f σ 1 , given in (29), is guaranteed to be finite by statement (iv) of Proposition 1. Using Formulae (43) and recalling the Riemann-Lebesgue theorem, we prove thatF(w) is a function of class C (p−1) , (p 1), since λ (p−1)f (λ) ∈ L 1 (−∞, +∞) for any Re λ −1/2 (see Proposition 1). Inverting Formulae (43), we have: which is the Laplace transform of the jump functionF(w), holomorphic for Re λ > −1/2. Finally, recalling thatf (σ + iν) belongs to L 2 (−∞, +∞) for any fixed value σ −1/2, we obtain the Plancherel equality (31), which for σ = − 1 2 reads: We can now state the following proposition (see also [8]).

Proposition 4. If in the trigonometrical series
the coefficients f satisfy the assumptions of Theorem 2, then: (1) the series converges to a continuous functionf (t), (t ∈ R), and the convergence is uniform on any compact subdomain of the real line; (2) the functionf (t) admits a holomorphic extension to the cut-domainİ Proof. Statement (1) follows by observing that, in view of the hypothesis on the coefficients f , we have: Next, in view of the Weierstrass theorem on the uniformly convergent series of continuous functions, the result follows. Statements (2) and (3) can be proved in a way completely analogous to the one followed for proving the corresponding statements of Theorem 2.
Suppose now that the partial-waves f satisfy the following bound, which is precisely the M.T. condition [2]: From inequality (49) it follows that the series (47) converges to a functionf (τ), analytic in the strip | Im τ| < ξ 0 , the convergence being uniform in any compact subdomain of this strip. With slight modifications of Theorem 2, it is also straightforward to prove thatf (τ) admits a holomorphic extension to the cut-domainİ Figure 3). The physical applications could also lead to assume that only the subset {a } ∞ = 0 , ( 0 > 0), satisfies condition (5). In this case it is easy to prove thatf (λ) is analytic in the half-plane Re λ > 0 − 1 2 , continuous at Re λ = 0 − 1 2 ; furthermore, bound (28) holds for σ 0 − 1 2 (instead of σ − 1 2 ). Analogous modifications must be introduced in connection with Formulae (30) and (31).

Holomorphic Extension Associated with the Legendre Series
In Appendix B we consider a real one-sheteed hyperboloid, denoted by X, and the associated meridian hyperbola, denoted byX + . Next, we introduce through a Radon transformation a fibration with basisX + . This Radon transformation allows us to pass from polar coordinates (θ, φ) to horocyclic coordinates (τ, ζ) (see Appendix B). In this transformation, the horocycles play the same role played by the planes in the ordinary Radon transform. Then, we replace the complex cos θ-plane (see the Introduction) by a complex hyperboloid X (c) , whose submanifolds are the real one-sheteed hyperboloid X and the Euclidean sphere S 2 = (iR × R 2 ) ∩ X (c) . Next, we extend the fibration with basisX + to the complex hyperboloid X (c) and to the related complex meridian hyperbolaX (c) . This fibration is realized by a Radon transform which allows us to pass from the complex polar coordinates (θ, φ), (θ = u + iv), to the complex horocyclic coordinates (τ, ζ), (τ = t + iw). Finally, in Appendix C, we study the inversion of the Radon transformation. By using this latter, we move back from the τ-plane to the θ-plane (see, in particular, representation (A53), which holds if the sequence a = ( + 1) p f satisfies the Hausdorff condition (5) with p 2; see Appendix C). Then, the holomorphic extension in the τ-plane, proved in the previous section, can now be re-formulated in the θ-plane. More precisely, we can state the following propositions.

Proof.
From the assumptions on f we can state thatf (τ) is a 2π-periodic function holomorphic inİ which follows from equality (A17) for the uniqueness of the analytic continuation (see Appendix B, Formulae (A41)). These properties imply thef (τ) is of the following form: τ }. Now, we adopt the following parametrization for γ θ : cos τ = 1 + α(cos θ − 1), 0 α 1, the r.h.s. of Formulae (50) can be rewritten as that represents an even, 2π-periodic function, holomorphic inİ (0) θ . Since we shall prove below that this function can be represented by the Fourier-Legendre expansion (1), it can be properly denoted by f (θ).
From formulae (50) we can compute the boundary values f ± (v), defined as: f η (v) = lim u→0 + f (ηu + iv), η = ±, on the semiaxis {θ = iv, v 0} in terms of the corresponding boundary valuesf ± (w), with γ iv = {τ = iw, 0 w v}, provided thatf ± (w) satisfies a C 1 -type regularity condition. This latter conditon is necessary in order to perform the inversion of the Radon-Abel transform at the boundary. The C 1 -continuity of the boundary values is indeed guaranteed by the assumption that the sequence a = ( + 1) p f satisfies the Hausdorff condition (5) with p 2. Thus, the jump function F(v) across the cuts in the θ-plane is related to the jump functionF(w) across the cuts in the τ-plane by the following Radon-Abel transformation: . We can then apply the inverse Radon-Abel transform operator, defined by Formulae (50), to the series on the r.h.s. of Formula (A19): and integrate term by term for the uniform convergence of this series, which, in turn, follows from the Hausdorff condition on the coefficients f . Next, it is convenient to introduce the following functions: which are connected with the Legendre polynomials P (cos u) as follows (see Formula (II.91) of Ref. [14,15]): Finally, recalling that a = −a − −1 ( ∈ Z), we re-obtain the original Legendre expansion (1): Remark 1. If the partial-waves f satisfy both conditions (5) and (49), then we can prove (following the arguments of Proposition 5) that the function f (θ), given in (50), is even, 2π-periodic and holomorphic inİ

Inversion of the Laplace-Abel Transformation by the Means of the Pollaczek Polynomials
From Formulae (30) and (A42) the following equality follows: By exchanging the order of integration on the r.h.s. of (58), we obtain for Re λ > − 1 2 : Next, recalling the integral representation of the Legendre function of the second kind Q λ (cosh v): Formula (59) can be rewritten in the following form: Next, we recall the following equality [16]: (where we have adopted a non-standard normalization of Q λ , which is appropriate to our joint consideration of P λ and Q λ , the discrepancy with the usual notation being merely a factor 1/π). Hence, we have the following equality which we will use later: Now, we can prove the following proposition.

Lemma 1.
Suppose that the sequence a = ( + 1) p f , ( ∈ N 0 ), satisfies the Hausdorff condition (5) with p 0; then, the functionF(w)e w/2 can be represented by the following expansion which converges in the L 2 -norm: where and P k (·) and L k (·) being the Pollaczeck and Laguerre polynomials, respectively.
Since the coefficients {a } ∞ =0 satisfy condition (5), then there exists a unique Carlsonian interpolant of the sequence { f } ∞ 0 , denotedf (λ), such thatf (−1/2 + iν) belongs to L 2 (−∞, +∞) (see Proposition 1). Therefore, the functionf (−1/2 + iν) can be expanded in terms of Pollaczek functions as follows: the convergence of this expansion being in the L 2 -norm. The coefficients d k are given by: Exploiting the asymptotic behaviour of the Euler gamma function, integral (79) can be evaluated by the contour integration method along the path shown in Figure 4. Note that the poles of the gamma function Γ(1/2 − iν) are located at ν = −i( + 1/2), and allow us to pick up correctly the input data. In fact, we have: [f (−1/2 + iν)] (ν=−i( +1/2)) =f ( ) = f . We obtain: Next, we observe that: where F denotes the Fourier integral operator. Let us note that the function exp(− exp(−w)) exp(−w/2) belongs to the Schwartz space S of the C ∞ (R) functions that decrease with all their derivatives for |w| tending to +∞, faster than any negative power of |w|. Therefore we can write: Inserting representation (82) of the Pollaczek functions into expansion (78), we obtain: Now, we apply the inverse operator F −1 to the r.h.s. of Formula (83) and exchange the integral operator F −1 with the sum, this operation being legitimate within the L 2 -norm convergence. We obtain: Finally, recalling Formula (74), we obtain the following expansion for the function F(w) exp(w/2):F the convergence being in the L 2 -norm. It can be easily verified that: where L k denotes the Laguerre polynomials [19]. It can be readily checked that the polynomials L k (w) are a set of polynomials orthonormal on the real line with weight function W(w) = exp(−w) exp(−2 exp(−w)). Consequently, the set of functions Φ k (w), defined by (73), represents an orthonormal basis in L 2 (−∞, +∞).
Finally, from Formula (85) we obtain: where c k = d k / √ 2π, and the functions Φ k (w) are given by (73).
The successive step of the inversion procedure consists in inverting the Abel transform. Consider Formula (A42) and put: cosh v = y, cosh w = x. The Abel transform can be written as a convolution product of the following form: From (53), we analogously obtain: We can now prove the following Lemma (see also [8]).

Lemma 2.
If the sequence a = ( + 1) p f ( ∈ N 0 ) satisfies the Hausdorff condition (5) with p 2, then the following limit holds: where and φ ∈ S ∞ (R) (S ∞ (R) is the Schwartz space of the C ∞ (R) functions φ(x) that, along with all their derivatives, decrease for |x| → ∞ faster than any negative power of |x|), and f , φ denotes the Lebesgue integral Proof. From (89) and (91) we have: From inequality (28) and Formula (A42) it follows that AF has a power-like behaviour in x. Next, for the Fubini theorem, we have: where Next, from Schwarz's inequality it follows: Remarks 1. (i) As mentioned at the end of Section 3, if only a subset of partial-waves { f l } ∞ 0 satisfy condition (49), thenf (λ) is analytic only in the half-plane Re λ > 0 − 1 2 . However, even in this case, the above inversion procedure can be applied with only minor technical modifications. (ii) So far, we have assumed that all the partial-waves f are known and, in addition, that are noiseless. A much more realistic analysis accounts that only a finite number of partial-waves can actually be known and that they are inevitably affected by noise. In this case, the reconstruction ofF(w) from this set of noisy data meets the difficulties related to the ill-posedness of the inverse problems. See Ref. [5] for a more detailed discussion of the algorithm of reconstruction from a finite set of noisy input data.

Born Approximation as a Laplace-Abel Transformation
Let us now consider the Born approximation. The partial-waves f , at the first-order Born approximation, are given (suitablly normalized) by where j (kr) are the spherical Bessel functions of order . Then, exploiting the following formula: +∞ 0 e −µr J 2 (kr) dr = 1 recalling that j (kr) = √ π/2krJ +1/2 (kr) and observing that V µ 0 (r) belongs to the Yukawian class (see representation (2)), we obtain: We put µ = √ t = k[2(cosh v − 1)] 1/2 and since on the cut we have θ equals iv, then Formula (98) can be rewritten as Therefore, we can see that the restriction off (λ) to ∈ N 0 coincides with the firstorder Born approximation of the partial-waves, after the following identification (see also [2]): Formula (100) allows us to recover the spectral density σ(µ) from the jump function F, but only in the domain µ 2 0 < t < 4µ 2 0 , that is, where F coincides with the jump function associated with the first-order Born approximation (see [2]). In fact, we can write an equation of the Lippman-Schwinger-type for the discontinuities ∆ f .
where (∆ f ) 0 represents the discontinuity function associated with the first-order Born approximation. The convolution product has as cycle of integration the compact oneshaped region determined by the support properties of the involved jump function [20]. Equation (101) has been investigated by many authors (see, e.g., [2] and the references quoted therein). Here, we limit ourselves to consider the results which are worth in our case. Equation (101) can be formally solved through the following Neumann-Born series: where These convolution products form an algebra, named Volterra algebra, which has found particular relevance in the study of the Bethe-Salpeter equation (see [4,20]). The domains on which ∆ f and (∆ f ) 0 are defined permit us to implement an iterative procedure that, in principle, can reconstruct the spectral density σ(µ) completely (see [2]). In fact, if the support of (∆ f ) 0 is given by the values t > µ 2 0 , then the support of (∆ f ) 0 (∆ f ) 0 will be t > 4µ 2 0 . Therefore, within the interval 4µ 2 0 < t < 9µ 2 0 , the spectral density σ(µ) is This procedure can be iterated: for instance, for values of t in 9µ 2 0 < t < 16µ 2 0 it involves a double convolution product. We have earlier shown that ∆ f and (∆ f ) 0 enjoy a power-like behavior in µ, but the entire iterative procedure is not guaranteed to maintain this type of behavior. Therefore, the existence of V µ 0 (r) cannot be proved rigorously. In other words, the behavior in µ of each multiple convolution product can be controlled, but when the amplitudes are restricted to be on the energy-shell, all the possible combination of the factors must be taken into consideration. This leads to an alternating series, whose global µ-behavior cannot be controlled.
Summarizing, we can prove the existence of a spectral density σ(µ) which possesses a good local behavior and, in addition, we can reconstruct this function from the amplitudes by using an iterative procedure. Unfortunately, to date we are not able to keep under control the global µ-behavior of the entire iterative procedure so to have guaranteed the convergence of integral (2) and, consequently, the existence of the potential V µ 0 (r).

Conclusions
We have shown that if the partial-waves f of the Fourier-Legendre expansion of the scattering amplitude satisfy a suitable Hausdorff condition and are exponentially bounded, then this expansion converges to an analytic function f (θ), which can be holomorphically extended to a slit domain in the cos θ-plane. The geometrical tool adopted to reach this result consists essentially in connecting the analysis on the complex cos θ-plane with that on a complex hyperboloid X (c) , which contains as submanifolds the support of the SO(3, R) harmonic analysis and the support of the cut. Through a suitable fibration, obtained by a Radon-type transformation, we can limit ourselves to study the harmonic analysis on a one dimensional complex hyperbola, thus allowing us to regard the Fourier-Legendre series as a trigonometrical series. The holomorphic extension of the latter can then be translated back to the slit cos θ-plane by inverting the Radon-type transform. The jump function describing the discontinuity on the cut is found to be related to the Laplace-Radon tranform of the analytic function interpolating uniquely the partial-waves f . The inversion of this transformation finally yields the spectral density σ(µ) and, under some additional restrictions, the potential V µ 0 (r). The key ingredient to achieve this result is the Haudorff condition (5) we imposed on the partial-waves. Such a condition, which firstly emerged in the theory of moments, calls for the analysis of the broader class of completely monotonic potentials which, in turn, find application in various problems, e.g., of extremal geometry. These issues will be the object of future investigations. Proof. Let us start from the integral representation (A1), and set the variable γ as: The following equalities are easily verified: Let us now regard e iγ as the integration variable; the integration path in integral (A1), corresponding to η ∈ [0, π], is the (oriented) linear segment δ 0 (u) starting at e iu and ending at e −iu (see Figure A1). As shown by equalities (A8) and (A9), the integrand is an analytic function of e iγ in the disk |e iγ | < 1 (since ∈ N). Then, the integration path δ 0 (u) can be replaced by the circular path δ + (u) = {e iγ : γ = t, u t 2π − u}. Moreover, by using the fact that (e iγ − e iu )(e iγ − e −iu ) −1/2 is positive for e iγ ∈ δ 0 (u) ∪ R, and therefore at e iγ = e iπ , we conclude from the left equality in (A6) that in the r.h.s. of Formula (A8) there holds the following specification (for γ = t, u t 2π − u): which gives wheref (t) is given by (A4).