Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
\addbibresource

references.bib

Gauge freedoms in the anisotropic elastic Dirichlet-to-Neumann map

Joonas Ilmavirta1, Hjørdis Schlüter1,∗ 1Department of Mathematics and Statistics, University of Jyväskylä, Finland Corresponding author joonas.ilmavirta@jyu.fi, hjordis.a.schluter@jyu.fi
(Date: August 22, 2024)
Abstract.

We address the inverse problem of recovering the stiffness tensor and density of mass from the Dirichlet-to-Neumann map. We study the invariance of the Euclidean and Riemannian elastic wave equation under coordinate transformations. Furthermore, we present gauge freedoms between the parameters that leave the elastic wave equations invariant. We use these results to present gauge freedoms in the Dirichlet-to-Neumann map associated to the Riemannian elastic wave equation.

Key words and phrases:
Dirichlet-to-Neumann map, inverse problems, elastic wave equation, Riemannian geometry, anisotropic stiffness tensor.
1991 Mathematics Subject Classification:
35R30

1. Introduction

We study the elastic wave equation (EWE) in the n𝑛nitalic_n-dimensional Euclidean space and on an n𝑛nitalic_n-dimensional Riemannian manifold. We address the inverse problem of recovering the anisotropic stiffness tensor and the density of mass from the Dirichlet-to-Neumann (DN) map. In this context it is essential for the reconstruction procedure whether the DN map determines the stiffness tensor and density uniquely. For that purpose we study invariance of the elastic wave equation: How the EWE behaves under coordinate transformations and what gauge freedoms there are between the parameters. Based on these results we present gauge freedoms for the DN map in the Riemannian setting and conjecture that these are the correct and full gauge groups. If this holds true then the Euclidean DN map determines the stiffness tensor and density uniquely in two (and possibly higher) dimensions. This is of particular interest as the three-dimensional Euclidean EWE is the natural setting in seismology, where the EWE models seismic waves.

Let Mn𝑀superscript𝑛M\subset\mathbb{R}^{n}italic_M ⊂ blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT be the closure of a smooth domain, and let g𝑔gitalic_g be a Riemannian metric on M𝑀Mitalic_M. On M𝑀Mitalic_M we use the Euclidean coordinates x𝑥xitalic_x. The material parameters are the stiffness tensor c𝑐citalic_c (contravariant of rank 4444) and the density ρ𝜌\rhoitalic_ρ. Following the symmetries of the stress and strain tensors and their relationship, the stiffness tensor has the minor and major symmetries

cijk=cjik=ckij,superscript𝑐𝑖𝑗𝑘superscript𝑐𝑗𝑖𝑘superscript𝑐𝑘𝑖𝑗c^{ijk\ell}=c^{jik\ell}=c^{k\ell ij},italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT = italic_c start_POSTSUPERSCRIPT italic_j italic_i italic_k roman_ℓ end_POSTSUPERSCRIPT = italic_c start_POSTSUPERSCRIPT italic_k roman_ℓ italic_i italic_j end_POSTSUPERSCRIPT , (1)

and is positive in the sense that111Here and throughout the paper we use the Einstein summation convention, also in the Euclidean setting where all indices are kept down. Whenever there is a non-repeated index, the corresponding equation holds for all values 1,,n1𝑛1,\dots,n1 , … , italic_n of it.

cijkAijAkδgjkgiAijAksuperscript𝑐𝑖𝑗𝑘subscript𝐴𝑖𝑗subscript𝐴𝑘𝛿superscript𝑔𝑗𝑘superscript𝑔𝑖subscript𝐴𝑖𝑗subscript𝐴𝑘c^{ijk\ell}A_{ij}A_{k\ell}\geq\delta g^{jk}g^{i\ell}A_{ij}A_{k\ell}italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_k roman_ℓ end_POSTSUBSCRIPT ≥ italic_δ italic_g start_POSTSUPERSCRIPT italic_j italic_k end_POSTSUPERSCRIPT italic_g start_POSTSUPERSCRIPT italic_i roman_ℓ end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_k roman_ℓ end_POSTSUBSCRIPT (2)

for some δ>0𝛿0\delta>0italic_δ > 0 and for all symmetric matrices A𝐴Aitalic_A. Additionally, the density of mass ρ𝜌\rhoitalic_ρ is positive:

ρK,𝜌𝐾\rho\geq K,italic_ρ ≥ italic_K ,

for some constant K>0𝐾0K>0italic_K > 0.

See [dehoop2023reconstruction] for a discussion on these structures. The density is assumed positive in the usual sense. The metric tensor is not a property of the material but of the space itself, and it is therefore forced to be Euclidean in practice — this will play an important role.

These material parameters give rise to the elastic Laplacian ΔcgsuperscriptsubscriptΔ𝑐𝑔\Delta_{c}^{g}\,roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT defined by

(Δcgu(t,x))i=1det(g)xj(det(g)cijkgmxkum(t,x)),superscriptsuperscriptsubscriptΔ𝑐𝑔𝑢𝑡𝑥𝑖1𝑔subscriptsuperscript𝑥𝑗𝑔superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑘superscript𝑢𝑚𝑡𝑥\displaystyle(\Delta_{c}^{g}\,u(t,x))^{i}=\frac{1}{\sqrt{\det(g)}}\,\partial_{% x^{j}}\left(\sqrt{\det(g)}c^{ijk\ell}g_{\ell m}\partial_{x^{k}}u^{m}(t,x)% \right),( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT italic_u ( italic_t , italic_x ) ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG square-root start_ARG roman_det ( italic_g ) end_ARG end_ARG ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( square-root start_ARG roman_det ( italic_g ) end_ARG italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_t , italic_x ) ) , (3)

which maps vector fields to vector fields. The elastic wave operator is

P(ρ,c,g)=t2ρ1Δcg.subscript𝑃𝜌𝑐𝑔superscriptsubscript𝑡2superscript𝜌1superscriptsubscriptΔ𝑐𝑔P_{(\rho,c,g)}=\partial_{t}^{2}-\rho^{-1}\Delta_{c}^{g}\,.italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT = ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT . (4)

The displacement vector u(t,x)n𝑢𝑡𝑥superscript𝑛u(t,x)\in\mathbb{R}^{n}italic_u ( italic_t , italic_x ) ∈ blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT satisfies the Riemannian version P(ρ,c,g)u=0subscript𝑃𝜌𝑐𝑔𝑢0P_{(\rho,c,g)}u=0italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT italic_u = 0 of the elastic wave equation.

We are therefore naturally led to the boundary value problem

{P(ρ,c,g)u(t,x)=0in (0,T)×M,u(t,x)=f(t,x)on (0,T)×M,u(0,x)=tu(0,x)=0in Mcasessubscript𝑃𝜌𝑐𝑔𝑢𝑡𝑥0in 0𝑇𝑀𝑢𝑡𝑥𝑓𝑡𝑥on 0𝑇𝑀𝑢0𝑥subscript𝑡𝑢0𝑥0in 𝑀\begin{cases}P_{(\rho,c,g)}u(t,x)=0&\text{in }(0,T)\times M,\\ u(t,x)=f(t,x)&\text{on }(0,T)\times\partial M,\\ u(0,x)=\partial_{t}u(0,x)=0&\text{in }M\end{cases}{ start_ROW start_CELL italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT italic_u ( italic_t , italic_x ) = 0 end_CELL start_CELL in ( 0 , italic_T ) × italic_M , end_CELL end_ROW start_ROW start_CELL italic_u ( italic_t , italic_x ) = italic_f ( italic_t , italic_x ) end_CELL start_CELL on ( 0 , italic_T ) × ∂ italic_M , end_CELL end_ROW start_ROW start_CELL italic_u ( 0 , italic_x ) = ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_u ( 0 , italic_x ) = 0 end_CELL start_CELL in italic_M end_CELL end_ROW (5)

in the spacetime (0,T)×M0𝑇𝑀(0,T)\times M( 0 , italic_T ) × italic_M and where f𝑓fitalic_f satisfies: f(x,0)=tf(x,0)=0𝑓𝑥0subscript𝑡𝑓𝑥00f(x,0)=\partial_{t}f(x,0)=0italic_f ( italic_x , 0 ) = ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_f ( italic_x , 0 ) = 0 for xM𝑥𝑀x\in\partial Mitalic_x ∈ ∂ italic_M.

For the special case g=gE𝑔subscript𝑔𝐸g=g_{E}italic_g = italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT the Riemannian EWE reduces to the familiar Euclidean EWE

{(P(ρ,c,gE)u(t,x))=0in (0,T)×Mu(t,x)=f(t,x)on (0,T)×M,u(0,x)=tu(0,x)=0in M.casessubscript𝑃𝜌𝑐subscript𝑔𝐸𝑢𝑡𝑥0in 0𝑇𝑀𝑢𝑡𝑥𝑓𝑡𝑥on 0𝑇𝑀𝑢0𝑥subscript𝑡𝑢0𝑥0in 𝑀\begin{cases}(P_{(\rho,c,g_{E})}u(t,x))=0&\text{in }(0,T)\times M\\ u(t,x)=f(t,x)&\text{on }(0,T)\times\partial M,\\ u(0,x)=\partial_{t}u(0,x)=0&\text{in }M.\end{cases}{ start_ROW start_CELL ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT italic_u ( italic_t , italic_x ) ) = 0 end_CELL start_CELL in ( 0 , italic_T ) × italic_M end_CELL end_ROW start_ROW start_CELL italic_u ( italic_t , italic_x ) = italic_f ( italic_t , italic_x ) end_CELL start_CELL on ( 0 , italic_T ) × ∂ italic_M , end_CELL end_ROW start_ROW start_CELL italic_u ( 0 , italic_x ) = ∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_u ( 0 , italic_x ) = 0 end_CELL start_CELL in italic_M . end_CELL end_ROW (6)

The Euclidean elastic Laplacian acts as

(ΔcgEu(t,x))i=xj(cijkxku(t,x)).subscriptsuperscriptsubscriptΔ𝑐subscript𝑔𝐸𝑢𝑡𝑥𝑖subscriptsubscript𝑥𝑗subscript𝑐𝑖𝑗𝑘subscriptsubscript𝑥𝑘subscript𝑢𝑡𝑥(\Delta_{c}^{g_{E}}\,u(t,x))_{i}=\partial_{x_{j}}\left(c_{ijk\ell}\,\partial_{% x_{k}}u_{\ell}(t,x)\right).( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_u ( italic_t , italic_x ) ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ∂ start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_c start_POSTSUBSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( italic_t , italic_x ) ) .

For the Euclidean elastic wave operator P(ρ,c,gE)subscript𝑃𝜌𝑐subscript𝑔𝐸P_{(\rho,c,g_{E})}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT we will use the short hand notation P(ρ,c,gE)=P(ρ,c)subscript𝑃𝜌𝑐subscript𝑔𝐸subscript𝑃𝜌𝑐P_{(\rho,c,g_{E})}=P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT = italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT in the following. Assuming that cijkL(M),ρL(M)formulae-sequencesubscript𝑐𝑖𝑗𝑘superscript𝐿𝑀𝜌superscript𝐿𝑀c_{ijk\ell}\in L^{\infty}(M),\rho\in L^{\infty}(M)italic_c start_POSTSUBSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUBSCRIPT ∈ italic_L start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( italic_M ) , italic_ρ ∈ italic_L start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( italic_M ) and with fL2(0,T;H1(M))H1(0,T;L2(M))𝑓superscript𝐿20𝑇superscript𝐻1𝑀superscript𝐻10𝑇superscript𝐿2𝑀f\in L^{2}(0,T;H^{1}(\partial M))\cap H^{1}(0,T;L^{2}(\partial M))italic_f ∈ italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( ∂ italic_M ) ) ∩ italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( ∂ italic_M ) ) the initial value boundary value problem (6) has the unique solution uL2(0,T;H1(M))H1(0,T;L2(M))𝑢superscript𝐿20𝑇superscript𝐻1𝑀superscript𝐻10𝑇superscript𝐿2𝑀u\in L^{2}(0,T;H^{1}(M))\cap H^{1}(0,T;L^{2}(M))italic_u ∈ italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_M ) ) ∩ italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_M ) ) and additionally, uC(0,T;H1(M))C1(0,T;L2(M))𝑢𝐶0𝑇superscript𝐻1𝑀superscript𝐶10𝑇superscript𝐿2𝑀u\in C(0,T;H^{1}(M))\cap C^{1}(0,T;L^{2}(M))italic_u ∈ italic_C ( 0 , italic_T ; italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_M ) ) ∩ italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_M ) ) [stolk00, Thm. 2.4.5],[BeLa02, Lemma 3.5]. Following the analysis by Stolk in [stolk00, Sec. 2.3] these well-posedness results can be extended to the Riemannian initial boundary value problem (5).

1.1. The Finsler metric arising from elasticity

Inverse problems related to elasticity are mainly concerned with recovering the stiffness tensor c𝑐citalic_c or the reduced stiffness tensor a=ρ1c𝑎superscript𝜌1𝑐a=\rho^{-1}citalic_a = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c everywhere from boundary data. One way to define elastic geometry is in terms of travel time distance between two points. Here one considers an elastic body that is modeled as a manifold and distance is measured by the shortest time it takes for an elastic wave to go from one point to the other. When the elastic material is elliptically anisotropic the elastic geometry is Riemannian, but this puts very stringent restrictions on the stiffness tensor.

Recent research is devoted to the fully anisotropic setting for the stiffness tensor, where only the physically necessary assumptions as in (1)–(2) are needed. The resulting elastic geometry is not Riemannian but Finslerian [HILSmay21, HILSaug21]. The Finsler metric can be derived from the principal behavior of the elastic wave operator P(ρ,c,g)subscript𝑃𝜌𝑐𝑔P_{(\rho,c,g)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT defined above. The components of its matrix-valued principal symbol σ(t,x,ω,p)𝜎𝑡𝑥𝜔𝑝\sigma(t,x,\omega,p)italic_σ ( italic_t , italic_x , italic_ω , italic_p ) are

σ(P(ρ,c,g))mi(t,x,ω,ξ)=δmiω2+ρ1(x)cijk(x)gm(x)ξjξk,𝜎superscriptsubscriptsubscript𝑃𝜌𝑐𝑔𝑚𝑖𝑡𝑥𝜔𝜉superscriptsubscript𝛿𝑚𝑖superscript𝜔2superscript𝜌1𝑥superscript𝑐𝑖𝑗𝑘𝑥subscript𝑔𝑚𝑥subscript𝜉𝑗subscript𝜉𝑘\sigma(P_{(\rho,c,g)})_{\,\,m}^{i}(t,x,\omega,\xi)=-\delta_{\,\,m}^{i}\omega^{% 2}+\rho^{-1}(x)c^{ijk\ell}(x)g_{\ell m}(x)\xi_{j}\xi_{k},italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( italic_t , italic_x , italic_ω , italic_ξ ) = - italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_x ) italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT ( italic_x ) italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( italic_x ) italic_ξ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ,

with (t,x,ω,ξ)T((0,T)×M)𝑡𝑥𝜔𝜉superscript𝑇0𝑇𝑀(t,x,\omega,\xi)\in T^{*}((0,T)\times M)( italic_t , italic_x , italic_ω , italic_ξ ) ∈ italic_T start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( ( 0 , italic_T ) × italic_M ). Introducing the slowness vector p=ω1ξ𝑝superscript𝜔1𝜉p=\omega^{-1}\xiitalic_p = italic_ω start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ξ and the Christoffel matrix Γ(x,p)Γ𝑥𝑝\Gamma(x,p)roman_Γ ( italic_x , italic_p ) defined by

Γmi(x,p):=ρ1(x)cijk(x)gm(x)pjpk,assignsuperscriptsubscriptΓ𝑚𝑖𝑥𝑝superscript𝜌1𝑥superscript𝑐𝑖𝑗𝑘𝑥subscript𝑔𝑚𝑥subscript𝑝𝑗subscript𝑝𝑘\Gamma_{\,\,m}^{i}(x,p):=\rho^{-1}(x)c^{ijk\ell}(x)g_{\ell m}(x)p_{j}p_{k},roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( italic_x , italic_p ) := italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_x ) italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT ( italic_x ) italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( italic_x ) italic_p start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ,

one can rewrite the principal symbol as

σ(P(ρ,c,g))(t,x,ω,p):=ω2[Γ(x,p)I].assign𝜎subscript𝑃𝜌𝑐𝑔𝑡𝑥𝜔𝑝superscript𝜔2delimited-[]Γ𝑥𝑝𝐼\sigma(P_{(\rho,c,g)})(t,x,\omega,p):=\omega^{2}[\Gamma(x,p)-I].italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT ) ( italic_t , italic_x , italic_ω , italic_p ) := italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ roman_Γ ( italic_x , italic_p ) - italic_I ] .

By symmetry and positivity of c𝑐citalic_c in (1)–(2) the Christoffel matrix ΓΓ\Gammaroman_Γ is symmetric and positive definite. Eigenvectors of ΓΓ\Gammaroman_Γ are called polarizations and the corresponding eigenvalues are related to the wave speeds. It turns out that the Finsler metric FaqPsuperscriptsubscript𝐹𝑎𝑞𝑃F_{a}^{qP}italic_F start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_q italic_P end_POSTSUPERSCRIPT associated to the qP𝑞𝑃qPitalic_q italic_P-polarization (the fastest waves) and the density-normalized stiffness tensor field a=ρ1c𝑎superscript𝜌1𝑐a=\rho^{-1}citalic_a = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c for the elastic geometry can be derived from the largest eigenvalue of ΓΓ\Gammaroman_Γ as described in [HILSaug21].

1.2. Dirichlet-to-Neumann maps and inverse problems

The (hyperbolic) Dirichlet-to-Neumann (DN) map is the map that sends the Dirichlet boundary condition f𝑓fitalic_f to its corresponding Neumann boundary condition. For the elastic wave operator P(ρ,c,g)subscript𝑃𝜌𝑐𝑔P_{(\rho,c,g)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT defined in (5) the associated DN operator Λ(ρ,c,g)subscriptΛ𝜌𝑐𝑔\Lambda_{(\rho,c,g)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT reads

(Λ(ρ,c,g)f)i=νjcijkgmxkum|(0,T)×M,superscriptsubscriptΛ𝜌𝑐𝑔𝑓𝑖evaluated-atsubscript𝜈𝑗superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑘superscript𝑢𝑚0𝑇𝑀(\Lambda_{(\rho,c,g)}f)^{i}=\,\nu_{j}\,c^{ijk\ell}g_{\ell m}\partial_{x^{k}}u^% {m}\Bigg{|}_{(0,T)\times\partial M},( roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT italic_f ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT = italic_ν start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT | start_POSTSUBSCRIPT ( 0 , italic_T ) × ∂ italic_M end_POSTSUBSCRIPT , (7)

where ν𝜈\nuitalic_ν is the unit outward normal to M𝑀\partial M∂ italic_M. The inverse problem is to recover ρ𝜌\rhoitalic_ρ, c𝑐citalic_c and g𝑔gitalic_g from Λ(ρ,c,g)subscriptΛ𝜌𝑐𝑔\Lambda_{(\rho,c,g)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT (or only ρ𝜌\rhoitalic_ρ and c𝑐citalic_c when g𝑔gitalic_g is assumed known).

In the Euclidean setting the DN operator Λ(ρ,c)subscriptΛ𝜌𝑐\Lambda_{(\rho,c)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT associated to P(ρ,c)subscript𝑃𝜌𝑐P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT defined in (6) simplifies to

(Λ(ρ,c)f)i=νjcijkxku|(0,T)×M.subscriptsubscriptΛ𝜌𝑐𝑓𝑖evaluated-atsubscript𝜈𝑗subscript𝑐𝑖𝑗𝑘subscriptsubscript𝑥𝑘subscript𝑢0𝑇𝑀(\Lambda_{(\rho,c)}f)_{i}=\nu_{j}\,c_{ijk\ell}\,\partial_{x_{k}}u_{\ell}\Bigg{% |}_{(0,T)\times\partial M}.( roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT italic_f ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_ν start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT | start_POSTSUBSCRIPT ( 0 , italic_T ) × ∂ italic_M end_POSTSUBSCRIPT . (8)

The inverse problem is to recover ρ𝜌\rhoitalic_ρ and c𝑐citalic_c from Λ(ρ,c)subscriptΛ𝜌𝑐\Lambda_{(\rho,c)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT. Following [HoopNakamuraZhai19] the operator Λ(ρ,c)subscriptΛ𝜌𝑐\Lambda_{(\rho,c)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT is defined as the map Λ(ρ,c):C2(0,T;H12(M))L2(0,T;H12(M)):subscriptΛ𝜌𝑐superscript𝐶20𝑇superscript𝐻12𝑀superscript𝐿20𝑇superscript𝐻12𝑀\Lambda_{(\rho,c)}:C^{2}\left(0,T;H^{\frac{1}{2}}(\partial M)\right)% \rightarrow L^{2}\left(0,T;H^{-\frac{1}{2}}(\partial M)\right)roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT : italic_C start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_H start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( ∂ italic_M ) ) → italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 0 , italic_T ; italic_H start_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( ∂ italic_M ) ). In the same lines as assessing well-posedness for the Riemmannian initial boundary value problem, the results for the Euclidean DN operator Λ(ρ,c)subscriptΛ𝜌𝑐\Lambda_{(\rho,c)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT can be extended to the Riemannian DN operator Λ(ρ,c,g)subscriptΛ𝜌𝑐𝑔\Lambda_{(\rho,c,g)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT.

This inverse problem has been addressed for the case when the stiffness tensor is isotropic [Rachele001, RACHELE002, HoopNakamuraZhai17, Stefanov2017] and for the case when the stiffness tensor is transversely isotropic or orthorhombic (with additional assumptions on the symmetry axis and symmetry planes)  [HoopNakamuraZhai19, deHoop2019]. Very recent results address the reconstruction procedure for anisotropic stiffness tensors [dehoop2023reconstruction].

For the DN maps corresponding to the Calderón problems for conductivities and Riemannian and Lorentzian metrics, see Appendix A.1.

1.3. Acknowledgements

This work was supported by the Research Council of Finland (Flagship of Advanced Mathematics for Sensing Imaging and Modelling grant 359208 and Centre of Excellence of Inverse Modelling and Imaging 353092 and other grants 351665, 351656, 358047). We thank the referee for useful suggestions and feedback.

2. Results

All material parameters, metric tensors, and diffeomorphisms are assumed to be Csuperscript𝐶C^{\infty}italic_C start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT throughout the paper. The only possibly non-smooth function is the displacement field u𝑢uitalic_u.

2.1. Dirichlet-to-Neumann maps

The point of this paper is to identify the correct gauge freedom in the inverse boundary value problems posed above. The natural gauges appear to be surprisingly different in Euclidean and Riemannian geometries, differing not only by changes of coordinates.

See section 2.2 for the definitions of the various pushforwards.

In the following we assume that the qP geometry is simple.

Theorem 2.1.

Let (M,g)𝑀𝑔(M,g)( italic_M , italic_g ) be a Riemannian manifold of any dimension n2𝑛2n\geq 2italic_n ≥ 2 with boundary, let c𝑐citalic_c be a stiffness tensor satisfying the symmetry and positivity conditions in (1)–(2), and let ρ>0𝜌0\rho>0italic_ρ > 0 be the density of mass. Let ϕ:MM:italic-ϕ𝑀𝑀\phi\colon M\to Mitalic_ϕ : italic_M → italic_M be a diffeomorphism fixing the boundary and let μ0𝜇0\mu\neq 0italic_μ ≠ 0 be a function so that μ|M=1evaluated-at𝜇𝑀1\mu|_{\partial M}=1italic_μ | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = 1. Then the DN map defined in (7) satisfies

Λ(μn2+nϕρ,μϕc,μ22+nϕg)=Λ(ρ,c,g).subscriptΛsuperscript𝜇𝑛2𝑛subscriptitalic-ϕ𝜌𝜇subscriptitalic-ϕ𝑐superscript𝜇22𝑛subscriptitalic-ϕ𝑔subscriptΛ𝜌𝑐𝑔\Lambda_{\left(\mu^{\frac{n}{2+n}}\,\phi_{*}\rho,\mu\,\phi_{*}c,\mu^{-\frac{2}% {2+n}}\,\phi_{*}g\right)}=\Lambda_{(\rho,c,g)}.roman_Λ start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_μ italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT . (9)

The previous observation concerns gauge freedoms in the DN map and on the full level of the PDE. The following recent observation concerns gauge freedoms on the principal level of the PDE:

Remark 2.2 ([IlmavirtaCocan]).

Let ΩdΩsuperscript𝑑\Omega\subset\mathbb{R}^{d}roman_Ω ⊂ blackboard_R start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT with d=2𝑑2d=2italic_d = 2 or d=3𝑑3d=3italic_d = 3 be a bounded, smooth, connected, and simply connected domain. There is an open and dense set W𝑊Witalic_W of density-normalized stiffness tensor fields on ΩΩ\Omegaroman_Ω so that if aW𝑎𝑊a\in Witalic_a ∈ italic_W and ϕFaqP=FbqPsuperscriptitalic-ϕsuperscriptsubscript𝐹𝑎𝑞𝑃superscriptsubscript𝐹𝑏𝑞𝑃\phi^{*}F_{a}^{qP}=F_{b}^{qP}italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_F start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_q italic_P end_POSTSUPERSCRIPT = italic_F start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_q italic_P end_POSTSUPERSCRIPT with a diffeomorphism ϕ:Ω¯Ω¯:italic-ϕ¯Ω¯Ω\phi\colon\bar{\Omega}\to\bar{\Omega}italic_ϕ : over¯ start_ARG roman_Ω end_ARG → over¯ start_ARG roman_Ω end_ARG fixing the boundary, then ϕ=iditalic-ϕid\phi=\operatorname{id}italic_ϕ = roman_id and b=a𝑏𝑎b=aitalic_b = italic_a.

The core of the proof of Remark 2.2 is to show that ϕitalic-ϕ\phiitalic_ϕ has to be conformal. Then it follows from fixing the boundary that it has to be the identity.

Note that the Finsler metric FaqPsuperscriptsubscript𝐹𝑎𝑞𝑃F_{a}^{qP}italic_F start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_q italic_P end_POSTSUPERSCRIPT is linked to the Christoffel matrix and thus the principal symbol of the elastic wave operator P(ρ,c)subscript𝑃𝜌𝑐P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT. Therefore this result suggests that the principal behaviour of P(ρ,c)subscript𝑃𝜌𝑐P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT determines a=ρ1c𝑎superscript𝜌1𝑐a=\rho^{-1}citalic_a = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c; cf. part (3) of Proposition 2.6. Information about the density is contained in lower order terms; cf. Lemma 2.7.

Using the fact that only conformal diffeomorphisms preserve the principal behavior of the elastic wave operator, we can derive the Euclidean version of Theorem 2.1: In the Euclidean setting g1=g2=gEsubscript𝑔1subscript𝑔2subscript𝑔𝐸g_{1}=g_{2}=g_{E}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT, so the metric component of the conclusion of Conjecture 2.2 becomes gE=μ2/(2+n)ϕgEsubscript𝑔𝐸superscript𝜇22𝑛subscriptitalic-ϕsubscript𝑔𝐸g_{E}=\mu^{-2/(2+n)}\phi_{*}g_{E}italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT = italic_μ start_POSTSUPERSCRIPT - 2 / ( 2 + italic_n ) end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT. This makes ϕitalic-ϕ\phiitalic_ϕ a conformal map of (Ω¯,gE)¯Ωsubscript𝑔𝐸(\bar{\Omega},g_{E})( over¯ start_ARG roman_Ω end_ARG , italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ) fixing the boundary, which in fact forces ϕitalic-ϕ\phiitalic_ϕ to be the identity map222Suppose there is a conformal map ϕ:Ω¯Ω¯:italic-ϕ¯Ω¯Ω\phi\colon\bar{\Omega}\to\bar{\Omega}italic_ϕ : over¯ start_ARG roman_Ω end_ARG → over¯ start_ARG roman_Ω end_ARG fixing the boundary. By the Riemann mapping theorem there is a conformal map η:DΩ:𝜂𝐷Ω\eta\colon D\to\Omegaitalic_η : italic_D → roman_Ω from the unit disc, and by smoothness of ΩΩ\Omegaroman_Ω it extends to a map η¯:D¯Ω¯:¯𝜂¯𝐷¯Ω\bar{\eta}\colon\bar{D}\to\bar{\Omega}over¯ start_ARG italic_η end_ARG : over¯ start_ARG italic_D end_ARG → over¯ start_ARG roman_Ω end_ARG. Thus η¯1ϕη¯:D¯D¯:superscript¯𝜂1italic-ϕ¯𝜂¯𝐷¯𝐷\bar{\eta}^{-1}\circ\phi\circ\bar{\eta}\colon\bar{D}\to\bar{D}over¯ start_ARG italic_η end_ARG start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ∘ italic_ϕ ∘ over¯ start_ARG italic_η end_ARG : over¯ start_ARG italic_D end_ARG → over¯ start_ARG italic_D end_ARG is a conformal map of the disc fixing its boundary. But the conformal self-maps of the disc are the Möbius transformations, and it is easily checked that only the trivial one fixes the boundary. Thus η¯1ϕη¯=idDsuperscript¯𝜂1italic-ϕ¯𝜂subscriptid𝐷\bar{\eta}^{-1}\circ\phi\circ\bar{\eta}=\operatorname{id}_{D}over¯ start_ARG italic_η end_ARG start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ∘ italic_ϕ ∘ over¯ start_ARG italic_η end_ARG = roman_id start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT and ϕ=idΩitalic-ϕsubscriptidΩ\phi=\operatorname{id}_{\Omega}italic_ϕ = roman_id start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT. and μ=1𝜇1\mu=1italic_μ = 1. This yields:

Corollary 2.3.

Let ΩdΩsuperscript𝑑\Omega\subset\mathbb{R}^{d}roman_Ω ⊂ blackboard_R start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT with d=2𝑑2d=2italic_d = 2 or d=3𝑑3d=3italic_d = 3 be a bounded, smooth, connected, and simply connected domain, let c𝑐citalic_c be a stiffness tensor field satisfying the symmetry and positivity conditions in (1)–(2), and let ρ>0𝜌0\rho>0italic_ρ > 0 be the density of mass. Then the gauge freedoms for the DN map Λ(ρ,c)subscriptΛ𝜌𝑐\Lambda_{\left(\rho,c\right)}roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT by a diffeomorphism ϕitalic-ϕ\phiitalic_ϕ or a function μ𝜇\muitalic_μ in (9) are eliminated.

As Remark 2.2 implies that the principal behavior of P(ρ,c)subscript𝑃𝜌𝑐P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT determines a=ρ1c𝑎superscript𝜌1𝑐a=\rho^{-1}citalic_a = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c uniquely, Corollary 2.3 implies that the most likely gauge freedom for the full behavior of P(ρ,c)subscript𝑃𝜌𝑐P_{(\rho,c)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT is eliminated and we usually expect more gauge freedom on the principal level than the full level (cf. Table 1), we find it likely that the following result holds true:

Conjecture 2.1.

Let ΩdΩsuperscript𝑑\Omega\subset\mathbb{R}^{d}roman_Ω ⊂ blackboard_R start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT with d=2𝑑2d=2italic_d = 2 or d=3𝑑3d=3italic_d = 3 be a bounded, smooth, connected, and simply connected domain. Let ρi>0subscript𝜌𝑖0\rho_{i}>0italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > 0 be density fields and cisubscript𝑐𝑖c_{i}italic_c start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT be generic stiffness tensor fields that satisfy the symmetry and positivity conditions in (1)–(2) for i=1,2𝑖12i=1,2italic_i = 1 , 2. If the DN maps defined in (8) satisfy Λ(ρ1,c1)=Λ(ρ2,c2)subscriptΛsubscript𝜌1subscript𝑐1subscriptΛsubscript𝜌2subscript𝑐2\Lambda_{(\rho_{1},c_{1})}=\Lambda_{(\rho_{2},c_{2})}roman_Λ start_POSTSUBSCRIPT ( italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT, then ρ1=ρ2subscript𝜌1subscript𝜌2\rho_{1}=\rho_{2}italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and c1=c2subscript𝑐1subscript𝑐2c_{1}=c_{2}italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

As a consequence we also find it likely that the gauge freedom found in Theorem 2.1 is indeed the whole gauge in the Riemannian setting:

Conjecture 2.2.

In the setting of Theorem 2.1 with generic stiffness tensor fields cisubscript𝑐𝑖c_{i}italic_c start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, density fields ρisubscript𝜌𝑖\rho_{i}italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and Riemannian metrics gisubscript𝑔𝑖g_{i}italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT the following holds: If Λ(ρ1,c1,g1)=Λ(ρ2,c2,g2)subscriptΛsubscript𝜌1subscript𝑐1subscript𝑔1subscriptΛsubscript𝜌2subscript𝑐2subscript𝑔2\Lambda_{(\rho_{1},c_{1},g_{1})}=\Lambda_{(\rho_{2},c_{2},g_{2})}roman_Λ start_POSTSUBSCRIPT ( italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT, then (ρ2,c2,g2)=(μn2+nϕρ1,μϕc1,μ22+nϕg1)subscript𝜌2subscript𝑐2subscript𝑔2superscript𝜇𝑛2𝑛subscriptitalic-ϕsubscript𝜌1𝜇subscriptitalic-ϕsubscript𝑐1superscript𝜇22𝑛subscriptitalic-ϕsubscript𝑔1(\rho_{2},c_{2},g_{2})=\left(\mu^{\frac{n}{2+n}}\,\phi_{*}\rho_{1},\mu\,\phi_{% *}c_{1},\mu^{-\frac{2}{2+n}}\,\phi_{*}g_{1}\right)( italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_μ italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) for some function μ𝜇\muitalic_μ and diffeomorphism ϕitalic-ϕ\phiitalic_ϕ.

Note how the distinction between principal and full behavior and recovery of the normalized stiffness tensor a=ρ1c𝑎superscript𝜌1𝑐a=\rho^{-1}citalic_a = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c and the full stiffness tensor should be compared with Rachele’s work [Rachele001, RACHELE002, Rachele03] in the isotropic333The stiffness tensor of an isotropic material enjoys more symmetry than we assumed, and it can be parametrized by the two Lamé parameters λ𝜆\lambdaitalic_λ and μ𝜇\muitalic_μ as cijk=λδijδk+μ(δikδj+δiδjk)subscript𝑐𝑖𝑗𝑘𝜆subscript𝛿𝑖𝑗subscript𝛿𝑘𝜇subscript𝛿𝑖𝑘subscript𝛿𝑗subscript𝛿𝑖subscript𝛿𝑗𝑘c_{ijk\ell}=\lambda\delta_{ij}\delta_{k\ell}+\mu(\delta_{ik}\delta_{j\ell}+% \delta_{i\ell}\delta_{jk})italic_c start_POSTSUBSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUBSCRIPT = italic_λ italic_δ start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_k roman_ℓ end_POSTSUBSCRIPT + italic_μ ( italic_δ start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_j roman_ℓ end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_i roman_ℓ end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_j italic_k end_POSTSUBSCRIPT ). The wave speeds for pressure and shear waves are cP=(λ+2μ)/ρsubscript𝑐𝑃𝜆2𝜇𝜌c_{P}=\sqrt{(\lambda+2\mu)/\rho}italic_c start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT = square-root start_ARG ( italic_λ + 2 italic_μ ) / italic_ρ end_ARG and cS=μ/ρsubscript𝑐𝑆𝜇𝜌c_{S}=\sqrt{\mu/\rho}italic_c start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = square-root start_ARG italic_μ / italic_ρ end_ARG. setting:

  1. (1)

    The DN map determines the travel times between boundary points for both wave speeds cPsubscript𝑐𝑃c_{P}italic_c start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT and cSsubscript𝑐𝑆c_{S}italic_c start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT.

  2. (2)

    If the manifolds (Ω¯,cP/S2gE)¯Ωsuperscriptsubscript𝑐𝑃𝑆2subscript𝑔𝐸(\bar{\Omega},c_{P/S}^{-2}g_{E})( over¯ start_ARG roman_Ω end_ARG , italic_c start_POSTSUBSCRIPT italic_P / italic_S end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ) are simple, then this geometric information determines the Riemannian metrics uniquely.

  3. (3)

    The metrics are known to be conformally Euclidean, which eliminates the coordinate gauge freedom intrinsic to a Riemannian manifold. Thus the DN map determines cPsubscript𝑐𝑃c_{P}italic_c start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT and cSsubscript𝑐𝑆c_{S}italic_c start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT.

  4. (4)

    The density ρ𝜌\rhoitalic_ρ is a lower order term and can be recovered from the DN map via a ray transform.

2.2. Various pushforwards

We define various different kinds of pushforwards, some of which are somewhat non-standard. Let U,Vn𝑈𝑉superscript𝑛U,V\subset\mathbb{R}^{n}italic_U , italic_V ⊂ blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT be any open sets and let ϕ:UV:italic-ϕ𝑈𝑉\phi\colon U\to Vitalic_ϕ : italic_U → italic_V be a diffeomorphism.

For the stiffness tensor cijksuperscript𝑐𝑖𝑗𝑘c^{ijk\ell}italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT, the metric tensor gijsubscript𝑔𝑖𝑗g_{ij}italic_g start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT and the density ρ𝜌\rhoitalic_ρ we define the pushforwards ϕcsubscriptitalic-ϕ𝑐\phi_{*}citalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c, ϕgsubscriptitalic-ϕ𝑔\phi_{*}gitalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g, and ϕρsubscriptitalic-ϕ𝜌\phi_{*}\rhoitalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ in the usual fashion:

(ϕc)i^j^k^^superscriptsubscriptitalic-ϕ𝑐^𝑖^𝑗^𝑘^\displaystyle(\phi_{*}c)^{\hat{i}\hat{j}\hat{k}\hat{\ell}}( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c ) start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG over^ start_ARG italic_j end_ARG over^ start_ARG italic_k end_ARG over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT =(cijk(Dϕ)ii^(Dϕ)jj^(Dϕ)kk^(Dϕ)^)ϕ1,absentsuperscript𝑐𝑖𝑗𝑘superscriptsubscript𝐷italic-ϕ𝑖^𝑖superscriptsubscript𝐷italic-ϕ𝑗^𝑗superscriptsubscript𝐷italic-ϕ𝑘^𝑘superscriptsubscript𝐷italic-ϕ^superscriptitalic-ϕ1\displaystyle=\left(c^{ijk\ell}(D\phi)_{\hskip 5.69054pti}^{\hat{i}}(D\phi)_{% \hskip 5.69054ptj}^{\hat{j}}(D\phi)_{\hskip 5.69054ptk}^{\hat{k}}(D\phi)_{% \hskip 5.69054pt\ell}^{\hat{\ell}}\right)\circ\phi^{-1},= ( italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_k end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , (10)
(ϕg)i^j^subscriptsubscriptitalic-ϕ𝑔^𝑖^𝑗\displaystyle(\phi_{*}g)_{\hat{i}\hat{j}}( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) start_POSTSUBSCRIPT over^ start_ARG italic_i end_ARG over^ start_ARG italic_j end_ARG end_POSTSUBSCRIPT =(gij(Dϕ1)i^i(Dϕ1)j^j)ϕ1,absentsubscript𝑔𝑖𝑗superscriptsubscript𝐷superscriptitalic-ϕ1^𝑖𝑖superscriptsubscript𝐷superscriptitalic-ϕ1^𝑗𝑗superscriptitalic-ϕ1\displaystyle=\left(g_{ij}(D\phi^{-1})_{\hskip 5.69054pt\hat{i}}^{i}(D\phi^{-1% })_{\hskip 5.69054pt\hat{j}}^{j}\right)\circ\phi^{-1},= ( italic_g start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_D italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT over^ start_ARG italic_i end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( italic_D italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , (11)
ϕρsubscriptitalic-ϕ𝜌\displaystyle\phi_{*}\rhoitalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ =ρϕ1.absent𝜌superscriptitalic-ϕ1\displaystyle=\rho\circ\phi^{-1}.= italic_ρ ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT . (12)

All the familiar formulas from differential geometry hold for these standard pushforwards — or, equivalently, pullbacks over ϕ1superscriptitalic-ϕ1\phi^{-1}italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT.

When the map ϕitalic-ϕ\phiitalic_ϕ is conformal with respect to the Euclidean metric, we define

(ϕc)i^j^k^^superscriptsubscriptitalic-ϕ𝑐^𝑖^𝑗^𝑘^\displaystyle(\phi_{\square}c)^{\hat{i}\hat{j}\hat{k}\hat{\ell}}( italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c ) start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG over^ start_ARG italic_j end_ARG over^ start_ARG italic_k end_ARG over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT =(det(Dϕ)12/ncijk(Dϕ)ii^(Dϕ)jj^(Dϕ)kk^(Dϕ)^)ϕ1,absentsuperscript𝐷italic-ϕ12𝑛superscript𝑐𝑖𝑗𝑘superscriptsubscript𝐷italic-ϕ𝑖^𝑖superscriptsubscript𝐷italic-ϕ𝑗^𝑗superscriptsubscript𝐷italic-ϕ𝑘^𝑘superscriptsubscript𝐷italic-ϕ^superscriptitalic-ϕ1\displaystyle=\left(\det(D\phi)^{-1-2/n}c^{ijk\ell}(D\phi)_{\hskip 5.69054pti}% ^{\hat{i}}(D\phi)_{\hskip 5.69054ptj}^{\hat{j}}(D\phi)_{\hskip 5.69054ptk}^{% \hat{k}}(D\phi)_{\hskip 5.69054pt\ell}^{\hat{\ell}}\right)\circ\phi^{-1},= ( roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 1 - 2 / italic_n end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_k end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , (13)
ϕρsubscriptitalic-ϕ𝜌\displaystyle\phi_{\square}\rhoitalic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_ρ =(det(Dϕ)1ρ)ϕ1absentsuperscript𝐷italic-ϕ1𝜌superscriptitalic-ϕ1\displaystyle=\left(\det(D\phi)^{-1}\rho\right)\circ\phi^{-1}= ( roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ρ ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (14)

and

(ϕc)i^j^k^^superscriptsubscriptitalic-ϕ𝑐^𝑖^𝑗^𝑘^\displaystyle(\phi_{\blacksquare}c)^{\hat{i}\hat{j}\hat{k}\hat{\ell}}( italic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT italic_c ) start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG over^ start_ARG italic_j end_ARG over^ start_ARG italic_k end_ARG over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT =(det(Dϕ)2/ncijk(Dϕ)ii^(Dϕ)jj^(Dϕ)kk^(Dϕ)^)ϕ1.absentsuperscript𝐷italic-ϕ2𝑛superscript𝑐𝑖𝑗𝑘superscriptsubscript𝐷italic-ϕ𝑖^𝑖superscriptsubscript𝐷italic-ϕ𝑗^𝑗superscriptsubscript𝐷italic-ϕ𝑘^𝑘superscriptsubscript𝐷italic-ϕ^superscriptitalic-ϕ1\displaystyle=\left(\det(D\phi)^{-2/n}c^{ijk\ell}(D\phi)_{\hskip 5.69054pti}^{% \hat{i}}(D\phi)_{\hskip 5.69054ptj}^{\hat{j}}(D\phi)_{\hskip 5.69054ptk}^{\hat% {k}}(D\phi)_{\hskip 5.69054pt\ell}^{\hat{\ell}}\right)\circ\phi^{-1}.= ( roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 2 / italic_n end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_k end_ARG end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT . (15)

It is important to note that the usual pushforward ϕsubscriptitalic-ϕ\phi_{*}italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is defined on the triplet (ρ,c,g)𝜌𝑐𝑔(\rho,c,g)( italic_ρ , italic_c , italic_g ), the pushforward ϕsubscriptitalic-ϕ\phi_{\square}italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT on the pair (ρ,c)𝜌𝑐(\rho,c)( italic_ρ , italic_c ) and pushforward ϕsubscriptitalic-ϕ\phi_{\blacksquare}italic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT only on c𝑐citalic_c.

Remark 2.4.

If we want to define a pushforward over a possibly non-conformal diffeomorphism ϕitalic-ϕ\phiitalic_ϕ, then the definitions of ϕ(c,ρ)subscriptitalic-ϕ𝑐𝜌\phi_{\square}(c,\rho)italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT ( italic_c , italic_ρ ) and ϕcsubscriptitalic-ϕ𝑐\phi_{\blacksquare}citalic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT italic_c need to be adjusted. The correct formulas are obtained by treating the stiffness tensor as the mixed rank object cijksubscriptsuperscript𝑐𝑖𝑗𝑘c^{ijk}_{\phantom{ijk}\ell}italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT and applying the usual formulas from differential geometry. Then the EWE is indeed invariant, but the issue is that the so obtained ϕcsubscriptitalic-ϕ𝑐\phi_{\square}citalic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c no longer satisfies the symmetry conditions. It is shown in the proof of Remark 2.2 in [IlmavirtaCocan] that only conformal distortions can give rise to a valid slowness polynomial that corresponds to a stiffness tensor with the required symmetries. For a conformal map ϕitalic-ϕ\phiitalic_ϕ this pushforward does preserve symmetry, and it gives rise exactly to our equation (13).

One way to see this is as follows: The usual Euclidean EWE is not coordinate invariant because it is tied to the Euclidean concept of distance. Conformal coordinate invariance only follows because Lemma 2.7 allows shifting scalar factors within the triplet (ρ,c,g)𝜌𝑐𝑔(\rho,c,g)( italic_ρ , italic_c , italic_g ).

Remark 2.5.

The elastic wave operator (3) on a Riemannian manifold can be written as

(Δcgu(t,x))i=xj(cijkgmxkum(t,x))superscriptsuperscriptsubscriptΔ𝑐𝑔𝑢𝑡𝑥𝑖subscriptsuperscript𝑥𝑗superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑘superscript𝑢𝑚𝑡𝑥(\Delta_{c}^{g}\,u(t,x))^{i}=\nabla_{x^{j}}\left(c^{ijk\ell}g_{\ell m}\nabla_{% x^{k}}u^{m}(t,x)\right)( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT italic_u ( italic_t , italic_x ) ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT = ∇ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_t , italic_x ) ) (16)

using covariant derivatives. Formula (3) and (16) are equivalent; the only difference is in the coordinate representation of the covariant divergence. Note that these two ways to define the elastic wave operator correspond to the following equivalent definition for the Laplace-Beltrami operator: Δgf=1detgxi(detggijxjf)subscriptΔ𝑔𝑓1𝑔subscriptsuperscript𝑥𝑖𝑔superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑓\Delta_{g}f=\frac{1}{\sqrt{\det g}}\partial_{x^{i}}\left(\sqrt{\det g}g^{ij}% \partial_{x^{j}}f\right)roman_Δ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_f = divide start_ARG 1 end_ARG start_ARG square-root start_ARG roman_det italic_g end_ARG end_ARG ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( square-root start_ARG roman_det italic_g end_ARG italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_f ) and Δgf=xi(gijxjf)subscriptΔ𝑔𝑓subscriptsuperscript𝑥𝑖superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑓\Delta_{g}f=\nabla_{x^{i}}\left(g^{ij}\nabla_{x^{j}}f\right)roman_Δ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_f = ∇ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_f ) respectively.

2.3. Invariance of the elastic wave equation

The following results concern local behavior of the EWE where the boundary is disregarded. These results address both the Euclidean and Riemannian setting and the global results can be derived from these.

By Remark 2.2, in most cases any diffeomorphism that preserves the Finsler metric arising from elasticity in dimension two and three has to be conformal. In the following proposition we write down explicitly how the Euclidean pushforward ϕcsubscriptitalic-ϕ𝑐\phi_{\blacksquare}citalic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT italic_c by a diffeomorphism ϕ:MM:italic-ϕ𝑀𝑀\phi\colon M\to Mitalic_ϕ : italic_M → italic_M is defined so that the principal behavior of the Euclidean EWE is preserved.

Proposition 2.6.

The Euclidean and Riemannian elastic wave equations have the following invariance properties:

  1. (1)

    P(ρ,c,g)=ϕP(ϕρ,ϕc,ϕg)ϕsubscript𝑃𝜌𝑐𝑔superscriptitalic-ϕsubscript𝑃subscriptitalic-ϕ𝜌subscriptitalic-ϕ𝑐subscriptitalic-ϕ𝑔subscriptitalic-ϕP_{(\rho,c,g)}=\phi^{*}P_{(\phi_{*}\rho,\phi_{*}c,\phi_{*}g)}\phi_{*}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT = italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT for any diffeomorphism ϕitalic-ϕ\phiitalic_ϕ.

  2. (2)

    P(ρ,c)=ϕP(ϕρ,ϕc)ϕsubscript𝑃𝜌𝑐superscriptitalic-ϕsubscript𝑃subscriptitalic-ϕ𝜌subscriptitalic-ϕ𝑐subscriptitalic-ϕP_{(\rho,c)}=\phi^{*}P_{(\phi_{\square}\rho,\phi_{\square}c)}\phi_{*}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT = italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_ρ , italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c ) end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT for any conformal map ϕitalic-ϕ\phiitalic_ϕ.

  3. (3)

    σ(P(ρ,c))=σ(P(1,ρ1c))=ϕσ(P(1,ϕ(ρ1c)))ϕ𝜎subscript𝑃𝜌𝑐𝜎subscript𝑃1superscript𝜌1𝑐superscriptitalic-ϕ𝜎subscript𝑃1subscriptitalic-ϕsuperscript𝜌1𝑐subscriptitalic-ϕ\sigma(P_{(\rho,c)})=\sigma(P_{(1,\rho^{-1}c)})=\phi^{*}\sigma(P_{(1,\phi_{% \blacksquare}(\rho^{-1}c))})\phi_{*}italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT ) = italic_σ ( italic_P start_POSTSUBSCRIPT ( 1 , italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c ) end_POSTSUBSCRIPT ) = italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_σ ( italic_P start_POSTSUBSCRIPT ( 1 , italic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT ( italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c ) ) end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT for any conformal map ϕitalic-ϕ\phiitalic_ϕ.

The pullbacks ϕsuperscriptitalic-ϕ\phi^{*}italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT and pushforwards ϕsubscriptitalic-ϕ\phi_{*}italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT surrounding the operators simply correspond to how the solution vector fields transform.

We only prove that in the Euclidean setting conformal diffeomorphisms preserve the desired structure of the stiffness tensor. For the other direction — that only conformal maps preserve such the symmetry — we refer to the theorems and conjectures of section 2.1.

Proposition 2.6 is largely based on the following lemma on scaling freedoms which may also be of independent interest:

Lemma 2.7.

The operator ΔcgsuperscriptsubscriptΔ𝑐𝑔\Delta_{c}^{g}\,roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT satisfies the following relationship for the scalar functions μ𝜇\muitalic_μ and λ𝜆\lambdaitalic_λ which satisfy μ0𝜇0\mu\neq 0italic_μ ≠ 0 and λ0𝜆0\lambda\neq 0italic_λ ≠ 0:

(λμρ)1(Δμcλg)mi=ρ1(Δcg)mi+Qmi,superscript𝜆𝜇𝜌1subscriptsuperscriptsuperscriptsubscriptΔ𝜇𝑐𝜆𝑔𝑖𝑚superscript𝜌1subscriptsuperscriptsuperscriptsubscriptΔ𝑐𝑔𝑖𝑚subscriptsuperscript𝑄𝑖𝑚(\lambda\mu\rho)^{-1}(\Delta_{\mu c}^{\lambda g}\,)^{i}_{\hskip 2.84526ptm}=% \rho^{-1}(\Delta_{c}^{g}\,)^{i}_{\hskip 2.84526ptm}+Q^{i}_{\hskip 2.84526ptm},( italic_λ italic_μ italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_μ italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_λ italic_g end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT + italic_Q start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ,

where

Qmi=λ2+n2μ1ρ1cijkgmxj(μλ2+n2)xk.subscriptsuperscript𝑄𝑖𝑚superscript𝜆2𝑛2superscript𝜇1superscript𝜌1superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑗𝜇superscript𝜆2𝑛2subscriptsuperscript𝑥𝑘Q^{i}_{\hskip 2.84526ptm}=\lambda^{-\frac{2+n}{2}}\mu^{-1}\rho^{-1}c^{ijk\ell}% \,g_{\ell m}\partial_{x^{j}}\left(\mu\lambda^{\frac{2+n}{2}}\right)\partial_{x% ^{k}}.italic_Q start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = italic_λ start_POSTSUPERSCRIPT - divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_μ italic_λ start_POSTSUPERSCRIPT divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ) ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT .

Table 1 lists the consequences of Lemma 2.7, where we distinguish between principal and full behavior of the PDE and between the Euclidean and the Riemannian setting. For the full behavior we list the case when μλ2+n2𝜇superscript𝜆2𝑛2\mu\lambda^{\frac{2+n}{2}}italic_μ italic_λ start_POSTSUPERSCRIPT divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT is a constant so that Qu=0𝑄𝑢0Qu=0italic_Q italic_u = 0.

Euclidean EWE Riemannian EWE
Principal behavior One conformal freedom μ𝜇\muitalic_μ: (ρ,c)(μρ,μc)similar-to𝜌𝑐𝜇𝜌𝜇𝑐(\rho,c)\sim(\mu\rho,\mu c)( italic_ρ , italic_c ) ∼ ( italic_μ italic_ρ , italic_μ italic_c ) Two conformal freedoms μ,λ𝜇𝜆\mu,\lambdaitalic_μ , italic_λ: (ρ,c,g)(λμρ,μc,λg)similar-to𝜌𝑐𝑔𝜆𝜇𝜌𝜇𝑐𝜆𝑔(\rho,c,g)\sim(\lambda\mu\rho,\mu c,\lambda g)( italic_ρ , italic_c , italic_g ) ∼ ( italic_λ italic_μ italic_ρ , italic_μ italic_c , italic_λ italic_g )
Full behavior No conformal freedom One conformal freedom μ𝜇\muitalic_μ: (ρ,c,g)(μn2+nρ,μc,μ22+ng)similar-to𝜌𝑐𝑔superscript𝜇𝑛2𝑛𝜌𝜇𝑐superscript𝜇22𝑛𝑔(\rho,c,g)\sim(\mu^{\frac{n}{2+n}}\rho,\mu c,\mu^{-\frac{2}{2+n}}g)( italic_ρ , italic_c , italic_g ) ∼ ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ , italic_μ italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g )
Table 1. Consequences of Lemma 2.7 for inverse problems. The relation similar-to\sim means that the two parameters give rise to the same boundary data on ΩΩ\partial\Omega∂ roman_Ω or M𝑀\partial M∂ italic_M. These are scaling freedoms, not diffeomorphism freedoms; the freedom to choose coordinates in the Euclidean setting is eliminated by the demand that boundary be fixed by the diffeomorphism.

3. Proofs

3.1. Proofs of auxiliary results

Proof of Lemma 2.7.

A calculation gives

(λμρ)1(Δμcλgu)isuperscript𝜆𝜇𝜌1superscriptsuperscriptsubscriptΔ𝜇𝑐𝜆𝑔𝑢𝑖\displaystyle(\lambda\mu\rho)^{-1}\left(\Delta_{\mu c}^{\lambda g}\,u\right)^{i}( italic_λ italic_μ italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_μ italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_λ italic_g end_POSTSUPERSCRIPT italic_u ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT
=λ1μ1ρ11det(λg)xj(det(λg)μcijkλgmxkum)absentsuperscript𝜆1superscript𝜇1superscript𝜌11𝜆𝑔subscriptsuperscript𝑥𝑗𝜆𝑔𝜇superscript𝑐𝑖𝑗𝑘𝜆subscript𝑔𝑚subscriptsuperscript𝑥𝑘superscript𝑢𝑚\displaystyle\hskip 17.07164pt=\lambda^{-1}\mu^{-1}\rho^{-1}\frac{1}{\sqrt{% \det(\lambda g)}}\partial_{x^{j}}\left(\sqrt{\det(\lambda g)}\,\mu\,c^{ijk\ell% }\,\lambda\,g_{\ell m}\,\partial_{x^{k}}u^{m}\right)= italic_λ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG square-root start_ARG roman_det ( italic_λ italic_g ) end_ARG end_ARG ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( square-root start_ARG roman_det ( italic_λ italic_g ) end_ARG italic_μ italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_λ italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT )
=ρ11det(g)j=1nxj(det(g)cijkgmxkum)=ρ1(Δcgu)iabsentsubscriptsuperscript𝜌11𝑔superscriptsubscript𝑗1𝑛subscriptsuperscript𝑥𝑗𝑔superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑘superscript𝑢𝑚absentsuperscript𝜌1superscriptsuperscriptsubscriptΔ𝑐𝑔𝑢𝑖\displaystyle\hskip 17.07164pt=\underbrace{\rho^{-1}\frac{1}{\sqrt{\det(g)}}\,% \sum_{j=1}^{n}\partial_{x^{j}}\left(\sqrt{\det(g)}\,c^{ijk\ell}\,g_{\ell m}\,% \partial_{x^{k}}u^{m}\right)}_{=\rho^{-1}(\Delta_{c}^{g}\,u)^{i}}= under⏟ start_ARG italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG square-root start_ARG roman_det ( italic_g ) end_ARG end_ARG ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( square-root start_ARG roman_det ( italic_g ) end_ARG italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ) end_ARG start_POSTSUBSCRIPT = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT italic_u ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT
+λ2+n2μ1ρ1j=1ncijkgmxj(μλ2+n2)xkum=(Qu)i.subscriptsuperscript𝜆2𝑛2superscript𝜇1superscript𝜌1superscriptsubscript𝑗1𝑛superscript𝑐𝑖𝑗𝑘subscript𝑔𝑚subscriptsuperscript𝑥𝑗𝜇superscript𝜆2𝑛2subscriptsuperscript𝑥𝑘superscript𝑢𝑚absentsuperscript𝑄𝑢𝑖\displaystyle\hskip 34.1433pt+\underbrace{\lambda^{-\frac{2+n}{2}}\mu^{-1}\rho% ^{-1}\sum_{j=1}^{n}c^{ijk\ell}\,g_{\ell m}\partial_{x^{j}}\left(\mu\lambda^{% \frac{2+n}{2}}\right)\,\partial_{x^{k}}u^{m}}_{=(Qu)^{i}}.+ under⏟ start_ARG italic_λ start_POSTSUPERSCRIPT - divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_μ italic_λ start_POSTSUPERSCRIPT divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ) ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT end_ARG start_POSTSUBSCRIPT = ( italic_Q italic_u ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT .

as claimed. ∎

Consequences of Lemma 2.7 as listed in Table 1.

The elliptic term (λμρ)1(Δμcλgu)isuperscript𝜆𝜇𝜌1superscriptsuperscriptsubscriptΔ𝜇𝑐𝜆𝑔𝑢𝑖(\lambda\mu\rho)^{-1}\left(\Delta_{\mu c}^{\lambda g}\,u\right)^{i}( italic_λ italic_μ italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_μ italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_λ italic_g end_POSTSUPERSCRIPT italic_u ) start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT corresponds to the operator P(λμρ,μc,λg)subscript𝑃𝜆𝜇𝜌𝜇𝑐𝜆𝑔P_{(\lambda\mu\rho,\mu c,\lambda g)}italic_P start_POSTSUBSCRIPT ( italic_λ italic_μ italic_ρ , italic_μ italic_c , italic_λ italic_g ) end_POSTSUBSCRIPT. As Qmisubscriptsuperscript𝑄𝑖𝑚Q^{i}_{\hskip 2.84526ptm}italic_Q start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT is a lower order term it follows directly that the principal symbols of the operators P(λμρ,μc,λg)subscript𝑃𝜆𝜇𝜌𝜇𝑐𝜆𝑔P_{(\lambda\mu\rho,\mu c,\lambda g)}italic_P start_POSTSUBSCRIPT ( italic_λ italic_μ italic_ρ , italic_μ italic_c , italic_λ italic_g ) end_POSTSUBSCRIPT and P(ρ,c,g)subscript𝑃𝜌𝑐𝑔P_{(\rho,c,g)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT are preserved:

σ(P(λμρ,μc,λg))=σ(P(ρ,c,g))𝜎subscript𝑃𝜆𝜇𝜌𝜇𝑐𝜆𝑔𝜎subscript𝑃𝜌𝑐𝑔\sigma(P_{(\lambda\mu\rho,\mu c,\lambda g)})=\sigma(P_{(\rho,c,g)})italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_λ italic_μ italic_ρ , italic_μ italic_c , italic_λ italic_g ) end_POSTSUBSCRIPT ) = italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT ) (17)

as highlighted in the upper right corner of Table 1. Conformal freedoms for the full Riemannian EWE can only be obtained in the case when Qu=0𝑄𝑢0Qu=0italic_Q italic_u = 0. Hence, when the term μλ2+n2𝜇superscript𝜆2𝑛2\mu\lambda^{\frac{2+n}{2}}italic_μ italic_λ start_POSTSUPERSCRIPT divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT is constant. Choosing λ=μ22+n𝜆superscript𝜇22𝑛\lambda=\mu^{-\frac{2}{2+n}}italic_λ = italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT implies μλ2+n2=1𝜇superscript𝜆2𝑛21\mu\lambda^{\frac{2+n}{2}}=1italic_μ italic_λ start_POSTSUPERSCRIPT divide start_ARG 2 + italic_n end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT = 1. Hence, in this case

(μn2+nρ)1Δμcμ22+ng=ρ1(Δcgu)superscriptsuperscript𝜇𝑛2𝑛𝜌1superscriptsubscriptΔ𝜇𝑐superscript𝜇22𝑛𝑔superscript𝜌1superscriptsubscriptΔ𝑐𝑔𝑢\left(\mu^{\frac{n}{2+n}}\rho\right)^{-1}\Delta_{\mu c}^{\mu^{-\frac{2}{2+n}}g% }\,=\rho^{-1}(\Delta_{c}^{g}\,u)( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_μ italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g end_POSTSUPERSCRIPT italic_u ) (18)

and thus

P(μn2+nρ,μc,μ22+ng)=P(ρ,c,g)subscript𝑃superscript𝜇𝑛2𝑛𝜌𝜇𝑐superscript𝜇22𝑛𝑔subscript𝑃𝜌𝑐𝑔P_{\left(\mu^{\frac{n}{2+n}}\rho,\mu c,\mu^{-\frac{2}{2+n}}g\right)}=P_{(\rho,% c,g)}italic_P start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ , italic_μ italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g ) end_POSTSUBSCRIPT = italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT (19)

as highlighted in the lower right corner of Table 1.

In the Euclidean case the metric is fixed; this eliminates the conformal freedom corresponding to the function λ𝜆\lambdaitalic_λ:

(μρ)1(ΔμcgE)i=ρ1(ΔcgE)i+Qi,superscript𝜇𝜌1subscriptsuperscriptsubscriptΔ𝜇𝑐subscript𝑔𝐸𝑖superscript𝜌1subscriptsuperscriptsubscriptΔ𝑐subscript𝑔𝐸𝑖subscript𝑄𝑖(\mu\rho)^{-1}(\Delta_{\mu c}^{g_{E}}\,)_{i\ell}=\rho^{-1}(\Delta_{c}^{g_{E}}% \,)_{i\ell}+Q_{i\ell},( italic_μ italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_μ italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT italic_i roman_ℓ end_POSTSUBSCRIPT = italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT italic_i roman_ℓ end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_i roman_ℓ end_POSTSUBSCRIPT ,

with

Qi=μ1ρ1cijkxj(μ)xk.subscript𝑄𝑖superscript𝜇1superscript𝜌1subscript𝑐𝑖𝑗𝑘subscriptsuperscript𝑥𝑗𝜇subscriptsuperscript𝑥𝑘Q_{i\ell}=\mu^{-1}\rho^{-1}c_{ijk\ell}\partial_{x^{j}}\left(\mu\right)\partial% _{x^{k}}.italic_Q start_POSTSUBSCRIPT italic_i roman_ℓ end_POSTSUBSCRIPT = italic_μ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_i italic_j italic_k roman_ℓ end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_μ ) ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT end_POSTSUBSCRIPT .

It follows that the principal behavior is preserved for the pairs (μρ,μc)𝜇𝜌𝜇𝑐(\mu\rho,\mu c)( italic_μ italic_ρ , italic_μ italic_c ) and (ρ,c)𝜌𝑐(\rho,c)( italic_ρ , italic_c ) so that

σ(P(μρ,μc))=σ(P(ρ,c)).𝜎subscript𝑃𝜇𝜌𝜇𝑐𝜎subscript𝑃𝜌𝑐\sigma(P_{(\mu\rho,\mu c)})=\sigma(P_{(\rho,c)}).italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_μ italic_ρ , italic_μ italic_c ) end_POSTSUBSCRIPT ) = italic_σ ( italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c ) end_POSTSUBSCRIPT ) . (20)

As the full Euclidean EWE can only be preserved when μ𝜇\muitalic_μ is constant it follows that there are no conformal freedoms in this case. These observations are highlighted in the upper and lower left corner respectively of Table 1. ∎

Proof of Proposition 2.6.

Part (1): All structures used in this claim are tensor fields and covariant derivatives in the sense of Riemannian geometry, and therefore they are invariant under a change of coordinates by design. Verification using equations (10), (11) and (12) by hand is also straightforward.

Part (2): Note that a conformal map ϕitalic-ϕ\phiitalic_ϕ with respect to the Euclidean geometry satisfies

ϕgE=det(Dϕ)2/ngE.subscriptitalic-ϕsubscript𝑔𝐸superscript𝐷italic-ϕ2𝑛subscript𝑔𝐸\phi_{*}g_{E}=\det(D\phi)^{-2/n}g_{E}.italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT = roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 2 / italic_n end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT . (21)

Such a map is conformal whenever the Jacobian Dϕ𝐷italic-ϕD\phiitalic_D italic_ϕ at each point is a scalar times a rotation matrix and the formula (21) can for instance be derived from (11) exploiting that fact. Combining this with part (1) gives

ϕ(ρ1ΔcgE)u=(ϕρ)1Δϕcdet(Dϕ)2/ngEϕusubscriptitalic-ϕsuperscript𝜌1superscriptsubscriptΔ𝑐subscript𝑔𝐸𝑢superscriptsubscriptitalic-ϕ𝜌1superscriptsubscriptΔsubscriptitalic-ϕ𝑐superscript𝐷italic-ϕ2𝑛subscript𝑔𝐸subscriptitalic-ϕ𝑢\phi_{*}\left(\rho^{-1}\Delta_{c}^{g_{E}}\,\right)u=(\phi_{*}\rho)^{-1}\Delta_% {\phi_{*}c}^{\det(D\phi)^{-2/n}g_{E}}\,\phi_{*}uitalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) italic_u = ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 2 / italic_n end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u (22)

for any u𝑢uitalic_u.

We then apply Lemma 2.7 with λ=det(Dϕ)2/n𝜆superscript𝐷italic-ϕ2𝑛\lambda=\det(D\phi)^{-2/n}italic_λ = roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 2 / italic_n end_POSTSUPERSCRIPT and μ=det(Dϕ)1+2/n𝜇superscript𝐷italic-ϕ12𝑛\mu=\det(D\phi)^{1+2/n}italic_μ = roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT 1 + 2 / italic_n end_POSTSUPERSCRIPT (so that μλn/2+1=1𝜇superscript𝜆𝑛211\mu\lambda^{n/2+1}=1italic_μ italic_λ start_POSTSUPERSCRIPT italic_n / 2 + 1 end_POSTSUPERSCRIPT = 1) to get

(ϕρ)1Δϕcdet(Dϕ)2/ngE=det(Dϕ)(ϕρ)1Δdet(Dϕ)12/nϕcgE.superscriptsubscriptitalic-ϕ𝜌1superscriptsubscriptΔsubscriptitalic-ϕ𝑐superscript𝐷italic-ϕ2𝑛subscript𝑔𝐸𝐷italic-ϕsuperscriptsubscriptitalic-ϕ𝜌1superscriptsubscriptΔsuperscript𝐷italic-ϕ12𝑛subscriptitalic-ϕ𝑐subscript𝑔𝐸(\phi_{*}\rho)^{-1}\Delta_{\phi_{*}c}^{\det(D\phi)^{-2/n}g_{E}}\,=\det(D\phi)(% \phi_{*}\rho)^{-1}\Delta_{\det(D\phi)^{-1-2/n}\phi_{*}c}^{g_{E}}\,.( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 2 / italic_n end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = roman_det ( italic_D italic_ϕ ) ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 1 - 2 / italic_n end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT . (23)

The definitions in (13) and (14) were set up so that ϕc=det(Dϕ)12/nϕcsubscriptitalic-ϕ𝑐superscript𝐷italic-ϕ12𝑛subscriptitalic-ϕ𝑐\phi_{\square}c=\det(D\phi)^{-1-2/n}\phi_{*}citalic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c = roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 1 - 2 / italic_n end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c and ϕρ=det(Dϕ)1ϕρsubscriptitalic-ϕ𝜌superscript𝐷italic-ϕ1subscriptitalic-ϕ𝜌\phi_{\square}\rho=\det(D\phi)^{-1}\phi_{*}\rhoitalic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_ρ = roman_det ( italic_D italic_ϕ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ. Therefore

ρ1ΔcgEu=ϕ((ϕρ)1ΔϕcgE)ϕusuperscript𝜌1superscriptsubscriptΔ𝑐subscript𝑔𝐸𝑢superscriptitalic-ϕsuperscriptsubscriptitalic-ϕ𝜌1superscriptsubscriptΔsubscriptitalic-ϕ𝑐subscript𝑔𝐸subscriptitalic-ϕ𝑢\rho^{-1}\Delta_{c}^{g_{E}}\,u=\phi^{*}\left((\phi_{\square}\rho)^{-1}\Delta_{% \phi_{\square}c}^{g_{E}}\,\right)\phi_{*}uitalic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_u = italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( ( italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_ρ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ start_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u (24)

from which the desired claim follows.

Part (3): This is similar to part (2) above, but we no longer need to satisfy μλn/2+1=1𝜇superscript𝜆𝑛211\mu\lambda^{n/2+1}=1italic_μ italic_λ start_POSTSUPERSCRIPT italic_n / 2 + 1 end_POSTSUPERSCRIPT = 1 as the subprincipal term of Lemma 2.7 is irrelevant. This corresponds to the principal symbol only depending on the density normalized stiffness tensor ρ1csuperscript𝜌1𝑐\rho^{-1}citalic_ρ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_c, which is evident. Our definition in (15) satisfies

ϕc=(ϕ1)1ϕcsubscriptitalic-ϕ𝑐superscriptsubscriptitalic-ϕ11subscriptitalic-ϕ𝑐\phi_{\blacksquare}c=(\phi_{\square}1)^{-1}\phi_{\square}citalic_ϕ start_POSTSUBSCRIPT ■ end_POSTSUBSCRIPT italic_c = ( italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT 1 ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT □ end_POSTSUBSCRIPT italic_c (25)

with the constant density 1111, so the claim follows from essentially the same calculation as above. ∎

3.2. Proof of the main result

Proof of Theorem 2.1.

As emphasized in part (1) of Proposition 2.6, the pushforwards of c𝑐citalic_c, g𝑔gitalic_g and ρ𝜌\rhoitalic_ρ defined in (10)–(12) and the pushforward ϕusubscriptitalic-ϕ𝑢\phi_{*}uitalic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u of the solution leave the Riemannian EWE invariant with respect to any diffeomorphism ϕitalic-ϕ\phiitalic_ϕ. Additionally we note that the DN map is invariant with respect to coordinate changes and with respect to the push forwards of ρ,c𝜌𝑐\rho,citalic_ρ , italic_c and g𝑔gitalic_g such a coordinate change can be expressed by:

(Λ(ϕρ,ϕc,ϕg)f)i^=(ϕν)j^(ϕc)i^j^k^^(ϕg)^m^xk^(ϕu)m^superscriptsubscriptΛsubscriptitalic-ϕ𝜌subscriptitalic-ϕ𝑐subscriptitalic-ϕ𝑔𝑓^𝑖subscriptsubscriptitalic-ϕ𝜈^𝑗superscriptsubscriptitalic-ϕ𝑐^𝑖^𝑗^𝑘^subscriptsubscriptitalic-ϕ𝑔^^𝑚subscriptsuperscript𝑥^𝑘superscriptsubscriptitalic-ϕ𝑢^𝑚(\Lambda_{\left(\phi_{*}\rho,\phi_{*}c,\phi_{*}g\right)}f)^{\hat{i}}=(\phi_{*}% \nu)_{\hat{j}}(\phi_{*}c)^{\hat{i}\hat{j}\hat{k}\hat{\ell}}(\phi_{*}g)_{\hat{% \ell}\hat{m}}\partial_{x^{\hat{k}}}(\phi_{*}u)^{\hat{m}}( roman_Λ start_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT italic_f ) start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG end_POSTSUPERSCRIPT = ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ν ) start_POSTSUBSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c ) start_POSTSUPERSCRIPT over^ start_ARG italic_i end_ARG over^ start_ARG italic_j end_ARG over^ start_ARG italic_k end_ARG over^ start_ARG roman_ℓ end_ARG end_POSTSUPERSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) start_POSTSUBSCRIPT over^ start_ARG roman_ℓ end_ARG over^ start_ARG italic_m end_ARG end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT over^ start_ARG italic_k end_ARG end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u ) start_POSTSUPERSCRIPT over^ start_ARG italic_m end_ARG end_POSTSUPERSCRIPT

where,

(ϕν)j^subscriptsubscriptitalic-ϕ𝜈^𝑗\displaystyle(\phi_{*}\nu)_{\hat{j}}( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ν ) start_POSTSUBSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUBSCRIPT =(νj(Dϕ1)j^j)ϕ1,absentsubscript𝜈𝑗superscriptsubscript𝐷superscriptitalic-ϕ1^𝑗𝑗superscriptitalic-ϕ1\displaystyle=\left(\nu_{j}(D\phi^{-1})_{\hskip 5.69054pt\hat{j}}^{j}\right)% \circ\phi^{-1},= ( italic_ν start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_D italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT over^ start_ARG italic_j end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ,
(ϕu)m^superscriptsubscriptitalic-ϕ𝑢^𝑚\displaystyle(\phi_{*}u)^{\hat{m}}( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u ) start_POSTSUPERSCRIPT over^ start_ARG italic_m end_ARG end_POSTSUPERSCRIPT =(um(Dϕ)mm^)ϕ1.absentsuperscript𝑢𝑚superscriptsubscript𝐷italic-ϕ𝑚^𝑚superscriptitalic-ϕ1\displaystyle=\left(u^{m}(D\phi)_{\hskip 5.69054ptm}^{\hat{m}}\right)\circ\phi% ^{-1}.= ( italic_u start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_D italic_ϕ ) start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over^ start_ARG italic_m end_ARG end_POSTSUPERSCRIPT ) ∘ italic_ϕ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT .

If ϕitalic-ϕ\phiitalic_ϕ is a diffeomorphism that fixes the boundary and thus is the identity in a neighborhood of M𝑀\partial M∂ italic_M one can assume Dϕ(x)=I𝐷italic-ϕ𝑥𝐼D\phi(x)=Iitalic_D italic_ϕ ( italic_x ) = italic_I for x𝑥xitalic_x in that neighborhood. This implies that ϕρ(x)=ρ(x),ϕc(x)=c(x),ϕg(x)=g(x),ϕν(x)=ν(x)formulae-sequencesubscriptitalic-ϕ𝜌𝑥𝜌𝑥formulae-sequencesubscriptitalic-ϕ𝑐𝑥𝑐𝑥formulae-sequencesubscriptitalic-ϕ𝑔𝑥𝑔𝑥subscriptitalic-ϕ𝜈𝑥𝜈𝑥\phi_{*}\rho(x)=\rho(x),\phi_{*}c(x)=c(x),\phi_{*}g(x)=g(x),\phi_{*}\nu(x)=\nu% (x)italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ ( italic_x ) = italic_ρ ( italic_x ) , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c ( italic_x ) = italic_c ( italic_x ) , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ( italic_x ) = italic_g ( italic_x ) , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ν ( italic_x ) = italic_ν ( italic_x ) and ϕu(x)=u(x)subscriptitalic-ϕ𝑢𝑥𝑢𝑥\phi_{*}u(x)=u(x)italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_u ( italic_x ) = italic_u ( italic_x ) for x𝑥xitalic_x in that neighborhood. Therefore the triplet (ϕρ,ϕc,ϕg)subscriptitalic-ϕ𝜌subscriptitalic-ϕ𝑐subscriptitalic-ϕ𝑔(\phi_{*}\rho,\phi_{*}c,\phi_{*}g)( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) gives rise to the same DN map as the triplet (ρ,c,g)𝜌𝑐𝑔(\rho,c,g)( italic_ρ , italic_c , italic_g ):

Λ(ϕρ,ϕc,ϕg)=Λ(ρ,c,g).subscriptΛsubscriptitalic-ϕ𝜌subscriptitalic-ϕ𝑐subscriptitalic-ϕ𝑔subscriptΛ𝜌𝑐𝑔\Lambda_{\left(\phi_{*}\rho,\phi_{*}c,\phi_{*}g\right)}=\Lambda_{(\rho,c,g)}.roman_Λ start_POSTSUBSCRIPT ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT . (26)

Additionally it was shown as a consequence of Lemma 2.7 in (18) that there is the following conformal freedom for the elastic wave operator P(ρ,c,g)subscript𝑃𝜌𝑐𝑔P_{(\rho,c,g)}italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT for any non-vanishing function μ𝜇\muitalic_μ:

P(μn2+nρ,μc,μ22+ng)=P(ρ,c,g).subscript𝑃superscript𝜇𝑛2𝑛𝜌𝜇𝑐superscript𝜇22𝑛𝑔subscript𝑃𝜌𝑐𝑔P_{\left(\mu^{\frac{n}{2+n}}\rho,\mu c,\mu^{-\frac{2}{2+n}}g\right)}=P_{(\rho,% c,g)}.italic_P start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ , italic_μ italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g ) end_POSTSUBSCRIPT = italic_P start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT . (27)

We note that

Λ(μn2+nρ,μc,μ22+ng)=μn2+nΛ(ρ,c,g).subscriptΛsuperscript𝜇𝑛2𝑛𝜌𝜇𝑐superscript𝜇22𝑛𝑔superscript𝜇𝑛2𝑛subscriptΛ𝜌𝑐𝑔\Lambda_{\left(\mu^{\frac{n}{2+n}}\rho,\mu c,\mu^{-\frac{2}{2+n}}g\right)}=\mu% ^{\frac{n}{2+n}}\Lambda_{(\rho,c,g)}.roman_Λ start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ , italic_μ italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g ) end_POSTSUBSCRIPT = italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT .

This implies for any function μ0𝜇0\mu\neq 0italic_μ ≠ 0 that satisfies μ|M=1evaluated-at𝜇𝑀1\mu|_{\partial M}=1italic_μ | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = 1:

Λ(μn2+nρ,μc,μ22+ng)=Λ(ρ,c,g).subscriptΛsuperscript𝜇𝑛2𝑛𝜌𝜇𝑐superscript𝜇22𝑛𝑔subscriptΛ𝜌𝑐𝑔\Lambda_{\left(\mu^{\frac{n}{2+n}}\rho,\mu c,\mu^{-\frac{2}{2+n}}g\right)}=% \Lambda_{(\rho,c,g)}.roman_Λ start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ρ , italic_μ italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_g ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT . (28)

Combining (26) and (28) yields

Λ(μn2+nϕρ,μϕc,μ22+nϕg)=Λ(ρ,c,g)subscriptΛsuperscript𝜇𝑛2𝑛subscriptitalic-ϕ𝜌𝜇subscriptitalic-ϕ𝑐superscript𝜇22𝑛subscriptitalic-ϕ𝑔subscriptΛ𝜌𝑐𝑔\Lambda_{\left(\mu^{\frac{n}{2+n}}\,\phi_{*}\rho,\mu\,\phi_{*}c,\mu^{-\frac{2}% {2+n}}\,\phi_{*}g\right)}=\Lambda_{(\rho,c,g)}roman_Λ start_POSTSUBSCRIPT ( italic_μ start_POSTSUPERSCRIPT divide start_ARG italic_n end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_ρ , italic_μ italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_c , italic_μ start_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 2 + italic_n end_ARG end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_ρ , italic_c , italic_g ) end_POSTSUBSCRIPT

as claimed. ∎

Appendix A Relation to Euclidean, Riemannian, and Lorentzian Calderón problems

Coordinate gauge freedom and conformal freedoms occur in some other inverse boundary value problems, and we discuss some examples to provide a point of comparison to our results.

The citations given in this appendix are only examples of the vast literature on these problems. For more details and background, we refer the reader to [30years, Belishev_2007, Lassas18, KatchalovKurylevLassas01] and references therein.

A.1. DN maps for Euclidean, Riemannian, and Lorentzian Calderón problems

The inverse problems we introduced in section 1.2 are closely related to the Calderón problem concerned with the conductivity equation describing electrostatics:

{Lγu(x):=(γ(x)u(x))=0in Ωu=fon Ω,casesassignsubscript𝐿𝛾𝑢𝑥𝛾𝑥𝑢𝑥0in Ω𝑢𝑓on Ω\begin{cases}L_{\gamma}u(x):=\nabla\cdot(\gamma(x)\nabla u(x))=0&\text{in }% \Omega\\ u=f&\text{on }\partial\Omega,\end{cases}{ start_ROW start_CELL italic_L start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_u ( italic_x ) := ∇ ⋅ ( italic_γ ( italic_x ) ∇ italic_u ( italic_x ) ) = 0 end_CELL start_CELL in roman_Ω end_CELL end_ROW start_ROW start_CELL italic_u = italic_f end_CELL start_CELL on ∂ roman_Ω , end_CELL end_ROW

where γ𝛾\gammaitalic_γ denotes an isotropic or anisotropic conductivity. The DN map corresponding to Lγsubscript𝐿𝛾L_{\gamma}italic_L start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT is defined by

Λγf=νiγijxju|Ω.subscriptΛ𝛾𝑓evaluated-atsubscript𝜈𝑖subscript𝛾𝑖𝑗subscriptsubscript𝑥𝑗𝑢Ω\Lambda_{\gamma}f=\nu_{i}\,\gamma_{ij}\,\partial_{x_{j}}u\Bigg{|}_{\partial% \Omega}.roman_Λ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT italic_f = italic_ν start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_u | start_POSTSUBSCRIPT ∂ roman_Ω end_POSTSUBSCRIPT . (29)

The Calderón problem is then concerned with recovering γ𝛾\gammaitalic_γ from ΛγsubscriptΛ𝛾\Lambda_{\gamma}roman_Λ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT. For an anisotropic conductivity γ𝛾\gammaitalic_γ it turns out that this inverse problem is closely related to a corresponding geometric inverse problem. Consider the Dirichlet problem associated to the Laplace–Beltrami operator on (Ω,g)Ω𝑔(\Omega,g)( roman_Ω , italic_g )

{Δgu=0in Ωu=fon Ω,casessubscriptΔ𝑔𝑢0in Ω𝑢𝑓on Ω\begin{cases}\Delta_{g}u=0&\text{in }\Omega\\ u=f&\text{on }\partial\Omega,\end{cases}{ start_ROW start_CELL roman_Δ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_u = 0 end_CELL start_CELL in roman_Ω end_CELL end_ROW start_ROW start_CELL italic_u = italic_f end_CELL start_CELL on ∂ roman_Ω , end_CELL end_ROW

where (cf. Remark 2.5 for an alternative formula)

Δgu=1det(g)xi(det(g)gijxju)subscriptΔ𝑔𝑢1𝑔subscriptsuperscript𝑥𝑖𝑔superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑢\Delta_{g}u=\frac{1}{\sqrt{\det(g)}}\partial_{x^{i}}\left(\sqrt{\det(g)}g^{ij}% \partial_{x^{j}}u\right)roman_Δ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_u = divide start_ARG 1 end_ARG start_ARG square-root start_ARG roman_det ( italic_g ) end_ARG end_ARG ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( square-root start_ARG roman_det ( italic_g ) end_ARG italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u ) (30)

and with corresponding DN map defined by

Λgf=νigijxju|Ω.subscriptΛ𝑔𝑓evaluated-atsubscript𝜈𝑖superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑢Ω\Lambda_{g}f=\nu_{i}\,g^{ij}\,\partial_{x^{j}}u\Bigg{|}_{\partial\Omega}.roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_f = italic_ν start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u | start_POSTSUBSCRIPT ∂ roman_Ω end_POSTSUBSCRIPT .

Then the geometric Calderón problem asks to recover g𝑔gitalic_g from ΛgsubscriptΛ𝑔\Lambda_{g}roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT. One can also consider the following geometric version of the conductivity equation on a Riemannian manifold (M,g)𝑀𝑔(M,g)( italic_M , italic_g ):

{divg(βgu)=0in Mu=fon M,casessubscriptdiv𝑔𝛽subscript𝑔𝑢0in 𝑀𝑢𝑓on 𝑀\begin{cases}\text{div}_{g}(\beta\nabla_{g}u)=0&\text{in }M\\ u=f&\text{on }\partial M,\end{cases}{ start_ROW start_CELL div start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_β ∇ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_u ) = 0 end_CELL start_CELL in italic_M end_CELL end_ROW start_ROW start_CELL italic_u = italic_f end_CELL start_CELL on ∂ italic_M , end_CELL end_ROW

where β𝛽\betaitalic_β is an isotropic conductivity. In this setting the DN map is defined as

Λ(β,g)f=νiβgijxju|M.subscriptΛ𝛽𝑔𝑓evaluated-atsubscript𝜈𝑖𝛽superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑢𝑀\Lambda_{(\beta,g)}f=\nu_{i}\,\beta\,g^{ij}\,\partial_{x^{j}}u\Bigg{|}_{% \partial M}.roman_Λ start_POSTSUBSCRIPT ( italic_β , italic_g ) end_POSTSUBSCRIPT italic_f = italic_ν start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_β italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT . (31)

This gives rise to the Calderón problem on a Riemannian manifold that asks to recover β𝛽\betaitalic_β and g𝑔gitalic_g from Λ(β,g)subscriptΛ𝛽𝑔\Lambda_{(\beta,g)}roman_Λ start_POSTSUBSCRIPT ( italic_β , italic_g ) end_POSTSUBSCRIPT.

There is also a similarity to the Lorentzian Calderón problem in that there is a conformal gauge freedom in all dimensions. An example is the non-linear wave equation

{gu(x)+a(x)u(x)4=0,on M,u(x)=f(x),on M,u(t,x)=0,t<0,casessubscript𝑔𝑢𝑥𝑎𝑥𝑢superscript𝑥40on 𝑀𝑢𝑥𝑓𝑥on 𝑀𝑢𝑡superscript𝑥0𝑡0\begin{cases}\square_{g}u(x)+a(x)u(x)^{4}=0,&\text{on }M,\\ u(x)=f(x),&\text{on }\partial M,\\ u(t,x^{\prime})=0,&t<0,\end{cases}{ start_ROW start_CELL □ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_u ( italic_x ) + italic_a ( italic_x ) italic_u ( italic_x ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT = 0 , end_CELL start_CELL on italic_M , end_CELL end_ROW start_ROW start_CELL italic_u ( italic_x ) = italic_f ( italic_x ) , end_CELL start_CELL on ∂ italic_M , end_CELL end_ROW start_ROW start_CELL italic_u ( italic_t , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = 0 , end_CELL start_CELL italic_t < 0 , end_CELL end_ROW

where gsubscript𝑔\square_{g}□ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT is the d’Alembert operator of g𝑔gitalic_g (the Laplace–Beltrami operator on a Lorentzian manifold) and M=×N𝑀𝑁M=\mathbb{R}\times Nitalic_M = blackboard_R × italic_N. The corresponding DN map is defined as

Λ(g,a)f=νjxju|(0,T)×NsubscriptΛ𝑔𝑎𝑓evaluated-atsuperscript𝜈𝑗subscriptsuperscript𝑥𝑗𝑢0𝑇𝑁\Lambda_{(g,a)}f=\nu^{j}\partial_{x^{j}}u|_{(0,T)\times\partial N}roman_Λ start_POSTSUBSCRIPT ( italic_g , italic_a ) end_POSTSUBSCRIPT italic_f = italic_ν start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u | start_POSTSUBSCRIPT ( 0 , italic_T ) × ∂ italic_N end_POSTSUBSCRIPT (32)

and the Lorentzian Calderón problem asks to recover the metric g𝑔gitalic_g and the function a𝑎aitalic_a from Λ(g,a)subscriptΛ𝑔𝑎\Lambda_{(g,a)}roman_Λ start_POSTSUBSCRIPT ( italic_g , italic_a ) end_POSTSUBSCRIPT.

A.2. Gauge freedom in the Euclidean and Riemannian Calderón problems

It was observed by L. Tartar that the map ΛγsubscriptΛ𝛾\Lambda_{\gamma}roman_Λ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT defined in (29) does not determine γ𝛾\gammaitalic_γ uniquely (see [KohnVog84] for an account). This is due to the fact that any Csuperscript𝐶C^{\infty}italic_C start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT diffeomorphism ϕ:Ω¯Ω¯:italic-ϕ¯Ω¯Ω\phi\colon\overline{\Omega}\to\overline{\Omega}italic_ϕ : over¯ start_ARG roman_Ω end_ARG → over¯ start_ARG roman_Ω end_ARG with ϕ|Ω=idevaluated-atitalic-ϕΩid\phi|_{\partial\Omega}=\operatorname{id}italic_ϕ | start_POSTSUBSCRIPT ∂ roman_Ω end_POSTSUBSCRIPT = roman_id and with the Euclidean pullback of the conductivity

ϕγ=1det(Dϕ)DϕγDϕt,superscriptitalic-ϕ𝛾1𝐷italic-ϕ𝐷italic-ϕ𝛾𝐷superscriptitalic-ϕ𝑡\phi^{\square}\gamma=\frac{1}{\det(D\phi)}D\phi\,\gamma\,D\phi^{t},italic_ϕ start_POSTSUPERSCRIPT □ end_POSTSUPERSCRIPT italic_γ = divide start_ARG 1 end_ARG start_ARG roman_det ( italic_D italic_ϕ ) end_ARG italic_D italic_ϕ italic_γ italic_D italic_ϕ start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT , (33)

produces the same DN map as γ𝛾\gammaitalic_γ:

Λϕγ=Λγ.subscriptΛsuperscriptitalic-ϕ𝛾subscriptΛ𝛾\Lambda_{\phi^{\square}\gamma}=\Lambda_{\gamma}.roman_Λ start_POSTSUBSCRIPT italic_ϕ start_POSTSUPERSCRIPT □ end_POSTSUPERSCRIPT italic_γ end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT . (34)

Notice that the pullback defined in (33) and the pushforward defined in (13) are not the usual pullbacks and pushforwards of tensor fields in differential geometry.

Similarly to the Euclidean case ΛgsubscriptΛ𝑔\Lambda_{g}roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT as defined in (A.1) does not determine g𝑔gitalic_g uniquely due to the coordinate gauge freedom

Λϕg=Λg,subscriptΛsubscriptitalic-ϕ𝑔subscriptΛ𝑔\Lambda_{\phi_{*}g}=\Lambda_{g},roman_Λ start_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , (35)

where ϕsubscriptitalic-ϕ\phi_{*}italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT denotes the classical pushforward of g𝑔gitalic_g by ϕitalic-ϕ\phiitalic_ϕ. It was shown in [LasUhl01] that this is the only gauge freedom for real analytic metrics in dimension n3𝑛3n\geq 3italic_n ≥ 3.

In dimension n=2𝑛2n=2italic_n = 2 the Laplace–Beltrami operator is conformally invariant which gives rise to an additional gauge for the inverse problem:

Λα(ϕg)=Λg.subscriptΛ𝛼subscriptitalic-ϕ𝑔subscriptΛ𝑔\Lambda_{\alpha(\phi_{*}g)}=\Lambda_{g}.roman_Λ start_POSTSUBSCRIPT italic_α ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT .

This is proven by [LasUhl01] to be the only obstruction for unique identifiability of g𝑔gitalic_g.

The inverse problems for the DN maps of (29) and (A.1) are equivalent when n3𝑛3n\geq 3italic_n ≥ 3. Namely, if we set

g=(detγ)1n2γ1or equivalentlyγ=(detg)12g1,formulae-sequence𝑔superscript𝛾1𝑛2superscript𝛾1or equivalently𝛾superscript𝑔12superscript𝑔1\begin{split}g&=(\det\gamma)^{\frac{1}{n-2}}\gamma^{-1}\quad\text{or % equivalently}\\ \gamma&=(\det\,g)^{\frac{1}{2}}g^{-1},\end{split}start_ROW start_CELL italic_g end_CELL start_CELL = ( roman_det italic_γ ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_n - 2 end_ARG end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT or equivalently end_CELL end_ROW start_ROW start_CELL italic_γ end_CELL start_CELL = ( roman_det italic_g ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT italic_g start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , end_CELL end_ROW

then (see [LeeUhlmann89])

Λg=Λγ,subscriptΛ𝑔subscriptΛ𝛾\Lambda_{g}=\Lambda_{\gamma},roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT ,

if we express ΛgsubscriptΛ𝑔\Lambda_{g}roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT with respect to the Euclidean surface measure:

Λgf=detgνigijxju|Ω.subscriptΛ𝑔𝑓evaluated-at𝑔subscript𝜈𝑖superscript𝑔𝑖𝑗subscriptsuperscript𝑥𝑗𝑢Ω\Lambda_{g}f=\sqrt{\det g}\,\nu_{i}\,g^{ij}\,\partial_{x^{j}}u\Bigg{|}_{% \partial\Omega}.roman_Λ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_f = square-root start_ARG roman_det italic_g end_ARG italic_ν start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_g start_POSTSUPERSCRIPT italic_i italic_j end_POSTSUPERSCRIPT ∂ start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u | start_POSTSUBSCRIPT ∂ roman_Ω end_POSTSUBSCRIPT .

A.3. Gauge freedom in the Calderón problem on manifolds

It was shown in [SunUhlmann03] that the conformal freedom of the Laplace–Beltrami operator in two dimensions yields the following conformal freedom for the DN operator defined in (31) in two dimensions for ϕ:MM:italic-ϕ𝑀𝑀\phi\colon M\to Mitalic_ϕ : italic_M → italic_M with ϕ|M=idevaluated-atitalic-ϕ𝑀id\phi|_{\partial M}=\operatorname{id}italic_ϕ | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = roman_id:

Λ(α(ϕg),ϕβ)=Λ(g,β),subscriptΛ𝛼subscriptitalic-ϕ𝑔subscriptitalic-ϕ𝛽subscriptΛ𝑔𝛽\Lambda_{(\alpha(\phi_{*}g),\phi_{*}\beta)}=\Lambda_{(g,\beta)},roman_Λ start_POSTSUBSCRIPT ( italic_α ( italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_g ) , italic_ϕ start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT italic_β ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_g , italic_β ) end_POSTSUBSCRIPT ,

for any positive scalar function α𝛼\alphaitalic_α.

A.4. Gauge freedom in the Lorentzian Calderón problem

Similarly to the previous Calderón problems there is a gauge freedom by a diffeomorphism for the DN map defined in (32). It was shown in [HintzUhlmannZhai22] that for ϕ:MM:italic-ϕ𝑀𝑀\phi\colon M\to Mitalic_ϕ : italic_M → italic_M with ϕ|M=idevaluated-atitalic-ϕ𝑀id\phi|_{\partial M}=\operatorname{id}italic_ϕ | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = roman_id we have

Λ(α(ϕg),ϕa)=Λ(g,a)subscriptΛ𝛼superscriptitalic-ϕ𝑔superscriptitalic-ϕ𝑎subscriptΛ𝑔𝑎\Lambda_{(\alpha(\phi^{*}g),\phi^{*}a)}=\Lambda_{(g,a)}roman_Λ start_POSTSUBSCRIPT ( italic_α ( italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_g ) , italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_a ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_g , italic_a ) end_POSTSUBSCRIPT

and that there is the conformal freedom

Λ(e2βg,eβa)=Λ(g,a)subscriptΛsuperscript𝑒2𝛽𝑔superscript𝑒𝛽𝑎subscriptΛ𝑔𝑎\Lambda_{(e^{-2\beta}g,e^{-\beta}a)}=\Lambda_{(g,a)}roman_Λ start_POSTSUBSCRIPT ( italic_e start_POSTSUPERSCRIPT - 2 italic_β end_POSTSUPERSCRIPT italic_g , italic_e start_POSTSUPERSCRIPT - italic_β end_POSTSUPERSCRIPT italic_a ) end_POSTSUBSCRIPT = roman_Λ start_POSTSUBSCRIPT ( italic_g , italic_a ) end_POSTSUBSCRIPT

for a smooth function β𝛽\betaitalic_β on M𝑀Mitalic_M such that

β|M=0,νβ|M=0,geβ=0.formulae-sequenceevaluated-at𝛽𝑀0formulae-sequenceevaluated-atsubscript𝜈𝛽𝑀0subscript𝑔superscript𝑒𝛽0\beta|_{\partial M}=0,\quad\partial_{\nu}\beta|_{\partial M}=0,\quad\square_{g% }e^{-\beta}=0.italic_β | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = 0 , ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_β | start_POSTSUBSCRIPT ∂ italic_M end_POSTSUBSCRIPT = 0 , □ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_β end_POSTSUPERSCRIPT = 0 .
\printbibliography