{'fraction_non_alphanumeric': 0.03810325139731955, 'fraction_numerical': 0.01775948073802512, 'mean_word_length': 4.316134185303515, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Approximation of the zero-index transmission eigenvalues with a conductive boundary and parameter estimation

In this paper, we present a Spectral-Galerkin Method to approximate the zero-index transmission eigenvalues with a conductive boundary condition. This is a new eigenvalue problem derived from the scalar inverse scattering problem for an isotropic media with a conductive boundary condition. In our analysis, we will consider the equivalent fourth-order eigenvalue problem where we establish the convergence when the approximation space is the span of finitely many Dirichlet eigenfunctions for the Laplacian. We establish the convergence rate of the spectral approximation by appealing to Weyl's law. Numerical examples for computing the eigenvalues and eigenfunctions for the unit disk and unit square are presented. Lastly, we provide a method for estimating the refractive index assuming the conductivity parameter is either sufficiently large or small but otherwise unknown.


Introduction
In this paper, we consider the numerical approximation of the zero-index transmission eigenvalues that are associated with the scalar scattering problem with a conductive boundary. In general, the transmission eigenvalues can be seen as the wave numbers where the associated far-field operator fails to be injective. The zero-index transmission eigenvalue problem is derived by mathematically imbedding the scattering object in a background with refractive index equalling zero in the interior of the scatterer. It can be shown that the resulting far-field operator fails to be injective with a dense range at the wave numbers corresponding to these eigenvalues see [7] for the case when the conductivity is zero. The zero-index transmission eigenvalue problem has two main advantages over the classical transmission eigenvalue problem. First, is that they avoid the assumption that the contrast must be either positive or negative definite in the scatterer. Next, is the fact that they are linear eigenvalue problems. This zero-index eigenvalue problem with a conductive boundary condition was introduced in [22] and was motivated by the work in [9,23] for the classical transmission eigenvalue problem with a conductive boundary and [7] for the scattering problem without a conductive boundary. We are also interested in the inverse spectral problem of estimating the refractive index with little a prior knowledge of the boundary conductivity parameter. There have been manuscripts written on the computation and application of transmission eigenvalue problems to parameter identification such as [11,13,21,30] to name a few. For the classical transmission eigenvalue problem we refer to [2,3] for the application of spectral-Galerkin methods to compute the eigenvalues. See for example [17,19] for some of the previous work for computing the classical transmission eigenvalues via the finite element method. Recently, the method of fundamental solutions for computing the classical transmission eigenvalues was studied and implemented in [25]. Due to the monotonicity property of the transmission eigenvalues one can estimate the refractive index from the knowledge of the eigenvalues(see for example [15,16]). The main contributions of this paper is the convergence analysis with error estimates of the spectral-Galerkin method with the Dirichlet eigenfunctions taken as the basis and the estimation of the refractive index from the zero-index transmission eigenvalues.
The zero-index transmission eigenvalue problem can be written as a fourth-order eigenvalue problem that depends on the refractive index and conductivity. We now derive the fourth-order formulation of the eigenvalue problem. To this end, we define the zero-index transmission eigenvalue problem from the scalar isotropic scattering problem as the values k ∈ C \ {0} such that there exists a nontrivial pair (u, u 0 ) ∈ H 1 (D) × H 1 (D) satisfying the system ∆u + k 2 nu = 0 and ∆u 0 = 0 in D (1) u − u 0 = 0 and ∂ ν u = ∂ ν u 0 + ηu 0 on ∂D. ( Here, we assume that D ⊂ R d (for d = 2, 3) is a simply connected open set where the boundary ∂D is either polygonal with no reentrant corners or class C 2 where ν is the outward unit normal vector. The eigenvalue k corresponds to the wave number for the associated scattering problem. Let the refractive index n ∈ L ∞ (D) and conductivity η ∈ L ∞ (∂D) where we assume that they are uniformly positive definite functions such that there exists positive constants n min ≤ n(x) ≤ n max a.e. x ∈ D and η min ≤ η(x) ≤ η max a.e. x ∈ ∂D.
Therefore, we can define the difference of the eigenfunctions w = u − u 0 . It is clear that w satisfies the equation Due to standard elliptic regularity results ( [18] page 334 for a C 2 boundary) we have that w ∈ H 2 (D) ∩ H 1 0 (D). Now, by appealing to the fact that u 0 is harmonic in D and the boundary condition (2) we can conclude that w satisfies the homogeneous boundary value problem In [22] it is shown that k ∈ C \ {0} is a zero-index transmission eigenvalue problem if and only if there is a nontrivial w ∈ H 2 (D) ∩ H 1 0 (D) satisfying (3). By studying the variational formulation of (3) it is shown that there exists infinitely many real zero-index transmission eigenvalues. This eigenvalue problem is derived by mathematically embedding scatterer D in a background where the refractive index is equal to zero in the interior of the object. This is done by studying the difference of the far-field operators for the standard scattering problem and the augmented far-field operator for the scattering problem where the refractive index is equal zero in D. In general, it is known that the transmission eigenvalues can be determined from the scattering data. In [24] it is shown that the classical transmission eigenvalues can be determined from the far-field data. While in [23] it is shown that the classical transmission eigenvalues with a conductive boundary can also be recovered from farfield data. This implies that these eigenvalues can be used as a target signature to determine the material properties.
The rest of the paper is ordered as follows. In the next section we will study the solution operator corresponding to the zero-index transmission eigenvalue problem with a conductive boundary (3). We will then consider the approximation of the eigenvalues via a Dirichlet spectral-Galerkin method where the approximation space is taken to be the span of finitely many Dirichlet eigenfunctions for the Laplacian. This method of representing the solution to a PDE by the eigenfunctions of an auxiliary eigenvalue problem is studied for physical applications in the manuscript [1]. We study the approximation properties of this space as well as prove convergence of the Dirichlet spectral-Galerkin method for computing the zero-index transmission eigenvalues and provide error estimates. We will then provide some numerical examples in two dimensions to show that the proposed spectral method is effective for computing the eigenvalues. Once we have a method to approximate the eigenvalues we will turn our attention to estimating the refractive index for either large or small valued conductivity parameters.

The Zero-Index Transmission Eigenvalues
This section focuses on the variational formulation of the zero-index transmission eigenvalue problem (3). In particular, we study the associated solution operator. The analysis of the solution operator will be used in the convergence analysis of the spectral method. We define the variational space for (3) as From [22] we have that the equivalent variational form for the zero-index transmission eigenvalue problem (3) is given by the values k ∈ C such that there is a nontrivial We will assume that the eigenfunctions are normalized with w L 2 (D) = 1. The bounded sesquilinear forms on are defined by Recall, that we assume that there exists positive constants n min ≤ n(x) ≤ n max a.e. x ∈ D and η min ≤ η(x) ≤ η max a.e. x ∈ ∂D.
We will study the variational formulation for the zero-index transmission eigenvalue problem in this section. Even though this is a linear eigenvalue problem for k 2 notice that the sesquilinear form b(· , ·) is not sign definite due to the opposing signs in the definition. Which does not give a semi-norm on the variational space which is usually the case for standard elliptic eigenvalue problems.
The well-posedness estimate for the Poisson problem with zero trace along with the H 2 elliptic regularity estimate gives have that for any H 2 (D) function with zero trace ∆ · L 2 (D) is equivalent to the · H 2 (D) . Therefore, we let Clearly, X(D) is a Hilbert space with the associated inner-product. This implies that a(· , ·) is a coercive and Hermitian sesquilinear form on X(D). This implies that k = 0 is not a zero-index transmission eigenvalue. Now, by the Lax-Milgram Lemma we can define the solution operator T : From the definition of T in (6) we have the following result.
Theorem 2.1. Let the operator T : X(D) → X(D) be as defined by (6). Then T is an a(· , ·) self-adjoint compact operator and satisfies the estimate Proof. Since a(· , ·) is a coercive and Hermitian sesquilinear form on X(D) it is an equivalent inner-product on X(D). Therefore, we have that for all f, ϕ ∈ X(D) a T f, ϕ = b(f, ϕ) = b(ϕ, f ) = a T ϕ, f = a f, T ϕ since n and η are real-valued. Proving that T is a(· , ·) self-adjoint on X(D). By the compact embedding of H 1/2 (∂D) into L 2 (∂D) and H 2 (D) into H 1 (D) the compactness of T will follow immediately from the estimate. To prove the estimate notice that by (5) and (6) we can conclude that where we have used the bounds on the coefficients. By appealing to the Trace Theorem and the continuous embedding of H 2 (D) into H 1 (D) we further have that proving the claim.
Notice that since T is a self-adjoint operator on the Hilbert space X(D) with the a(· , ·) inner-product the Hilbert-Schmidt Theorem implies that there exists infinitely many eigenvalues counting multiplicity µ ∈ R for the operator T such that Note that since T is not sign definite there can be complex transmission eigenvalues k that are purely imaginary. In [22] it has been shown that there are infinitely many zero-index transmission eigenvalues k ∈ R. Also, note that again by the Hilbert-Schmidt Theorem we have that there are infinity many eigenfunctions w that form an a(· , ·) orthonormal basis of X(D).

The Dirichlet Spectral-Galerkin Approximation
With the results given in the previous section we can prove the convergence and error estimates of the Dirichlet spectral-Galerkin approximation method of the zero-index transmission eigenvalue problem. We will use the approximation space that is the span of finitely many Dirichlet eigenfunctions for the Laplacian in the domain D. To prove the convergence and error estimates we must show the approximation properties of this space and use the variational formulation (4) to show the convergence of the eigenvalues and eigenfunctions. Even though we focus on the approximation space of Dirichlet eigenfunctions similar analysis as in Section 3.2 will work for any conforming approximation space such as the Legendre-Galerkin approximation which is used for the fourth order formulation of the classical transmission eigenvalue problem in [3].

Approximation Space
Here we will define the approximation space of Dirichlet eigenfunctions and study the approximation properties of the space. To begin, we let φ j ∈ H 1 0 (D) be the jth Dirichlet eigenfunction for the Laplacian and the corresponding eigenvalue λ j ∈ R + for the domain D. The Dirichlet eigenpair satisfy By again appealing to elliptic regularity we have that φ j ∈ X(D). From [22] we have the following result.
The eigenvalues are assumed to be arranged in non-decreasing order such that 0 < λ j ≤ λ j+1 for all j ∈ N. It is well known that the eigenfunctions {φ j } ∞ j=1 form an orthonormal basis of L 2 (D) which implies that for all f ∈ L 2 (D) as a convergent series in the L 2 (D) norm. This series representation will be used to show the approximation rates for this set of basis functions. To do so, we will show the convergence of this series in the X(D) norm.
Proof. The result follows by Lemma 3.1 which gives that ψ j = φ j /λ j is an orthonormal basis for X(D) by equation (7) such that as a X(D) convergent series. Then appealing to Green's 2nd Theorem we have that Therefore, we obtain that is a convergent series in the X(D) norm. Then applying the Laplacian term by term to the series representation proves the claim.
By the series representation (8) for any f ∈ L 2 (D) along with (7) we can define the powers of the Laplacian ∆ m for m ∈ N by Note that this is often done to define (fractional) powers of an elliptic operator(see for example [26]). We will denote the domain of the mth power of the Laplacian in the set L 2 (D) as (9) is an L 2 (D) convergent series .
Therefore, we have that D ∆ m is a Hilbert space with the associated norm By Theorem 3.1 we have that X(D) ⊆ D(∆). For our spectral approximation method we will take the conforming computational subspace of X(D) to be given by for some fixed N ∈ N.
A key ingredient to determining the approximation rate for this set of basis functions is Weyl's law for the Dirichlet eigenvalues(see [5]). Weyl's law states that there exists two constants c 1 , c 2 > 0 independent of j such that where again the dimension d = 2, 3. We now consider the L 2 (D) projection onto the approximation space X N (D) which we denote Π N : X(D) → X N (D) and is given by It is clear that we have the point-wise convergence for any f ∈ X(D) by Theorem 3.1. We now give a convergence estimate for any function such that equation (10) is satisfied.
Proof. To prove the claim we estimate where we have used the series representation in Theorem 3.1. Now appealing to Weyl's law and by (10) we have that proving the estimate.

Convergence and Error Estimates
In this section, we will establish the Dirichlet spectral-Galerkin method for the zeroindex transmission eigenvalue problem. In our analysis we will use the approximation space X N (D) defined as the span of the first N Dirichlet eigenfunctions for the Laplacian in D. Similar results can be established by using other conforming approximation subspaces such as a finite element approximation space of piecewise polynomials. Using the convergence analysis for the approximation space in the previous section we are now ready to prove the convergence on the spectral approximation.
To begin, we will first consider the approximation for the operator T defined in (6) by the L 2 (D) projection of T onto the space X N (D). We will show that this approximation converges in the operator norm. This fact will be used to prove convergence and error estimates for the approximation of the eigenvalues. Theorem 3.3. Let the operator T : X(D) → X(D) be as defined by (6) and Π N : Proof. Since the operator T defined in equation (6) is compact by Theorem 2.1 we have that the point-wise convergence of Π N to the Identity operator on X(D) implies that Proving the claim.
We now define the spectral approximation of the zero-index transmission eigenvalue problem as find the values k N ∈ C such that there is a nontrivial w N ∈ X N (D) satisfying where the sesquilinear forms a(· , ·) and b(· , ·) are defined by (5). Here we again assume that the eigenfunctions are normalized such that w N L 2 (D) = 1. Therefore, just as in the continuous case we can define the spectral approximation of the solution operator as operator T N : Since the dimension of the range of T N is finite we have that it is compact. It is also clear that T N restricted to X N (D) is an a(· , ·) self-adjoint operator on X N (D). This gives that T N has N eigenvalues counting multiplicity. Arguing similarly as in Section 2 we have that the eigenpair (k N , w N ) ∈ C × X N (D) satisfying (11) is the eigenpair for the spectral approximation of the solution operator such that In order to attain the convergence as well as an error estimate we will study the convergence of the spectral approximation of the solution operator as N → ∞. Therefore, by appealing to Galerkin orthogonality and Cea's Lemma( [6] page 372) we have that and for all f ∈ X(D). From the above estimate we conclude that where again Π N is the L 2 (D) projection onto the approximation space X N (D). This analysis gives the following result. Proof. This result follows from the fact that and then appealing to the results in [27]. Now that we have established the convergence of the spectral approximation our next step it to determine the convergence rate. To this end, we will argue similarly to Theorem 3.2 along with the using variational formulations (4) and (11). Simple calculations give that for (k, w) and (k N , w N ) being the jth egienpair for (4) and (11) then for any N . The equality (13) will be used to establish the convergence rate for the eigenvalues. So we need to establish that the sequence b w N , w N is bounded below. Therefore, we present the following result.
Theorem 3.5. Let (k N , w N ) be the jth egienpair for (11). Then there is a constant Proof. To prove the claim we proceed by way of contradiction. To this end, assume no such β exists, then we can extract a subsequence still denoted with N such that b w N , w N → 0 as N → ∞.
By the continuity of the sesquilinear form b(· , ·) and Theorem 3.4 we obtain where w is an eigenfunction corresponding to (4). The variational formulation implies that a(w, w) = 0 which contradicts the fact that w L 2 (D) = 1 since a(· , ·) defines an inner-product on X(D). Proving the claim.
From Theorem 3.5 we can conclude that for (k, w) and (k N , w N ) being the jth egienpair for (4) and (11) respectively then there is a C > 0 independent of N where Note that we have used (13) and the boundedness of sesquilinear forms a(· , ·) and b(· , ·). In order to obtain the error estimate for the spectral approximation of the zero-index transmission eigenvalues we must estimate the error in the Galerkin approximation in the approximation space X N (D) on the eigenspace corresponding to k. To this end, we will denote the eigenspace corresponding to the zero-index transmission eigenvalue k by E(k). It is clear that E(k) ⊂ X(D) is finite dimensional and any u ∈ E(k) satisfies T u = k −2 u. With this we can now prove the error estimate.
Theorem 3.6. Let k and k N be the jth eigenvalues for (4) and (11) respectively. If the corresponding eigenspace E(k) ⊂ D ∆ m for m ∈ N then there is a C > 0 independent of N such that Proof. Notice that according to the analysis in [8], we have that Therefore, by the definition of the norm · E(k) →X(D) we can estimate We further have that Where we have used the fact that u ∈ E(k). Now, since we have assumed that E(k) ⊂ D ∆ m we further estimate just as in Theorem 3.2 which we obtain by appealing to Weyl's law. This estimate give the convergence rate proving the claim.
We end this section by noting that the convergence rate for the eigenfunctions are the square root of the convergence rate for the eigenvalues. This is clear from the proof of Theorem 3.6 since estimates for w N − w 2 X(D) are used to derive the estimates for the eigenvalue convergence rate.
Theorem 3.7. Let w and w N be the jth eigenfunctions for (4) and (11) respectively. If the corresponding eigenspace E(k) ⊂ D ∆ m for m ∈ N then there is a C > 0 independent of N such that Proof. The result follows from the analysis given in the proof of Theorem 3.6.

Numerical Examples
In this section, we provided some numerical examples of computing the zero-index transmission eigenvalues via the Dirichlet spectral-Galerkin method studied in the previous sections. For simplicity, will assume that the domain D is the unit ball in R 2 . We refer to [4] for the approximation of the classical transmission eigenvalues via a Spectral Element Method for a spherically stratified media. We will check the accuracy of our Dirichlet spectral-Galerkin method by comparing to the eigenvalues computed by separation of variables for constant refractive index n and conductivity η. The computations in this section are motivated by the work in [20] where the authors studied the convergence of the spectral-Galerkin method for computing the classical transmission eigenvalues where the basis functions in the approximation space are the eigenfunctions for the bilaplacian with zero clamped plate boundary conditions. We will also consider estimating the refractive index n for η either small or large but unknown. To do so, we will use the convergence results given in [22] where it is shown that as the conductivity tends to zero or infinity one obtains an eigenvalue problem that only depends on the refractive index. The estimation of the refractive index from the scattering data has been studied in [12,29] to name a few examples. In these papers the authors show that the far and near field data can be used to recover the classical transmission eigenvalues and use monotonicity results to recover a constant refractive index or estimate a variable refractive index. Here we numerically study this problem for the zero-index transmission eigenvalues.

Computing the Eigenvalues
We are now ready to compute the zero-index transmission eigenvalues using the Dirichlet spectral-Galerkin approximation presented in Section 3. All of our experiments are done with the software MATLAB 2018a on an iMac with a 4.2 GHz Intel Core i7 processor with 8GB of memory. To begin, we will describe an effective method for computing the approximated eigenvalues that satisfy (11). Since the domain D is given by the unit circle we can determine the basis functions from separation of variables. Therefore, the basis functions are taken to be Here square root of the eigenvalues λ p,q corresponds to the qth positive root of the pth first kind Bessel function denoted J p for all p ∈ N ∪ {0} and q ∈ N.
In the following numerical examples we take 24 basis functions where 0 ≤ p ≤ 5 and 1 ≤ q ≤ 4 which will give that the approximation space is defined as For our Dirichlet spectral-Galerkin approximation we will solve (11) for w N in the aforementioned approximation space. This gives that the approximation of the eigenfunctions will have the form Therefore, the spectral approximation of the eigenvalues k N satisfying (11) correspond to the eigenvalues for the matrix equation We have that the 24 × 24 Galerkin mass and stiffness matrices in the Spectral approximation (11) are given by The integrals can be simplified by using (7) to evaluate the sesquilinear forms a(· , ·) and b(· , ·) giving that Notice that we have used that for the Dirichlet eigenfunctions φ i and φ j we have the following integral identities In our calculations, we use the fact that the Dirichlet eigenfunctions are orthogonal in L 2 (D) which gives that the volume integral in b(φ i , φ j ) corresponds to a diagonal matrix.
To compute the Galerkin matrices we implement a 2d Gaussian quadrature method. In the numerical examples the integrals are written in polar coordinates where 12 quadrature points are used to evaluate the radial and angular parts of the integrals. The discretized eigenvalue problem (14) is then solved using the 'polyeig' command in MATLAB since (14) is a quadratic eigenvalue problem for the parameter k N . From [22] we have that for n and η constant the eigenvalues k satisfy the transcendental equation Using the 'fzero' command in MATLAB we can determine the exact zero-index transmission eigenvalues. In Figures 1 and 2 we plot the function d m (k) for the values of    The plot of the relative error for the first eigenvalue k 1 is given where we let n vary in the interval [3,5] for η = 25 or η = 1/10. We use d 0 (k) to compute the exact first zero-index transmission eigenvalues. In Figure 3 we see that the relative error for η = 25 is on the order of 10 −4 where as the relative error for η = 1/10 is on the order of 10 −2 . This seems to imply that the Dirichlet spectral-Galerkin approximation method is better suited for problems with larger conductivity. We now consider computing the eigenvalues for variable coefficients. Here we take a smooth and a piece-wise constant refractive index defined as n 1 = 4 + exp(−r 2 ) and n 2 = 4 * 1 (0.25≤r<1) + 2 * 1 (r<0. 25) for a spherically stratified media. Here 1 I denotes the indicator function on the interval I. In Table 3 and 4 we report the first three real zero-index transmission eigenvalues computed via our approximation method for various conductivities. Here we take three different parameters η. Two of the conductivities are constants taken to be 25 and 1/10 while the third is a variable conductivity parameter η = 1/ 10 + sin 2 (2θ) .
Recall, that the analysis in Section 3.2 also gives the convergence of the eigenfunctions in Theorems 3.4 and 3.7. So, presented with Tables 3 and 4 are the plots for the first two corresponding eigenfunctions for the spherically stratified refractive indices n 1 and n 2 defined above.  Table 3: The first three real zero-index transmission eigenvalues for n 1 . We also plot the eigenfunctions corresponding to the eigenvalues k 1 and k 2 for η = 1/10.
Notice that the computed zero-index transmission eigenvalues are monotonically decreasing with respect to the coefficients which is predicted by the theory in [22]. We can see that for the first three real zero-index transmission eigenvalues that for η = 25 and 1/10 we have k j (n 1 , η) ≤ k j (4, η) ≤ k j (n 2 , η).  Table 4: The first three real zero-index transmission eigenvalues for n 2 . We also plot the eigenfunctions corresponding to the eigenvalues k 1 and k 2 for η = 25.
Similarly comparing the reported eigenvalues we see the monotonicity with respect to the conductivity η for various refractive indices.

Estimating the Refractive Index
In this section, we present numerical examples for estimating the refractive index from the knowledge of the zero-index transmission eigenvalues. It has been shown in [9,23] that the classical transmission eigenvalues with a conductive boundary condition can be recovered from the scattering data via the Linear Sampling Method and the Inside-Out Duality, see [14,24] for details of these methods to recover the transmission eigenvalues. Therefore, we will assume that the zero-index transmission eigenvalues can be recovered from the scattering data and we wish to estimate n.
In order to estimate the refractive index from the zero-index transmission eigenvalues k(n, η) we will restrict ourselves to the case where η is either sufficiently large or small. The case for η known was considered numerically in [22]. The limiting behavior of the zero-index transmission eigenvalues was studied in [22] as η tends to zero or infinity. It has been shown that k(n, η) → τ (n) as η → ∞ where τ is a 'modified' Dirichlet eigenvalue satisfying that there exist a nontrivial v such that or as η → 0 where τ is a 'modified' plate buckling eigenvalue satisfying that there exist a nontrivial v such that This limiting behavior will allow use to estimate the refractive index without the knowledge of η on the boundary. This gives that if it is known a prior that η 1 or η 1 then we can estimate the refractive index by finding the constant n approx such that k 1 (n, η) = τ 1 (n approx ) where τ 1 is the first 'modified' Dirichlet eigenvalue for η 1 or the first 'Modified' plate buckling eigenvalue for η 1. It is known that τ 1 depends monotonically on n by the max-min principle [30]. Since D is assumed to be known we can compute τ 1 for any constant refractive index n via the methods from [10,28]. To numerically approximate τ 1 we use separation of variables since D is the unit circle. Therefore, we have that the 'modified' Dirichlet eigenvalues for a constant n satisfies J |m| τ √ n = 0 for all m ∈ Z and the 'modified' plate buckling eigenvalues for a constant n satisfies To determine the approximate refractive index we first find the polynomial interpolation for τ 1 (n) for constant n ∈ [2, 7] via the 'polyfit' command in MATLAB. Then we solve for the constant n approx such that k 1 (n, η) = τ 1 (n approx ) via the 'fzero' command in MATLAB. Since τ 1 is a deceasing function of n the above equation has a unique solution n approx . The results are reported in Tables 5 and 6 for the spherically stratified refractive indices used in our previous calculations with variable coefficient conductivity parameters.   Table 6: Estimation of the refractive index n for η = 25(2 + sin 4 (θ)).
Simple calculus gives that the average value of n = 4 + exp(−r 2 ) in the unit disk to be 4 + (1 − exp(−1)) ≈ 4.6321205 which is fairly close to the approximation in Table 6. We can also compute the average value for the piece-wise constant refractive index n = 4 * 1 (0.25≤r<1) + 2 * 1 (r<0. 25) in the unit disk which is 3.875. Also notice that due to the monotonicity of τ 1 we have that n min ≤ n approx ≤ n max . In the case of the classical transmission eigenvalues it has been numerically documented that estimating the refractive index by a constant leads to determining it's average value [29]. Table 6 seems to suggest that for η 1 that estimating the refractive index by a constant may also lead to determining the average value.

A Numerical Example for the Unit Square
For completeness we provide numerical examples for the unit square D = (0, 1) 2 . This is given to show that this method also works for polygonal domains with no reentrant corners. Here we wish to show the accuracy of the approximation for this domain. To this end, we will compute the zero-index transmission eigenvalues for a constant η. To establish that the approximation is accurately computing the eigenvalues we will show test convergence as η → ∞ as well as the monotonicity. We will also estimate the refractive index n(x 1 , x 2 ) assuming η 1 just as we did for the unit sphere. Therefore, the zero-index transmission eigenvalues k should converge to the 'modified' Dirichlet eigenvalues for the unit square.
For the approximate we have that the basis functions are taken to be In the numerical examples we take 25 basis functions where 1 ≤ p, q ≤ 5 which gives the spectral approximation space as Span φ j(p,q) (x 1 , x 2 ) 5 p,q=1 ⊂ X(D).
To compute the zero-index transmission eigenvalues we proceed just as in the previous section. The spectral approximation of the eigenvalues k N satisfy the corresponding matrix eigenvalue problem with the appropriate mass and stiffness matrices. Again the matrix eigenvalue problem is solved by using 'polyeig' command in MATLAB. η = 1 η = 10 η = 100 η = 1000 k 1 2.42379135332 2.23942914304 2.22322075768 2.22161920572 Table 7: The first zero-index transmission eigenvalue for n = 4 and η = 10 m with m = 0, 1, 2, 3 of the unit square. Here we see that k 1 is converging to the first 'modified' Dirichlet eigenvalues for the unit square π/ √ 2 ≈ 2.2214414. η = 1 η = 10 η = 100 η = 1000 k 1 2.23128387981 2.00960194537 1.96995925736 1.96292112038 Table 8: The first zero-index transmission eigenvalue for n = Here we see in Table 7 the convergence of the first zero-index transmission eigenvalue k 1 to the first 'modified' Dirichlet eigenvalue. Also, notice that in Tables 7 and 8 the monotonicity of the transmission eigenvalue k 1 (n, η) with respect to n and η is verified by the calculations. Now from the approximated transmission eigenvalue we can again estimate the refractive index. Using the convergence as η → ∞ we proceed just as in the previous section. That is we find n approx such that k 1 (n, η) = τ 1 (n approx ) where the conductivity parameter η 1. To estimate the refractive index we can solve the above equation exactly since the 'modified' Dirichlet eigenvalues are known analytically. In Tables 9 we present the approximation of two refractive indices from the first zero-index transmission eigenvalue for the unit square.  Table 9: Estimation of the refractive index n for η = 10.

Summary and Conclusions
In conclusion, we have provided a numerical method to compute the zero-index transmission eigenvalues via the Dirichlet spectral-Galerkin approximation method. Our approximation space is taken to be the span of the first N Dirichlet eigenfunctions. Even though our numerical examples are only presented in for the unit disk in R 2 the analysis is also valid for any domain where the boundary ∂D is either polygonal with no reentrant corners or class C 2 in R d for d = 2, 3. In order to apply this method one needs the Dirichlet eigenpairs a prior for the domain of interest. In theory this can be done by pre-calculating a fixed number of Dirichlet eigenpairs for the domain D via BEM [28] or FEM [30]. We have also given numerical examples to validate the theoretical results as well as investigated estimating the refractive index from the first zero-index transmission eigenvalue. It seems that for η sufficiently large one can estimate the average value of n which can be used for nondestructive testing. Possible future work would consist of applying this method to compute 'classical' transmission eigenvalues with a conductive boundary and considering the inverse problem of recovering a variable coefficient refractive index from the eigenvalues.