Inverse Diffraction Theory and Computation of Minimum Source Regions of Far Fields

A methodology based on the multipole expansion is developed to estimate the minimum source region of a given far field. The support of any source that produces the given far field must contain this minimum source region. The results are derived in the framework of the scalar Helmholtz equation in two-dimensional free space, which is relevant to transverse magnetic electromagnetic waves.The proposed approach consists of two steps. First we address, via an exterior inverse diffraction framework, the estimation of theminimum convex source region, which is the convex hull of theminimum source region. Next we compute, via a complementary interior inverse diffraction approach, nonconvex bounds for the minimum source region. This allows, in theory, the estimation of the minimum source region which can be nonconvex. The derived approach is illustrated with analytical and numerical examples relevant to inverse source and scattering problems.


Introduction
We consider in the framework of the scalar Helmholtz equation the inverse problem of estimating the minimum source region of a given far field.The support of any source that produces the given far field must contain this minimum source region.This region is also sufficiently big so that there exists a source supported in it that radiates the given far field.Thus this region is optimal in the sense that it is the smallest support of any source that can produce the given far field.Previous works tackling this fundamental inverse problem have emphasized the computation of the convex hull of the minimum source region, which has been termed the convex scattering support [1] or the minimum convex source region [2] of the given far field.The purpose of the present paper is to further extend the scope of the available theory and algorithms to compute the actual minimum source region including nonconvex minimum source regions that cannot be estimated with the existing methods.We also derive new closed form formulas for the minimum source region that are based on the multipole expansion [3] and that complement the counterpart plane-wave-theorybased, closed form method of a previous paper [2].The derived formulas facilitate analytical demonstration of the computation of the minimum source region for far fields that are known in closed form, which is of interest in antenna synthesis and as a theoretical tool to handle certain inverse scattering problems analytically.We also address the respective numerical implementation and illustrate it with examples.The results are quite encouraging, rendering good bounds of source support for the particular sources considered.On the other hand, we do not specify in the following the class of sources for which the methods of this paper (the nonconvex bounds, in particular) are useful in practice.This challenging question is related to theoretical issues raised in [4] which go beyond the scope of our presentation and is left for the future.
Next we outline the key background and ideas and sketch the plan of this paper.In this introduction, we outline the problem statement implicitly assuming the physical threedimensional (3D) space.However, to facilitate the subsequent exposition as well as the tractability of illustrative examples, in the following sections of the paper we restrict attention to the simpler two-dimensional (2D) formulation.
It has been shown in a number of investigations that, given the far field associated with an unknown radiating source, it is possible to uniquely define and compute the so-called convex scattering support [1] or minimum convex source region [2] associated with the given far field.The theoretical principles can be traced back to the treatise of Müller [5,Theorems 26,27,29] who used the spherical wave or multipole expansion of the radiated field to show that, for a fixed origin of coordinates, say , there is a minimum spherical source region of radius  min ,  min () = {r ∈ R 3 :  ≡ |r| ≤  min ()}, such that, in order for a source to be able to produce the given far field radiation pattern, it must lie in the interior of a spherical volume centered about the same origin and having a radius that is at least as large as  min .Furthermore, the radiated field outside  min is an analytic function that can be analytically continued up to the boundary of a minimum source region  min located inside  min ; see [6, page 143].Figure 1 illustrates the minimum source region and the respective minimum spherical source region corresponding to a given origin.From the point of view of propagation or diffraction, then the far field can be inverse-diffracted up to the boundary of  min ⊆  min .In addition, clearly if one considers another origin   then one can also find another minimum spherical volume  min (  ) centered about this origin, and so on.Then if one considers a set of, say,  origins and  associated minimum spherical volumes, the intersection of all such source regions, say B min , establishes a sharper bound on the minimum source region  min ⊆ B min which lies in the interior of such intersection.By means of this procedure, which has been elaborated in rigor in [1,7], it is theoretically possible to uniquely define and compute the convex scattering support or minimum convex source region, say  min,conv ⊇  min , which is a subset of the convex hull of any source producing the given far field.This minimum convex source region is the smallest convex support of any source producing that field and represents important source localization information that can be inverted uniquely from the far field data despite the general nonuniqueness (without further priors) of the full inverse source problem of reconstructing the actual source that generated the given far field; see, for example, [8,9] for discussions of the inverse source problem and the relevant nonradiating sources corresponding to the null space of the forward source-to-field mapping.In this connection, it is worth commenting that for zero far fields the minimum source region does not exist (empty set).Any nonradiating source can produce that trivial far field, and of course no source at all also produces that same far field.However, for nontrivial far fields associated with fields obeying the radiation condition the minimum source region is a nonempty set since, as is well known, the only solution to the homogeneous Helmholtz equation obeying the radiation condition is the zero field.Consequently, here and henceforth we implicitly assume nontrivial far fields, for which one can always define corresponding minimum source regions.
The theory behind the convex scattering support or minimum convex source region has been detailed in the key papers [1,4,7,10,11].Kusiak and Sylvester [1] used the circular Paley-Wiener theorem to demonstrate the existence and computation of the convex scattering support.The original focus was, as in [2], on wave radiation in free space.Potthast et al. [10] proposed the "range test" to compute in practice the convex scattering support.Kusiak and Sylvester [7] demonstrated that the convex scattering support can be defined via a Picard test.This is an important generalization which allows the definition of the convex support (of the object function) for a broad class of inverse problems, including inverse source and scattering problems in nonhomogeneous media, and problems with limited data.Sylvester [4] demonstrated the existence of a unique union of well-supported convex sets (UWSC sets) corresponding to the given far field.A source having that support exists that radiates the given far field.Also, this unique UWSC sets is a subset of the UWSC support of any source that produces the same far field.Haddar et al. [11] addressed the related problem of computing the minimum convex support of a scatterer from back-scattering data.Most of these studies were clearly motivated by the relevance of the convex scattering support in inverse scattering problems, with linkages to the sampling and probe methods [12].
The computation of the minimum convex source region has also been addressed by means of a simpler methodology based on the plane wave expansion, rendering an alternative approach to estimate the minimum convex source region of a given far field ( [2]; see also [6, pages 143-146]).Unlike the multipole expansion or circular Paley-Wiener theorem approaches, or other range or Picard tests proposed to date, the plane wave expansion method in [2] does not depend on the position of the origin of the coordinate system.On the other hand, this method is useful only if the far field is known in closed form.For real data and computations, one needs to adopt alternative approaches such as those detailed in [1,7] and other related studies [10,11].Nevertheless, from both theoretical and practical points of view, there is a lot of interest in the analytical approach in [2] since it allows handling of problems where the far field is known in closed form, which cannot be studied in rigor with computational approaches which employ discrete data.
The method in [2] applies to the minimum convex source region.For geometrical reasons this plane wave expansion approach cannot be used to estimate nonconvex minimum source regions.In this paper we use the multipole expansion to obtain the respective closed form formula to estimate the full minimum source region including nonconvex minimum source regions that cannot be derived with the existing methods.To achieve this goal, first we derive, in 2D space, the respective multipole theory counterpart of the key formula (6) in [2], for estimation of the convex support.Later we derive the respective generalization to tackle the pending nonconvex support, which completes the computation of the minimum source region.In addition to studying the analytical computation of minimum source regions, we also develop and illustrate with examples the corresponding numerical version of the approach.
We wish to point out that part of our motivation to generalize the estimation of the minimum convex source region to the actual minimum source region is provided by the generality of the Picard test as shown in [7] and the related paper [10].We show that, in theory, just as we can estimate from the far field the respective minimum convex source region via a Picard test based on an exterior inverse problem from far fields to sources in the near field region, we can also formulate the estimation of the possibly nonconvex minimum source region via a similar Picard test based on an interior inverse problem involving near fields.In the exterior problem one goes from the field data in a large spherical surface to a source contained inside a spherical volume of smaller radius.In the interior problem one goes from the near field data in a small spherical surface to a source supported outside a spherical volume of larger radius.In problems involving real data, the near field can be estimated from the far field only up to a finite-dimensional approximation due to the ill-posedness of the near-to-far field map.However, it may be possible, even in that realistic computational framework, to make an estimate of the possibly nonconvex minimum source region from the far field data.The evanescent near field components decay a short distance (a wavelength or two) away from the source, so that the inverse diffraction of the far field back to the near zone remains reliable up to a short distance away from the source, and consequently the inference of some nonconvex support information is possible as we show in the paper.For far fields that are known in closed form, the near field can be computed uniquely up to the minimum convex source region and the minimum source region can be estimated subsequently via the proposed approach.
We wish to further comment on the fundamental and practical relevance of minimum source regions.The last decade has seen the development of many methods for imaging and shape reconstruction of scatterers from scattering data.The efforts on minimum source regions are an important part of the associated state-of-the-art for a number of reasons.First, the definition of the minimum source region of a far field is rigorous.It sheds light on the information about the support of a source or scatterer that is contained in far field data.Even though the inverse source problem exhibits general nonuniqueness, the minimum source region is a piece of information about the sought-after source that can be deduced uniquely from the far field.In relation to alternative imaging methods, this approach has the particular appeal that it applies in theory to the extraction of information about the support of a scatterer using a single probing field.In numerical scattering experiments, this method is quite robust [7,10,11], even though it defines only a bound for the true scatterer support.Finally, thanks to the generalization in [7], perhaps the main appeal of this theory is its universality, suggesting its implementation as both a theoretical tool and an algorithm in a broad class of inverse problems.
The rest of the paper is organized as follows.Section 2 presents, in the framework of the Helmholtz equation in 2D space, the key theoretical results for estimating bounds for the minimum convex source region as well as the possibly nonconvex minimum source region.First we address the estimation of convex bounds, in theory up to the minimum convex source region.This is done via an exterior inverse diffraction approach.Next we refine the estimate by allowing nonconvex bounds, in theory up to the actual minimum source region which can be nonconvex.This is done by means of a complementary, interior inverse diffraction approach.In Section 3 we consider examples motivated by inverse scattering that can be handled analytically.Of particular interest is the support estimation of cylindrical objects, within exact scattering theory, from the scattering amplitude of a single scattering experiment (for a single probing field).In Section 4 we discuss the practical computational implementation aspects of the inverse support methods derived in the paper and present the numerical counterpart of the analytical results obtained in Section 3. The obtained computational results show the practical robustness of minimum source region computations in inverse scattering.In addition, the provided examples nicely illustrate the UWSC sets theorized in earlier work [4].Section 5 provides concluding remarks and suggestions for future research.

Theory
2.1.Forward Problem.We consider electromagnetic radiation and scattering problems involving transverse magnetic  (TM  ) polarization for which radiation is governed by the Helmholtz equation in 2D free space, where   is the  component of the electric field,   is the  component of the current distribution that produces the field,  is the angular oscillation frequency,  is the free space permeability, and  is the wavenumber, defined by where  is the free space permittivity.In Cartesian coordinates r = (,) while in cylindrical coordinates r = (,) where  ≡ |r| and  is the azimuthal angle.This context is relevant, for example, in wave modeling of long cylindrical objects (relative to the wavelength  = 2/) that essentially depend only on the transverse coordinates.The source   can represent the transverse dependence of an active source (antenna) or an induced source such as a scatterer excited by incident fields.The radiated field obeying the radiation condition at infinity is well known to be given by [13, page 598] where  (1)  0 is the Hankel function of the first kind and order 0. From the addition theorem for Hankel functions we have  (1)  0 where  < = min(,   ) and  > = max(,   ) where  is the distance from r to the origin of coordinates and   is the distance from r  to the origin.For a source whose support  is confined inside a disk or 2-ball of radius  centered about the origin  (see Figure 2) we have from ( 3) and ( 4) where the multipole moments Using the large argument approximation for the Hankel function we get from ( 14) the far zone behavior where the far field radiation pattern so that in view of the orthogonality of the complex exponentials for integer It can be shown from ( 8) and ( 9), Faraday's law, and Poynting's theorem that the average radiated power  is given by which must be finite, implying that the multipole moments are square-summable.
In addition, for nontrivial far fields, there exists a source inside the minimum source region  min ⊆  (see Figure 2) that generates the given far field, corresponding to the given multipole moments   .By applying the above formulation to this source we obtain where  min is the radius of the smallest ball that is centered about the same origin and contains the minimum source region  min ⊆  of the given far field (see Figure 2).
On the other hand, if the origin of coordinates  is such that  is contained outside a circle of radius  int centered at the origin (see Figure 3), then the field in the interior of this circle of radius  int is given from ( 3) and ( 4) by where Furthermore, there is a source inside  min ⊆  that generates the same field outside  so that where  int min is the radius of the largest ball centered about the same origin that is disjoint to the minimum source region  min and is tangential to it in a way that defines a concave boundary of  min .This is illustrated in Figure 3.

Minimum Convex Source Region.
Let us derive next convex and nonconvex boundaries for the minimum source region of a given far field based on the results above.We tackle first the convex support information by means of a Picard test applied to an exterior inverse diffraction problem.
Equations ( 8) and ( 9) define the far field while the field for  =  >  is defined by (5); that is, Next, we inverse-diffract this boundary-value field to a more interior (smaller) circle of radius  =  < , requiring that the resulting field for  =  <  be square-integrable ( 2 ) over the unit circle, which is the Picard condition of this linear inverse problem.The reconstructed field in the circle of radius  =  <  that obeys the homogeneous form of (1) must be of the form The smallest circle onto which the field can be inversediffracted via ( 16) and ( 17) defines the sought-after minimum circular region  min associated with the given far field () or equivalently the multipole moments   and is defined as where inf means infimum.The following result follows from (19) and the asymptotic behavior (see [5, page 77]) Main Result 1.We have Here sup means supremum.The proof is given in Appendix A.
Main Result 1 (21) holds for any origin .The intersection of the minimum spherical source regions  min = {r ∈ R 2 :  ≤  min } for different origins defines a convex region where the minimum source region resides.It is theoretically possible for sufficiently many such origins to actually define the minimum convex source region  min,conv of the given far field.The method can be summarized as follows.12) can be used for any of the test origins   , it follows that we can uniquely associate, with the given far field, a valid radiated field that is well defined everywhere outside the respective convex region B min that contains the minimum convex source region  min,conv .We discuss next how to use this near field information to further refine our estimate of the minimum source region, in theory up to the possibly nonconvex minimum source region.

Inference of Nonconvex Support Information.
In order to sharpen the estimate of the minimum source region, including nonconvex regions, we consider an interior inverse diffraction problem in the near field, associated with a complementary geometry where the field data are known in the interior of a given circle (interior data) and the field is inverse-diffracted into outer concentric circles.If closed form expressions for the far field or the multipole moments are available, then it is theoretically possible to estimate also the associated near field.If, on the other hand, the data are experimental or computational and thus subject to noise and perturbations, then the near field can be estimated only up to an essentially finite-dimensional representation, for example, the well known truncation   ≃ 0 for  > , where  is the radius of the smallest circle that bounds the source support .We assume next that the near field outside the convex region B min can be estimated analytically and leave the practical question of estimating the near field computationally for other sections of this paper which address the computational implementation and examples.
Once we compute the convex region B min that contains the minimum convex source region  min,conv as well as the radiated field that it produces in its exterior, it becomes possible to estimate another, possibly nonconvex region that is contained inside B min and that contains the minimum source region.We assume that the field in the interior of a ball of radius  centered at the origin can be estimated, for example, as part of the computations giving B min .We use this interior field information to compute the radius  int min of the largest circle up to which said interior field can be inversediffracted.The geometry of this interior inverse diffraction can pierce or intersect concave portions of  min,conv , allowing  the extraction of the pending nonconvex support information of  min that is not contained in  min,conv .The concept is illustrated in Figure 4.
In particular, consider for a given origin , chosen outside the derived bounding region for the minimum convex source region, a ball centered at that origin and residing outside the bounding region for the minimum source region.It follows from ( 13) and ( 14) with the substitution of  by the derived bounding region that in the interior of the circular region in question the field obeys an expansion of the form (13).The field at a circle of radius  =  inside this region is then Inverse diffraction to a more exterior circle of radius  =  >  is done via which is meaningful so long as The largest circle of radius  int min onto which the field can be inverse-diffracted this way defines a bounding region for the minimum source region associated with the given far field () or equivalently the multipole moments   .From (24) The minimum source region  min ⊆  The proof uses the same general approach used to derive Main Result 1 in Appendix A and does not need to be (essentially) repeated.
The procedure to estimate the minimum source region can now be summarized as follows.
(1) Applying the exterior inverse diffraction approach of the previous subsection, compute the convex region B min where  min ⊆  min,conv ⊆ B min .For sufficiently many computational origins it is possible that the estimate B min =  min,conv .
(2) Consider next reference origins outside B min , say   ,  = 1, 2, . . ., , and evaluate the field (via the exterior inverse diffraction approach) over circles centered about those origins and having radius as large as possible but not intersecting the convex region B  .Let us assume that such computational circles correspond to radii   , with centers   .
(3) Using the value of the field on such circles compute the respective minimum radius  int min ().(5) Compute the intersection of  int and B min , which defines the revised estimate  min,est ⊇  min for the minimum source region which can be nonconvex, in particular, (6) The information used in the previous steps is the field outside the convex region B min .However, as the estimate  min,est becomes closer to  min than B min we need to repeat steps 2-4 for computational circles located outside the best current estimate  min,est .Computation of the field outside the best current estimate relies on the interior inverse diffraction approach explained in this section.The counterpart of step 5 associated with the second and following iterations is  min,est (current) =  min,est (previous) ∩  int .By continuing this process (see Figure 4) it is theoretically possible to obtain  min,est =  min , as desired.
As an example of the implementation of this methodology, we discuss in Appendix B the computation of the minimum convex source region and minimum source region associated with the far field of two point sources and generalize to  sources.In the following section we consider more examples, involving canonical scattering objects that can be addressed analytically.The complementary numerical counterpart of the same examples is presented in another section.
Mathematical Problems in Engineering 7

Plane Wave Scattering by a Dielectric Circular Cylinder.
To demonstrate the application of Main Result 1, we consider scattering by canonical objects for which analytical results are available.First we consider electromagnetic wave scattering by a lossless dielectric circular cylinder of radius  centered about the origin.It is assumed that the cylinder's length is so much larger than the wavelength that it can be modelled as being infinitely long along the -axis, which facilitates description of the corresponding wave radiation, propagation, and scattering via the 2D Helmholtz equation model.The relative permittivity and permeability of this cylinder are   and   , respectively.The free space wavenumber  is defined in (2).On the other hand, the dielectric wavenumber  1 =  √     .It is well known that the field scattered in free space by this dielectric circular cylinder upon excitation by a TM  plane wave of amplitude  0 traveling in the  direction can be written as ( [14], page 666) which corresponds to the far field radiation pattern defined in (9) with where    () ≡ (/)  () and  (1)  () ≡ (/) (1)   ().We compute next the minimum convex source region of the far field radiation pattern defined in ( 9) and (30).For this purpose we use Main Result 1.Using (20) and (26) and Stirling's formula we get from (30) that lim where It follows from these results and Main Result 1 that the minimum source region of the far field scattering amplitude associated with this scattered field is a point, in particular, the center of the circular cylinder (corresponding to  = 0) that scatters the incident plane wave.This is a somewhat surprising result.It indicates that the center of the scatterer can be deduced unequivocally from the far field, but the information about the size of the scatterer is very limited since a very small source could produce the same far field, at least in theory (disregarding practical physical limitations of small antennas).Before we discuss its interpretation, which will naturally lead us to a more practical numerical computation counterpart of the theory and algorithms considered up to this point, we wish to briefly consider a variant of the scattering cylinder above and a generalization of the results above to near field incident wave excitation, which will prove to be quite insightful.

Scattering by a Conducting Circular Cylinder.
In particular, for the case of a circular cylinder made of a perfect electric conductor (PEC) we have that the scattered field   is given by (29) with the multipole moments given by ( [14], page 604) and by manipulations analogous to those leading to (32) we again obtain lim which means that, as in the dielectric case, the minimum source region of the far field scattering amplitude associated with the field scattered by a PEC circular cylinder upon plane wave incidence is simply the center of the cylinder.

Line Source Scattering by a PEC Circular Cylinder.
In the previous two examples we assumed that the incident wave is a plane wave.Next, we generalize the results by considering the tractable case of scattering of the field due to an infinitely long antenna (in the respective 3D picture) or a point source located outside the PEC circle (in the equivalent 2D framework).In particular, the field scattered by the PEC cylinder upon incidence of the field, due to an infinite line of constant current   located at position r  given in cylindrical coordinates as r  = (  ,   ), is given by ([14], section 11.5.5) where (1)   () (1)   (  )    (37) so that lim In view of Main Result 1, this means that the smallest ball circumscribing the minimum source region is the ball  min = {r ∈ R 2 :  ≤  2 /  }.In other words, the largest dimension  of the minimum source region must be  = 2 2 /  .Note that if   → ∞ as is the case for the plane wave excitation of the previous section, we obtain  = 0 corresponding to a multipole source localized in the origin as we found in the previous section.

Example of the Union of Well-Separated Convex Sets.
The result above gives a first estimate of the minimum PEC cylinder Source Figure 5: Estimated UWSC sets (shown in black), for the total far field due to a source (infinite line source in 3D space or equivalently point source in the 2D formulation) and a circular cylinder (PEC or dielectric) in its vicinity.The estimate for the minimum convex source region of the scatterer source depends on the proximity to the radiating source.The closer it is, then the larger the minimum source region needs to be in order to be able to radiate the corresponding scattering field.
convex source region of the scatterer.Furthermore, the field due to the combined line source and scatterer is the linear superposition of the fields of the line source and the source induced in the scatterer.In 2D space the minimum source region of the line source is the single point r  = (  ,   ), where the source is located.The field due to the point source cannot be annulled by that of the scatterer source since it does not fully cover (enclose) the point source.In particular, it is not hard to show via Green's theorem that this field and the one due to the scatterer are linearly independent.Therefore the minimum source region of the field due to both is the union of the minimum source regions of the point source and the scatterer source individually.It follows that we have an estimate for the minimum source region of this total field as given by the union of the convex set  min = {r ∈ R 2 :  ≤  2 /  } and the single point r  = (  ,   ).This is an exciting example of the existence and analytical computation of the UWSC sets of a far field whose existence was predicted in [4].This result is illustrated in Figure 5.

Infinite-versus
Finite-Dimensional Approaches.The analytical inverse support results above establish a fundamental constraint on the support of any source that produces the given scattered field.In the case of the plane wave excitation, which corresponds also to the limit   → ∞, we found that the analytical minimum source region is the center of the scattering cylinder, a result that is quite surprising.On the other hand, for excitation by a line source at (  ,   ) we found that, for both PEC and dielectric cylinders, the minimum spherical source region has the finite radius  2 /  .The larger minimum source region can be thought of as representing greater information about the source support that is accessible in the far field when the scatterer is excited by a near field source, which intuitively makes sense.These results establish bounds for the true scatterer support.In the present case they define for far field excitation the center of the scatterer and for near field probing a larger region which must be in the interior of the actual scatterer responsible for the given scattered field data.It is also important to emphasize that these results apply to the infinite-dimensional multipole domain and are thus relevant only for perfect far field data.For a source whose radius is as large as the predicted minimum radius, the terms in the infinite series expansion for the source energy decay sufficiently fast at infinity, so that the resulting source energy is not infinite; however, it can be prohibitively large.Thus strictly speaking, the source energy is not infinite; however, it can be quite large for a source in the predicted minimum region.Furthermore, it can be shown that regardless of   the number of essential dimensions of the multipole domain is governed by the ruleof-thumb formula   ≃ 0,  > , so that, by plotting   versus  and estimating the knee of the plot, one can approximately estimate via this rule-of-thumb formula the circular cylinder's radius.We experimented with this finitedimensional approach and found it to be quite robust.

Numerical Illustration
The numerical counterpart of Main Results 1 and 2 employs finite-dimensional field data and a finite number of test origins.For a reference origin, for example, we compute the far field radiation pattern and the corresponding multipole moments via (10).We keep only the nonnegligible multipole moments and take the rest to be zero.We found that, in the examples to be discussed in the following, the number of significant multipole moments corresponds approximately to the truncation   ≃ 0,  > , where  is the radius of the smallest circle that circumscribes the original source support and is centered at the same origin.This step makes the far field pattern and multipole moment data finite-dimensional for the remaining computations.Main Results 1 and 2 are implemented requiring that the functional energy of the reconstructed field, as defined in ( 18) and (24), exhibit a reasonable finite value.We experimented with several ways of implementing this cutoff and settled in an approach resembling the familiar L-curve method.In particular, we automatically estimated the peak curvature point of the "log(energy) versus radius" curve and used this peak curvature point as the numerical estimate of the minimum radius  min in the convex support estimation (Main Result 1) or of  int min in the complementary nonconvex support estimation (Main Result 2).Next we discuss this computational implementation with two examples.The first illustrates the use of support information in an inverse source problem.The second sheds light on the numerical counterpart of the scattering PEC cylinder case discussed in Section 3.

Inverse Source Problem with Support
Estimation.As a first numerical example, we consider the inverse source problem of estimating an unknown radiating source from knowledge of the far field radiation pattern it generates.As part of the inversion, we incorporate estimation of the minimum source region consistent with the given far field.The support information is used as prior for the reconstruction of the minimum energy source that resides in that support and gives the desired far field data.To estimate the minimum source region, we focus first on the estimation of a convex bound for the minimum source region.Later we refine the estimate of the minimum source region up to a nonconvex region, via the methodology presented in the paper.The final result is that the reconstructed source is quite close to the original source generating the far field data, as desired.
Figure 6 shows the source used in the simulations.It is formed by two 0.2×0.2squares, where  is the wavelength.The total source is contained inside a  ×  square and thus has maximal dimension √ 2.In the following all dimensions are given as multiples of .The source distribution inside the squares is uniform and has unit amplitude.Figure 7 shows the polar diagram of the amplitude of the far field radiation pattern generated by this source.First, some insight about the source to be reconstructed and its near field can be gained by backpropagating the far field to the near zone.Figure 8 illustrates this backpropagated field.The peaks in this plot offer qualitative insight about the location and size of the two square sources composing the total source.On the other hand, the near field associated with this backpropagation is not the correct near field (as we show next).Furthermore, the value of the source cannot be deduced from this qualitative insight alone.In the following, we implement the full reconstruction of the source by means of a three-step procedure based on the methods developed in this paper.First we estimate the minimum convex source region and compute the minimum energy source associated with this estimated region.Next we compute the near field based on the reconstructed minimum energy source and use this information to obtain a tighter, nonconvex estimate of the minimum source region.Finally, we compute the minimum energy source associated with this nonconvex estimate.The thus reconstructed source is very similar to the original source used in generating the far field data, as desired.
To estimate the minimum convex source region, we considered 6 test origins: (−10, 0), (10, 0), (0, −10), (0, 10), (−10, 10), and (10, −10).Figure 9 shows the computed estimate of the minimum convex source region obtained using these test origins.Also shown in the figure is the true source support, which is contained inside the convex source The high amplitude region of this source (value near 1) coincides very well with the correct source support of the original unitamplitude source, as desired.region estimate, as expected.Using this prior source support information, we computed next the minimum energy source contained in this source region that is consistent with the far field data.This is done by standard matrix inversion giving the solution of minimum 2 norm to the inverse problem of finding a source in the given support that produces the given data.The contour plot of the reconstructed source is given in Figure 10.The source regions where the amplitude is significant (close to 1) are seen to coincide quite well with the true source support.Furthermore, this step also allows us to estimate the near field, by computing the radiation associated with the reconstructed minimum energy source.Figure 11 illustrates the near field associated with this source.
The figure shows both the computed near field and the correct near field based on the original source, and the two coincide very well, as desired.This estimate of the near field is used in the following step, which seeks to refine the estimate of the minimum source region, up to its possibly nonconvex nature.Using the estimated near field, we implemented the method presented in the paper to extract further information about the possibly nonconvex minimum source region.Using two test origins, (1.2, −.3) and (−.3, 1.2), we obtained the better, nonconvex estimate of the minimum source region shown in Figure 12.As a final step in the reconstruction, we did the standard linear inversion leading to the minimum energy source that is consistent with the far field data and is confined within this better estimate of the minimum source region.The reconstructed source is shown in Figure 13.The source strength is clearly focused on the vicinity of the true source support, and the value of the source is near unity in this region, as desired.Thus this three-step methodology, exploiting the extracted source support information via the multipole inverse theory established in the paper, is quite practical, as well as insightful about the possibilities and limitations of inverse source problems with far field data.
For comparison purposes, we present in Figure 14 a plot of the minimum energy source that would have been obtained in the same computational grid if no prior information about the source support would have been sought during the inversion process.Note that this source is quite similar to the backpropagated field in Figure 8 (which is not surprising in view of the multipole expansion).The reconstruction in Figure 13, which seeks and exploits support information from far field data, is clearly superior.

Numerical Scattering Example.
As a second numerical case, we considered the numerical computations associated with the scattering PEC cylinder discussed in Section 3. First we addressed the estimation, from the far field, of the minimum convex source region associated with the source that is induced in a PEC cylinder due to different excitations.We found that very accurate estimates of the scatterer support can be obtained via the computational counterpart of Main Result 1 for scatterers whose radius  is greater than about a couple of wavelengths.For smaller cylinders, the method gives excellent reconstruction of the section of the cylinder that is illuminated but gives a more ambiguous boundary for the nonilluminated portion of the cylinder, which cuts part of that half of the cylinder.A more interesting problem consists of estimating bounds for the total source support, including both the excitation point source and the scattering PEC cylinder.Next we focus on this problem, illustrating the application of the computational counterparts of Main Results 1 and 2. In the following calculations, we begin with the convex estimate of the minimum source region.Later we refine the bounding region for the minimum source region using the complementary nonconvex boundaries.We use a PEC cylinder of radius  = 5 (where, as before, all dimensions are given as multiples of ).We consider two particular cases: (1) an excitation source at position (  ,   ) = (−10, 0) and (2) two excitation sources, at positions (−15, −5) and (15, −5).
For the case of a single excitation source at (−10, 0), we employed the following 4 test origins to obtain a first convex boundary for the total source (the point source plus the scatterer): (−50, 0), (50, 0), (10,50), and (10, −50).Figure 15 shows the computational sequence corresponding to these test origins.The plots show the correct source support, the estimated boundary of the minimum circle of the current for the nonconvex part.We tried other origins but concluded that only marginal support information was available beyond this point.The nonconvex part, in particular, could not be used to isolate the supports of the point source and the scatterer.

Conclusion
We developed a theoretical and computational framework to exactly or approximately estimate the minimum source region of a far field.The methodology constitutes the multipole theory counterpart of the plane wave expansion approach of a previous paper [2].The proposed method allows handling of problems for which the far field is given as a multipole series expansion (e.g., cylindrical scatterers), for which the method in [2] cannot be used since it requires knowledge of the far field pattern versus the observation direction in closed form.We found that, for canonical cylinders and under plane wave excitation, the analytical results define only the center of the scatterer.On the other hand, under near field excitation, the boundary of the minimum source region associated with the induced source in the scatterer is a circle of finite radius inside the scatterer.This can be thought of as rigorous quantification of the greater information about the scatterer available with near field excitation.Furthermore, due to reciprocity, the role of the (near field) excitation and the (far field) receiver can be interchanged, so that the same results also apply for plane wave excitation (far field excitation) and sensing of the scattered field in the near zone.Thus the same results, with this role reversal, coincidentally show that the same developments derived in the paper can be applied to the companion problem of estimating minimum source regions from near field data.This, in turn, further demonstrates the broad applicability of the ideas about minimum source regions considered in this paper and in previous work.
The analytical results can be thought of as a quantification of the scatterer support information that is inferrable from the far scattered field.They show that this information is, in theory, critically affected by the proximity of the excitation Mathematical Problems in Engineering source to the scatterer.The provided analytical results also illustrate the UWSC sets proposed in a previous paper.The analytical results are fundamental (they are strict bounds to the minimum source regions), but at least in the particular case of the PEC and dielectric cylindrical scatterers considered in the paper they are also quite restrictive if the excitation source is far from the scatterer (the minimum source region becomes very small as   becomes larger, relative to the support of the true scatterer responsible for the data).This limitation naturally led us to the more practical computational implementation of the ideas derived in the paper, where the field data are essentially finite-dimensional and the realizability constraint becomes a bounding of the functional energy.This regularization led to results such as those shown in Section 4, which are more conservative but also more stable.
The provided analytical and numerical inverse support results based on the methods derived in the paper were quite encouraging.Thus the proposed methods are useful for the particular sources we considered, and we suspect that they work for a broad class of sources.We acknowledge, however, that the nonconvex bounds may not be correct for some sources.We have not specified in this paper the class of sources for which these methods (the nonconvex bounds, in particular) are expected to be applicable in practice.This pending question is related to issues raised in [4] and is left for the future.
We also wish to comment that, since the estimation of the near field from the far field is ill-posed, the nonconvex support estimation component of the proposed methodology (which relies on such near field estimation) is clearly ill-posed as well.Yet we have established the basic inverse-diffractionbased construction procedure which shows that it is theoretically possible to compute the possibly nonconvex minimum source region.The proposed methodology complements the previous efforts in this area which had emphasized only the minimum nonconvex source region.Furthermore, being aware of the practical challenge of the required near field estimation, we have examined the practical applicability of the proposed technique in the context of two numerical examples.In both cases we were able to extract some nonconvex support information, as desired.The usefulness of the extra nonconvex support information was illustrated in the numerical example related to the inverse source problem, where we were able to obtain a solid reconstruction of the source thanks to the nonconvex support estimation component of the proposed methodology.
We are currently working on a purely computational variant of the methods derived in this paper, which employs manifolds of realizable field data, corresponding to bounded functional energy sources.This approach offers insight onto the observed robustness of minimum source regions in inverse scattering formulations (as shown in the numerical simulations of this paper and those of other studies).We are also exploring the possibility of extending the derived methods to characterize and compute minimum source regions in inverse problems with partially correlated fields which arise in optical coherence theory.Finally, we have considered here the scalar formulation.Extension to the full vector electromagnetic fields, based on the corresponding vector multipole fields, appears at first sight equally tractable and provides another important avenue for future research.We plan to report our research findings in these and related areas elsewhere.

2 MathematicalFigure 1 :
Figure 1: Minimum source region  min and the minimum ball  min () for a given origin .

Figure 2 :Figure 3 :
Figure 2: Relations between  min and  and  min and  for a given origin .
(1) Consider a number  of origins of coordinates,   ,  = 1, 2, . . ., , within the region of interest.(2) For each origin, compute via the test (21) the associated minimum radius  min (  ) and associated minimum circular region  min (  ) = {r ∈ R 2 : |r −   | ≤  min (  )}.(3) The intersection of the regions  min (  ),  = 1, 2, . . ., , defines a convex region B min bounding the minimum convex source region  min,conv in the sense that  min,conv ⊆ B min .The estimate B min becomes closer to  min,conv as one takes more sample origins and their corresponding minimum circular regions.Now, since (

Figure 4 :
Figure 4: Successive application of the diffraction-based estimation of the critical sphere of radius  intmin , based on field information available from previous iterations.This method allows the estimation of concave boundaries of the minimum source region.Thus it allows, in theory, the estimation of the minimum source region  min which can be nonconvex.

Figure 6 :
Figure 6: Contour plot of the unit-amplitude source used in the simulations.

Figure 7 :Figure 8 :
Figure 7: Polar diagram of the radiation pattern amplitude.

Figure 11 :
Figure 11: (a) Contour plot of the magnitude of the near field radiated by the minimum energy source.The support of this source is the minimum convex source region estimate corresponding to test origins (−10, 0), (10, 0), (0, −10), (0, 10), (−10, 10), and (10, −10) (shown in black color).(b) Contour plot of the magnitude of the near field radiated by the original source.The two plots are almost identical, as desired.

Figure 12 :Figure 13 :
Figure 12: Estimated nonconvex region bounding the minimum source region.It was obtained using test origins (1.2, −.3) and (−.3, 1.2).Also shown are the true source support and the boundary of the minimum source region associated with test origin (−.3, 1.2).

Figure 14 :
Figure 14: Amplitude of the minimum energy solution to the inverse source problem without support information.

Figure 18 :Figure 19 :
Figure 18: Cummulative estimated nonconvex support based on 8 test origins for the convex part and 4 test origins for the nonconvex part.

Figure 20 :
Figure 20: Cummulative estimated nonconvex support based on 6 test origins for the convex part and 5 test origins for the nonconvex part.