Reconstruction via interior transmission eigenfunctions

# Reconstruction via the intrinsic geometric structures of interior transmission eigenfunctions

Jingzhi Li Department of Mathematics, Southern University of Science and Technology, Shenzhen, China. Xiaofei Li Department of Mathematics, Inha University, Incheon, South Korea.  and  Hongyu Liu Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Hong Kong SAR. and HKBU Institute of Research and Continuing Education, Virtual University Park, Shenzhen, P. R. China.
###### Abstract.

We are concerned with the inverse scattering problem of extracting the geometric structures of an unknown/inaccessible inhomogeneous medium by using the corresponding acoustic far-field measurement. Using the intrinsic geometric properties of the so-called interior transmission eigenfunctions, we develop a novel inverse scattering scheme. The proposed method can efficiently capture the cusp singularities of the support of the inhomogeneous medium. If further a priori information is available on the support of the medium, say, it is a convex polyhedron, then one can actually recover its shape. Both theoretical analysis and numerical experiments are provided. Our reconstruction method is new to the literature and opens up a new direction in the study of inverse scattering problems.

Key words. Inverse medium scattering, non-destructive reconstruction, polyhedral object, interior transmission eigenfunctions

AMS subject classifications. 35P25, 78A46, 35R30, 65N25

## 1. Introduction

This work is concerned with the inverse medium scattering problem of extracting the geometric structures of a penetrable scatter by the corresponding acoustic wave detection. Let , be a bounded Lipschitz domain with a connected complement . Assume that the refractive index of the background space is , and the refractive index of is , where is the sound speed in , and is the sound speed in . We assume that the refractive index is always positive, and let .

Given an incident field , the presence of the scatter will give rise to a scattered field . We take to be a time-harmonic plane wave, where , is the incident direction, and is the wave number. We define to be the total field, which satisfies the following Helmholtz system

 (1.1)

The last limit in (1.1) holds uniformly in all directions and is known as the Sommerfeld radiation condition.

The system (1.1) is well understood, and it is known that there exists a unique solution . Moreover, is (real) analytic and the asymptotic behaviour of the scattered field is governed by

 us(x)=eik|x||x|(l−1)/2u∞(^x)+O(1|x|(l+1)/2),  |x|→∞,

uniformly in all directions . The analytic function is defined on the unit sphere and called the far-field pattern, which is given by

 u∞(^x)=k24π∫Rle−ik^x⋅yq(y)u(y)dy, (1.2)

with known as the observation angle/direction (see ). In what follows, we write to specify its dependence on the observation direction , the incident direction, and the wave number . In the current article, we are concerned with a geometrical inverse problem of reconstructing the shape/support of the scatterer, namely , disregarding its content , by knowledge of the measurement of . This problem has been playing an indispensable role in many areas of sciences and technology such as radar and sonar, medical imaging, geophysical exploration, and nondestructive testing (see, e.g., ).

A number of numerical reconstruction methods have been developed for the aforementioned inverse scattering problem in various scenarios, including the linear sampling method, factorization method, MUSIC-type methods, time reversal method, single-shot method and topological-optimization-type method; we refer the readers to [1, 2, 3, 4, 5, 6, 7, 12, 16, 19, 13, 28, 31, 37] and the references therein for these methods and some other related developments. In particular, a type of so-called interior transmission eigenvalue problem arises in the study of linear sampling and factorization methods. These methods succeed only at wave numbers which are not the transmission eigenvalues. The study of the transmission eigenvalue problem is mathematically interesting and challenging since it is a type of non-elliptic and non self-adjoint problem. Recently, the transmission eigenvalue problems were discovered to be useful in inverse medium scattering problem. In , a lower bound on the index of refraction is obtained by knowledge of the smallest real transmission eigenvalue corresponding to the medium. In , transmission eigenvalues are used to the nondestructive testing of dielectrics. The aforementioned studies as well as the other existing ones, to name a few [23, 35, 36], focus on the intrinsic properties of transmission eigenvalues. In the recent papers [8, 9, 10], the intrinsic geometric structure of interior transmission eigenfunctions are studied. Specifically, the vanishing and localizing properties of the interior transmission eigenfunctions at cusps on the support of the refractive index are investigated both theoretically and numerically. Finally, we would like to mention in passing that the study of interior transmission eigenvalue problems was also connected to that of invisibility cloaking in the recent papers [27, 30, 32].

Motivated by the aforementioned works on the intrinsic geometric properties of the interior transmission eigenfunctions, we shall develop a novel scheme for the inverse scattering problem of extracting the geometric structures of an unknown and inaccessible inhomogeneous medium by using the corresponding acoustic far-field measurement. The rationale of the proposed method can be briefly described as follows. Consider the following interior transmission eigenvalue problem associated to ,

 ⎧⎪⎨⎪⎩Δu+k2(1+q)u=0 in D,Δv+k2v=0 in D,u−v∈H20(D), (1.3)

where it is noted that if is a Lipschitz domain

 H20(D)={v∈H2(D):v=0, ∂νv=0 on ∂D},

with signifying the unit normal vector directed into the exterior of .

###### Definition 1.1.

A number for which the transmission problem (1.3) has nontrivial solutions such that is called an interior transmission eigenvalue associated with . The nontrivial solutions are called the corresponding interior transmission eigenfunctions.

The existence, discreteness and infiniteness of transmission eigenvalues for (1.3) can be found [20, 23, 18, 17] and the references therein. However, there are few results concerning the intrinsic properties of the transmission eigenfunctions. The intrinsic geometric structures of interior transmission eigenfunctions were recently investigated in [8, 9, 10]. Specifically, it is shown that if there are cusp singularities on in (1.3), then the interior transmission eigenfunctions would reveal a certain intrinsic quantitative behaviours near the cusps. Here, by a cusp we mean a point where the tangential vector field of the boundary is discontinuous. Corner singularities are very typical cusp singularities. To be more precise, in [8, 10], it is shown that if there is a cusp on the support of , then the transmission eigenfunction vanishes near the cusp if its interior angle is less than , whereas the transmission eigenfunction localizes near the cusp if its interior angle is bigger than . Furthermore, it is shown that the vanishing and blowup orders are inversely proportional to the interior angle of the cusp: the sharper the angle, the higher the convergence order.

The proposed reconstruction method in this paper is first to make use of the far-field data to determine the interior transmission eigenvalue as well as the corresponding transmission eigenfunctions. We are aware of some existing study by making use of the so-called inside-outside duality to determine the transmission eigenvalues . But for our need, we shall also need to determine the corresponding eigenfunctions. The determination of transmission eigenvalue and corresponding eigenfunctions is based on the far-field regularization techniques. To our best knowledge, the study in this aspect is also new to the literature. After the determination of the transmission eigenvalue, we seek the Herglotz wave function in a certain domain which is the approximation of the corresponding transmission eigenfunction. The places where the Herglotz wave function is vanishing or localizing are the locations of those cusp singularities of the support of the medium scatterer. If further a priori information is available on the support of the medium, say, it is a convex polyhedron, then one can actually recover its shape by simply joining the cusp singularities by line. If not much a priori information is given, at least we can recover the locations of the cusp singularities. Extensive numerical experiments show that our reconstruction method is efficient and effective. Our study is first of its type in the literature and opens up a new direction in the study of inverse scattering problems.

The rest of the paper is organized as follows. In Section 2, we review some theoretic results on the properties of transmission eigenvalues and eigenfunctions. In Section 3, we focus on the property of the kernel of far-field operator and the Herglotz approximation to the interior transmission eigenfunctions. In Section 4, we present the recovery scheme in detail. Section 5 is devoted to numerical experiments to validate the applicability and effectiveness of the proposed method. The paper is concluded in Section 6 with some discussion.

## 2. Preliminaries on transmission eigenvalue problem

In this section, we first collect some theoretical results for the interior transmission eigenvalue problem (1.3). Throughout this paper, we assume that the refractive index is real-valued. Denote by

 n∗=infx∈Dn(x),n∗=supx∈Dn(x).

The following theorem in  gives the existence of interior transmission eigenvalues.

###### Theorem 2.1.

Let satisfy either one of the following assumptions:

1. ,

for some constants and . Then there exists an infinite set of interior transmission eigenvalues with as the only accumulation point.

Recall from  that the lower bound of the first transmission eigenvalue has the following estimation.

###### Theorem 2.2.

Let be the radius of the smallest ball containing . Denote by and the first positive transmission eigenvalue corresponding to the ball of radius with index of refraction and , respectively. Then the first transmission eigenvalue corresponding to and the given index of refraction has the following estimations:

1. If for some constant , then

 k1≥max⎛⎝k1,n∗r,√λ1(D)n∗⎞⎠; (2.1)
2. If for some constant , then

 k1≥max(k1,n∗r,√λ1(D)), (2.2)

where is the first Dirichlet eigenvalue for in .

The numerical experiments in  indicate the existence of complex transmission eigenvalues, but this has not been established theoretically in general. The following theorem shows the non-existence of purely imaginary transmission eigenvalues .

###### Theorem 2.3.

If or almost everywhere in , then there exists no purely imaginary transmission eigenvalues.

Next, we give a more definite description of the vanishing and localizing of transmission eigenfunctions from [8, 10].

###### Definition 2.1.

Assume that is a transmission eigenvalue, then there exist such that

 ⎧⎪⎨⎪⎩Δu+k2(1+q)u=0in D,Δv+k2v=0in D,u−v∈H20(D), ∥v∥L2(D)=1.

Let be a point and be a ball of radius centered at . Set . Assume that

 ∥q∥L∞(Dr(p))≥ϵ0,ϵ0∈R+.

Then we say that vanishing occurs near if

 limr→+01√|Dr(P)|∥v(x)∥L2(Dr(P))=0,

whereas we say that localizing occurs near if

 limr→+01√|Dr(P)|∥v(x)∥L2(Dr(P))=+∞,

where signifies the area or volume of the region , respectively, in two and three dimensions.

The vanishing and localizing property of transmission eigenfunctions near a corner point comes from the corner scattering study in [9, 11, 33, 26, 24, 25]. If is the vertex of a corner with an interior angle less than , the vanishing of the transmission eigenfunctions has been rigorously verified in . An important consequence of the study in [9, 11, 33, 24] is the fact that the transmission eigenfunction cannot be analytically extended across a corner point to form an entire solution to the Helmholtz equation, in . However, we note that due to the interior regularity, the transmission eigenfunction is always analytic away from the corner point. Hence, heuristically, the failure of the analytic extension may indicate that either vanishes or blows up when approaching the corner point. Clearly, the failure of the analytical extension should also hold across any irregular point on the support of , and hence in , it is numerically shown that the vanishing or localizing behaviours of the transmission eigenfunctions would occur near any cusps on the support of the underlying. Theoretical proof of such a conjecture is fraught with significant difficulties. Indeed, the proof in  of the vanishing of the transmission eigenfunction in the special case near a corner with an interior angle less than already involves much technical analysis and advanced tools. In , the order of vanishing and localizing convergence rate has been investigate numerically, and show that it is related to the angle of the corner. In the three dimensional case, it turns out that edges also posses the vanishing phenomena.

## 3. Property of the kernel of far-field operator and the Herglotz approximation to the interior transmission eigenfunction

In this paper, we are concerned with the inverse medium problem of determining the support of unknown domain with cusp singularities based on the vanishing and localizing property of transmission eigenfunctions. Consider the scattering problem (1.1). We first introduce the far-field operator : defined by

 Fk(g)(^x):=∫Sl−1u∞(^x,d,k)g(d)ds(d), ^x∈Sl−1.

The far field operator is injective and has dense range if and only if there does not exist a Dirichlet eigenfunction for which is a Herglotz wave function (see, for example, ).

###### Definition 3.1.

A Herglotz wave function is a function of the form

 Hk(g)(x)=∫Sl−1eikx⋅dg(d)ds(d),  x∈Rl, (3.1)

where . The function is called the Herglotz kernel of .

We recall from  the following result concerning the Herglotz approximation .

###### Theorem 3.1.

Let denote the space of all Herglotz wave functions of the form (3.1). Define, respectively,

 Wk(D):={u∈C∞(D):(Δ+k2)u=0},

and

 Wk(D):={u|D:u∈Wk}.

Then is dense in with respect to the topology induced by the -norm.

If happens to be the interior transmission eigenvalue to the following interior transmission problem

 ⎧⎪⎨⎪⎩Δu+k2(1+q)u=0in D,Δv+k2v=0in D,u−v∈H20(D). (3.2)

By the denseness property of Herglotz wave functions, see Theorem 3.1, there exists Herglotz wave function

 vg=∫Sl−1eikx⋅dg(d)ds(d),

where the kernel , such that for any sufficiently small , there is

 ∥vg−v∥H1(D)≤ϵ. (3.3)

We say that is normalized if . Then extends to the whole as the Herglotz wave function . It is shown that transmission eigenfunctions cannot be extended analytically to a neighbourhood of a corner (see, for example, ). Therefore, the transmission eigenfunctions must vanish near the corner point.

In this paper, we aim at finding the support of when it is supported in a polyhedral domain. If domain is not of polyhedral shape but possesses cusp point, we aim at locating cusp singularities. The main idea is to find the vanishing and localizing points of the Herglotz wave extension function of transmission eigenfunctions. To do so, we first need to retrieve the transmission eigenvalue from a knowledge of far-field pattern .

Consider the scattering problem (1.1). Let be the simple incident plane wave, where is the incident direction. The far field pattern of the scattered field can be collected at circular boundary which enclose .

By superposition we note that

 Fk(g)(^x)=∫SN−1u∞(^x,d,k)g(d)ds(d),  ^x∈SN−1 (3.4)

is the far-field pattern of the scattered field

 vs(x)=∫SN−1us(x,d)g(d)ds(d)

corresponding to the incident field (see, for example, ). For , it has the Fourier expansion

 g(x)=∞∑n=0n∑m=−namnYmn(^x), (3.5)

where denotes a normalized spherical harmonic, and the series is uniformly and absolutely convergent. The coefficients are given by

 amn=∫Sl−1g(d)¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯Ymn(d)ds(d).

By (1.2) and the reciprocity relation , see, for example , we have

 Fk(g)(^x)=∫Sl−1u∞(−d,−^x)g(d)ds(d)=k24π∫Sl−1∫Beikd⋅yq(y)u(y,−^x)dyg(d)ds(d), (3.6)

where is a bounded open ball with radius enclosing scatter . Inserting (3.5) into (3.6), then (3.6) becomes

 Fk(g)(^x)=k24π∫Sl−1∫Beikd⋅yq(y)u(y,−^x)dy∞∑n=0n∑m=−namnYmn(d)ds(d)=k24π∫Bq(y)u(y,−^x)∞∑n=0n∑m=−namn∫Sl−1eikd⋅yYmn(d)ds(d)dy. (3.7)

By the Funk-Hecke formula (see for example ), (3.7) becomes

 Fk(g)(^x)=k24π∫Bq(y)u(y,−^x)∞∑n=0n∑m=−namn4πinjn(k|y|)Ymn(y/|y|)dy, (3.8)

where denotes a spherical Bessel function.

With the aid of Stirling’s formula , we obtain

 jn(kr)=O(ekr2n)n, n→∞ (3.9)

uniformly on . Then we have for any , there holds

 Fk(g)(^x)≈0

for sufficiently large , by letting , .

In order to distinguish the far-field operator between two cases that is a transmission eigenvalue and is not a transmission eigenvalue, we rigorously prove the following three crucial theorems in detail.

Define

 gN:=N∑n=0n∑m=−namnYmn(^x), (3.10)

and

 Fk,N(g):=∫Sl−1u∞(^x,d,k)gN(d)ds(d).

We prove the following theorem.

###### Theorem 3.2.

If is sufficiently large, then the following holds

 ∥Fk(g)−Fk,N(g)∥L2(Sl−1)≤O(ek2N)N,

where .

###### Proof.

Let be a bounded open set enclosing scatter . Without loss of generality, we assume that is a ball centered at with radius , i.e., .

By (1.2) and (3.10), can be expanded as

 Fk,N(g)(^x)=k24π∫Bq(y)u(y,−^x)N∑n=0n∑m=−namn4πinjn(k|y|)Ymn(y/|y|)dy. (3.11)

From (3.8) and (3.11) we have

 Fk(g)(^x)−Fk,N(g)(^x)=k24π∫Bq(y)u(y,−^x)∞∑n=N+1n∑m=−namn4πinjn(k|y|)Ymn(y/|y|)dy.

Then

 ∥Fk(g)−Fk,N(g)∥2L2(Sl−1)=k416π2∫Sl−1∣∣ ∣∣∫Bq(y)u(y,−^x)∞∑n=N+1n∑m=−namn4πinjn(k|y|)Ymn(y/|y|)dy∣∣ ∣∣2ds(^x)≤k416π2∫Sl−1supy∈B|q(y)u(y,−^x)|∞∑n=N+1n∑m=−n16π2|amn|2∫R0rl−1|jn(kr)|2drds(^x)≤k4M|Sl−1|∞∑n=N+1n∑m=−n|amn|2∫R0rl−1|jn(kr)|2dr,

where , is the area of . Due to (3.9), we have

 ∥Fk(g)−Fk,N(g)∥2L2(Sl−1)≤k4M|Sl−1∥g∥2L2(Sl−1)supn>N∫R0rl−1|jn(kr)|2dr≤k4M|Sl−1|supn>N∫R0rl−1|jn(kr)|2dr=k4M|Sl−1|supn>NO(ek2n)2n∫R0r2n+l−1dr≤k4M|Sl−1|O(ek2N)2NR2N+l2N+l=O(ek2N)2N (3.12)

for sufficiently large. The proof is done.

Let denote the set of all the interior transmission eigenvalues. We have the following theorem.

###### Theorem 3.3.

If with respect to the transmission eigenvalue problem (1.3) , then for and a sufficiently large , there holds

 ∥Fk,N(g)∥L2(Sl−1)≤Cϵ+O(ek2N+1)N+1, (3.13)

where .

###### Proof.

Let , and are the corresponding eigenfunctions. For any sufficiently small , by the denseness property of Herglotz functions, see Theorem 3.1, there exists Herglotz wave function (3.3) with , such that for any sufficiently small , there is

 ∥vg−v∥H1(D)≤ϵ. (3.14)

Denote by , and define

 Hk,N(g)(x)=∫Sl−1eikx⋅dgN(d)ds(d),  x∈Rl.

By (3.5), (3.10) and the Funk-Hecke formula, and can be expanded as

 Hk(g)(x)=4π∞∑n=0n∑m=−namninjn(k|x|)Ymn(^x) (3.15)

and

 Hk,N(g)(x)=4πN∑n=0n∑m=−namninjn(k|x|)Ymn(^x),

for all . Then

 Hk(g)(x)−Hk,N(g)(x)=4π∞∑n=Nn∑m=−namninjn(k|x|)Ymn(^x).

Similarly to Theorem 3.2, let be a bounded open set enclosing scatter . Without loss of generality, we assume that is a ball centered at with radius , i.e., . Then

 ∥Hk(g)−Hk,N(g)∥2L2(B)=∞∑n=Nn∑m=−n16π2|amn|2∫R0rl−1|jn(kr)|2dr≤16π2|amn|2maxn≥N∫R0rl−1|jn(kr)|2dr. (3.16)

Since

 ∫Rrl−1|jn(kr)|2dr=O(ek2N+1)2N+2R2N+l2N+l→0,

as for any fixed , then we have

 ∥Hk(g)−Hk,N(g)∥L2(B)≤O(ek2N+1)N+1, (3.17)

for large .

Let be the incident wave on and be the corresponding scattered wave. The total wave satisfies

 {Δvk+k2(1+q)vk=0 in Rl,limr→∞r(l−1)/2(∂vs∂ν−ikvs)=0, r=|x|, (3.18)

where in . Clearly we can see that is the far-field pattern of the scattered field corresponding to the incident wave .

Set

 us={u−v in D,0 in Rl∖¯D. (3.19)

Then we have

 {Δus+k2(1+q)us=k2qv, in Rl,limr→∞r(l−1)/2(∂us∂ν−ikus)=0, r=|x|. (3.20)

On the other hand,

 {Δvs+k2(1+q)vs=k2qHk,N(g) in Rl,limr→∞r(l−1)/2(∂vs∂ν−ikvs)=0, r=|x|. (3.21)

Subtracting (3.21) from (3.20), we have

 ⎧⎨⎩Δ(us−vs)+k2(1+q)(us−vs)=k2q(v−Hk,N(g)) in Rl,limr→∞r(l−1)/2(∂(us−vs)∂ν−ik(us−vs))=0, r=|x|. (3.22)

By the well-posedness of the scattering problem (3.22), see for example , we have the following estimate

 ∥us−vs∥H1(Rl)≤C∥v−Hk,N(g)∥H1(Rl),

where is positive constant depending on . It follows from (3.17) that

 ∥us−vs∥H1(Rl)≤Cϵ+O(ek2N+1)N+1.

Therefore, by (3.19) and the well-posedness of the scattering problem, one readily has (3.13). ∎

If does not belong to the interior transmission eigenvalue class , we have the following theorem.

###### Theorem 3.4.

If , i.e., , for some , then for , , there exists large enough and some constant such that

 ∥Fk,N(g)∥L2(Sl−1)≥C,

where the constant depends on .

###### Proof.

By the expansion (3.15), we have

 ∥Hk,N(g)∥2L2(B)=∫R0rl−1∫Sl−1|Hk,N(g)(rd)|2ds(d)dr=N∑n=0n∑m=−n16π2|amn|2∫R0rl−1|jn(kr)|2dr.

By Proposition 3.1 in , and the fact that we see that

 jn(kr)≥(1−ϵ)en+12kn√2(2n+1)n+1rn

for any , and , for some constant depending on . Then

 ∫R0rl−1|jn(kr)|2dr≥(1−ϵ)2e2n+1k2n2(2n+1)2n+2R2n+l2n+l:=Cn=O(ek2n+1)2n+2,

and

 ∫R0rl−1|jn(kr)|2dr→0,

as for any fixed . Thus

 ∥Hk,N(g)∥2L2(B)≥16π2∥g∥2L2(Sl−1)min0≤n≤N∫R0rl−1|jn(kr)|2dr=16π2min0≤n≤NCn. (3.23)

Note that we consider for only finite order , then is a positive number.

We prove that has a lower bound using the contradiction argument. Suppose that . Since is the far-field pattern corresponding to the incident wave , then it follows that is either zero, or otherwise is a Herglotz wave function approximation to the transmission eigenfunction corresponding to transmission eigenvalue . This contradicts to (3.23) and . The proof is done. ∎

## 4. Recovery Scheme

Based on our study in the previous section, we shall present a recovery scheme of locating scatters which possesses cusp singularities and reconstructing the shape of penetrable scatter when it is of polyhedral shape. There is no restriction on the size of the scatter in both situations. Our scheme for the reconstruction is based on the intrinsic geometric properties of the transmission eigenfunctions. If the unknown scattering obstacle has a transmission eigenfunction which is approximated by an entire Herglotz wave function

 vg(x)=∫Sl−1eikx⋅dg(d)ds(d)

with proper chosen kernel function , where is the wave number. Then by finding the vanishing and localizing points of the Herglotz wave function we can locate all the cusp singularities of the scatter.

We shall first concentrate on the problem of retrieving the transmission eigenvalues from far field data over a range of test frequency region. Theoretically, we can retrieve transmission eigenvalues from far field data under no a prior information on . In fact, without any a prior assumption on , the range of test frequency could be wide, therefore the computation is a relatively huge work. We recall the fact that the first transmission eigenvalue, the first Dirichlet eigenvalue and the refractive index of scatter have the estimation in Theorem 2.2 in preliminary. Therefore, if we know a priori assumption on the size of domain and on the upper and lower bounds of refractive index , using the estimation in Theorem 2.2, we can significantly narrow the searching region of frequency, further reduce the computing cost.

Given proper searching frequency region , the forward problem (1.1) is first solved under the incident wave , where . We collect the far field pattern on the circular boundary which enclose the desired polyhedron for incident directions and observation directions, where are positive integers. All the incident directions for each are uniformly distributed on the unit circle or unit sphere, and all the observation points are uniformly distributed on the closed circle (two-dimensional) or surface (three-dimensional) enclosing the scatterer.

From the discussion in the previous section, if is a transmission eigenvalue, then there should exist integer sufficiently large and satisfying , such that

 Fk0,N(g)(^x)≈0, (4.1)

for all incident directions and observation directions . By reciprocity relation, (4.1) is equivalent to

 ∫Sl−1u∞(^x,k0,d)gN(^x)ds(^x)≈0,

for all incident directions . Then, the work is reduced to be a nonlinear minimization problem. We define the cost functional of the minimization problem as

 F(k,g)=min∥g∥L2(Sl−1)=1∑d∣∣Fk,N(g)∣∣. (4.2)

for sufficiently large. The main issue here is to find proper and such that the cost functional is minimized to sufficiently small number. From Theorem 3.3 and 3.4, we see that if the above can be found, then it must be a transmission eigenvalue. Equipped with the finding transmission eigenvalue and Herglotz kernel , the Herglotz wave function can be calculated

 vg(x)=∫Sl−1eik0x⋅dgN(d)ds(d),  x∈Rl, (4.3)

where we know it is the extension of corresponding transmission eigenfunction. The locating of cusp singularities then can be achieved by finding the vanishing and localizing points of Herglotz wave function (4.3).

Summarizing the above discussion, we present the scheme of locating cusp singularities scatter as follows:

Recovery Scheme.

Step 1.  Fix the searching frequency region via estimation in Theorem 2.2, discretize by step size as , and set j=0.

Step 2.  Set . Collect the far field pattern on the circular boundary of radius which enclose the desired polyhedron for uniformly distributed observation directions and uniformly distributed incident directions , where the incident directions are uniformly distributed as , the observation directions are uniformly distributed on the unit circle or unit sphere.

Step 3.  Solve the minimization problem (4.2) for sufficiently large. If the cost functional can be minimized to zero, then go to the next step; Otherwise, go to step 2.

Step 4.  Define

 vg(k)=∫Sl−1eikx⋅dg(d)ds(d) (4.4)

to be the Herglotz function (which is the anlytic extension of transmission eigenfunction). Find the vanishing and localizing points of . Those points are the cusp singularities of the desired scatter.

If further a priori information is available on the support of the medium, say, it is a convex polyhedron, then we can recover its shape by simply joining the cusp singularities by line. The recovery scheme works in a very general and practical setting. It is no need to give any a priori knowledge of the scatter. There is also no restriction to the size of the scatter.

## 5. Numerical experiments and discussions

In this section, we present some numerical tests to verify the applicability and effectiveness of the proposed recovery scheme.

we first concentrate on the problem of retrieving the transmission eigenvalues and eigenfunctions from far field data over a range of frequencies. We collect far field data by solving the forward equation (1.1) under incident wave , for each with step and each observation direction . The forward problem is solved by using the quadratic finite element discretization on a truncated circular (two-dimensional) or spherical (three-dimensional) domain enclosed by a perfectly matched layer (PML). The forward equation is solved on a sequence of successively refined meshes till the relative error of two successive finite element solutions between the two adjacent meshes is below . The scattered data are transformed into the far-field data by employing the Kirchhoff integral formula on a closed circle (two-dimensional) or surface (three-dimensional) enclosing the scatterer. We collect the far field pattern on the circular boundary which enclose the desired polyhedral scatter for incident directions and observation directions (precisely, directions for each uniformly distributed on the unit circle in 2D or unit sphere in 3D, observation points uniformly distributed on the inner part of PML) and a range of wave numbers . Solve the minimization problem (4.2) and find the desired and . The minimizer of the minimization problem (4.2) is obtained by employing a derivative-free trust region method via a local quadratic surrogate model-based search algorithm. Finally, the vanishing and localizing points of Herglotz wave equation (4.4) give us the location of cusp singularities.

The following numerical experiments consist two parts, one is to locate penetrable scatter of general shape with cusp points, the other is to reconstruct the support of polyhedral type shape.

### 5.1. Locating of general shaped scatters with cusps

#### 5.1.1. Example: rain drop shape

The true scatter is of rain drop shape with boundary illustrated in Figure 1,

and described by rotation of the parametric representation

 x(t)=(15sint2,−110sint), 0≤t≤2π.

The location of the corner is . The refractive index of is constant which is not known a prior. From the estimation (2.1), we define the searching frequency region to be and the searching step is . We find that when , the minimization problem (4.2) achieve its minimum. We retrieve the corresponding density function . The corresponding Herglotz function is shown in Figure 2 (a) and it has zero values at shown in Figure 2 (b).