\documentclass[twoside]{article}
\usepackage{amssymb} % font used for R in Real numbers
\pagestyle{myheadings}
\markboth{\hfil Riemann-Lebesgue properties of Green's functions 
 \hfil EJDE--2000/07}
{EJDE--2000/07\hfil Richard Ford \hfil}
\begin{document}
\title{\vspace{-1in}\parbox{\linewidth}{\footnotesize\noindent
{\sc  Electronic Journal of Differential Equations},
Vol.~{\bf 2000}(2000), No.~07, pp.~1--19. \newline
ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
\newline ftp  ejde.math.swt.edu \quad ftp ejde.math.unt.edu (login: ftp)}
 \vspace{\bigskipamount} \\
%
 Riemann-Lebesgue properties of Green's  functions with applications 
to inverse scattering  
\thanks{ {\em 1991 Mathematics Subject Classifications:} 
35J10, 35P25, 35R30, 81U40.
\hfil\break\indent
{\em Key words and phrases:} inverse scattering, Riemann-Lebesgue properties.
\hfil\break\indent
\copyright 2000 Southwest Texas State University  and University of
North Texas. \hfil\break\indent
Submitted August 24, 1999. Published January 21, 2000.} }
\date{}
%
\author{Richard Ford}
\maketitle

\begin{abstract}  
 Sait\={o}'s method has been applied successfully for measuring 
 potentials with compact support in three dimensions. Also potentials
 have been reconstructed in the sense of distributions using a 
 weak version of the method.
 Sait\={o}'s method does not depend on the decay of the boundary value
 of the resolvent operator, but instead on certain Reimann-Lebesgue type 
 properties of convolutions of the kernel of the unperturbed resolvent. 
 In this paper these properties are extended from three to higher 
 dimensions. We also provide an important application to inverse 
 scattering by extending reconstruction results to measure potentials 
 with unbounded support. 
\end{abstract} 

\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\renewcommand{\theequation}{{\arabic{section}}.{\arabic{equation}}}

\section{Introduction}

The field of inverse potential scattering has matured over the 
last decade resulting in the emergence of a variety of effective 
techniques.  Early work by Faddeev \cite{faddeev} uses a high 
energy limit of the scattering amplitude which he shows converges 
to the Fourier Transform of the potential for a certain class of 
potentials.   A different method that also utilizes the high 
energy data from the scattering amplitude has been developed by 
Sait\={o} and applied to general short-range potentials in 
${\mathbb R}^3$  \cite{saito1}, \cite{saito2} and subsequently to 
${\mathbb R}^n$, $n\geq 2$ \cite{saito3}.  Newton's method 
\cite{newtonbook} exploits the reciprocity relations to reduce the 
inverse problem in ${\mathbb R}^3$ to the Marchenko equation, 
which is subsequently shown to be uniquely solvable under a variety 
of conditions on the potential.   Newton's method has subsequently 
been generalized from $R^3$ to $R^2$ by Cheney \cite{cheney} and 
to ${\mathbb R}^n$ by Weder \cite{weder}.    Other authors have 
been successful in generalizing the class of admissible potentials 
to include certain types of singularities (see P\"{a}iv\"{a}rinta, 
Serov, and Somersalo \cite{serov2}, Serov \cite{serov1}, 
\cite{serov3}). Our interest lies specifically in Schr\"{o}dinger 
operators associated with measure potentials.   Delta potentials 
and other highly singular measure potentials have captured the 
interest of many authors, (e.g.  \cite{albeverio}, \cite{brasche},  
\cite{shabani} and the references therein).  The direct scattering 
problem has been solved in ${\mathbb R}^n$ for a broad class of 
measure potentials \cite{ford1}, \cite{ford3}, \cite{ford5}.    
Inverse scattering for measure potentials is not so well developed 
and few results are currently available.   One approach that has 
produced results is that of Sait\={o}, but modified to accommodate 
reconstruction of the measure potential in the sense of 
distributions \cite{ford6}.  The results only apply in 3 
dimensions and to potentials with bounded support. Most of the 
inversion methods we have mentioned above depend on the decay 
properties of the boundary value of the resolvent operator in 
operator norm when viewed as an operator from weighted to 
unweighted $L^p$ spaces (usually $p=2$).   In stark contrast, the 
modified Sait\={o} method relies on certain Riemann-Lebesgue like 
properties of convolutions of the kernel of the unperturbed 
resolvent operator.   Essentially,  the vibrations of the 
convolutions annihilate the Born remainder shedding additional 
light on the importance of the Born approximation.   It has long 
been known that the Born approximation to the scattering amplitude 
carries all of the essential information for reconstruction for a 
wide variety of potentials.  The kernel convolutions and their 
properties provide evidence that this variety may be substantially 
wider than results to date indicate.  This paper shall address 
extending the results of \cite{ford6} regarding the kernel 
convolutions from 3 to $n \geq 3$ dimensions.  Applications to the 
inverse scattering problem shall be illustrated by obtaining new 
reconstruction results for measure potentials with unbounded 
support. An explicit example is provided showing that even highly 
singular delta potentials can be recovered through the modified 
Sait\={o}'s method.  

\section{Results} \setcounter{equation}{0}

  Let $H$ be the standard self-adjoint realization of the 
Laplacian on ${\mathbb R}^n$ and let $R(z)$ denote the associated 
resolvent operator, $(z-H)^{-1}$.  It is well known that $R(z)$ is 
an integral operator with kernel given by  $G(x,y; \kappa) = G(|x- 
y|,\kappa)$   where $ G(r,\kappa) = \frac{i}{4} (\frac{\kappa}{2 
\pi r})^{\frac{n-1}{2}} H^{1}_{\frac{n-1}{2}}(r \kappa )$, 
$\kappa^2 = z$, $\mathop{\rm im}\kappa > 0$ and $H^{1}_{\nu}$ is 
the Hankel function of the first kind.  Now let 
\begin{equation} J(r)  = \int_{S^{n-1}} e^{i r \omega \cdot \omega' 
}\,d\omega'\,. \label{introeq1}
 \end{equation} 
By Alsholm and Schmidt \cite{alsholm} we have that 
$ \frac{4\pi}{i} \left( \frac{2\pi}{k} \right)^{n-2} 
(G(x,k) - G(x,- k)) = J(k|x|)  $ 
where $k >0$. The primary result of this work is the following.

\begin{theorem} Let $\psi(x) \in C_{0}^{\infty}({\mathbb R}^{n})$.  Then 
for all $x \neq y$ we have
\begin{equation}   \lim_{k \to \infty} 
k^{n-1}\int_{{\mathbb R}^{n}}J(k|x-\xi|) J(k|\xi - y|)  
\psi(\xi)\,d\xi = 0 \,.
\end{equation}
\label{maintheorem}
\end{theorem}

 The impact of this theorem on the inverse scattering problem 
lies in 
the following observations.  In the case of ordinary Schr\"{o}dinger 
Operator scattering with short range potentials and valid generalized 
eigenfunctions, $\phi_{\pm}(x,k\omega)$, the scattering matrix, $S(k)$, 
is given by
\begin{equation} [S(k)f](\omega) = f(\omega) - \int \int 
e^{ik\omega\cdot x} 
\phi_{+}(x,k\omega') f(\omega')V(x)\,dx\,d\omega'\,.
\end{equation} 
When one views the the potential, $|V(x)|$ as a weight 
defining a Hilbert space, ${\bf K} = L^{2}({\mathbb R}^n,|V(x)|)$, 
this operator takes the form, 
$$S(k) = 1- \gamma(k) (1-Q^{+}(k^2))^{-1} \gamma(k)^{*} \,,$$
where $\gamma(k) $ is a mapping from ${\bf K}$ to $L^{2}(S^{n-1})$ 
where $S^{n-1}$ denotes the unit sphere in ${\mathbb R}^n$.  This mapping is 
sometimes referred to as a {\em trace} operator and is given by 
\begin{equation} 
\gamma(k)f(\omega) = \int_{{\mathbb R}^n} e^{-ik\omega\cdot x} 
f(x) V(x)\,dx \label{gammaeq} \end{equation} 
and $Q^{+}(k^2)$ is the 
boundary value of the modified unperturbed resolvent operator mapping 
${\bf K}$ to ${\bf K}$ given by
\begin{equation} 
Q^{+}(k^2)f(x) = \int_{{\mathbb R}^n} G(x-y;k) f(y) V(y)\,dy 
\label{qeq} \end{equation} Under various limiting absorbtion principles 
the perturbed resolvent admits a boundary value of its associated Green's 
function, $G_{1}^{+}(x,y:k)$.  The associated modified resolvent 
operator is then given by
\begin{equation} Q_{1}^{+}(k^2)f(x) = \int_{{\mathbb R}^n} G_{1}^{+}(x,y;k) f(y) 
V(y)\,dy. \label{qmodeq} \end{equation}  
The Marchenko-Newton, Sait\={o}, and other inverse scattering methods 
essentially require that this modified perturbed resolvent vanish in 
some sense as $k \to \infty$.  
The Riemann-Lebesgue properties established in Theorem 
\ref{maintheorem} can weakly annihilate the Born remainder terms in 
the inner product between the scattered and unscattered plane waves. 
This remainder is given in weak formulation by
\begin{equation} -\pi\left(\frac{k}{2\pi}\right)^{n-1} \int_{{\mathbb R}^n} \, 
\langle [\gamma(k)Q_{1}^{+}(k^2)\gamma(k)^{*}e^{-ik(\cdot)\cdot 
x}](\omega), e^{-ik\omega\cdot x}\rangle _{\omega} \psi(x)\,dx \,,
\end{equation}
with $\psi(x)$ in  $C_{0}^{\infty}({\mathbb R}^n)$, and  
where the integration $\langle\cdot,\cdot\rangle_{\omega}$ is taken over 
the surface of the unit sphere.  As we will show in our application to 
measure potentials, this remainder will reduce to  
\begin{equation} 
\mbox{const.} \times k^{n-1} \int\int G_{1}(x,y;k)J(k|\xi-
y|)J(k|x-y|) \psi(x)V(d\xi)V(dy)\,dx 
\end{equation} 
which will vanish as $k \to \infty$ under appropriate 
conditions on $G_{1}$ due to the Riemann-Lebesgue properties of $J$ (which is 
independent of the potential) 
rather than the decay of modified {\em perturbed} resolvent.  The potential can 
therefore be reconstructed weakly through the knowledge of the high energy 
scattering data with the Sait\={o} method.  We will first carry out the details 
of this process in the next section.  The proof of Theorem 
\ref{maintheorem} will then be carried out in section 4.

\section{Application to Inverse Scattering} \setcounter{equation}{0}

 Theorem \ref{maintheorem} shall now be applied to the inverse 
scattering problem for a broad class of measure potentials. This section 
will conclude with an explicit example, the delta function on a sphere. 
 Let us consider the Schr\"{o}dinger equation with real potential in 
${\mathbb R}^n \, (n>2)$.  Throughout this section our 
measure potentials, $V(dx)$, will satisfy the following:

\paragraph{Assumption 1.}
\begin{eqnarray}
&\sup_{y>M} \int_{|x-y| > 1 } |x-y|^{\frac{1-n}{2}} 
|V|(\!dx\!) \to 0 \quad\mbox{as }M \to \infty\,,&
 \label{p8app1} \label{ass1eq1} \\
&\sup_{y>M} \int_{|x-y| \leq 1 } |x-y|^{2-n} |V|(\!dx\!) 
\to 0 \quad \mbox{as }M \to \infty\,, &\label{p8app2} 
\label{ass1eq2} \\
&\sup_{y} \int_{|x-y|<\delta} |x-y|^{2-n} |V|(\!dx\!) 
\to 0 \quad\mbox{as } \delta \to 0\,.&
\label{ass1eq3} \end{eqnarray}

 Our conditions in Assumption 1 will allow for the possibility of 
highly singular measure potentials, $V(dx)$, such as the delta 
function on a manifold. We denote by ${\bf  K}$ the space of square
integrable functions under the measure, $|V|(\!dx\!)$ with norm and 
inner product given  respectively by 
 \begin{equation} \| f \|_{K} = \left\{\int 
|f(x)|^{2}|V|(\!dx\!)\right\}^{1/2},\quad\langle f,g\rangle _{K}=
\int f(x) \overline{g(x)} |V|(\!dx\!)\,. \label{results2a} 
\end{equation} 

 Let ${\cal S}$ be the Schwartz space of rapidly decreasing 
functions on ${\mathbb R}^n$ and ${\cal S}'$ the corresponding set 
of tempered distributions.   Let ${\cal K}$ be the space of 
$C_{0}^{\infty}({\bf R}^n)$-functions with the usual ``test 
function" topology, i.e.   $f_{k}\to 0$ in ${\cal K}$ if and only 
if some bounded set contains all of their supports and the 
$\infty$-norm of $f_{k}$ and all of its derivatives vanish.  Let 
${\cal K}'$ be the set of continuous linear functionals on ${\cal 
K}$. We shall define $\mathop{\rm A}$ as the identification operator, 
$\mathop{A}f(x)=f(x)$, from ${\cal K}$ to ${\bf K}$. It is shown in 
\cite{ford1} that under Assumption 1,  $AR(z)$ extends to a 
bounded operator from $L^{2}({\mathbb R}^{n})$ to ${\bf K}$.  An 
important operator on the space, ${\bf  K}$ is the modified 
resolvent given by $Q(z) = A[AR(\overline{z})]^{*}$ so that for $z 
= \kappa^2$, $\mathop{\rm im}\kappa > 0$, 
\begin{equation} 
Q(z)f(x) = \int  G(x,y;\kappa) f(y) V(\!dy\!)\,.
\label{Qdef} \end{equation}  
It is also shown in \cite{ford1} that 
there is a self-adjoint operator $H_{1}$ satisfying 
\begin{equation} (H_{1}u,v) = (u,Hv) + \int u(x)\overline{v(x)} 
V(\!dx\!)\quad u\in D(H_{1}),\,\,v\in D(H)\,.  \label{opequat} 
\end{equation}  
We denote by $R_{1}(z)$ the associated resolvent 
operator, $[z-H_{1}]^{-1}$ and let $Q_{1}(z) = 
A[AR_{1}(\overline{z})]^{*}$. Under Assumption~1 the wave 
operators 
$$W_{\pm}(H_{1},H) = \mathop{\rm s-lim}_{t\to\pm{\infty}} 
e^{itH_{1}}e^{-itH}$$ 
exist and are strongly complete and the scattering operator, $S = 
W_{+}^{*}W_{-}$ exists and is unitary. It is shown in \cite{ford5} 
that there exists a family of unitary operators, $\{S(k)\}_{k>0}$ 
on $L^{2}(S^{n-1})$ such that 
        \begin{equation} [S(k)f(k,\cdot)](\omega) = \hat{S}f(k\omega) 
\quad f\in {\cal K} \label{results4} \end{equation}
where $\hat{S} = {\cal F}S{\cal F}^{*}$, ${\cal F}$ is the unitary 
Fourier 
Transform and $k\omega$ are polar coordinates in ${\mathbb R}^n$.  Modified 
{\em trace} operators, $\gamma(k)$ and $\gamma_{\pm}(k)$ are defined 
by:
        \begin{eqnarray}
\gamma(k)f(\omega) & = & \int e^{-ik\omega\cdot 
x} f(x) V(\!dx\!) \label{results5} \\ 
\gamma_{\pm}(k)f(\omega) & = &  \int \phi_{\pm}(x,k\omega) f(x) 
V(\!dx\!)\,, \label{results6} \end{eqnarray}
 where $\phi_{\pm}$ are 
the generalized eigenfunctions associated with $H_{1}$.  It is 
known \cite{ford5} that these are bounded operators from ${\bf K}$ 
to $L^{2}(S^{n-1})$ under the conditions in Assumption~1.  
Furthermore, we have the following representation of the 
scattering matrix.
        \begin{equation} 1-S(k) = 
\frac{i}{4\pi}\left(\frac{k}{2\pi}\right)^{n-
2} \gamma_{+}(k)\gamma(k)^{*}. \label{results7} \end{equation} 
We will denote by $V$ the linear functional associated with 
$V(\!dx\!)$,
        \begin{equation} (V,u) = \int u(x)  V(\!dx\!), \quad \, u \in 
{\cal S}.\label{results10} \end{equation}  
We denote by  $\nu $ the linear functional defined by ,
\begin{equation}  (\nu, f) = \int \int |x-y|^{1-n} f(x) V(dy)\,dx, 
\quad u \in {\cal S}.
\label{p8app4} \end{equation} 
Finally, we define the function, $h(k,x)$ on 
$({\mathbb R}^{+} \times {\mathbb R}^n)  $ by 
        \begin{equation} h(k,x) = 2\pi i k \langle (1-S(k))e^{-ik\omega\cdot 
x},e^{-ik\omega\cdot x}\rangle _{\omega} \label{results8} \end{equation} 
where $\langle\cdot,\cdot\rangle_{\omega}$ is the inner product in $L^2(S^{n-1})$. 
Noting that $h(k,x)$ is bounded for each fixed $k$, we consider $h(k)$ 
as the associated distribution in ${\cal S}'$ given by 
        \begin{equation} (h(k),u) = \int \,h(k,x) \overline{u(x)} \,dx  
\label{results9} \end{equation} 

\paragraph{Assumption 2A.} There exists 
some positive $k_{0}$ such that $Q_{1}(k^2+ i\epsilon)$ admits a 
boundary value, $Q^{+}_{1}(k^2) = \lim_{\epsilon\searrow 
0}Q_{1}(k^2+ i\epsilon)$ with integral kernel, $G^{+}_{1}(x,y;k)$ 
for all $k > k_{0}$.  Furthermore, there exists a function 
$F(x,y)$ on ${\bf R}^{2n}$ with $F(x,y) > |G^{+}_{1}(x,y;k)|$ for 
all $k > k_{0}$ and $\int \int F(x,y) |V|(\!dx\!) |V|(\!dy\!) < 
\infty$.  $G_{1}^{+}(x,y;k)$ is the Green's function and is known 
to exist under a variety of conditions and is described as the 
response at $x$ to a point source at $y$ (see \cite{cheney2} and 
the references provided).  As an alternative to Assumption 2A, we 
can offer 

\paragraph{Assumption 2B.} There exists a positive number,  $\alpha < 1$, 
some $k_{0}$, and a function, $F(x,y)$ on ${\mathbb R}^{2n}$  such that 
for all $k > k_{0}$,  $F(x,y) \geq |G(|x-y|;k)|$ and 
\begin{equation} \sup_{y} \int F(x,y) |V|(dx)\leq \alpha
\label{assumption2Beq} \end{equation} 
We define $\Omega_{V} = \int |V|(dx)$ and  assume it is finite.
We now provide;

\begin{theorem} 
If $V$ satisfies Assumption~1 above,  then the following holds:

\noindent{\bf 1.}  $V$ and $\nu$ as given above are both continuous linear functionals 
on 
${\cal S}$, hence in ${\cal S}'$.

\noindent{\bf 2.}   The linear functional, $\Lambda\nu$, defined on ${\cal K}$ by 
\begin{equation}(\Lambda \nu, \psi) = \left(\nu, {\cal F}^{*}(|\xi|{\cal F}\psi) 
\right) 
\label{Lambdanueq} \end{equation}
 extends to a well-defined element of ${\cal S}'$.

 If, in addition $V$ satisfies either of assumptions 2A or 2B 
then 
we also have

\noindent{\bf 3.}  $lim_{k\to\infty} h(k)  = -2\pi \nu$  in ${\cal S}'$ and 

\noindent{\bf 4.}  $V$ can be recovered in the sense of distributions through:
\begin{equation} \langle V,\psi\rangle  = (2\pi)^{-\frac{n}{2}} \alpha^{-
1} (\Lambda \nu, \psi)  \quad \psi \in {\cal S}\,, \label{item3} 
\end{equation} 
where $\alpha = 2^{\frac{2-n}{2}} \sqrt{\pi}\Gamma(\frac{n-1}{2})^{-1}$.
\label{p8appthm} 
\end{theorem}

\paragraph{Proof:} To show that $V$ is in ${\cal S}'$ we let 
$\psi_{k}(x)$ be a sequence of functions vanishing in ${\cal S}$.  We 
therefore have for any $\alpha>0$ a vanishing sequence $C^{\alpha}_{k}$ 
such that 
\begin{equation} |x^{\alpha} \psi_{k}(x)| \leq C^{\alpha}_{k} 
\quad \forall \,x\in{\mathbb R}^n \end{equation}
  We now have that
\begin{eqnarray} |(V, \psi_{k})| & \leq & \int_{|x| \geq 1} 
|\psi_{k}(x)| 
|V|(dx) + \int_{|x|<1} |\psi_{k}(x) |V|(dx) \nonumber \\
& \leq & C_{k}^{\alpha} \int_{|x|\geq 1} |x|^{-\alpha} |V|(dx) + 
C_{k}^{0} \int_{|x|<1} |V|(dx). \nonumber \end{eqnarray}
Taking $\alpha = (n-1)/2$ and applying (\ref{ass1eq1}) and 
(\ref{ass1eq2}) we see this vanishes as $k \to \infty$ showing 
that 
$V$ is in ${\cal S}'$.   Turning our attention to $\nu$ we write
\begin{eqnarray*} 
(\nu, \psi_{k}) & = &  \int \int _{|x-y|\leq 1} |x-y|^{1-n} 
\psi_{k}(x)\,dx V(dy) \\ 
&&+\int \int _{|x-y|> 1} |x-y|^{1-n} \psi_{k}(x)\,dx V(dy) \\ 
& = & I_{1}(k) + I_{2}(k) 
\end{eqnarray*}
 Since $\psi_{k}$ vanishes in ${\cal S}$ we 
also have for each $\alpha \geq 0$ vanishing sequences of constants, 
$c^{\alpha}_{k}$ such that 
\begin{equation} |\psi_{k}(x)| \leq c^{\alpha}_{k}(2+|x|)^{\alpha} 
\quad \forall x \in {\mathbb R}^n . 
\end{equation} 
Thus we have that
\begin{eqnarray} |I_{1}| & \leq  & \int  \left\{ \sup_{x: |x-y|\leq 1} 
c_{k}^{\alpha}(2+|x|)^{-\alpha} \right \} \int_{|x-y|\leq 1} |x-y|^{1-
n}\,dx 
|V|(dy) \nonumber \\ & \leq & c_{k}^{\alpha} \Omega^{n} \int (1+|y|)^{-
\alpha} |V|(dy)\,, \nonumber \end{eqnarray}
 where $\Omega^{n}$ is the surface area of $S^{n-1}$.  Again by 
(\ref{ass1eq1}) and (\ref{ass1eq2}) we see by taking $\alpha \geq 
\frac{n- 1}{2}$ that $I_{1}(k)$ vanishes as $k \to 
\infty$.  The verification that $I_{2}(k)$ vanishes is similar and 
follows from the observation that for $|x-y|>1$, we have 
$|x-y|^{1-n} \leq |x-y|^{(1-n)/2}$. This completes the proof 
of item 1 of Theorem \ref{p8appthm}.  Item 2 requires that we 
show that $\Lambda \nu \in {\cal S}'$.  Here we will borrow some 
ideas from Sait\={o} \cite{saito3}.  We define for $s\geq 0$ a 
norm $\|\cdot\|_{s}$ on ${\cal S}$ by 
\begin{equation} \| \psi \|_{s} = \sum_{|\beta|\leq n} 
\sum_{\alpha\leq\beta} \int |\xi|^{s-|\beta|+|\alpha|} \left| 
D^{\alpha}\psi(\xi) \right|\,d\xi \label{normdef} \end{equation} 
where $\alpha$ and $\beta$ are multi-indices.  As noted by 
Sait\={o} the topology induced by this norm is weaker than the 
proper topology on ${\cal S}$.  Furthermore,  it follows from the 
proof of lemma A.1 in \cite{saito3} that for any $\psi \in {\cal 
S}$ 
\begin{equation}|{\cal F}^{*}(|\xi| \psi)(x)| \leq C_{s} (2+|x|)^{-n} 
\| 
\psi \|_{s}  \end{equation}  where $C_{s}$ is a constant depending 
on 
$s$.  With $\psi_{k}$ as before, we have
\begin{equation} |(\Lambda\nu, \psi_{k})|  \leq C_{s} \| {\cal 
F}\psi_{k} \|_{s} \left| \int \! \int \, |x-y|^{1-n} (2+|x|)^{-n} 
V(dy)\,dx \right|
\end{equation}  
By the method of proof applied to $V$ and $\nu$ one can 
verify the integral remaining is finite, hence the right hand 
side will vanish as $k \to \infty$ proving item 2 of the 
theorem. \smallskip

 Moving on to item 3, we let $\psi(x) \in {\cal K}$. By the 
definition of $h(k,x)$ above and (\ref{results7}) we have that
\begin{equation}  (h(k), \psi)  =    -\pi  
\left(\frac{k}{2\pi}\right)^{n-1} \int  
 \langle  \gamma_{+}(k)\gamma(k)^{*} e^{-ik\omega\cdot x}, e^{-
ik\omega\cdot x}\rangle _{\omega} \overline{\psi(x)}\,dx\,. \label{heq1} 
\end{equation}
By Assumption 1 and the results of \cite{ford1} we know that $1-
Q(\kappa^2)$ has a bounded inverse given by $1+Q_{1}(\kappa^2)$ 
for im  $\kappa^2$ sufficiently large and where $Q_{1}(\kappa^2)$ 
is defined by, 
\begin{equation} Q_{1}(\kappa^2)f(x) = \int G_{1}(x,y;\kappa) f(y) 
V(dy)\,, 
\end{equation} where $G_{1}(x,y;\kappa)$ is the kernel of the perturbed 
resolvent operator.  By Assumption 2A, analytic continuation, and 
continuity we have that 
\begin{equation}  \left[1-Q^{+}(k^{2})\right]^{-1} 
=1+Q^{+}_{1}(k^{2}) \mbox{ for $k > k_{0}$ }. \end{equation}
  By \cite{ford5} (Eqn.~4.23) we have that $ \gamma_{+}(k) = 
\gamma(k)(1+Q^{+}_{1}(k^2))$.  Inserting this into (\ref{heq1}) we 
obtain, 
\begin{equation} (h(k), \psi(x)) = (h_{1}(k), \psi(x)) + (h_{2}(k), 
\psi(x))
\end{equation} where 
\begin{equation} (h_{1}(k), \psi(x)) = -\pi  
\left(\frac{k}{2\pi}\right)^{n-1} 
\int  \langle  \gamma(k)\gamma(k)^{*} e^{-ik\omega\cdot x}, e^{-
ik\omega\cdot x}\rangle _{\omega} \overline{\psi(x)}\,dx \end{equation} 
and
\begin{equation} (h_{2}(k), \psi(x)) = -\pi  
\left(\frac{k}{2\pi}\right)^{n-1} 
\int  \langle  \gamma(k) Q_{1}^{+}(k^2) \gamma(k)^{*} e^{-ik\omega\cdot 
x}, e^{-ik\omega\cdot x}\rangle _{\omega} \overline{\psi(x)}\,dx 
\end{equation}  
The first term, $h_{1}$, corresponds to the Born approximation while 
the second term, $h_{2}$ is the Born remainder.  We will look at $h_{2}$ 
first to  illustrate how Theorem \ref{maintheorem} applies.  By writing all 
of the integrals involved including the integral form of $Q_{1}(k)$,  and 
employing Assumption 2A and Fubini's theorem we arrive at 
\begin{equation} (h_{2}(k), \psi(x)) = C k^{n-1} \int G_{1}(x,y;k) 
J^{\psi}_{k}(x,y) V(dx) V(dy)\,, \end{equation} 
 where $C$ is a constant and 
\begin{equation} J^{\psi}_{k}(x,y) = \int J(k|x-\xi|)J(k|\xi-y|) 
\overline{\psi(\xi)}\,d\xi. \end{equation} 
By Assumption 2A we can apply the Lebesque dominated convergence theorem 
and Theorem \ref{maintheorem} to conclude this term vanishes as 
$k \to \infty$. 
If Assumption 2B holds instead of Assumption 2A, then the Neumann 
series for $[1-Q^{+}(k^2)]^{-1}$ will converge uniformly for 
$k>k_{0}$.  In this case, we have the expansion
\begin{eqnarray}
\lefteqn{ (h_{2}(k),\psi(x)) }  \label{h2ass2Beq} \\
&=& -\pi \left(\frac{k}{2\pi}\right)^{n-1}\sum_{N=1}^{\infty} \int 
\langle \gamma(k)Q^{+}(k)^{N}\gamma(k)^{*}e^{-ik(\cdot)\cdot 
x}(\omega), e^{-ik\omega\cdot x}\rangle _{\omega} 
\overline{\psi(x)}\,dx\,. \nonumber
\end{eqnarray}
For each fixed $N$, we can construct the iterated kernel, 
$Q^{N}(x,y,k)$ for $Q^{+}(k)^{N}$ given by
\begin{equation} Q^{0}(x,y,k) = G(|x-y|;k) ,\quad  Q^{N}(x,y;k) = \int 
Q^{N-1}(x,\xi,k) G(|\xi-y|;k) V(d\xi)\,.
\label{Qneq} \end{equation} 
 We next define $F^{N}(x,y)$ similarly by replacing 
$G$ with $F$ in (\ref{Qneq}). By induction we can see that 
$F^{N}(x,y) \geq |Q^{N}(x,y;k)| $ for all 
$k$ and furthermore,
\begin{equation}   \int F^{N}(x,y) |V|(dx)|V|(dy) \leq \Omega_{V} 
\alpha^{N+1}. \end{equation}  
We therefore have
\begin{eqnarray*}
\lefteqn{ \left|-\pi \left(\frac{k}{2\pi}\right)^{n-1} 
\int \langle \gamma(k)Q^{+}(k)^{N}\gamma(k)^{*}(k),e^{-
ik(\cdot)\cdot x}(\omega),e^{-ik\omega\cdot x}\rangle _{\omega} 
\overline{\psi(x)}\,dx \right| } \\ 
& = & \left|\pi \left(\frac{k}{2\pi} \right)^{n-1}  \int 
Q^{N}(x,y,k)J^{\psi}_{k}(x,y) V(dx) V(dy)\right| \\ 
& \leq & \pi \left(\frac{k}{2\pi} \right)^{n-1}  \int 
\left|F^{N}(x,y)J^{\psi}_{k}(x,y)\right| |V|(dx) |V|(dy)\,. 
\hspace{2.5cm}
\end{eqnarray*} 
The last integral will vanish as $k \to \infty$ by  dominated 
convergence and the fact that $ J^{\psi}_{k}(x,y) $ vanishes for all 
$x \neq y$ by Theorem \ref{maintheorem}.  
Furthermore, by the  uniform convergence of the Neumann series for 
$k>k_{0}$  we see  that the infinite sum yielding $h_{2}(k)$ will 
vanish as well. 
 We now turn our attention to the Born approximation, $h_{1}(k)$.  By 
again writing out the integrals involved and employing Fubini's 
theorem we obtain,
\begin{equation}  (h_{1}(k), \psi(x)) =-
\pi\left(\frac{k}{2\pi}\right)^{n- 1}\int \int 
\left|J(k|x-y|)\right|^{2} \overline{\psi(x)} V(dy)\,dx .  
\label{h1eq} \end{equation} 
Setting $\Delta(k)=J(k) - 2(\frac{k}{2 
\pi})^{\frac{1-n}{2}}\cos(k- \frac{(n- 1)\pi}{4})$ we have by 
\cite{saito3} equation (2.17) 
\begin{equation} 
 |\Delta(k)| \leq C \min\{k^{- \frac{(n-1)}{2}}, k^{- 
\frac{(n+1)}{2}}\}\,,   \label{Delta} 
\end{equation}  
where $C$ is a constant.  We can now write,
\begin{eqnarray*} 
\lefteqn{ - (h_{1}(k), \psi(x)) }\\ 
&=&\pi\left(\frac{k}{2\pi}\right)^{n-1}\int \! \int \big| 
2\left(\frac{k}{2 \pi}\right)^{\frac{1-n}{2}} \cos\left(k-\frac{(n-
1)\pi}{4}\right) + \Delta(k|x-y|) \big|^{2} \\
&&\hspace{3cm} \times \overline{\psi(x) }V(dy) \,dx \\ 
& = & \pi\left(\frac{k}{2\pi}\right)^{n-1}\int\!\int |\Delta(k|x-y|)|^2 
\overline{\psi(x) }V(dy)\,dx  \\ 
&  &  + 4\pi \left(\frac{k}{2 \pi}\right)^{\frac{n-1}{2}} \mathop{\rm Re} 
 \int\!\int \Delta(k|x-
y|) \cos\left(k-\frac{(n-1)\pi}{4}\right) \overline{\psi(x) }V(dy)\,dx \\ 
&  & + 4\pi \int\!\int |x-y|^{1-n} \cos\left(k-\frac{(n-1)\pi}{4}\right) 
\overline{\psi(x) }V(dy)\,dx  \,.
\end{eqnarray*}

  We set these three terms equal to $ L_{1}(k)$, 
$L_{2}(k)$,  and $L_{3}(k)$ respectively.  Our proof of item 3 of 
Theorem \ref{p8appthm} will be completed by showing that 
$L_{1}(k)$  and $L_{2}(k)$ vanish as $k \to \infty $ and 
that $\lim_{k\to\infty} L_{3}(k) =  2\pi (\nu,\psi)$.  
Regarding $L_{2}(k)$, by using (\ref{ass1eq3}), (\ref{Delta}) and 
the fact that for $|x-y|<1$, $ |x-y|^{\frac{1-n}{2}}|\leq |x-y|^{2-n}  
$ we see that
\begin{eqnarray} \lefteqn{ 4\pi \left(\frac{k}{2 \pi}\right)^{\frac{n-
1}{2}}  \int\!\int_{|x-y|\leq\delta}  \left| \Delta(k|x-y|) 
\cos\left(k-
\frac{(n-1)\pi}{4}\right) \overline{\psi(x) } \right|  |V|(dy)\,dx 
}\nonumber \\
& \leq & \mbox{constant} \times  \int\!\int_{|x-y|<\delta} |x-y|^{2-n} 
|\psi(x)| 
|V|(dy)\,dx \to 0\quad  \mbox{as $\delta \to 0$} 
\nonumber \end{eqnarray} 
Furthermore, for fixed ${\delta>0}$ we also have that 
\begin{eqnarray*} \lefteqn{ 4\pi \left(\frac{k}{2 \pi}\right)^{\frac{n-
1}{2}}  \int\!\int_{|x-y|>\delta}  \left| \Delta(k|x-y|) \cos\left(k-
\frac{(n-1)\pi}{4}\right) \overline{\psi(x) } \right| |V|(dy)\,dx 
} \\
& \leq & constant \times  k^{-1}  \int\!\int_{|x-y|>\delta} |x-
y|^{\frac{1-n}{2}} |\psi(x)| |V|(dy)\,dx \nonumber \\ & \leq & 
constant \times k^{-1} ||\psi||_{L^2} \left\{ \sup_{x} \int_{|x-
y|>\delta} |x-y|^{\frac{1-n}{2}} |\psi(x)| |V|(dy) 
\right\}\hspace{1.8cm}
\end{eqnarray*}
which vanishes as $k \to \infty$.  This shows that 
$\lim_{k\to\infty} L_{2}(k)= 0 $. A similar approach can 
be applied to show that $\lim_{k\to\infty} L_{1}(k)= 0 $.  
Considering $L_{3}(k)$ we first note that applying an elementary 
trig identity we have
\begin{eqnarray}  L_{3}(k) & =  & 2\pi \int \int |x-y|^{1-n}\left(1+ 
\cos(2k|x-y|-\frac{(n-1)\pi}{2})\right ) \overline{\psi(x)} V(dy)\,dx 
\nonumber  \\ & = & 2\pi(\nu, \psi) + L_{4}(k) + 
L_{5}(k)\,, \nonumber 
\end{eqnarray} where
\begin{eqnarray*} L_{4}(k) & = & 2\pi \int \int_{|x-y|\leq 1} |x-
y|^{1-n} \cos\left(2k|x-y|-\frac{(n-1)\pi}{2}\right)\overline{\psi(x)} 
V(dy)\,dx\\ L_{5}(k) & =& 2\pi \int\!\int_{|x-y|>1}  |x-y|^{1-n} 
\cos\left(2k|x-y|-\frac{(n-1)\pi}{2}\right) \overline{\psi(x)} V(dy)\,dx 
\end{eqnarray*}
Concerning $L_{4}(k)$ we set
\begin{equation} F_{1}(k,y) =  \int_{|x-y|\leq 1} |x-y|^{1-n} 
\cos(2k|x-y|-\frac{(n-1)\pi}{2}) \psi(x)\,dx
\label{eqF1} \end{equation}
It is easy to see that there is a bounded set containing the supports 
of 
$F_{1}(k,y)$ for all $k$, for if the support of $\psi$ is contained in 
$B_{R}=\{x\in {{\mathbb R}^n} : |x| < R\}$ then for any $k$,  $F_{1}(k,y) = 
0$ for $|y|> R+1$  It is also easy to see that $F_{1}(k,y)$ is 
uniformly bounded on ${\mathbb R}^{+} \times {\mathbb R}^n $. Furthermore, by 
the Riemann-Lebesgue lemma, for each fixed $y$, $F_{1}(k,y)$ 
vanishes as $k\to\infty$.  Applying the bounded 
convergence theorem we have $\lim_{k\to\infty} L_{4}(k) 
= 0$.    Turning our attention to $L_{5}(k)$  we have 
\begin{eqnarray} \left| L_{5}(k) \right| & \leq & \int \! \int_{|x-
y|>1} 
|x-y|^{1-n} |\psi(x)|  |V|(dy) \,dx \nonumber \\ & \leq & \left\{ 
\sup_{x} \int_{|x-y|>1} |x-y|^{\frac{1-n}{2}} |V|(dy)  \right\} 
||\psi(x) || _{L^{1}}\,, \nonumber  
\end{eqnarray} 
which is finite by (\ref{ass1eq1}). We can conclude 
by Fubini-Tonelli that if
$ F(y) =  \int_{|x-y|>1} |x-y|^{1-n} |\psi(x)|  \,dx $ then $F(y)$ is in 
$ L^{1}({\mathbb R}^n;V(dy))$.  Since $\psi$ is in ${\cal K}$ we  can 
conclude that for all $y$,  $|x-y|^{1-n} \psi(x)$ is in $L^{1}({\mathbb R}^n; 
dx)$. Setting 
\begin{equation} F_{2}(k,y) = \int_{|x-y|>1} |x-y|^{1-n}\cos(2k|x-
y|-\frac{(n-1)\pi}{2}) \overline{\psi(x)}\,dx
\label{eqF2} \end{equation} we can conclude by the Riemann-Lebesque 
lemma that $F_{2}(k,y)$ vanishes as $k\to \infty$ for almost 
all $y$.  
Since $F(y) \geq F(k,y)$ we can apply dominated convergence to conclude 
\begin{equation} |L_{5}(k)| = 2\pi \left| \int F(k,y) V(dy) \right| \, 
\to 0 \quad \mbox{as } k \to \infty\,.
\end{equation}
Since ${\cal K}$ is dense in ${\cal S}$, this completes the 
verification that $h_{1}(k) \to -2\pi\nu$ proving item 3 of 
Theorem \ref{p8appthm}.

 To verify item 4 of the theorem we let $\psi \in {\cal S}$ 
and set $\varphi (\xi) = \alpha^{-1}(2\pi)^{-\frac{n}{2}} |\xi| {\cal 
F}\psi (\xi)$   where $\alpha = 2^{\frac{2-n}{2}} 
\sqrt{\pi}\Gamma(\frac{n-1}{2})^{-1}$.  We claim that 
\begin{equation} \left[|x|^{1-n}* {\cal F}^{*}\varphi\right](x) = 
\psi(x)  \quad  \forall x \in {\mathbb R}^n. \label{claimeq} 
\end{equation} By using the fact that ${\cal F}(|x|^{1-n})(\xi)  =  
\alpha |\xi|^{-1}$ (e.g. see \cite{shilov}) we have for any $\phi 
\in {\cal K}$, 
\begin{eqnarray} \langle \phi, |x|^{1-n}* {\cal F}^{*}\varphi\rangle  & 
= & \alpha^{-1} (2\pi)^{-\frac{n}{2}} \langle |x|^{1-n}*\phi, {\cal 
F}^{*} |\xi| {\cal F}\psi (\xi)\rangle  \nonumber \\ & = & \alpha^{-1} 
(2\pi)^{-\frac{n}{2}} \langle {\cal F}[|x|^{1-n}*\phi], |\xi| {\cal 
F}\psi 
(\xi)\rangle  \nonumber \\ & = & \langle |\xi|^{-1}{\cal F}\phi, |\xi| 
{\cal 
F}\psi (\xi)\rangle  \nonumber \\ & =&  \langle \phi, \psi \rangle  
\nonumber \end{eqnarray}
This together with the fact that both sides of (\ref{claimeq}) are 
continuous proves our claim.  We now have;
\begin{eqnarray} \langle  \nu,{\cal F}^{*}\varphi \rangle  & = & \langle  
|x|^{1-n}* V, {\cal F}^{*}\varphi\rangle   \nonumber \\ & = & 
\langle   V, |x|^{1-n}*{\cal F}^{*}\varphi \rangle   \nonumber \\ & = 
&  \langle   V , \psi \rangle   \nonumber \end{eqnarray} Thus in the 
sense of distributions, $\alpha^{-1}(2\pi)^{- \frac{n}{2}} {\cal 
F}^{*} |\xi| {\cal F}\nu = V$ completing the proof of item 4.  
Note that once $V$ is recovered, it is a simple matter to 
reconstruct the associated measure potential, $V(dx)$ (see 
\cite{ford6} for details). 



\subsection*{A Singular Example} 

 We now provide a concrete example illustrating the use of 
these results.  We will let $n = 3$ and define $V(dx)$ as the 
measure, $\beta \delta(|x|-R)$, the delta function of the sphere of 
radius $R$ and strength parameter $\beta$.  For any $\psi \in {\cal 
K}$ we have
\begin{equation} \langle  V, \psi \rangle  = \int_{S_{R}}\beta \psi(x) 
\,d\omega(x)\,, \label{exeq1} \end{equation} 
where $S_{R}$ is the 
sphere of radius $R$ and $d\omega(x)$ is the inherited surface 
measure. The verification that Assumption 1 of Theorem
\ref{p8appthm} holds is straightforward.  In this explicit example 
$G(x,k) = -\frac{1}{4\pi}\frac{e^{i\kappa|x|}}{|x|}$ and so $Q(z)$ 
is given by;
\begin{equation} Q(\kappa^{2})f(x) = -\frac{\beta}{4\pi} \int_{S 
_{R}} \frac{e^{i\kappa|x-y|}}{|x-y|} f(y)\,d\omega(x)\,. 
\label{exQeq} \end{equation} Taking $F(x,y) = \frac{1}{4\pi|x-
y|}= |G(x-y;k)|$ we note that the $\sup_{y} \int F(x,y) |V|(dx)$ 
occurs when $|y| = R$, and therefore
\begin{equation} \sup_{y} \int F(x,y) |V|(dx)  =  \frac{\beta}{4\pi}
\int_{S_{R}} \frac{1}{|x-y|}\,d\omega(x) = \beta \, R\,.  
\nonumber \end{equation} In addition we have $\int |V|(dx) = 
\int_{S^{2}}\,d\omega = \Omega^{3}$. Thus we see that 
Assumption 2B holds whenever $\beta < R$.  Theorem 
\ref{p8appthm} then provides us the following:
\begin{equation}
\lim_{k\to\infty} \left( -2\pi i k \langle (1-
S(k))e^{-ik\omega\cdot x},e^{-ik\omega\cdot x}\rangle ,\psi(x) \right) 
= \int_{S_{R}} \frac{\beta \overline{\psi(y)}}{|x-y|}\,\, d\omega(y)
\label{example1} \end{equation} 
Furthermore,  (\ref{item3}) holds as well giving
\begin{eqnarray}
\lefteqn{ \lim_{k\to\infty} \left( -(2\pi)^{-
\frac{n}{2}} i k \alpha^{-1} \langle (1-S(k))e^{-ik\omega\cdot 
x},e^{-ik\omega\cdot x}\rangle _{\omega} ,{\cal F}^{*}(|\xi|{\cal 
F}\psi(x) \right) }\nonumber\\
& =& \int_{S_{R}} \beta \overline{\psi(x)} \,d\omega(x)\,. 
\hspace{7cm}
  \label{example2} \end{eqnarray}
We summarize this 
example by stating the following:

\begin{theorem} Fix $R>0$ and let $0 < \beta < R$.  Let $H$ be the self-adjoint 
realization of the Laplacian in ${\mathbb R}^{3}$.  Then the following 
hold:

\noindent{\bf 1.} There exists a self-adjoint operator, $H_{1}$, satisfying
$$ (H_{1} u , v) = (u, Hv) + \beta \int_{S_{R}} u(x) 
\overline{v(x)}\,d\omega(x)
$$ 
for all $u$ in $D(H_{1})$ and $v$ in $D(H)$. 

\noindent{\bf 2.} The associated wave operators, $W_{\pm}(H_{1},H)$ 
exist and are strongly complete.

\noindent{\bf 3.} There exist generalized eigenfunctions, $\phi_{\pm}(x,\xi)$ 
satisfying
$$ \phi_{\pm}(x,\xi) = e^{-ix\cdot\xi}  - 4\pi \beta 
\int_{S_{R}} \frac{\phi_{\pm}(y,\xi) e^{\pm i|\xi||x-y|}}{|x-y|} 
d\omega(y). $$ 

\noindent{\bf 4.} The scattering matrix, $S(k)$, exists as a unitary 
operator on $L^{2}(S^{2})$ satisfying 
$$ [S(k)f(k,\cdot)](\omega) = \left[{\cal F} 
W_{+}^{*}W_{-}{\cal F}^{*}f \right](k\omega) \quad \forall\,\,f \in 
C_{0}^{\infty}({\mathbb R}^{3})\,. $$

\noindent{\bf 5.} The scattering matrix admits the representation,
$$ S(k) = 1 - \frac{ik}{8\pi^{2}} \gamma_{+}(k)\gamma(k)\,,
$$ 
 where $\gamma$ and $\gamma_{+}$ are given by (\ref{results5}) 
and (\ref{results6}).

\noindent{\bf 6.} The inverse scattering results, equations (\ref{example1}) and 
(\ref{example2}) hold.
\end{theorem} 

 The remainder of this work is devoted to the proof of main 
theorem. 

\section{Proof of Theorem 2.1} \setcounter{equation}{0}

  We will prove the theorem through a series of lemmas.  

\begin{lemma}  Fix  $z \neq 0$ in ${\mathbb R}^n$. Let $ (\rho, 
\theta_{1}, \theta_{2}, \ldots , \theta_{n-1} ) $ be spherical 
coordinates with $\theta_{1}$ the polar angle measured from the 
positive $z$ direction.  Let $\Omega$ be an open $n$-rectangle in  
${\mathbb R}^n$, i.e. 
\[ \Omega = \{ (\rho, \theta_{1}, \theta_{2}, \ldots , \theta_{n-1} )| 
\rho\in (a_{0},b_{0}), \theta_{1}\in (a_{1},b_{1}), \ldots, 
\theta_{n-1} \in (a_{n-1}, b_{n-1}) \}\,, \]
where $0\leq a_{0} < b_{0}\leq \infty$, $0\leq a_{n-1}< b_{n-1} 
\leq 2\pi$, and $0\leq a_{i}<b_{i} \leq  \pi$, $i=1,2,\ldots, n-2$.  Let 
${\bf X}_{\Omega}(x)$ be the characteristic function  for $\Omega 
\subset {\mathbb R}^{n}$.  Set 
\begin{equation} c(k,r) = \cos(k r-\frac{(n-1)\pi}{4} ) \quad k, r \in 
{\mathbb R}, 
\label{eqforc} \end{equation}
and set
\begin{equation} I(k,z) = \int_{{\mathbb R}^{n}} \frac{c(k,|x|)}{|x|^{n-
2}} \frac{ c( k, |z-x| )}{|z-x|} (\sin \theta_{1})^{3-n} {\bf 
X}_{\Omega}(x)\,dx \,.
\label{lem1eq1} \end{equation}
Then for all $z \neq 0$
\begin{equation}
\lim_{k\to\infty} I(k,z) = 0\,. 
\label{lem1eq2} \end{equation}  \label{lem1} \end{lemma}

\paragraph{Proof:}   Noting that the integrand in (\ref{lem1eq1}) is 
dependent only on $\rho$ and $\theta_{1}$, and calculating the 
integral in spherical coordinates, we have that
\begin{equation} I(k,z) = C \int_{ a_{0}}^{b_{0}} 
\int_{a_{1}}^{b_{1}}  \frac{ c(k, \rho) }{\rho^{n-2}} \frac{c( k, u) 
}{u}  \rho^{n-1}  \sin \theta_{1} \, d\theta_{1} \,d\rho \,,
\label{lem1eq3} \end{equation} 
where $C$ is the constant obtained by integrating over the angular 
variables, $\theta_{i}$, $i  \geq 2$ and where $u = |z-x|$.  We make 
a change of coordinates in the $\theta_{1}$-variable by setting 
$$
u=  |z-x| =  ( |z|^2 + |x|^2 - 2 |z| |x| \cos \theta_{1})^{1/2}
$$ 
to obtain
\begin{eqnarray}
I(k,z) &=& C \int_{a_{0}}^{b_{0}} \int_{u_{a}}^{u_{b}}  \frac{ c( 
k, \rho ) c(k, u )}{|z| }  \, d\theta_{1}\, d\rho \nonumber\\
& =&  \frac{C}{k |z|}  \int_{a_{0}}^{b_{0}}  c(k, \rho) 
\left( s(k, u_{b}) - s(k, u_{a}) \right)\, d\rho \,,
\label{lem1eq4} \end{eqnarray}  
where $u_{a}$ and $u_{b}$ are the appropriate new limits of 
integration in the $u$-variable.  And where 
$$ s(k, \rho) = \sin(k \rho -\frac{(n-1)\pi}{4})   
$$ 
 We now can see that the final integral in 
(\ref{lem1eq4}) is uniformly bounded in $k$, hence 
$|I(k,z)| \leq \mbox{constant }/k$ which vanishes as $k\to\infty$. 
\hfill $ \Box$ \smallskip

 The following corollary is immediate.

\begin{corollary} \label{cor1}
Let $\phi(x)$ be any step function on ${\mathbb R}^{n}$ of the form 
$\sum_{j=1}^{m}  a_{j} {\bf X}_{\Omega_{j}}(x)$ where each 
$\Omega_{j}(x)$ is an open $n$-rectangle as in lemma 
(\ref{lem1}).   Then for all $z \neq 0$
\begin{equation}
\int_{{\mathbb R}^{n}} \frac{c( k,|x|)}{|x|^{n-2}} 
\frac{c(k , |z-x| ) }{|z-x|} (\sin \theta_{1})^{3-n} \phi(x)\,dx 
\to 0 \quad \mbox{as } k\to\infty \,.\label{lem1eq5} 
\end{equation} 
 \end{corollary}

\begin{lemma}
Let $z \neq 0$ be fixed in ${\mathbb R}^{n}$ and set  $\Psi(x,k) =  
\frac{c(k,|x|)}{|x|^{n-2}} 
\frac{c(k ,|z-x|)}{|z-x|} (\sin \theta_{1})^{3-n}$ where 
$\theta_{1}$ and $c(k,r)$ are as in Lemma \ref{lem1}. Then for $0< p < \frac{n}{n-2}$, and any bounded, 
measurable set, $\Omega 
\subset {\mathbb R}^{n}$,  $\Psi(x,k) \in L^{p}(\Omega)$ for each $k$.  
Furthermore, there 
exists a constant, $C(\Omega,z,p)$ depending on $\Omega$, $z$,  
and $p$ but independent 
of $k$ such that $||\Psi(x,k)||_{p} < C(\Omega,z,p)$.
\label{lem2} \end{lemma}

\paragraph{Proof:}   If we fix, $\Omega \subset {\mathbb R}^{n}$, 
and let $M = \sup \{|x|\, : \, x\in\Omega\} $.  We have 
\begin{equation} 
 \int_{\Omega} | \Psi(x,k) |^{p}\,dx =  I_{1} +  I_{2}  +  I_{3}\,, 
\end{equation} 
where 
\begin{eqnarray*}
I_{1} & = &   \int_{\Omega_1}  | \Psi(x,k) |^{p}\,dx\,,  \quad 
\Omega_1 = \Omega \cap \left\{ x \, : \, |x| < \frac{|z|}{2}\,, 
\right\}  \\
I_{2} & = &  \int_{\Omega_2} | \Psi(x) |^{p}\,dx \,,  
\quad \Omega_2 = \Omega  \cap \left\{x \, : \, |x-z|< 
\frac{|z|}{2} \right\}\,,
\end{eqnarray*}  
and $I_{3}$ is the integral over the remaining region.  Considering 
$I_{1}$ we switch to spherical coordinates taking the $z$ direction 
as the  polar axis and letting $\theta_{1}$ be as before.  We obtain
\begin{eqnarray}
|I_{1}|  & \leq&  \int_{|x|<\frac{|z|}{2}} |x|^{(2-n)p} \left( 
\frac{2}{|z|} \right)^{p} (\sin 
\theta_{1})^{(3-n)p}\,dx  \label{lem2eq2} \\  & \leq &
C  \left( \frac{2}{|z|} \right)^{p} \int_{\rho=0}^{\frac{|z|}{2}} 
\int_{\theta_{1}=0}^{\pi} \rho^{(2-
n)p+n-1}  (\sin \theta_{1})^{(3-n)p+n-2}\, d\theta_{1}\, d\rho\,, \nonumber
 \end{eqnarray} 
 where $C$ is as in Lemma \ref{lem1} with $a_{i}=0$ 
for each $i = 1,2,\ldots, n-1$ and 
$b_{i} = \pi$ for $i = 1, 2, \ldots, n-2$ and $b_{n-1} = 2 \pi$.  The 
singularities at $\rho = 0$ and $\theta_{1} = 0$ are integrable 
provided that 
\begin{equation}
(2-n)p+n-1 > -1  \mbox{ and } (3-n)p+n-2> -1
\label{lem2eq3} \end{equation}
both of which are satisfied when $p < \frac{n}{n-2}$.  Note that for 
such $p$ we have that $|I_{1}|$ is bounded by a constant depending 
only on $z$.  Looking at $I_{2}$ we make a change of variable and let 
$w = z-x$.  We will again switch 
to spherical coordinates taking  the $z$ direction as the polar axis, 
but setting 
$\tilde{\theta_{1}}$ to be the angle between $z$ and $w$.   By the 
law of sines we have 
that $\sin \theta_{1} = \frac{|w|}{|x|} \sin \tilde{\theta_{1}}$ and so 
we obtain
\begin{eqnarray}
|I_{2}| & \leq &  \int_{|w|<\frac{|z|}{2}} |x|^{p(2-n)} |w|^{-p}  
(\frac{|w|}{|x|} \sin 
\tilde{\theta_{1}})^{(3-n)p}\,dx \nonumber \\ 
 & = &\int_{|w|<\frac{|z|}{2}} |x|^{-1}  |w|^{(2-n)p} (\sin 
\theta_{1})^{(3-n)p}\,dx \label{lem2eq4} \\ 
&\leq & C  \int_{\rho=0}^{\frac{|z|}{2}} \int_{\tilde{\theta_{1}}=0}^{\pi}  
\left(  \frac{2}{|z|} \right)   \rho^{(2-n)p+n-1}  (\sin 
\tilde{\theta_{1}})^{(3-n)p+n-2}\, d\tilde{\theta_{1}}\, 
d\rho\,.\nonumber
 \end{eqnarray} 
Note that the integrability conditions on $p$ are identical 
with those found in  $I_{1}$.  
Thus for $p< \frac{n}{n-2}$, $|I_{2}|$ is also bounded by a 
constant depending only on 
$z$.

Finally, considering $I_{3}$ we note that over the 
remaining region we have
\begin{equation} 
|\Psi(x,k)|^{p} \leq \left(\frac{2}{|z|}\right)^{(n-1)p}
\label{lem2eq5} \end{equation} 
and so $|I_{3}|$ is bounded by this constant times the 
volume of the sphere of radius $M$ in ${\mathbb R}^{n}$  Combining 
these results for 
$I_{1}$, $I_{2}$, and $I_{3}$ we have that $\int_{\Omega} 
|\Psi(x,k)|^{p}\,dx$ is bounded 
by a constant depending only on $z$ and $\Omega$.   \hfill   
$\Box$

\begin{lemma}
Let  $z$ be fixed in ${\mathbb R}^{n}$, $n\geq 3$ and let $\psi(x) \in 
C_{0}^{\infty} ({\mathbb R}^{n})$.  Let $\theta_{1}$ be the angle 
between $z$ and $x$ as in Lemma \ref{lem2}. Then 
$F(x,\psi) = (\sin \theta_{1})^{n-3} |z-x|^{\frac{3-n}{2} } 
|x|^{\frac{n-3}{2} } \psi(x)$ is a 
bounded function with compact support and hence is in $L^{p}({\bf 
R}^{n})$ for all $p\geq 1$.
\label{lem3} \end{lemma}

\paragraph{Proof:}   To prove the lemma it is enough to verify 
that $|z-x|^{\frac{3-n}{2}}  (\sin \theta_{1})^{n-3} $ is bounded.  
To see this we first set for $r>0$,
\[ \Theta(r) = \sup \{ \theta_{1} : |z-x| = r \}\,.  \]

 By elementary geometric considerations we have for any 
$x$ such that $|z-x|<|z|$  that $|z-
x|= |z| \sin \Theta(|z-x|) $.   This fact in turn implies that 
\begin{equation}
|z-x| \geq |z| \sin \theta_{1}\quad \mbox{whenever } |z-x|<|z|\,.
\label{lem3eq1} \end{equation}
From this we see that for $x$ such that $|z-x|<|z|$ we 
have,
\begin{equation}
|z-x|^{\frac{3-n}{2} }  |\sin \theta_{1}|^{n-3}   \leq |z|^{\frac{3-
n}{2} } |\sin \theta_{1}|^{\frac{n-
3}{2}} \leq |z|^{\frac{3-n}{2} }
\label{lem3eq2} \end{equation}
and trivially for $x$ such that  $|z-x| \geq |z|$, we have
\begin{equation} |z-x|^{\frac{3-n}{2}} |\sin \theta_{1}|^{n-3}  \leq 
|z|^{\frac{3-n}{2}}. 
\label{lem3eq3} \end{equation}  \hfill  $\Box$


\begin{lemma}
Let $\psi(x) \in C_{0}^{\infty}({\mathbb R}^{n})$  and $c(k,r)$ be  
defined as in Lemma \ref{lem1} and set 
\begin{equation} 
J(x,y,k,\psi) = \int_{{\mathbb R}^{n}} \frac{ c(k, |x - \xi| )}{|x - 
\xi|^{\frac{n-1}{2}} } \frac{ c(k, |\xi - y|)}{ |\xi - y|^{\frac{n-
1}{2} } } \psi(\xi)\,d\xi \,.
\label{lem4eq1} \end{equation}
Then for all $x \neq y$,  $\lim_{k\to\infty} J(x,y,k,\psi) = 
0$.
\label{lem4} \end{lemma}

\paragraph{Proof:}   Setting $w = x - \xi$ and $z = x-y$ we 
obtain
\begin{eqnarray}
J(x,y,k,\psi) & = & \int_{{\mathbb R}^{n}} \frac{ c(k, |w|)}{|w|^{\frac{n-
1}{2}} } \frac{ c(k,|z-w|)}{ |z-w|^{\frac{n-1}{2} } } \psi_{x}(w) 
dw \nonumber \\
& = & \int_{{\mathbb R}^{n}} \Psi(k,w) F(w,\psi_{x})\, dw\,, 
\label{lem4eq2}
 \end{eqnarray}
where $\Psi$ is as defined in Lemma \ref{lem2} and $ 
F(w,\psi_{x})$ is as defined in Lemma \ref{lem3} with 
$\psi_{x}(w) = \psi(x-w)$.  We now fix $p \in (1, \frac{n}{n-2})$  
and let $p'$ satisfy $\frac{1}{p} +\frac{1}{p'} = 1 $.  Let 
$\epsilon>0$ be arbitrary but fixed.  Let $C(\Omega_{x},z,p)$ be 
the uniform bound from Lemma \ref{lem2} on the $p$-norm of 
$\Psi(k,w)$ and a step function, $\phi(w)$ such that
\[   ||F(w,\psi_{x})-\phi(w)||_{L^{p}(\Omega_{x})} < \epsilon\,.
 \]
Then by Lemmas \ref{lem2} and \ref{lem3}, we have
\begin{eqnarray} 
\lefteqn{ \lim_{k\to\infty} \left|J(x,y,k,\psi)\right| }\nonumber\\
& \leq & \lim_{k\to\infty}  \left|  \int_{{\mathbb R}^{n}} \Psi(k,w) 
  F(w,\psi_{x})\, dw  \right|  \nonumber \\ 
& \leq & \lim_{k\to\infty}  \left|  \int_{{\mathbb R}^{n}} 
 \Psi(k,w) \phi(w) dw  \right|
 +  \left|  \int_{{\mathbb R}^{n}} \Psi(k,w) (F(w,\psi_{x}) -\phi(w))dw  
 \right|  \nonumber \\
& \leq & C(\Omega_{x},z,p)   
 ||F(w,\psi_{x}) -\phi(w)||_{L^{p'}(\Omega_{x}) } \label{lem4eq3} \\  
& \leq & C(\Omega_{x},z,p)   \epsilon \,. \nonumber
 \end{eqnarray}
The present proof follows from the fact that $\epsilon$ was arbitrary.  \hfill  
$\Box$ 

\paragraph{Proof of Theorem \ref{maintheorem}}
 With the aid of this last lemma we now now provide the 
proof of the main theorem.  We begin by setting 
\begin{eqnarray*}
&C(r)  =   2 \left( \frac{2 \pi}{r} 
 \right)^{\frac{n-1}{2}} \cos \left( r - \frac{(n-1) \pi}{4}  \right)\,, &\\ 
&\Delta(r)  =  J(r) - C(r)\,, &
\end{eqnarray*} 
where $J(r)$ is defined by 
(\ref{introeq1}). The asymptotic behavior of $J(r)$ is given (see 
Sait\={o} \cite{saito3}) by 
$$ |\Delta(r)| = |J(r) - C(r)|  \leq {\bf A} r^{- 
\frac{(n+1)}{2}}\,,
$$ 
where ${\bf A}$ is a constant. With this notation, we have that
\begin{eqnarray*}
\lefteqn{k^{n-1}\int_{{\mathbb R}^{n}}J(k|x-\xi|) J(k|\xi - y|)  \psi(\xi) 
d\xi }  \\ 
& = & k^{n-1} \int_{{\mathbb R}^{n}} \left( C(k|x-\xi|)
  + \Delta(k|x-\xi|) \right)   \left( C(k|\xi - y|) + \Delta(k|\xi-
   y|) \right) \psi(\xi)\,d\xi   \\ 
&=&  k^{n-1} \int_{{\mathbb R}^{n}}  C(k|x-\xi) C(k|\xi - y|) \psi(\xi) 
   \,d\xi  \\
&&+ k^{n-1} \int_{{\mathbb R}^{n}}  \Delta(k|x-\xi|) C(k|\xi - y|) 
  \psi(\xi)\,d\xi \\ 
&& +  k^{n-1} \int_{{\mathbb R}^{n}}  
  C(k|x-\xi|) \Delta(k|\xi-y|) \psi(\xi)\,d\xi \\
&& + k^{n-1} \int_{{\mathbb R}^{n}}  \Delta(k|x-\xi|) \Delta(k|\xi-y|)
  \psi(\xi)\,d\xi  \\ 
&=& I_{1} + I_{2} + I_{3} + I_{4}\,. 
\end{eqnarray*}
We have that $I_{1}$ vanishes by Lemma \ref{lem4}.  Looking at 
$I_{2}$ we have
\begin{eqnarray*} 
|I_{2}|  & \leq & k^{n-1} \int_{{\mathbb R}^{n}}  
{\bf A} (k|x-\xi|)^{-\frac{(n+1)}{2}} |C(k|\xi - y|)| |\psi(\xi)|\,d\xi \\
&\leq & k^{n-1} \int_{{\mathbb R}^{n}}  {\bf A} 
 (k|x-\xi|)^{- \frac{(n+1)}{2}} 2\left(\frac{2\pi}{k|\xi - y|}
 \right)^{\frac{n-1}{2}} |\psi(\xi)|\,d\xi  \\ 
&\leq & \frac{{\bf A}}{k}  \int_{{\mathbb R}^{n}}\frac{|\psi(\xi)|}
{|x-\xi|^{\frac{n+1}{2}}{|\xi-y|^{\frac{n-1}{2}}}}\,d\xi 
\end{eqnarray*}
Clearly the last integral is finite and independent of $k$, hence 
$I_{2}$ vanishes as $k \to \infty$.
$I_{3}$ behaves just like $I_{2}$, hence it vanishes as $k 
\to \infty$.   For the final term, we have
\begin{eqnarray*}
|I_{4}| &\leq&  k^{n-1} \left| \int_{{\mathbb R}^{n}}  \Delta(k|x-\xi|) 
  \Delta(k|\xi-y|) \psi(\xi)\,d\xi\right|\\
& \leq& \frac{{\bf A}^{2}}{ k^{2}} 
\int_{{\mathbb R}^{n}}  |x-\xi|^{- \frac{(n+1)}{2}} |\xi-y|^{- 
\frac{(n+1)}{2}}  |\psi(\xi)|\,d\xi 
\end{eqnarray*}
 This integral is also finite and independent of $k$, hence $I_{4}$ 
vanishes as $k \to \infty$.  This completes the proof of our theorem.  

 \begin{thebibliography}{99}

\bibitem{albeverio} S. Albeverio, F. Gesztesy, R. H\/{o}egh-
Krohn, H. Holden, {\em Solvable Models in Quantum Mechanics}, 
Springer-Verlag: Berlin - Heidelberg, New York (1988)

\bibitem{alsholm} P. Alsholm \& G. Schmidt,  Spectral and 
Scattering Theory for Schr\"{o}dinger Operators, {\em Arch. 
Rational Mech. Anal}
{\bf 40} (1971) 281-311

\bibitem{brasche} J. F. Brasche, P. Exner, Y. A. Kuperin \& P. 
\v{S}eba, Schr\"{o}dinger Operators with Singular Interactions, 
Preprint Nr. 132  Nov. 1991

\bibitem{cheney} M. Cheney,  Inverse Scattering in Dimension 
Two, {\em J. Math. Phys.}  {\bf 25}  (1984)  94-107

\bibitem{cheney2} M. Cheney, J. Rose, Generalization of the 
Fourier Transform: Implications for Inverse Scattering Theory, {\em 
Phys. Rev. Lett.} {\bf 60} 1221-1224 (1988)

\bibitem{shabani}  L. Dabrowski, J. Shabani,  Finitely Many 
Sphere interactions in Quantum Mechanics: Nonseparated Boundary 
Conditions  {\em J. Math. Phys.}  {\bf 29}   (1988) 2241-2244

\bibitem{faddeev}  L. D. Faddeev, The Uniqueness of Solutions 
for the Scattering Inverse Problem, {\em Vestnik Leningrad Univ.} 
{\bf 7} (1956) 126-130 

\bibitem{ford1} R. Ford, Generalized Potentials and Obstacle 
Scattering, {\em Trans. AMS}, {\bf 329}   (1992) 415-431

\bibitem{ford3} R. Ford, The Scattering Matrix for Singular 
Schr\"{o}dinger Operators,  {\em J. Math Anal.Appl.}, {\bf 175}   
(1993)  588-605

\bibitem{ford5} R. Ford, Scattering Results for Measure 
Potentials with Unbounded Support, {\em Math. Meth. in the Appl. 
Sci.} {\bf 17}   (1994)  1213-1230

\bibitem{ford6} R. Ford, An Inverse Scattering Result for 
Measure Potentials, {\em Inv. Prob.} {\bf 11}   (1995)  939-948

\bibitem{shilov} I. M. Gel'fand, G. E. Shilov  {\em Generalized 
Functions: Vol 1}, Academic Press, New York and London  (1964)


\bibitem{hormander} L. H\"{o}rmander, {\em The Analysis of 
Linear Partial Differential Operators I}, Springer-Verlag, Berlin 
(1983)

\bibitem{newtonbook} R. G. Newton, {\em Inverse 
Schr\"{o}dinger Scattering in Three Dimensions}, Springer-Verlag, 
Berlin (1989) 

\bibitem{serov2} L. P\"{a}iv\"{a}rinta, V. S. Serov, and E. 
Somersalo,  Reconstruction of Singularities of a Scattering 
Potential in Two Dimensions, {\em Adv. Appl. Math.}, {\bf15}  (1994), 
97-113  

\bibitem{paivarinta1} L. P\"{a}iv\"{a}rinta,  and E. 
Somersalo,  
Inversion of Discontinuities for the Schr\"{o}dinger Equation in 
Three Dimensions, {\em Siam J. Math. Anal.}, {\bf22}  (1991), 480-
499 

\bibitem{saito1} Y. Sait\={o}, Some Properties of the 
Scattering Amplitute and the Inverse Scattering Problem, {\em 
Osaka J. Math.} {\bf 19} (1982), 57-78

\bibitem{saito2} Y. Sait\={o},  An Inverse Problem in Potential 
Theory and the Inverse Scattering Problem,  {\em J. Math Kyoto 
Univ.} {\bf 22-2} (1982) 307-321

\bibitem{saito3} Y. Sait\={o},  An Asymptotic Behavior of the 
S-Matrix and the Inverse Scattering Problem, {\em J. Math. Phys.} 
{\bf 25} (1984) 3105-3109

\bibitem{serov1} V. S. Serov,  On Recovering the Potential in 
a 
Three-Dimensional Problem of Scattering Theory {\em Soviet Math. 
Dokl.}, {\bf 43} (1991)  471-475 

\bibitem{serov3} V. S. Serov,  A Scattering Problem for a 
Schr\"{o}dinger Operator with Singular Potential in the Two 
Dimensional Case II, {\em Diff. Eqns.} {\bf 26} (1990) 616-623

\bibitem{weder} R. Weder, Multidimension Inverse Scattering: 
The Reconstruction Problem, {\em Inverse Problems} {\bf 6} (1990)                       
        \end{thebibliography}                         
\medskip


\noindent {\sc Richard Ford } \\ 
Department of Mathematics \\ 
California State University \\ 
Chico,  CA 95929, USA \\
e-mail: rford@csuchico.edu 

\end{document}
