Consider the transmission eigenvalue problem for and :
where Ω is a ball in , . If σ and n are both radially symmetric, namely they are functions of the radial parameter r only, we show that there exists a sequence of transmission eigenfunctions associated with as such that the -energies of ’s are concentrated around . If σ and n are both constant, we show the existence of transmission eigenfunctions such that both and are localized around . Our results extend the recent studies in (SIAM J. Imaging Sci.14 (2021), 946–975; Chow et al.). Through numerics, we also discuss the effects of the medium parameters, namely σ and n, on the geometric patterns of the transmission eigenfunctions.
In this paper, we study the geometric patterns of transmission eigenfunctions. To begin with, we briefly discuss the physical origin of the transmission eigenvalue problem.
Let us consider the time-harmonic wave scattering caused by the interaction of an incident wave field and an inhomogeneous medium. Let denote the incident field, which is an entire solution to in , . Here, signifies the (normalized) angular frequency of the wave propagation. Let signify the inhomogeneous medium. Here, Ω denotes the support of the inhomogeneity of the medium, which is a bounded Lipschitz domain such that . σ and n signify the medium parameters. It is assumed that σ and n are functions and both are bounded below by a positive constant. We also let in . Let u and , respectively, denote the total and scattered wave fields. The wave scattering is governed by the following Helmholtz system:
where , for , and with . The limit in (1.1) is called the radiation condition and characterizes the outgoing nature of the scattered wave field . In the physical setup, when , (1.1) describes the transverse electromagnetic scattering, where and specify the electric permittivity and magnetic permeability of the optical medium [28]; whereas when , (1.1) describes the acoustic scattering, where and are respectively the density and modulus of the acoustic medium [30]. We refer to [30,35] for the well-posedness of the scattering problem (1.1) with a unique solution . It holds that:
In (1.2), is known as the far-field pattern, which encodes the scattering information caused by the perturbation of the incident field due to the scatterer .
Associated with the scattering problem describe above, a practical inverse scattering problem of industrial importance is to recover by knowledge of . The inverse problem can be abstractly recast as the following operator equation:
where is defined via the scattering problem (1.1). For (1.3), one peculiar case is that . In such a case, the scatterer produces no scattering information to the outside observation, namely it is invisible/transparent with respect to the wave probing. Noting that readily yields in by Rellich’s Theorem, one can directly derive that and fulfil that:
where is the exterior unit normal vector to . (1.4) is referred to as the transmission eigenvalue problem. It is clear that are trivial solutions to (1.4). If there exists a nontrivial pair of solutions , k is called a the transmission eigenvalue and u, v are the associated transmission eigenfunctions. Clearly, according to our discussion above, the transmission eigenfunctions depict the wave propagation inside the scatterer when invisibility/transparency occurs.
The spectral theory of transmission eigenvalue problem has received considerable interest in the literature. We refer to [29] for survey and review on the spectral properties of transmission eigenvalues. Recently, several intrinsic local and global geometric patterns of the transmission eigenfunctions have been revealed. Roughly and heuristically speaking, the transmission eigenfunctions tend to (globally) localize/concentrate on while (locally) vanish around singular/high-curvature points on . Here, by localization/concentration, we mean that the -energies of the eigenfunctions in Ω are localized/concentrated around ; and by a singular point, we mean the boundary point on at which the normal vector ν is no longer differentiable. A singular point can be regarded as having (extrinsic) curvature being infinity. The local geometric property was first discovered and investigated in [6,8], and was further studied in [3–5,10,16,19,20,22,25] for different geometric and physical setups. We also refer to [2,3,5,7,9–13,19,21–24,26,31–33,36,37] for related studies in characterizing non-scattering waves (locally) around singular/high-curvature points and the related applications in inverse problems. The global geometric property was first discovered and investigated in [14], and was further studied in [15,17,18] for different geometric and physical setups. Those geometric patterns are physically interpretable. In fact, in order to achieve invisibility/transparency, the wave propagates in a “smart” way which slides over the boundary surface of the scattering object while avoids the singular/highly-curved places to avoid being trapped. More intriguingly, the geometric properties have been used to produce super-resolution imaging schemes for inverse acoustic and electromagnetic scattering problems [14,27], artificial mirage [18], and pseudo surface plasmon resonance [18].
In this paper, we further study the boundary concentration of the transmission eigenfunctions and extend the related studies in [14,15] to more general setups. Specifically, in [14,15], the theoretical justifications are mainly concerned with the case that , n is constant and Ω is smooth and convex, though the general case with variable medium parameter n and non-convex and non-smooth Ω is numerically investigated. We shall include both (being possibly variable) σ and n into the current study and rigorously justify the boundary concentration phenomenon for the transmission eigenfunctions. Moreover, we shall present some novel numerical observations which strengthen the medium effect on the geometric patterns of the transmission eigenfunctions. More detailed discussion about the main results shall be given in Section 2, and the corresponding proofs are provided in Sections 3 and 4.
Statement of the main results and discussion
Following [14], we first provide a quantitative description of surface/boundary localization of a function . In what follows, for , we define
where dist signifies the Euclidean distance in , . Clearly, defines an ε-neighbourhood of .
A function is said to be boundary-localized (or, surface-localized) if there exists such that
According to Definition 2.1, the -energy of the function φ is mainly localized in a small neighbourhood of , namely . In what follows, the asymptotic parameters involved in Definition 2.1 shall become more rigorous. In the case that Ω is a ball in , , by scaling and translation if necessary, we can assume without loss of generality that Ω is the unit ball, namely . In such a case, we also set , to signify the ball of radius τ.
Our main results can be stated as follows.
Consider the transmission eigenvalue problem (
3.1
). Let Ω be the unit ball in,. Assume thatn, σ are functions of the radial parameter r only which fulfil Assumption
A
in what follows. Then for any given, there exists a sequence of eigenfunctionsassociated to eigenvaluesassuch that:
According to Theorem 2.1, if one takes τ to be sufficiently close to 1, namely , it is clear that in (2.7) with m sufficiently large are all boundary-localized according to Definition 2.1. The Assumption A in Theorem 2.1 is stated as follows.
Let , , , be positive constants, and the radial functions , satisfy the following properties:
It is remarked that in our subsequent analysis, we actually can combine conditions (A4) and (A5) by requiring a slightly less restrictive condition:
However, the condition (2.5) involves the eigenvalue , and we split it into the two conditions (A4) and (A5) in Assumption A. It is directly verified that when both σ and n are constant, the assumptions in (2.4) yield that
Nevertheless, if both σ and n are constant, we can prove a stronger boundary-localization result.
Consider the same setup as Theorem
2.1
and assume thatn, σ are both positive constants satisfying. Then for any given, there exists a sequence of eigenfunctionsassociated to eigenvaluesassuch that
That is, if σ and n are both constant, there exists a sequence of transmission eigenfunctions where both the u-parts and v-parts are localized around .
In Theorem 2.2, the condition is required. We would like to point that the condition can also be replaced to be . In fact, in the latter case, for to (1.4), we can set
It is directly verified that
Since , one clearly has . Hence, Theorem 2.2 applied to (2.9) readily yields the existence of a sequence of boundary-localized transmission eigenfunctions associated with as . Then by using the relations in (2.8), we have the existence of a sequence of boundary-localized transmission eigenfunctions associated with as . Finally, we would like to point out that if one takes , the result in Theorem 2.2 recovers those in [14,17].
So far, we have mainly considered the radially symmetric cases. In particular, in Theorem 2.1, we can only show the boundary-localization of the -part, though we believe that there exist infinitely many transmission eigenfunctions such that both u- and v-parts are boundary-localized. In [14], extensive numerical examples show that the surface/boundary-localization is a generic phenomenon occurring for transmission eigenfunctions, even associated with variable medium parameters and general domains. In particular, we note that for the case that Ω is divided into two connected subdomains, with and , if the medium parameters in and in are both constant, then there exist boundary-localized transmission eigenmodes generically even when (which corresponds to variable medium parameters). It is remarked that in the numerical examples in [14], it is always assumed that . Nevertheless, we would like to confirm that the same numerical conclusion holds when and are other constants; see Fig. 1 for typical illustration and comparison. It can be seen that the boundary-localization phenomenon is still every evident though less sharper than the case with . Next, we present two more numerical examples which were not considered in [14], and show that variable medium parameters can make the geometric patterns of the transmission eigenfunctions more intriguing; see Fig. 2, where and respectively signify the left-subdomain and right-subdomain of the domain Ω. It can be observed that in addition to the boundary localization, it may happen that the transmission eigenfunctions are localized around the material interface or even at two “exceptional” points. The numerical observations in Fig. 2 partly corroborates the necessity of introducing Assumption A in our current study. On the other hand, they are highly interesting spectral phenomena that are worth further investigation in our forthcoming work.
First row: (left) and (right) with and . Second row: (left) and (right) with and .
First row: (left) and (right) with and . Second row: (left) and (right) with and .
In this section, we present the proof of Theorem 2.1. That is, the transmission eigenvalue problem is given by:
Throughout the rest of this section, we assume that Ω is the unit ball in , . We shall divide our analysis into two parts, respectively, for the two and three dimensions.
Two-dimensional case
In two dimensions, we let , , denote the polar coordinate. In the sequel, , , signifies the mth order Bessel function [34].
First, we know that the solutions to (3.1) have the following Fourier expansions [34]:
where , . Set
then we have
where the differentiations are with respect to the variable r. Furthermore, we assume that
and then
which implies
For the subsequent use, we let and denote the sth positive root of and , respectively, that are arranged according to the magnitudes [34].
Under Assumption
A
, for anythe functionin (
3.5
) possesses at least one zero point infor any.
The condition (A1) guarantees the non-degeneracy of u in (3.1). Let
and it follows from (A2)-(A5) that
Consider the following differential equations
Then for each fixed , the solutions of the are denoted by , whose roots are given by . It follows from the Singular Strum theorem in [1] and (3.7) that the solution of has at least one zero in for any , denoted by ; see the dashed line in Fig. 3. □
Schematic demonstration of the curve and the point .
Since the solutions of and are continuous with respect to the parameter k, there exists a number such that
which is the intersection of the dashed line and in Fig. 3.
Let in (3.3). Using the transmission condition, it holds that
Set . Using the recursive formula of Bessel functions [34], we have
Next, we find the roots of on the interval .
Under Assumption
A
, for any given, there exists, depending onnand σ, such that when, the functionin (
3.11
) possesses at least one root in.
In the following, we let denote the sth negative zero of the Airy function [38]:
where satisfies the following estimate:
By (1.2) in [38], it holds that
where we recall that denotes the sth positive root of . For each fixed , when is large enough, we have
Consider the interval , there exist at least two consecutive zeros and of , as well as . This result together with the monotonicity of in the interval readily yields that
By virtue of (3.13) and (3.11), it can be directly verified that
Since , are two consecutive roots, and have opposite signs. Hence, we have
By applying Rolle’s theorem to (3.15), we immediately see that there exists at least one zero point of in .
The proof is complete. □
Lemma 3.2 shows the existence of transmission eigenvalues to (3.1) associated with spherically stratified media. In what follows, for a fixed , we let the transmission eigenvalue be denoted by
where is sufficiently large.
Consider the same setup as Theorem
2.1
inand assume that Assumption
A
holds. Letbe the pair of eigenfunctions associated within (
3.16
). We have
Let in (3.2). By direct calculations, we have
which in particular gives that
Combining
and the estimate (3.12), we know that for any , there exists such that when , the following estimate holds
Similar to the arguments in the proof of Theorem 2.6 in [17], for any , there exists such that
which readily implies (3.17).
The proof is complete. □
Three-dimensional result
The proof of Theorem 2.1 in three dimensions follows a similar argument to that of the two-dimensional case in Theorem 3.1. In what follows, we only sketch the necessary modifications in what follows.
In three dimensions, we let denote the polar coordinate, where , and . Let be the spherical harmonic function of order m and degree l, and be the spherical Bessel function [34]. The solutions to (3.1) in have the following Fourier expansions:
where . Set
Then we have
where the derivatives are associated to the variable r. Furthermore, we assume that
Then
which implies
Under Assumption
A
, for anythe functionin (
3.21
) possesses at least one zero point infor any.
The condition (A1) guarantees the non-degeneracy of u in equation (3.1). Let
and it follows from (A2)–(A5) that
Consider the following differential equations:
Then for each fixed , the solutions of the are denoted by , whose roots are given by . It follows from the Singular Strum theorem in [1] and (3.7) that the solutions of has at least one zero in for any , denoted by . □
Since the solutions of and are continuous with respect to the parameter k, there exists a number , such that
Let in (3.19). For the transmission condition, it holds that
Set . Using the recursive formula of the Bessel functions [34], we have
Next, we find the roots of within the interval .
Under Assumption
A
, for any given, there exists, depending onnand σ, such that when, the functionin (
3.30
) possesses at least one root in.
For any given , we have for large enough that
Consider the interval . There exist at least two consecutive zeros and of , as well as . Using such a fact, together with the monotonicity of within the interval , one can show that
By virtue of (3.32) and (3.30), one can show that
Since , are two consecutive roots, it is clear that and have opposite signs. Hence,
By applying Rolle’s theorem to (3.34), one sees that there exists at least one root of in .
The proof is complete. □
In the following, for a given , we let , , denote the transmission eigenvalues determined in Lemma 3.4, where is chosen to be sufficiently large.
Consider the same setup as Theorem
2.1
inand assume thatandis fixed. Letbe the pair of eigenfunctions associated withdefined above. Then it holds that
Let in (3.2). By direct calculations, one has that
which in particular gives that
Combining
and the estimate (3.12), we know that for any , there exists such that when , the following estimate holds
Similar to the arguments in the proof of Theorem 2.6 in [17], for any , there exists such that
which readily gives (3.35).
This section is devoted to the proof of Theorem 2.2. In what follows, we assume that σ and n are both constant. We shall construct a sequence of transmission eigenvalues and prove that the corresponding eigenfunctions are both boundary-localized.
First, we know that the solutions to (3.1) have the following Fourier representations:
Set
Let in (4.2), for the transmission condition, then set
Suppose that. For any given, there exists, depending onnand σ, such that when,in (
3.30
) possesses at least one root inandin (
3.30
) possesses at least one root in.
First, we have by direct calculations that
For a fixed , we have for large enough that
Note that the zeros of are interlaced with those of for (cf. [34]). One clearly has
which readily yields the desired results in the statement of the lemma by Rolle’s theorem. □
Assume thatandis fixed. Letbe the pair of eigenfunctions associated within (
4.2
). We have
For , it is similar to the nonconstant case in the previous section. Hence, we only prove that is boundary-localized in the two dimensions and the three-dimensional case can be proved by following similar arguments.
Since , one has
Next, we show that for m sufficiently large, it holds that . In fact, one can first deduce that
which gives that
Hence, for , there exists a such that when , we have
That is, . Next similar to the arguments in the proof of Theorem 2.6 in [17], we have
which readily proves the first limit in (4.7).
The proof is complete. □
Footnotes
Acknowledgements
The research of H. Liu was supported by the Hong Kong RGC General Research Funds (projects 11300821, 12301420 and 12302919), NSFC/RGC Joint Research Scheme (project N_CityU101/21) and ANR/RGC Joint Research Scheme, A-HKBU203/19. The research of J. Zhang was supported by the Natural Science Foundation of Jiangsu Province (grant no. BK20210540), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (grant no. 21KJB110015). The research of Y. Jiang and K. Zhang was supported in part by China Natural National Science Foundation (grant no. 11871245 and 11971198), the National Key R&D Program of China (grant no. 2020YFA0713601), and by the Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, Jilin University, China.
References
1.
D.Aharonov and U.Elias, Singular Sturm comparison theorems, J. Math. Anal. Appl.371 (2010), 759–763. doi:10.1016/j.jmaa.2010.05.071.
2.
H.Ammari, Y.T.Chow and H.Liu, Localized sensitivity analysis at high-curvature boundary points of reconstructing inclusions in transmission problems, SIAM J. Math. Anal.54(2) (2022), 1543–1592. doi:10.1137/20M1323576.
3.
E.Blåsten, Nonradiating sources and transmission eigenfunctions vanish at corners and edges, SIAM J. Math. Anal.50(6) (2018), 6255–6270. doi:10.1137/18M1182048.
4.
E.Blåsten, X.Li, H.Liu and Y.Wang, On vanishing and localizing of transmission eigenfunctions near singular points: A numerical study, Inverse Problems33 (2017), 105001. doi:10.1088/1361-6420/aa8826.
5.
E.Blåsten and Y.-H.Lin, Radiating and non-radiating sources in elasticity, Inverse Problems35(1) (2019), 015005. doi:10.1088/1361-6420/aae99e.
6.
E.Blåsten and H.Liu, On vanishing near corners of transmission eigenfunctions, J. Funct. Anal.273 (2017), 3616–3632, Addendum:arXiv:1710.08089. doi:10.1016/j.jfa.2017.08.023.
7.
E.Blåsten and H.Liu, Recovering piecewise-constant refractive indices by a single far-field pattern, Inverse Problems36(8) (2020), 085005. doi:10.1088/1361-6420/ab958f.
8.
E.Blåsten and H.Liu, Scattering by curvatures, radiationless sources, transmission eigenfunctions and inverse scattering problems, SIAM J. Math. Anal.53(4) (2021), 3801–3837. doi:10.1137/20M1384002.
9.
E.Blåsten and H.Liu, On corners scattering stably and stable shape determination by a single far-field pattern, Indiana Univ. Math. J.70 (2021), 907–947. doi:10.1512/iumj.2021.70.8411.
10.
E.Blåsten, H.Liu and J.Xiao, On an electromagnetic problem in a corner and its applications, Anal. PDE14(7) (2021), 2207–2224. doi:10.2140/apde.2021.14.2207.
X.Cao, H.Diao and H.Liu, Determining a piecewise conductive medium body by a single far-field measurement, CSIAM Trans. Appl. Math.1 (2020), 740–765. doi:10.4208/csiam-am.2020-0020.
13.
X.Cao, H.Diao, H.Liu and J.Zou, On nodal and generalized singular structures of Laplacian eigenfunctions and applications to inverse scattering problems, J. Math. Pures Appl. (9)143 (2020), 116–161. doi:10.1016/j.matpur.2020.09.011.
14.
Y.-T.Chow, Y.Deng, Y.He, H.Liu and X.Wang, Surface-localized transmission eigenstates, super-resolution imaging and pseudo surface plasmon modes, SIAM J. Imaging Sci.14(3) (2021), 946–975. doi:10.1137/20M1388498.
15.
Y.-T.Chow, Y.Deng, H.Liu and M.Sunkula, Surface concentration of transmission eigenfunctions, arXiv:2109.14361.
16.
Y.Deng, C.Duan and H.Liu, On vanishing near corners of conductive transmission eigenfunctions, Res. Math. Sci.9(1) (2022), 2. doi:10.1007/s40687-021-00299-8.
17.
Y.Deng, Y.Jiang, H.Liu and K.Zhang, On new surface-localized transmission eigenmodes, Inverse Probl. Imaging16(3) (2022), 595–611. doi:10.3934/ipi.2021063.
18.
Y.Deng, H.Liu, X.Wang and W.Wu, On geometrical properties of electromagnetic transmission eigenfunctions and artificial mirage, SIAM J. Appl. Math.82(1) (2022), 1–24. doi:10.1137/21M1413547.
19.
H.Diao, X.Cao and H.Liu, On the geometric structures of transmission eigenfunctions with a conductive boundary condition and applications, Comm. Partial Differential Equations46 (2021), 630–679. doi:10.1080/03605302.2020.1857397.
20.
H.Diao, X.Fei and H.Liu, Local geometric properties of conductive transmission eigenfunctions and applications, arXiv:2206.01933.
21.
H.Diao, X.Fei, H.Liu and K.Yang, Visibility, invisibility and unique recovery of inverse electromagnetic problems with conical singularities, arXiv:2204.02835.
22.
H.Diao, H.Liu and B.Sun, On a local geometric property of the generalized elastic transmission eigenfunctions and application, Inverse Problems37 (2021), 105015. doi:10.1088/1361-6420/ac23c2.
23.
H.Diao, H.Liu and L.Wang, On generalized Holmgren’s principle to the Lamé operator with applications to inverse elastic problems, Calc. Var. Partial Differential Equations59(5) (2020), 179. doi:10.1007/s00526-020-01830-5.
24.
H.Diao, H.Liu and L.Wang, Further results on generalized Holmgren’s principle to the Lamé operator and applications, J. Differential Equations309 (2022), 841–882. doi:10.1016/j.jde.2021.11.039.
25.
H.Diao, H.Liu, X.Wang and K.Yang, On vanishing and localizing around corners of electromagnetic transmission resonance, Partial Differ. Equ. Appl.2(6) (2021), 78. doi:10.1007/s42985-021-00131-6.
26.
H.Diao, H.Liu, L.Zhang and J.Zou, Unique continuation from a generalized impedance edge-corner for Maxwell’s system and applications to inverse problems, Inverse Problems37(3) (2021), 035004. doi:10.1088/1361-6420/abdb42.
27.
Y.He, H.Liu and X.Wang, Invisibility enables super-visibility in electromagnetic imaging, arXiv:2112.07896.
28.
H.Li and H.Liu, On anomalous localized resonance and plasmonic cloaking beyond the quasistatic limit, Proceedings of the Royal Society A474 (2018), 20180165. doi:10.1098/rspa.2018.0165.
29.
H.Liu, On local and global structures of transmission eigenfunctions and beyond, J. Inverse Ill-Posed Probl.30(2) (2022), 287–305. doi:10.1515/jiip-2020-0099.
30.
H.Liu, Z.Shang, H.Sun and J.Zou, On singular perturbation of the reduced wave equation and scattering from an embedded obstacle, J. Dynamics and Differential Equations24 (2012), 803–821. doi:10.1007/s10884-012-9270-5.
31.
H.Liu and C.-H.Tsou, Stable determination of polygonal inclusions in Calderón’s problem by a single partial boundary measurement, Inverse Problems36(8) (2020), 085010. doi:10.1088/1361-6420/ab9d6b.
32.
H.Liu and C.-H.Tsou, Stable determination by a single measurement, scattering bound and regularity of transmission eigenfunctions, Calc. Var. Partial Differential Equations61(3) (2022), 91. doi:10.1007/s00526-022-02211-w.
33.
H.Liu and J.Xiao, On electromagnetic scattering from a penetrable corner, SIAM J. Math. Anal.49(6) (2017), 5207–5241. doi:10.1137/16M110753X.
34.
H.Liu and J.Zou, Zeros of the Bessel and spherical Bessel functions and their applications for uniqueness in inverse acoustic obstacle scattering, IMA J. Appl. Math.72(6) (2007), 817–831. doi:10.1093/imamat/hxm013.
35.
W.McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge Univ. Press, 2000.
M.Salo and H.Shahgholian, Free boundary methods and non-scattering phenomena, arXiv:2106.15154.
38.
R.Wong and C.K.Qu, Best possible upper and lower bounds for the zeros of the Bessel function , Trans. Amer. Math. Soc.351 (1999), 2833–2859. doi:10.1090/S0002-9947-99-02165-0.