1. Introduction
We begin by introducing the geometric framework of our problem. We fix once and for all a natural number
\[ n \in \mathbb{N} {\setminus} \{0,1\} \]
that will be the dimension of the Euclidean space 
$\mathbb {R}^n$ we are going to work in. We also fix a parameter
\[ \alpha \in ]0,1[, \]
which we use to define the regularity of our sets and functions. In order to introduce the domains where our problem is defined, we take two sets 
$\Omega ^o$ and 
$\Omega ^i$ that satisfy the following conditions:
\begin{align}
&\Omega^o, \Omega^i\ \mbox{are bounded open connected
subsets of } \mathbb{R}^n \mbox{ of class } C^{1,\alpha},
\nonumber\\ &\mbox{with exteriors } \mathbb{R}^n{\setminus}
\overline{\Omega^o} \mbox{ and } \mathbb{R}^n{\setminus}
\overline{\Omega^i} \mbox{ connected and }
\overline{\Omega^i}\subset \Omega^o
\end{align}
(see figure 1). Here the superscript ‘
$o$’ stands for ‘outer domain’ whereas the superscript ‘
$i$’ stands for ‘inner domain’. We first want to introduce a transmission problem in the pair of domains consisting of 
$\Omega ^o {\setminus} \overline {\Omega ^i}$ and 
$\Omega ^i$. Therefore, to define the boundary conditions, we fix three functions
\begin{equation} F_1 \in C^0(\partial \Omega ^i \times \mathbb{R} \times \mathbb{R}),\quad F_2 \in C^0(\partial \Omega^i \times \mathbb{R} \times \mathbb{R}),\quad f^o \in C^{0,\alpha}(\partial\Omega^o). \end{equation}
Figure 1. The domains 
$\Omega ^o$ and 
$\Omega ^i$ (
$n=2$).
The functions 
$F_1$ and 
$F_2$ determine the transmission conditions on the inner boundary 
$\partial \Omega ^i$. Instead, 
$f^o$ plays the role of the Neumann datum on the outer boundary 
$\partial \Omega ^o$. We consider the following nonlinear transmission boundary value problem for a pair 
$(u^o,u^i) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i) \times C^{1,\alpha }(\overline {\Omega ^i})$:
\begin{equation} \begin{cases} \Delta u^o = 0 & \mbox{in } \Omega^o {\setminus} \overline{\Omega^i}, \\ \Delta u^i = 0 & \mbox{in } \Omega^i, \\ \nu_{\Omega^o}(x) \cdot \nabla u^o(x)=f^o(x) & \forall x \in \partial \Omega^o, \\ \nu_{\Omega^i}(x) \cdot \nabla u^o (x) = F_1(x,u^o(x),u^i(x)) & \forall x \in \partial \Omega^i, \\ \nu_{\Omega^i}(x) \cdot \nabla u^i (x) = F_2(x,u^o(x),u^i(x)) & \forall x \in \partial \Omega^i, \end{cases} \end{equation}
where 
$\nu _{\Omega ^o}$ and 
$\nu _{\Omega ^i}$ denote the outward unit normal vector field to 
$\partial \Omega ^o$ and to 
$\partial \Omega ^i$, respectively. We note that, a priori, it is not clear why problem (1.3) should admit a classical solution. As a first result, we prove that under suitable conditions on 
$F_1$ and 
$F_2$, problem (1.3) has at least one solution 
$(u^o,u^i) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i) \times C^{1,\alpha }(\overline {\Omega ^i})$. We notice that a problem similar to (1.3) has been studied in Dalla Riva and Mishuris [Reference Dalla Riva and Mishuris7]. More precisely, in [Reference Dalla Riva and Mishuris7] the authors consider a nonlinear transmission problem with a Dirichlet boundary condition on 
$\partial \Omega ^o$ and a jump type condition for the normal derivative across the interface 
$\partial \Omega ^i$. Then we study the existence and the analytic dependence of the solutions of the transmission problem (1.3) upon domain perturbation of the inclusion, i.e. of the inner set 
$\Omega ^i$. Hence, we introduce a ‘perturbed’ version of problem (1.3): we fix the external domain 
$\Omega ^o$ and we assume that the boundary of the internal domain is of the form 
$\phi (\partial \Omega ^i)$, where 
$\phi$ is a diffeomorphism of 
$\partial \Omega ^i$ into a subset of 
$\mathbb {R}^n$ that belongs to the class
\begin{equation} \begin{split} \mathcal{A}_{\partial\Omega^i} \equiv \Big\{ \phi & \in C^{1,\alpha}(\partial\Omega^i, \mathbb{R}^n): \, \phi \text{ is injective and}\\ & \text{the differential} \, d\phi(y) \text{ is injective for all } y \in \partial\Omega^i \Big\}. \end{split} \end{equation}
Clearly, the identity function of 
$\partial \Omega ^i$ belongs to the class 
$\mathcal {A}_{\partial \Omega ^i}$, and, for convenience, we set
\begin{equation} \phi_0 \equiv \text{id}_{\partial\Omega^i}. \end{equation}
Then by the Jordan Leray Separation Theorem (cf., e.g., Deimling [Reference Deimling9, theorem 5.2] and [Reference Dalla Riva, Lanza de Cristoforis and Musolino6, § A.4]), 
$\mathbb {R}^n {\setminus} \phi (\partial \Omega ^i)$ has exactly two open connected components for all 
$\phi \in \mathcal {A}_{\partial \Omega ^i}$, and we define 
$\Omega ^i[\phi ]$ to be the unique bounded open connected component of 
$\mathbb {R}^n {\setminus} \phi (\partial \Omega ^i)$. We set
\[ \mathcal{A}^{\Omega^o}_{\partial\Omega^i} \equiv \left\{ \phi \in \mathcal{A}_{\partial\Omega^i} : \overline{\Omega^i[\phi]} \subset \Omega^o \right\}. \]
By assumption (1), 
$\phi _0 \in \mathcal {A}^{\Omega ^o}_{\partial \Omega ^i}$. Now let 
$\phi \in \mathcal {A}^{\Omega ^o}_{\partial \Omega ^i}$. We wish to consider the following nonlinear transmission boundary value problem for a pair of functions 
$(u^o,u^i) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i[\phi ]) \times C^{1,\alpha }(\overline {\Omega ^i[\phi ]})$:
\begin{equation} \begin{cases} \Delta u^o = 0 & \mbox{in } \Omega^o {\setminus} \overline{\Omega^i[\phi]}, \\ \Delta u^i = 0 & \mbox{in } \Omega^i[\phi], \\ \nu_{\Omega^o}(x) \cdot \nabla u^o(x)=f^o(x) & \forall x \in \partial \Omega^o, \\ \nu_{\Omega^i[\phi]}(x) \cdot \nabla u^o (x) = F_1(\phi^{({-}1)}(x),u^o(x),u^i(x)) & \forall x \in \phi(\partial\Omega^i), \\ \nu_{\Omega^i[\phi]}(x) \cdot \nabla u^i (x) = F_2(\phi^{({-}1)}(x),u^o(x),u^i(x)) & \forall x \in \phi(\partial\Omega^i), \end{cases} \end{equation}
where 
$\nu _{\Omega ^i[\phi ]}$ denotes the outward unit normal vector field to 
$\Omega ^i[\phi ]$. We prove that, under suitable conditions, problem (1.6) admits a family of solutions 
$\{(u^o_\phi ,u^i_\phi )\}_{\phi \in Q_0}$, where 
$Q_0$ is a neighbourhood of 
$\phi _0$ in 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i}$ and 
$(u^o_\phi ,u^i_\phi ) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i[\phi ]) \times C^{1,\alpha }(\overline {\Omega ^i[\phi ]})$ for every 
$\phi \in Q_0$. In literature, the existence of solutions of nonlinear boundary value problems has been largely investigated by means of variational techniques (see, e.g., the monographs of Nečas [Reference Nečas26] and of Roubíček [Reference Roubíček28] and the references therein). Moreover, potential theoretic techniques have been widely exploited to study nonlinear boundary value problems with transmission conditions by Bergeret al. [Reference Berger, Warnecke and Wendland3], by Costabel and Stephan [Reference Costabel and Stephan5], by Gatica and Hsiao [Reference Gatica and Hsiao11], and by Barrenechea and Gatica [Reference Barrenechea and Gatica2]. Boundary integral methods have been applied also by Mityushev and Rogosin for the analysis of transmission problems in the plane (cf. [Reference Mityushev and Rogosin24, chapter 5]). Several authors have investigated the dependence upon domain perturbation of the solutions to boundary value problems and it is impossible to provide a complete list of contributions. Here we mention, for example, Henrot and Pierre [Reference Henrot and Pierre13], Henry [Reference Henry14], Keldysh [Reference Keldysh15], Novotny and Sokołowski [Reference Novotny and Sokołowski27], and Sokołowski and Zolésio [Reference Sokolowski and Zolésio31]. Most of the contributions on this topic deal with first- or second-order shape derivability of functionals associated to the solutions of linear boundary value problems. In the present paper, instead, we are interested into higher order regularity properties (namely real analiticity) of the solutions of a nonlinear problem. To do so, we choose to adopt the Functional Analytic Approach, which has revealed to be a powerful tool to analyse perturbed linear and nonlinear boundary value problems. This method has been first applied to investigate regular and singular domain perturbation problems for elliptic equations and systems with the aim of proving real analytic dependence upon the perturbation parameter (cf. Lanza de Cristoforis [Reference Lanza de Cristoforis16, Reference Lanza de Cristoforis18, Reference Lanza de Cristoforis19]). An application to the study of the behaviour of the effective conductivity of a periodic two-phase composite upon perturbations of the inclusion can be found in Luzzini and Musolino [Reference Luzzini and Musolino22]. The key point of the strategy of the method is the transformation of the perturbed boundary value problem into an equivalent functional equation that can be studied by the Implicit Function Theorem. Typically, such a transformation is achieved by exploiting classical results of potential theory, for example, integral representation of harmonic functions in terms of layer potentials. Nonlinear transmission problems in perturbed domains have been studied by Lanza de Cristoforis in [Reference Lanza de Cristoforis19] and by the authors of the present paper in [Reference Dalla Riva, Molinarolo and Musolino8, Reference Molinarolo25], where they have investigated the behaviour of the solution of a nonlinear transmission problem for the Laplace equation in a domain with a small inclusion shrinking to a point.
The paper is organized as follows. In § 2 we define some of the symbols used later on. In § 3 we introduce some classical results of potential theory that we need. Section 4 is devoted to the study of problem (1.3). We first prove a representation result for harmonic functions in 
$\overline {\Omega ^o}{\setminus} \Omega$ and 
$\overline {\Omega }$ (where 
$\Omega$ is an open bounded connected subset of class 
$C^{1,\alpha }$ contained in 
$\Omega ^o$) in terms of single-layer potentials with appropriate densities and constant functions (cf. lemma 4.1). Then we prove an uniqueness result in 
$C^{1,\alpha }(\Omega ^o{\setminus} \overline {\Omega ^i}) \times C^{1,\alpha }(\overline {\Omega ^i})$ for an homogeneous linear transmission problem in the pair of domains 
$\Omega ^o{\setminus} \overline {\Omega ^i}$ and 
$\Omega ^i$ and we analyse an auxiliary boundary operator arising from the integral formulation of that problem (cf. lemma 4.2 and proposition 4.3). In proposition 4.4 we provide a formulation of problem (1.3) in terms of integral equations. The obtained integral system is solved by means of a fixed-point theorem, namely the Leray-Schauder Theorem (cf. propositions 4.5 and 4.7). Finally, under suitable conditions on the functions 
$F_1$ and 
$F_2$, we obtain an existence results in 
$C^{1,\alpha }(\Omega ^o{\setminus} \overline {\Omega ^i}) \times C^{1,\alpha }(\overline {\Omega ^i})$ for problem (1.3) (cf. proposition 4.8). Section 5 is devoted to the study of problem (1.6). We provide a formulation of problem (1.6) in terms of integral equations depending on the diffeomorphism 
$\phi$ which we rewrite into an equation of the type 
$M[\phi ,\mu ] = 0$ for an auxiliary map 
$M: \mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times X \to Y$ (with 
$X$ and 
$Y$ suitable Banach spaces), where the variable 
$\mu$ is related to the densities of the integral representation of the solution (cf. proposition 5.1). Then, by analiticity results for the dependence of single- and double-layer potentials upon the perturbation of the support, we prove that 
$M$ is real analytic (cf. proposition 5.2) and the differential of 
$M$ with respect to the variable 
$\mu \in X$ is an isomorphism (cf. proposition 5.3). Hence, by the Implicit Function Theorem, we show the existence of a family of solutions 
$\{(u^o_\phi ,u^i_\phi )\}_{\phi \in Q_0}$ of (1.6) (cf. theorem 5.6) and we prove that it can be represented in terms of real analytic functions (cf. theorem 5.7).
2. Notation
We denote by 
$\mathbb {N}$ the set of natural numbers including 
$0$. We denote the norm of a real normed space 
$X$ by 
$\| \cdot \| _X$. We denote by 
$I_X$ the identity operator from 
$X$ to itself and we omit the subscript 
$X$ where no ambiguity can occur. If 
$X$ and 
$Y$ are normed spaces we consider on the product space 
$X \times Y$ the norm defined by 
$\| (x,y) \|_{X \times Y} \equiv \|x\|_X + \|y\|_Y$ for all 
$(x,y) \in X \times Y$, while we use the Euclidean norm for 
$\mathbb {R}^d$, 
$d\in \mathbb {N}{\setminus} \{0,1\}$. If 
$U$ is an open subset of 
$X$, and 
$F:U \to Y$ is a Fréchet-differentiable map in 
$U$, we denote the differential of 
$F$ by 
$\textrm {d}F$. The inverse function of an invertible function 
$f$ is denoted by 
$f^{(-1)}$, while the reciprocal of a non-zero scalar function 
$g$ or the inverse of an invertible matrix 
$A$ are denoted by 
$g^{-1}$ and 
$A^{-1}$ respectively. Let 
$\Omega \subseteq \mathbb {R}^n$. Then 
$\overline {\Omega }$ denotes the closure of 
$\Omega$ in 
$\mathbb {R}^n$, 
$\partial \Omega$ denotes the boundary of 
$\Omega$, and 
$\nu _\Omega$ denotes the outward unit normal to 
$\partial \Omega$. For 
$x \in \mathbb {R}^d$, 
$x_j$ denotes the 
$j$-th coordinate of 
$x$, 
$|x|$ denotes the Euclidean modulus of 
$x$ in 
$\mathbb {R}^d$. If 
$x \in \mathbb {R}^d$ and 
$r>0$, we denote by 
$B_d(x,r)$ the open ball of centre 
$x$ and radius 
$r$. Let 
$\Omega$ be an open subset of 
$\mathbb {R}^n$ and 
$m \in \mathbb {N} {\setminus} \{0\}$. The space of 
$m$ times continuously differentiable real-valued function on 
$\Omega$ is denoted by 
$C^m(\Omega ,\mathbb {R})$ or more simply by 
$C^m(\Omega )$. Let 
$r \in \mathbb {N} {\setminus} \{0\}$, 
$f \in (C^m(\Omega ))^r$. The 
$s$-th component of 
$f$ is denoted by 
$f_s$ and the gradient of 
$f_s$ is denoted by 
$\nabla f_s$. Let 
$\eta =(\eta _1, \dots ,\eta _n) \in \mathbb {N}^n$ and 
$|\eta |=\eta _1+ \dots +\eta _n$. Then 
$D^\eta f \equiv {\partial ^{|\eta |}f}/{\partial x^{\eta _1}_1, \dots , \partial x^{\eta _n}_n}$. We retain the standard notation for the space 
$C^{\infty }(\Omega )$ and its subspace 
$C^{\infty }_c(\Omega )$ of functions with compact support. The subspace of 
$C^m(\Omega )$ of those functions 
$f$ such that 
$f$ and its derivatives 
$D^\eta f$ of order 
$|\eta |\le m$ can be extended with continuity to 
$\overline {\Omega }$ is denoted 
$C^m(\overline {\Omega })$. We denote by 
$C^m_b(\overline {\Omega })$ the space of functions of 
$C^m(\overline {\Omega })$ such that 
$D^{\eta } f$ is bounded for 
$|\eta |\leq m$. Then the space 
$C^m_b(\overline {\Omega })$ equipped with the usual norm 
$\|f\|_{C^m_b(\overline {\Omega })} \equiv \sum _{|\eta |\leq m} \sup _{\overline {\Omega }} |D^\eta f|$ is well known to be a Banach space. Let 
$f \in C^0(\overline {\Omega })$. Then we define its Hölder constantas
\[ |f : \Omega|_\alpha\equiv \mbox{sup} \left\{\frac{|f(x)-f(y)|}{|x-y|^\alpha} : x,y \in \overline{\Omega}, x \neq y \right\}. \]
We define the subspace of 
$C^0(\overline {\Omega })$ of Hölder continuous functions with exponent 
$\alpha \in ]0,1[$ by 
$C^{0,\alpha }(\overline {\Omega }) \equiv \{f \in C^0(\overline {\Omega }) : \, |f : \Omega |_\alpha < \infty \}$. Similarly, the subspace of 
$C^m(\overline {\Omega })$ whose functions have 
$m$-th order derivatives that are Hölder continuous with exponent 
$\alpha \in ]0,1[$ is denoted 
$C^{m,\alpha }(\overline {\Omega })$. Then the space 
$C^{m,\alpha }_b(\overline {\Omega }) \equiv C^{m,\alpha }(\overline {\Omega }) \cap C^m_b(\overline {\Omega })$, equipped with its usual norm 
$\|f\|_{C^{m,\alpha }_b(\overline {\Omega })} \equiv \|f\|_{C^{m}_b(\overline {\Omega })} + \sum _{|\eta |=m}{|D^\eta f : \Omega |_\alpha }$, is a Banach space. If 
$\Omega$ is bounded, then 
$C^{m,\alpha }_b(\overline {\Omega }) = C^{m,\alpha }(\overline {\Omega })$, and we omit the subscript 
$b$. We denote by 
$C^{m,\alpha }_{\mathrm {loc}}(\mathbb {R}^n {\setminus} \Omega )$ the space of functions on 
$\mathbb {R}^n {\setminus} \Omega$ whose restriction to
$\overline {U}$belongs to 
$C^{m,\alpha }(\overline {U})$ for all open bounded subsets 
$U$ of 
$\mathbb {R}^n {\setminus} \Omega$. On 
$C^{m,\alpha }_{\mathrm {loc}}(\mathbb {R}^n {\setminus} \Omega )$ we consider the natural structure of Fréchet space. Finally if 
$\Omega$ is bounded, we set
\begin{align*}
C^{m,\alpha}_{\mathrm{h}}(\overline{\Omega}) &\equiv \{ u
\in C^{m,\alpha}(\overline{\Omega}) \cap C^2(\Omega):
\Delta u = 0 \text{ in } \Omega \},\\
C^{m,\alpha}_{\mathrm{h}}(\mathbb{R}^n {\setminus} \Omega)
&\equiv \{ u \in C^{m,\alpha}(\mathbb{R}^n {\setminus}
\Omega)\cap C^2(\mathbb{R}^n {\setminus} \overline{\Omega}):
\Delta u = 0 \text{ in } \mathbb{R}^n {\setminus}
\overline{\Omega}, \\ &\qquad |u(x)| = O(|x|^{2-n}) \text{
as } x \to +\infty \}.
\end{align*}
The condition 
$|u(x)| = O(|x|^{2-n})$ as 
$x \to +\infty$ in the above definition is equivalent for an harmonic function to the so-called harmonicity at infinity (see Folland [Reference Folland10, proposition (2.74), p. 112]). We say that a bounded open subset of 
$\mathbb {R}^n$ is of class 
$C^{m,\alpha }$ if it is a manifold with boundary imbedded in 
$\mathbb {R}^n$ of class 
$C^{m,\alpha }$. In particular if 
$\Omega$ is a 
$C^{1,\alpha }$ subset of 
$\mathbb {R}^n$, then 
$\partial \Omega$ is a 
$C^{1,\alpha }$ sub-manifold of 
$\mathbb {R}^n$ of co-dimension 
$1$. If 
$M$ is a 
$C^{m,\alpha }$ sub-manifold of 
$\mathbb {R}^n$ of dimension 
$d\ge 1$, we define the space 
$C^{m,\alpha }(M)$ by exploiting a finite local parametrization. We retain the standard definition of the Lebesgue spaces 
$L^p$, 
$p\ge 1$. If 
$\Omega$ is of class 
$C^{1,\alpha }$, we denote by 
$d\sigma$ the area element on 
$\partial \Omega$. If 
$Z$ is a subspace of 
$L^1(\partial \Omega )$, we set
\[ Z_0 \equiv \left\{ f \in Z : \int_{\partial\Omega} f \,{\rm d}\sigma = 0 \right\}. \]
Then we introduce a notation for superposition operators: if 
$H$ is a function from 
$\partial \Omega ^i \times {\mathbb{R}} \times \mathbb {R}$ to 
$\mathbb {R}$, then we denote by 
$\mathcal {N}_{H}$ the nonlinear nonautonomous superposition operator that take a pair 
$(h^1,h^2)$ of functions from 
$\partial \Omega ^i$ to 
$\mathbb {R}$ to the function 
$\mathcal {N}_{H}(h^1,h^2)$ defined by
\[ \mathcal{N}_{H}(h^1,h^2)(x) \equiv H(x,h^1(x),h^2(x)) \quad\forall x \in \partial\Omega^i. \]
Here the letter ‘
$\mathcal {N}$’ stands for ‘Nemytskii operator’. Finally, we have the following by Lanza de Cristoforis and Rossi [Reference Lanza de Cristoforis and Rossi21, lemma 3.3, proposition 3.13].
Lemma 2.1 Let 
$\Omega ^i$ be as in (1) and let 
$\mathcal {A}_{\partial \Omega ^i}$ be as in (1.4). Let 
$\phi \in \mathcal {A}_{\partial \Omega ^i}$. Then there exists a unique function 
$\tilde {\sigma }_n[\phi ] \in C^{0,\alpha }(\partial \Omega ^i)$ such that
\[ \int_{\phi(\partial\Omega^i)} f(y) \,{\rm d}\sigma_y = \int_{\partial\Omega^i} f(\phi(s)) \, \tilde{\sigma}_n[\phi](s) \,{\rm d}\sigma_s \quad\forall f \in L^1(\phi(\partial\Omega^i)). \]
Moreover the map from 
$\mathcal {A}_{\partial \Omega ^i}$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$\phi$ to 
$\tilde {\sigma }_n[\phi ]$ and the map from 
$\mathcal {A}_{\partial \Omega ^i}$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$\phi$ to 
$\nu _{\Omega ^i[\phi ]}(\phi (\cdot ))$ are real analytic.
3. Some preliminaries of potential theory
As we have mentioned, a key point of the Functional Analytic Approach is the reformulation of the boundary value problem in terms of an equivalent integral equation. To this aim, we exploit representation formulas for harmonic functions in terms of layer potentials. In this section, we collect some classical results of potential theory. We do not present proofs that can be found, for example, in Folland [Reference Folland10, chapter 3], in Gilbarg and Trudinger [Reference Gilbarg and Trudinger12, § 2].
Definition 3.1 We denote by 
$S_n$ the function from 
$\mathbb {R}^n {\setminus} \{0\}$ to 
$\mathbb {R}$ defined by
\[ S_n(x) \equiv
\begin{cases} \dfrac{1}{s_n} \log |x| & \forall x \in
\mathbb{R}^n {\setminus} \{0\} \quad \mbox{if } n=2 \\
\dfrac{1}{(2-n) s_n} |x|^{2-n} & \forall x \in \mathbb{R}^n
{\setminus} \{0\} \quad \mbox{if } n>2 \end{cases}
\]
where 
$s_n$ denotes the 
$(n-1)$-dimensional measure of 
$\partial B_n(0,1)$.
$S_n$ is well known to be a fundamental solution of the Laplace operator 
$\Delta =\sum _{j=1}^n\partial ^2_{x_j}$. We now assume that
\[\Omega\ \text{is an open bounded subset of } \mathbb{R}^n \text{ of class } C^{1,\alpha}. \]
In the following definition we introduce the single-layer potential, which we use to transform our problems into integral equations.
Definition 3.2 We denote by 
$v_{\Omega }[\mu ]$ the single-layer potential with density 
$\mu$, i.e. the function defined by
\[ v_{\Omega}[\mu](x) \equiv \int_{\partial \Omega}{S_n(x-y) \mu(y) \,{\rm d}\sigma_y} \quad \forall x \in \mathbb{R}^n , \forall \mu \in L^2(\partial\Omega). \]
It is well known that if 
$\mu \in C^{0,\alpha }(\partial \Omega )$, then 
$v_{\Omega }[\mu ] \in C^0(\mathbb {R}^n)$. We set
\[ v^+_{\Omega}[\mu]\equiv v_{\Omega}[\mu]_{| \overline{\Omega}}, \quad v^-_{\Omega}[\mu] \equiv v_{\Omega}[\mu]_{| \mathbb{R}^n {\setminus} \Omega}. \]
Then we define the boundary integral operators associated to the trace of the single-layer potential and its normal derivative.
Definition 3.3 We denote by 
$V_{\partial \Omega }$ the operator from 
$L^2(\partial \Omega )$ to itself that takes 
$\mu$ to the function 
$V_{\partial \Omega }[\mu ]$ defined in the trace sense by
\[ V_{\partial\Omega}[\mu]\equiv v_{\Omega}[\mu]_{|\partial\Omega}. \]
We denote by 
$W_{\partial \Omega }$ the integral operator from 
$L^2(\partial \Omega )$ to itself defined by
\[ W_{\partial\Omega}[\mu](x) \equiv{-} \int_{\partial \Omega} {\nu_\Omega(y) \cdot \nabla S_n(x-y) \mu(y) \,{\rm d}\sigma_y} \quad \text{for a.e. } x \in \partial \Omega, \forall \mu \in L^2(\partial\Omega). \]
We denote by 
$W^\ast _{\partial \Omega }$ the integral operator from 
$L^2(\partial \Omega )$ to itself which is the transpose of 
$W_{\partial \Omega }$ and that is defined by
\[ W^\ast_{\partial\Omega}[\mu](x) \equiv \int_{\partial \Omega} {\nu_\Omega(x) \cdot \nabla S_n(x-y) \mu(y) \,d\sigma_y} \quad \mbox{for a.e. } x \in \partial \Omega, \forall \mu \in L^2(\partial\Omega). \]
As it is well known, since 
$\Omega$ is of class 
$C^{1,\alpha }$, 
$W_{\partial \Omega }$ and 
$W^\ast _{\partial \Omega }$ are compact operators from 
$L^2(\partial \Omega )$ to itself (both display a weak singularity). In particular 
$( \pm \frac {1}{2} I + W^\ast _{\partial \Omega } )$ are Fredholm operators of index 
$0$ from 
$L^2(\partial \Omega )$ to itself. Moreover, one verifies that 
$W_{\partial \Omega }: C^{1,\alpha }(\partial \Omega ) \to C^{1,\alpha }(\partial \Omega )$ and 
$W^\ast _{\partial \Omega }: C^{0,\alpha }(\partial \Omega ) \to C^{0,\alpha }(\partial \Omega )$ are transpose to one another with respect to the duality of 
$C^{1,\alpha }(\partial \Omega ) \times C^{0,\alpha }(\partial \Omega )$ induced by the inner product of 
$L^2(\partial \Omega )$. We collect some well-known properties of the single-layer potential in the theorem below. In particular, we note that the operator of statement (iv) is an isomorphism both in the case of dimension 
$n=2$ and 
$n \geq 3$ (see, e.g., [Reference Dalla Riva, Lanza de Cristoforis and Musolino6, theorem 6.47]).
Theorem 3.4 Properties of the single-layer potential
The following statements hold.
(i) For all
$\mu \in L^2(\partial \Omega )$, the function $v_{\Omega }[\mu ]$ is harmonic in $\mathbb {R}^n{\setminus} \partial \Omega$. If $n\geq 3$ or if $n=2$ and $\int _{\partial \Omega }\mu \, \textrm {d}\sigma =0$ then $v_{\Omega }[\mu ]$ is also harmonic at infinity.(ii) If
$\mu \in C^{0,\alpha }(\partial \Omega )$, then $v^+_{\Omega }[\mu ] \in C^{1,\alpha }(\overline {\Omega })$ and the map from $C^{0,\alpha }(\partial \Omega )$ to $C^{1,\alpha }(\overline {\Omega })$ that takes $\mu$ to $v^+_{\Omega }[\mu ]$ is linear and continuous. Moreover, $v^-_{\Omega }[\mu ] \in C^{1,\alpha }_{\mathrm {loc}}(\mathbb {R}^n {\setminus} \Omega )$ and the map from $C^{0,\alpha }(\partial \Omega )$ to $C^{1,\alpha }_{\mathrm {loc}}(\mathbb {R}^n {\setminus} \Omega )$ that takes $\mu$ to $v^-_{\Omega }[\mu ]$ is linear and continuous.(iii) If
$\mu \in C^{0,\alpha }(\partial \Omega )$, then we have following jump relations
\[ \nu_\Omega(x) \cdot \nabla v^\pm_{\Omega}[\mu] (x) = \left({\mp} \frac{1}{2} I + W^\ast_{\partial\Omega} \right)[\mu](x) \quad \forall x \in \partial \Omega. \]
(iv) The map from
$C^{0,\alpha }(\partial \Omega )_0 \times \mathbb {R}$ to $C^{0,\alpha }(\partial \Omega )$ that takes a pair $(\mu ,\rho )$ to $V_{\partial \Omega }[\mu ] + \rho$ is an isomorphism.
Since 
$\Omega$ is of class 
$C^{1,\alpha }$, the following classical compactness result holds (cf. Schauder [Reference Schauder29, Reference Schauder30]).
Theorem 3.5 The map that takes 
$\mu$ to 
$W^\ast _{\partial \Omega }[\mu ]$ is compact from 
$C^{0,\alpha }(\partial \Omega )$ to itself.
Theorem 3.5 implies that 
$( \pm \frac {1}{2} I + W^\ast _{\partial \Omega } )$ are Fredholm operators of index 
$0$ from 
$C^{0,\alpha }(\partial \Omega )$ into itself. We now collect some regularity results for integral operators. We first introduce the following (see Folland [Reference Folland10, chapter 3 § B]).
Definition 3.6 Let 
$K$ be a measurable function from 
$\partial \Omega \times \partial \Omega$ to 
$\mathbb {R}$ and let 
$0 \leq \beta < n-1$. We say that 
$K$ is a continuous kernel of order 
$\beta$ if
\[ K(x,y) = k(x,y) |x-y|^{-\beta} \quad \forall(x,y)\in \partial\Omega \times \partial\Omega, \]
for some continuous function 
$k$ on 
$\partial \Omega \times \partial \Omega$.
If 
$K$ is a continuous kernel of order 
$\beta$, we denote by 
$\mathcal {K}_K$ the integral operator from 
$L^2(\partial \Omega )$ to itself defined by
\[ \mathcal{K}_K [\mu] (x) \equiv \int_{\partial\Omega} K(x,y) \mu (y) \,{\rm d}\sigma_y \quad \text{for a.e. } x \in \partial \Omega, \forall \mu \in L^2(\partial\Omega). \]
We observe that the functions 
$K_1(x,y) \equiv S_n(x-y)$ and 
$K_2(x,y) \equiv \nu _\Omega (y) \cdot \nabla S_n(x-y)$ of 
$(x,y) \in \partial \Omega \times \partial \Omega$, 
$x \neq y$, are continuous kernels of order 
$n-2$ (cf. Folland [Reference Folland10, proposition 3.17]). Clearly, we can extend the notion of integral operator with a continuous kernel to the vectorial case just applying the definition above component-wise. Then we present a vectorial version of a classical regularity result (see, for example, Folland [Reference Folland10, proposition 3.13]).
Theorem 3.7 Let 
$0\leq \beta < n-1$. Let 
$K_i^j$ with 
$i,j \in \{1,2\}$ be continuous kernels of order 
$\beta$. Let 
$\mathcal {K}=(\mathcal {K}_1,\mathcal {K}_2)$ be the operator from 
$(L^2(\partial \Omega ))^2$ to itself defined by
\[ \mathcal{K}_1 [\mu_1,\mu_2] = \mathcal{K}_{K_1^1}[\mu_1] + \mathcal{K}_{K_2^1}[\mu_2],\quad \mathcal{K}_2 [\mu_1,\mu_2] = \mathcal{K}_{K_1^2}[\mu_1] + \mathcal{K}_{K_2^2}[\mu_2], \]
for all 
$(\mu _1,\mu _2) \in (L^2(\partial \Omega ))^2$. If 
$(I + \mathcal {K})[\mu _1,\mu _2] \in (C^{0}(\partial \Omega ))^2$, then 
$(\mu _1,\mu _2) \in (C^{0}(\partial \Omega ))^2$.
Finally, we present in theorem 3.8 a regularity result that will be widely used in what follows. The proof exploits a standard argument on iterated kernels and can be found, e.g., in Dalla Riva and Mishuris [Reference Dalla Riva and Mishuris7, lemma 3.3].
Theorem 3.8 Let 
$\mu \in L^2(\partial \Omega )$. Let 
$\beta \in [0,\alpha ]$. If 
$( \frac {1}{2} I + W^\ast _{\partial \Omega } )[\mu ]$ or 
$( -\frac {1}{2} I + W^\ast _{\partial \Omega } )[\mu ]$ belongs to 
$C^{0,\beta }(\partial \Omega )$, then 
$\mu \in C^{0,\beta }(\partial \Omega )$.
4. Existence result for problem (1.3)
The aim of this section is to prove an existence result for problem (1.3). We start with the following representation result for harmonic functions in 
$\Omega ^o {\setminus} \overline {\Omega }$ and in 
$\Omega$ in terms of single-layer potentials plus constant functions. The set 
$\Omega$ in the lemma 4.1 will be later replaced by the set 
$\Omega ^i$ and by the perturbed set 
$\Omega ^i[\phi ]$.
Lemma 4.1 Let 
$\Omega$ be an open bounded connected subset of 
$\mathbb {R}^n$ of class 
$C^{1,\alpha }$, such that 
$\mathbb {R}^n{\setminus} \overline {\Omega }$ is connected and 
$\overline {\Omega }\subset \Omega ^o$. Then the map from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times C^{0,\alpha }(\partial \Omega )_0 \times \mathbb {R}^2$ to 
$C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega ^o} {\setminus} \Omega ) \times C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega })$ that takes a quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the pair of functions 
$(U^o_\Omega [\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i], U^i_\Omega [\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i])$ defined by
\begin{equation} \begin{aligned} & U^o_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] \equiv (v^+_{\Omega^o} [\mu^o] + v^-_{\Omega}[\mu^i] + \rho^o)_{| \overline{\Omega^o} {\setminus} \Omega} \\ & U^i_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] \equiv v^+_{\Omega}[\eta^i] + \rho^i \end{aligned} \end{equation}is bijective.
Proof. The map is well defined. Indeed, by the harmonicity and regularity properties of single-layer potentials (cf. theorem 3.4(i)–(ii)), we know that
\begin{align*} &\Delta U^o_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i]=0 \quad \text{on } \Omega^o {\setminus} \overline{\Omega}, \\ &\Delta U^i_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i]=0 \quad \text{on } \Omega, \\ &(U^o_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i],U^i_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i]) \in C^{1,\alpha}(\overline{\Omega^o} {\setminus} \Omega) \times C^{1,\alpha}(\overline{\Omega}), \end{align*}
for all 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times C^{0,\alpha }(\partial \Omega )_0 \times \mathbb {R}^2$. We now show that it is bijective. So, we take a pair of functions 
$(h^o,h^i) \in C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega ^o} {\setminus} \Omega ) \times C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega })$ and we prove that there exists a unique quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times C^{0,\alpha }(\partial \Omega )_0 \times \mathbb {R}^2$ such that
\begin{equation} (U^o_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i],U^i_\Omega[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i])=(h^o,h^i). \end{equation}By the uniqueness of the classical solution of the Dirichlet boundary value problem, the second equation in (4.2) is equivalent to
\begin{equation} V_{\partial\Omega}[\eta^i] + \rho^i = h^i_{| \partial \Omega} \end{equation}
(notice that, since 
$h^i$ is an element of 
$C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega })$, we have 
$h^i_{|\partial \Omega } \in C^{1,\alpha } (\partial \Omega ) \subseteq C^{0,\alpha } (\partial \Omega )$). By theorem 3.4(iv), there exists a unique pair 
$(\eta ^i,\rho ^i) \in C^{0,\alpha }(\partial \Omega )_0 \times \mathbb {R}$ such that (4.3) holds. Then it remains to show that there exists a unique triple 
$(\mu ^o,\mu ^i,\rho ^o) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times \mathbb {R}$ such that
\begin{equation} (v^+_{\Omega^o} [\mu^o] + v^-_{\Omega^i}[\mu^i] + \rho^o)_{| \overline{\Omega^o} {\setminus} \Omega} = h^o. \end{equation}By the jump relations for the single-layer potential (cf. theorem 3.4(iii)) and by the uniqueness of the classical solution of the Neumann-Dirichlet mixed boundary value problem, equation (4.4) is equivalent to the following system of integral equations:
\begin{equation} \begin{aligned} & V_{\partial\Omega^o} [\mu^o] + v^-_{\Omega}[\mu^i]_{|\partial\Omega^o} + \rho^o = h^o_{|\partial\Omega^o}, \\ & \left( \frac{1}{2} I + W^\ast_{\partial\Omega} \right) [\mu^i] + \nu_{\Omega} \cdot \nabla v^+_{\Omega^o}[\mu^o]_{|\partial\Omega} = \nu_{\Omega} \cdot \nabla h^o_{|\partial\Omega} \end{aligned} \end{equation}
(notice that, by 
$h^o \in C^{1,\alpha }_{\mathrm {h}}(\overline {\Omega })$, we get 
$h^o_{|\partial \Omega } \in C^{1,\alpha } (\partial \Omega ) \subseteq C^{0,\alpha } (\partial \Omega )$ and 
$\nu _{\Omega } \cdot \nabla h^o_{|\partial \Omega } \in C^{0,\alpha } (\partial \Omega )$). Then we observe that by theorem 3.4(iv), the map from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times \mathbb {R}$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times C^{0,\alpha }(\partial \Omega )$ that takes a triple 
$(\mu ^o,\mu ^i,\rho ^o)$ to the pair of functions 
$(V_{\partial \Omega ^o} [\mu ^o] + \rho ^o, \frac {1}{2} \mu ^i )$ is an isomorphism. Moreover, by the properties of integral operators with real analytic kernel and no singularities (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]) and by theorem 3.5, the map from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times \mathbb {R}$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times C^{0,\alpha }(\partial \Omega )$ that takes a triple 
$(\mu ^o,\mu ^i,\rho ^o)$ to the pair of functions 
$(v^-_{\Omega }[\mu ^i]_{|\partial \Omega ^o}, W^\ast _{\partial \Omega }[\mu ^i] + \nu _{\Omega } \cdot \nabla v^+_{\Omega ^o}[\mu ^o]_{|\partial \Omega })$ is compact. Hence, the map from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ) \times \mathbb {R}$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times C^{0,\alpha }(\partial \Omega )$ that takes a triple 
$(\mu ^o,\mu ^i,\rho ^o)$ to the pair of functions
\[ \left(V_{\partial\Omega^o} [\mu^o] + v^-_{\Omega}[\mu^i]_{|\partial\Omega^o} + \rho^o, \left( \frac{1}{2} I + W^\ast_{\partial\Omega} \right) [\mu^i] + \nu_{\Omega} \cdot \nabla v^+_{\Omega^o}[\mu^o]_{|\partial\Omega} \right) \]
is a compact perturbation of an isomorphism and therefore it is a Fredholm operator of index 0. Thus, to complete the proof, it suffices to show that (4.5) with 
$(h^o_{|\partial \Omega ^o}, \nu _{\Omega } \cdot \nabla h^o_{|\partial \Omega }) = (0,0)$ implies 
$(\mu ^o,\mu ^i,\rho ^o)=(0,0,0)$. If
\begin{align} \left(V_{\partial\Omega^o} [\mu^o] + v^-_{\Omega}[\mu^i]_{|\partial\Omega^o} + \rho^o, \left( \frac{1}{2} I + W^\ast_{\partial\Omega} \right) [\mu^i] + \nu_{\Omega} \cdot \nabla v^+_{\Omega^o}[\mu^o]_{|\partial\Omega} \right) =(0,0), \end{align}
then by the jump relations for the single-layer potential (cf. theorem 3.4(iii)) and by the uniqueness of the classical solution of Neumann-Dirichlet mixed boundary value problem, one deduces that 
$(v^+_{\Omega ^o} [\mu ^o] + v^-_{\Omega }[\mu ^i] + \rho ^o)_{| \overline {\Omega ^o} {\setminus} \Omega }=0$. Moreover, by the continuity of 
$v_{\Omega }[\mu ^i]$ in 
$\mathbb {R}^n$, we have that 
$(v^+_{\Omega ^o} [\mu ^o] + v^-_{\Omega }[\mu ^i] + \rho ^o)_{| \partial \Omega } = (v^+_{\Omega ^o} [\mu ^o] + v^+_{\Omega }[\mu ^i] + \rho ^o)_{|\partial \Omega } = 0$. Then by the uniqueness of the classical solution of Dirichlet boundary value problem in 
$\Omega$ we deduce that
\begin{equation} (v^+_{\Omega^o} [\mu^o] + v^+_{\Omega}[\mu^i] + \rho^o)_{| \overline{\Omega} }=0. \end{equation}
By the jump relations for the single-layer potential (cf. theorem 3.4(iii)), adding and subtracting the term 
$\nu _{\Omega } \cdot \nabla ( v^+_{\Omega ^o}[\mu ^o] +\rho ^o)_{|\partial \Omega }$ and taking into account (4.7), we get
\begin{align*} \mu^i &= \nu_{\Omega} \cdot \nabla v^-_{\Omega}[\mu^i]_{|\partial\Omega} - \nu_{\Omega} \cdot \nabla v^+_{\Omega}[\mu^i]_{|\partial\Omega} \\ &= \nu_{\Omega} \cdot \nabla ( v^+_{\Omega^o}[\mu^o] + v^-_{\Omega}[\mu^i] +\rho^o)_{|\partial\Omega} - \nu_{\Omega} \cdot \nabla ( v^+_{\Omega^o}[\mu^o] + v^+_{\Omega}[\mu^i] +\rho^o)_{|\partial\Omega} = 0. \end{align*}
Thus, by (4.6), we obtain 
$V_{\Omega ^o} [\mu ^o] + \rho ^o = 0$ on 
$\partial \Omega ^o$, which implies 
$(\mu ^o,\rho ^o)=(0,0)$ (cf. theorem 3.4(iv)). Hence 
$(\mu ^o,\mu ^i,\rho ^o)=(0,0,0)$ and the proof is complete.
To represent the boundary conditions of a linearized version of problem (1.3), we find convenient to introduce a matrix function
\[ A({\cdot}) = \begin{pmatrix} A_{11}({\cdot}) & A_{12}({\cdot}) \\ A_{21}({\cdot}) & A_{22}({\cdot}) \end{pmatrix} : \partial \Omega^i \to M_2(\mathbb{R}). \]
Here above, the symbol 
$M_2(\mathbb {R})$ denotes the set of 
$2\times 2$ matrices with real entries. We set
\[ \tilde{A}({\cdot})\equiv \begin{pmatrix} A_{11} ({\cdot}) & A_{12}({\cdot}) \\ -A_{21} ({\cdot}) & -A_{22} ({\cdot}) \end{pmatrix}. \]
We will assume the following conditions on the matrix 
$A$:
\begin{equation} \begin{split} & \bullet \, A_{j,k} \in C^{0,\alpha} (\partial\Omega^i) \mbox{ for all } j,k \in \{1,2\}; \\ & \bullet \,\mbox{For every } (\xi_1,\xi_2) \in \mathbb{R}^2, (\xi_1,\xi_2) \tilde{A} (\xi_1,\xi_2)^T \geq 0 \mbox{ on } \partial\Omega^i; \\ & \bullet \, \mbox{If } (c_1,c_2) \in \mathbb{R}^2 \mbox{ and } A(x)(c_1,c_2)^T = 0 \mbox{ for all } x \in \partial \Omega^i, \mbox{ then } (c_1,c_2)=(0,0). \end{split}\end{equation}
We remark that in literature the third condition in (4.8) is often replaced by a condition on the invertibility of the matrix 
$A$, namely
\begin{equation} \bullet \mbox{ There exists a point } x \in \partial\Omega^i \mbox{ such that } A(x) \mbox{ is invertible}. \end{equation}
We point out that, for instance, the matrix 
$A(x) = \left (\begin {smallmatrix} x_1^2 & x_1\\ -x_1 & -1 \end {smallmatrix}\right )$ with 
$x = (x_1,\dots ,x_n) \in \partial \Omega ^i$ satisfies the third condition in (4.8) but not condition (4.9). Then by a standard energy argument we deduce the following result on the uniqueness of the solution of a transmission problem.
Lemma 4.2 Let 
$A$ be as in (4.8). Then the unique solution in 
$C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i) \times C^{1,\alpha }(\overline {\Omega ^i})$ of problem
\begin{equation} \begin{cases} \Delta u^o = 0 & \mbox{in } \Omega^o {\setminus} \overline{\Omega^i}, \\ \Delta u^i = 0 & \mbox{in } \Omega^i, \\ \nu_{\Omega^o}(x) \cdot \nabla u^o(x)= 0 & \forall x \in \partial \Omega^o, \\ \nu_{\Omega^i}(x) \cdot \nabla u^o (x) - A_{11}(x) u^o(x) - A_{12}(x) u^i(x) = 0 & \forall x \in \partial \Omega^i, \\ \nu_{\Omega^i}(x) \cdot \nabla u^i (x) - A_{21}(x) u^o(x) - A_{22}(x) u^i(x) = 0 & \forall x \in \partial \Omega^i, \end{cases} \end{equation}
is 
$(u^o,u^i)=(0,0)$.
In the following proposition, we investigate the properties of an auxiliary boundary operator, 
$J_A$, which we will exploit in the integral formulation of our problem, in order to recast a fixed point equation. More precisely, we prove that 
$J_A$ is an isomorphism in 
$L^2$, in 
$C^0$, and in 
$C^{0,\alpha }$. All those three frameworks will be important: the first setting is suitable to use Fredholm theory and to directly prove the isomorphic property of 
$J_A$, the second setting will be used in order to apply Leray-Shauder Theorem to the aforementioned obtained fixed point equation (see propositions 4.5 and 4.7 below) and the third setting will be central to deduce that the solution of problem (1.3) we built is actually a classical solution, in particular of class 
$C^{1,\alpha }$ (cf. propositions 4.4 and 4.8)
Proposition 4.3 Let 
$A$ be as in (4.8). Let 
$J_A$ be the map from 
$L^2(\partial \Omega ^o)_0 \times L^2(\partial \Omega ^i) \times L^2(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2$ that takes a quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the triple 
$J_A[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i]$ defined by
\begin{align}
J_{A,1}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &
\equiv \left( -\frac{1}{2} I + W^\ast_{\partial\Omega^o}
\right) [\mu^o] + \nu_{\Omega^o} \cdot \nabla
v^-_{\Omega^i}[\mu^i]_{|\partial\Omega^o} \qquad \mbox{on }
\partial\Omega^o, \notag\\
J_{A,2}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & \equiv \left(
\frac{1}{2} I + W^\ast_{\partial\Omega^i} \right) [\mu^i] +
\nu_{\Omega^i} \cdot \nabla
v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} -
(A_{11},A_{12})\notag\\ & \quad \cdot
(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} +
V_{\partial\Omega^i}[\mu^i] + \rho^o
,V_{\partial\Omega^i}[\eta^i] + \rho^i ) \quad\mbox{on }
\partial\Omega^i, \notag\\
J_{A,3}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & \equiv \left(
-\frac{1}{2} I + W^\ast_{\partial\Omega^i} \right) [\eta^i]
- (A_{21},A_{22})\notag\\ & \quad \cdot
(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} +
V_{\partial\Omega^i}[\mu^i] + \rho^o
,V_{\partial\Omega^i}[\eta^i] + \rho^i ) \quad \mbox{on }\partial\Omega^i. \end{align}Then the following statements hold.
(i)
$J_A$ is a linear isomorphism from $L^2(\partial \Omega ^o)_0 \times L^2(\partial \Omega ^i) \times L^2(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to $L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2$.(ii)
$J_A$ is a linear isomorphism from $C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to $C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$.(iii)
$J_A$ is a linear isomorphism from $C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to $C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$.
Proof. We first prove (i). We write 
$J_A$ in the form 
$J_A = \tilde {J}^{+}_A \circ \tilde {J}_A \circ \tilde {J}^{-}_A$, where 
$\tilde {J}^{-}_A$ is the inclusion of 
$L^2(\partial \Omega ^o)_0 \times L^2(\partial \Omega ^i) \times L^2(\partial \Omega ^i)_0 \times \mathbb {R}^2$ into 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2 \times \mathbb {R}^2$, 
$\tilde {J}_A$ is the map from 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2 \times \mathbb {R}^2$ into itself that takes 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the quintuple 
$\tilde {J}_A[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i]$ defined by
\begin{align*} \tilde{J}_{A,1}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &\equiv \left( -\frac{1}{2} I + W^\ast_{\partial\Omega^o} \right) [\mu^o] + \nu_{\Omega^o} \cdot \nabla v^-_{\Omega^i}[\mu^i]_{|\partial\Omega^o} \quad \mbox{on } \partial\Omega^o, \\ \tilde{J}_{A,2}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &\equiv \left( \frac{1}{2} I + W^\ast_{\partial\Omega^i} \right) [\mu^i] + \nu_{\Omega^i} \cdot \nabla v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i}\\ &\quad - (A_{11},A_{12}) \cdot (v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] ,V_{\partial\Omega^i}[\eta^i] ) \quad \mbox{on } \partial\Omega^i, \\ \tilde{J}_{A,3}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &\equiv \left( -\frac{1}{2} I + W^\ast_{\partial\Omega^i} \right) [\eta^i]\\ &\quad - (A_{21},A_{22}) \cdot (v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i],V_{\partial\Omega^i}[\eta^i]) \quad \mbox{on } \partial\Omega^i, \\ \tilde{J}_{A,4}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &\equiv \rho^o, \\ \tilde{J}_{A,5}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &\equiv \rho^i, \end{align*}
and 
$\tilde {J}^{+}_A$ is the map from 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2 \times \mathbb {R}^2$ into 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2$ that takes a quintuple 
$(f,g_1,g_2,c_1,c_2)$ to the triple 
$\tilde {J}^{+}_A[f,g_1,g_2,c_1,c_2]$ defined by
\[ \tilde{J}^{+}_A[f,g_1,g_2,c_1,c_2] \equiv (f, g_1 - (A_{11},A_{12}) \cdot (c_1,c_2),g_2 - (A_{21},A_{22}) \cdot (c_1,c_2) ) . \]
Then we observe that 
$\tilde {J}^{+}_A$ is a Fredholm operator of index 
$2$, because 
$\mathrm {Coker}\, \tilde {J}^{+}_A = \{0\}$ and 
$\mathrm {Ker}\, \tilde {J}^{+}_A = \mathrm {Span} \{(0,A_{11},A_{21},1,0), (0,A_{12},A_{22},0,1)\}$, and that 
$\tilde {J}^{-}_A$ is Fredholm of index 
$-2$, because 
$\mathrm {Ker}\, \tilde {J}^{-}_A = \{0\}$ and 
$\mathrm {Coker}\, \tilde {J}^{-}_A = \mathrm {Span} \{(1,0,0,0,0), (0,0,1,0,0)\}$. Next, we observe that the map from 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2 \times \mathbb {R}^2$ into itself that takes a quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the quintuple 
$(-\frac {1}{2}\mu ^o,\frac {1}{2}\mu ^i,-\frac {1}{2}\eta ^i,\rho ^o,\rho ^i)$ is a linear isomorphism. Moreover, by the mapping properties of the integral operators with real analytic kernel and no singularity (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]), by the compactness of the operators 
$W^\ast _{\partial \Omega ^o}$ and 
$W^\ast _{\partial \Omega ^i}$ from 
$L^2(\partial \Omega ^o)$ to itself and from 
$L^2(\partial \Omega ^i)$ to itself, respectively (see comments below definition 3.3), by the compactness of the operator 
$V_{\partial \Omega ^i}$ from 
$L^2(\partial \Omega ^i)$ into itself (see Costabel [Reference Costabel4, theorem 1]), and by the bilinearity and continuity of the product from 
$C^{0,\alpha }(\partial \Omega ^i) \times L^2(\partial \Omega ^i)$ to 
$L^2(\partial \Omega ^i)$, we deduce that the map from 
$L^2(\partial \Omega ^o) \times (L^2(\partial \Omega ^i))^2 \times \mathbb {R}^2$ into itself that takes a quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the quintuple 
$\tilde {J}^C_A[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i]$ defined by
\begin{align*} \tilde{J}^C_{A,1}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] &= W^\ast_{\partial\Omega^o} [\mu^o] + \nu_{\Omega^o} \cdot \nabla v^-_{\Omega^i}[\mu^i]_{|\partial\Omega^o} \quad \mbox{on }\partial\Omega^o, \\ \tilde{J}^C_{A,2}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & = W^\ast_{\partial\Omega^i} [\mu^i] + \nu_{\Omega^i} \cdot \nabla v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i}\\ &\quad - (A_{11},A_{12}) \cdot (v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] ,V_{\partial\Omega^i}[\eta^i] ) \quad \mbox{on } \partial\Omega^i, \\ \tilde{J}^C_{A,3}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & = W^\ast_{\partial\Omega^i} [\eta^i] \\ &\quad - (A_{21},A_{22}) \cdot (v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i],V_{\partial\Omega^i}[\eta^i]) \quad \mbox{on } \partial\Omega^i, \\ \tilde{J}^C_{A,4}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & = 0, \\ \tilde{J}^C_{A,5}[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] & = 0, \end{align*}
is compact. Hence, we conclude that 
$\tilde {J}_A$ is a compact perturbation of an isomorphism and therefore it is Fredholm of index 0. Since the index of a composition of Fredholm operators is the sum of the indexes of the components, we deduce that 
$J_A$ is a Fredholm operator of index 
$0$. Therefore, in order to complete the proof of point 
$(i)$, it suffices to prove that 
$J_A$ is injective. Thus, we now assume that 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in L^2(\partial \Omega ^o)_0 \times L^2(\partial \Omega ^i) \times L^2(\partial \Omega ^i)_0 \times \mathbb {R}^2$ and that
\begin{equation} J_A[\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] = (0,0,0). \end{equation}
We first verify that 
$(\mu ^o,\mu ^i,\eta ^i)$ is actually in 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0$. In fact, by the definition of 
$J_{A,1}$ in (4.11), by the fact that 
$\nu _{\Omega ^o} \cdot \nabla v^-_{\Omega ^i}[\mu ^i]_{|\partial \Omega ^o} \in C^{0}(\partial \Omega ^o)$ (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]), and by theorem 3.8, we deduce that 
$\mu ^o \in C^{0}(\partial \Omega ^o)$. Let 
$\mathcal {K}\equiv (\mathcal {K}_1, \mathcal {K}_2)$ be the map from 
$L^2(\partial \Omega ^i)\times L^2(\partial \Omega ^i)_0$ to itself that takes a pair 
$(\mu ^i,\eta ^i)\in L^2(\partial \Omega ^i)\times L^2(\partial \Omega ^i)_0$ to
\begin{align*} \mathcal{K}_1[\mu^i,\eta^i] &\equiv 2W^\ast_{\partial\Omega^i}[\mu^i] - 2A_{11} V_{\partial\Omega^i}[\mu^i] - 2A_{12}V_{\partial\Omega^i}[\eta^i] \quad\mbox{on } \partial\Omega^i, \\ \mathcal{K}_2[\mu^i,\eta^i] &\equiv{-}2W^\ast_{\partial\Omega^i} [\eta^i] +2 A_{21} V_{\partial\Omega^i}[\mu^i] +2 A_{22} V_{\partial\Omega^i}[\eta^i] \quad \mbox{on } \partial\Omega^i. \end{align*}
Notice that each component of 
$\mathcal {K}$ is a linear combination of integral operators with a continuous kernel of order 
$n-2$ (see definition 3.6 and comments below). By the fact that 
$\mu ^o \in C^0(\partial \Omega ^i)$ and by the first condition in (4.8), we know that 
$2(A_{11},A_{12}) \cdot (v^+_{\Omega ^o}[\mu ^o]_{|\partial \Omega ^i} + \rho ^o , \rho ^i ), -2(A_{21},A_{22}) \cdot (v^+_{\Omega ^o}[\mu ^o]_{|\partial \Omega ^i} + \rho ^o ,\rho ^i )$ belong to 
$C^{0,\alpha }(\partial \Omega ^i) \subseteq C^{0}(\partial \Omega ^i)$. Then (4.12) and the definition of the operator 
$\mathcal {K}$ imply that 
$(I + \mathcal {K})[\mu ^i,\eta ^i] \in (C^{0}(\partial \Omega ^i))^2$. Hence, by theorem 3.7 we conclude that 
$(\mu ^i,\eta ^i) \in C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0$. Then by mapping properties of integral operators with real analytic kernel and no singularity (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]) and by classical known results in potential theory (cf. Miranda [Reference Miranda23, chapter II, § 14]), we know that 
$\nu _{\Omega ^o} \cdot \nabla v^-_{\Omega ^i}[\mu ^i]_{|\partial \Omega ^o} \in C^{0,\alpha }(\partial \Omega ^o)$ and 
$v^+_{\Omega ^o}[\mu ^o]_{|\partial \Omega ^i}$, 
$\nu _{\Omega ^i} \cdot \nabla v^+_{\Omega ^o}[\mu ^o]_{|\partial \Omega ^i}$, 
$V_{\partial \Omega ^i}[\eta ^i]$, 
$V_{\partial \Omega ^i}[\mu ^i]\in C^{0,\alpha }(\partial \Omega ^i)$. Hence, by (4.12) and by the membership of 
$A \in M_2(C^{0,\alpha } (\partial \Omega ^i))$ (cf. first condition in (4.8)), we obtain that 
$( -\frac {1}{2} I + W^\ast _{\partial \Omega ^o} ) [\mu ^o] \in C^{0,\alpha }(\partial \Omega ^o)$ and 
$( \frac {1}{2} I + W^\ast _{\partial \Omega ^i} ) [\mu ^i], ( -\frac {1}{2} I + W^\ast _{\partial \Omega ^i} ) [\eta ^i] \in C^{0,\alpha }(\partial \Omega ^i)$. Then theorem 3.8 implies 
$(\mu ^o,\mu ^i,\eta ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0$. By the jump relations (cf. theorem 3.4(iii)), by lemma 4.1, and by (4.12), we deduce that the pair 
$(U^o_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i],U^i_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i])$ defined by (4.1) is a solution of the boundary value problem (4.10). Then by lemma 4.2, we have that 
$(U^o_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i],U^i_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i])=(0,0)$, which implies 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)=(0,0,0,0,0)$, by the uniqueness of the representation provided by lemma 4.1. We now prove statement (ii). First we note that the integral operators that appear in the definition (4.11) of 
$J_A$ have either a weakly singular or a real analytic kernel. It follows that 
$J_A$ is continuous from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1] for the properties of integral operators with real analytic kernels). Then we observe that, by theorems 3.7 and 3.8, if we have 
$J_A[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i] \in C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ for some 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in L^2(\partial \Omega ^o)_0 \times L^2(\partial \Omega ^i) \times L^2(\partial \Omega ^i)_0 \times \mathbb {R}^2$, then 
$(\mu ^o,\mu ^i,\eta ^i) \in C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0$ (see also the argument used after (4.12) to prove that 
$(\mu ^o,\mu ^i,\eta ^i)$ belongs to 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0$). Then, by statement (i) we deduce that 
$J_A$ is a bijective continuous linear map from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$. By the Open Mapping Theorem it follows that 
$J_A$ is a linear homeomorphism from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$. The proof of statement (iii) is similar to that of statement (ii) and we leave it to the zealous reader (see also the argument used after (4.12) to prove that 
$(\mu ^o,\mu ^i,\eta ^i)$ belongs to 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0$).
We are now ready to convert (1.3) into a system of integral equations.
Proposition 4.4 Let 
$A$ be as in (4.8). Let 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$. Let 
$(U^o_{\Omega ^i}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ],U^i_{\Omega ^i}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ])$ be defined by (4.1). Let 
$J_A$ be as in proposition 4.3. Then 
$(U^o_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i],U^i_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i])$ is a solution of (1.3) if and only if
\begin{align} \begin{pmatrix} \mu^o \\ \mu^i \\ \eta^i \\ \rho^o \\ \rho^i \end{pmatrix} &= J_A^{({-}1)} \left[ \begin{pmatrix} f^o \\ \mathcal{N}_{F_1}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \\ \mathcal{N}_{F_2}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \end{pmatrix} \right. \nonumber\\ & \quad \left. - \begin{pmatrix} 0 & 0 & 0 \\ 0 & A_{11} & A_{12} \\ 0 & A_{21} & A_{22} \end{pmatrix} \begin{pmatrix} 0 \\ v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o \\ V_{\partial\Omega^i}[\eta^i] +\rho^i \end{pmatrix} \right]. \end{align}Proof. By lemma 4.1 and by the jump relations of theorem 3.4, we know that if 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ then the pair 
$(U^o_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i],U^i_{\Omega ^i}[\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i])$ defined by (4.1) is a solution of problem (1.3) if and only if
\begin{align} &
\begin{pmatrix} \left( -\frac{1}{2} I +
W^\ast_{\partial\Omega^o} \right) [\mu^o] + \nu_{\Omega^o}
\cdot \nabla v^-_{\Omega^i}[\mu^i]_{|\partial\Omega^o} \\
\left( \frac{1}{2} I + W^\ast_{\partial\Omega^i} \right)
[\mu^i] + \nu_{\Omega^i} \cdot \nabla
v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} \\ \left(
-\frac{1}{2} I + W^\ast_{\partial\Omega^i} \right) [\eta^i]
\end{pmatrix}\nonumber\\ &\quad = \begin{pmatrix} f^o \\
\mathcal{N}_{F_1}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i}
+ V_{\partial\Omega^i}[\mu^i] +\rho^o
,V_{\partial\Omega^i}[\eta^i] +\rho^i) \\
\mathcal{N}_{F_2}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i}
+ V_{\partial\Omega^i}[\mu^i] +\rho^o
,V_{\partial\Omega^i}[\eta^i] +\rho^i)
\end{pmatrix}.
\end{align}Then, by subtracting in both sides of (4.14) the term
\[ \begin{pmatrix} 0 & 0 & 0 \\ 0 & A_{11} & A_{12} \\ 0 & A_{21} & A_{22} \end{pmatrix} \begin{pmatrix} 0 \\ v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o \\ V_{\partial\Omega^i}[\eta^i] +\rho^i \end{pmatrix} \in C^{0,\alpha}(\partial\Omega^o) \times (C^{0,\alpha}(\partial\Omega^i))^2 \]
and by the invertibility of 
$J_A$ from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$ provided by proposition 4.3(iii), the validity of the statement follows.
We now introduce an auxiliary map. If 
$A$ is as in (4.8) and 
$J_A$ is as in proposition 4.3, we denote by 
$T_A$ the map from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ defined by
\begin{align} &T_A(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i) \nonumber\\ &\quad \equiv J_A^{({-}1)} \left[ \begin{pmatrix} f^o \\ \mathcal{N}_{F_1}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \\ \mathcal{N}_{F_2}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \end{pmatrix} \right.\nonumber\\ &\qquad \left. - \begin{pmatrix} 0 & 0 & 0 \\ 0 & A_{11} & A_{12} \\ 0 & A_{21} & A_{22} \end{pmatrix} \begin{pmatrix} 0 \\ v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o \\ V_{\partial\Omega^i}[\eta^i] +\rho^i \end{pmatrix} \right]. \end{align}
We study the continuity and compactness of 
$T_A$ in the following proposition.
Proposition 4.5 Let 
$A$ be as in (4.8). Let 
$T_A$ be as in (4.15). Then 
$T_A$ is a continuous (nonlinear) operator from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ and is compact.
Proof. By the properties of integral operators with real analytic kernel and no singularities (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]) and by the compactness of the embedding of 
$C^{0,\alpha }(\partial \Omega ^i)$ into 
$C^0(\partial \Omega ^i)$, 
$v^+_{\Omega ^o}[\cdot ]_{|\partial \Omega ^i}$ is compact from 
$C^0(\partial \Omega ^o)_0$ into 
$C^0(\partial \Omega ^i)$. By mapping properties of the single-layer potential (cf. Miranda [Reference Miranda23, chapter II, § 14, III]) and by the compactness of the embedding of 
$C^{0,\alpha }(\partial \Omega ^i)$ into 
$C^0(\partial \Omega ^i)$, 
$V_{\partial \Omega ^i}$ is compact from 
$C^0(\partial \Omega ^i)$ into itself. Hence, by the bilinearity and continuity of the product of continuous functions, the map from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ that takes the quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the triple given by
\[ \begin{pmatrix} 0 & 0 & 0 \\ 0 & A_{11} & A_{12} \\ 0 & A_{21} & A_{22} \end{pmatrix} \begin{pmatrix} 0 \\ v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o \\ V_{\partial\Omega^i}[\eta^i] +\rho^i \end{pmatrix} \]
is continuous and maps bounded sets into sets with compact closure, i.e., is compact. Moreover, by assumption (1.2), one readily verifies that the operators 
$\mathcal {N}_{F_1}$ and 
$\mathcal {N}_{F_2}$ are continuous from 
$(C^0(\partial \Omega ^i))^2$ into 
$C^0(\partial \Omega ^i)$. Hence the map from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ that takes the quintuple 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the triple
\[ \begin{pmatrix} f^0 \\ \mathcal{N}_{F_1}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \\ \mathcal{N}_{F_2}(v^+_{\Omega^o}[\mu^o]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i] +\rho^o ,V_{\partial\Omega^i}[\eta^i] +\rho^i) \end{pmatrix} \]
is compact. Finally, by proposition 4.3(ii), 
$J_A$ is a linear isomorphism from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ and, accordingly, 
$T_A$ is compact.
In what follows we will assume the following growth condition on the pair 
$(F_1,F_2)$ with respect to the matrix function 
$A$ defined as in (4.8):
\begin{align} \bullet\, &\mbox{There exist two constants } C_F \in ]0,+\infty[ \mbox{ and } \delta \in ]0,1[ \mbox{ such that} \nonumber\\ & \qquad \left|\begin{pmatrix} F_1(x, \zeta_1,\zeta_2) \\ F_2(x, \zeta_1,\zeta_2) \end{pmatrix} - A(x) \begin{pmatrix} \zeta_1\\ \zeta_2 \end{pmatrix} \right| \leq C_F (1 + |\zeta_1| + |\zeta_2|)^\delta\nonumber\\ & \mbox{for all } (x,\zeta_1,\zeta_2) \in \partial\Omega^i \times \mathbb{R}^2. \end{align}
In proposition 4.7 below we prove the existence of a solution in 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ of (4.13). Our argument exploits the Leray-Schauder Fixed-Point Theorem (cf. Gilbarg and Trudinger [Reference Gilbarg and Trudinger12, theorem 11.3]).
Theorem 4.6 Leray-Schauder Theorem
Let 
$X$ be a Banach space. Let 
$T$ be a continuous operator from 
$X$ into itself. If 
$T$ is compact and there exists a constant 
$M \in ]0,+\infty [$ such that 
$\|x\|_{X} \leq M$ for all 
$(x,\lambda ) \in X \times [0,1]$ satisfying 
$x=\lambda T(x)$, then 
$T$ has at least one fixed point 
$x \in X$ such that 
$\|x\|_{X} \leq M$.
Then we have the following.
Proposition 4.7 Let 
$A$ be as in (4.8). Let assumption (4.16) holds. Let 
$J_A$ be as in proposition 4.3. Then the nonlinear system (4.13) has at least one solution 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$.
Proof. We plan to apply the Leray-Schauder Theorem 4.6 to the operator 
$T_A$ defined by (4.15) in the Banach space 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$. By proposition 4.5 we already know that 
$T_A$ is a continuous operator from 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ and maps bounded sets into sets with compact closure. So in order to apply the Leray-Schauder Theorem 4.6, we are left to show that if 
$\lambda \in ]0,1[$ and if
\begin{equation} (\mu^o,\mu^i,\eta^i,\rho^o,\rho^i) = \lambda T_A(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i) \end{equation}
with 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$, then there exists a constant 
$C \in ]0,+\infty [$ (which does not depend on 
$\lambda$ and 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$), such that
\begin{equation} \|(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)\|_{C^0(\partial\Omega^o) \times (C^0(\partial\Omega^i))^2\times \mathbb{R}^2} \leq C. \end{equation}
By (4.17) and by 
$|\lambda |<1$, we readily deduce that
\begin{align} &\|(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)\|_{C^0(\partial\Omega^o) \times (C^0(\partial\Omega^i))^2 \times \mathbb{R}^2} \nonumber\\ &\quad\leq \|T_A(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)\|_{C^0(\partial\Omega^o) \times (C^0(\partial\Omega^i))^2}. \end{align}By the growth condition (4.16), we can show that
\begin{align} &\left\| \begin{pmatrix} \mathcal{N}_{F_1}(h^i_1,h^i_2) \\ \mathcal{N}_{F_2}(h^i_1,h^i_2) \end{pmatrix} - A \begin{pmatrix} h^i_1 \\ h^i_2 \end{pmatrix} \right\|_{(C^0(\partial\Omega^i))^2} \nonumber\\ & \quad\leq C_F (1+\|h^i_1\|_{C^0(\partial\Omega^i)}+\|h^i_2\|_{C^0(\partial\Omega^i)})^\delta \end{align}
for all pair of functions 
$(h^i_1,h^i_2) \in (C^0(\partial \Omega ^i))^2$. Hence, by (4.19) and by the definition of 
$T_A$ in (4.15), we deduce that there exist two constants 
$C_1,C_2 \in ]0,+\infty [$, which depend only on the operator norm of 
$J_A^{(-1)}$ from 
$C^0(\partial \Omega ^o) \times (C^0(\partial \Omega ^i))^2$ to 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ (cf. theorem 4.3(ii)), on 
$\|f^o\|_{\partial \Omega ^o}$, on the constant 
$C_F \in ]0,+\infty [$ provided by the growth condition (4.16) (cf. (4.20)), on the norm of the bounded linear operator 
$v^+_{\Omega ^o}[\cdot ]_{|\partial \Omega ^i}$ from 
$C^0(\partial \Omega ^o)$ to 
$C^0(\partial \Omega ^i)$, and on the norm of the bounded linear operator 
$V_{\partial \Omega ^i}$ from 
$C^0(\partial \Omega ^i)$ into itself, such that
\begin{align*} &\|(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)\|_{C^0(\partial\Omega^o) \times (C^0(\partial\Omega^i))^2 \times \mathbb{R}^2}\\ & \quad \leq C_1 (C_2 + \|(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)\|_{C^0(\partial\Omega^o) \times (C^0(\partial\Omega^i))^2 \times \mathbb{R}^2} )^\delta. \end{align*}
Then, by a straightforward calculation, we can show the existence of a constant 
$C>0$ such that inequality (4.18) holds true (cf. Lanza de Cristoforis [Reference Lanza de Cristoforis17, proof of theorem 7.2]). Hence, by the Leray-Schauder Theorem 4.6 there exists at least one solution 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ of 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) = T_A(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$. Finally, by the definition of 
$T_A$ (cf. (4.15)), we conclude that 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0)$ is a solution in 
$C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ of the nonlinear system (4.13).
In what follows we will exploit a continuity condition on the superposition operators generated by 
$F_1$ and 
$F_2$, namely
\begin{align} \bullet & \, \mbox{The superposition operators } \mathcal{N}_{F_1} \mbox{ and } \mathcal{N}_{F_2} \mbox{ are continuous from } \nonumber\\ &\quad (C^{0,\alpha}(\partial\Omega^i))^2 \mbox{ into } C^{0,\alpha}(\partial\Omega^i). \end{align}
For conditions on 
$F_1$ and 
$F_2$ which imply the validity of assumption (4.21), we refer to Appell and Zabrejko [Reference Appell and Zabrejko1, chapter 8] and to Valent [Reference Valent32, chapter II]. Then we can prove a regularity result for the fixed point provided by proposition 4.7, and, thus, an existence result for problem (1.3).
Proposition 4.8 Let A be as in (4.8). Let assumptions (4.16) and (4.21) hold. Then the nonlinear system (4.13) has at least one solution 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$. In particular, problem (1.3) has at least one solution 
$(u^o_0,u^i_0) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i) \times C^{1,\alpha }(\overline {\Omega ^i})$ given by
\begin{equation} (u^o_0,u^i_0) \equiv (U^o_{\Omega^i}[\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0],U^i_{\Omega^i}[\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0]) \end{equation}
where the pair 
$(U^o_{\Omega ^i}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ],U^i_{\Omega ^i}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ])$ is defined by (4.1).
Proof. Let 
$T_A$ be as in (4.15). By proposition 4.7, we deduce the existence of a quintuple 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^0(\partial \Omega ^o)_0 \times C^0(\partial \Omega ^i) \times C^0(\partial \Omega ^i)_0 \times \mathbb {R}^2$ such that 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) = T_A(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0)$. By the mapping properties of integral operators with real analytic kernel and no singularities (cf. Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]), 
$v^+_{\Omega ^o}[\mu ^o_0]_{|\partial \Omega ^i}$ belongs to 
$C^{0,\alpha }(\partial \Omega ^i)$. By classical results in potential theory (cf. Miranda [Reference Miranda23, chapter II, § 14, III]), 
$V_{\Omega ^i}[\mu ^i_0]$ and 
$V_{\Omega ^i}[\eta ^i_0]$ belong to 
$C^{0,\alpha }(\partial \Omega ^i)$. Then, by condition (4.21) and by the membership of 
$A \in M_2(C^{0,\alpha } (\partial \Omega ^i))$ and of 
$f^o \in C^{0,\alpha } (\partial \Omega ^o)$, we obtain that
\begin{align*} &\begin{pmatrix} f^o \\ \mathcal{N}_{F_1}(v^+_{\Omega^o}[\mu^o_0]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i_0] +\rho^o_0 ,V_{\partial\Omega^i}[\eta^i_0] +\rho^i_0) \\ \mathcal{N}_{F_2}(v^+_{\Omega^o}[\mu^o_0]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i_0] +\rho^o_0 ,V_{\partial\Omega^i}[\eta^i_0] +\rho^i_0) \end{pmatrix} \\ &\quad - \begin{pmatrix} 0 & 0 & 0 \\ 0 & A_{11} & A_{12} \\ 0 & A_{21} & A_{22} \end{pmatrix} \begin{pmatrix} 0 \\ v^+_{\Omega^o}[\mu^o_0]_{|\partial\Omega^i} + V_{\partial\Omega^i}[\mu^i_0] +\rho^o_0 \\ V_{\partial\Omega^i}[\eta^i_0] +\rho^i_0 \end{pmatrix} \end{align*}
belongs to the product space 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$. Finally, by the invertibility of the operator 
$J_A$ from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$, we obtain that 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$. In particular, by proposition 4.4 we deduce that the pair given by (4.22) is a solution of (1.3) (cf. (4.15)).
5. The perturbed transmission problem (1.6)
This section is devoted to the study of the perturbed transmission problem (1.6). We introduce the map 
$M=(M_1,M_2,M_3)$ from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$ defined by
\begin{align} & M_1[\phi,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i](x) \nonumber\\ &\quad \equiv \left( -\frac{1}{2} I + W^\ast_{\partial\Omega^o} \right) [\mu^o] (x) + \nu_{\Omega^o}(x) \cdot \nabla v^-_{\Omega^i[\phi]}[\mu^i \circ \phi^{({-}1)}] (x) - f^o(x) \quad \forall x \in \partial\Omega^o\nonumber\\ & M_2[\phi,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] (t) \nonumber\\ &\quad \equiv \left( \frac{1}{2} I + W^\ast_{\partial\Omega^i[\phi]} \right) [\mu^i \circ \phi^{({-}1)}] (\phi(t)) + \nu_{\Omega^i[\phi]}(\phi(t)) \cdot \nabla v^+_{\Omega^o}[\mu^o](\phi(t))\nonumber\\ & \qquad - F_1\bigg(t,v^+_{\Omega^o}[\mu^o](\phi(t)) + V_{\partial\Omega^i[\phi]}[\mu^i\circ \phi^{({-}1)}](\phi(t)) +\rho^o , \nonumber\\ & \qquad \qquad V_{\partial\Omega^i[\phi]}[\eta^i\circ \phi^{({-}1)}](\phi(t)) +\rho^i \bigg) \quad \forall t \in \partial\Omega^i\nonumber\\ & M_3[\phi,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] (t) \nonumber\\ &\quad \equiv \left( -\frac{1}{2} I + W^\ast_{\partial\Omega^i[\phi]} \right) [\eta^i \circ \phi^{({-}1)}] (\phi(t)) \nonumber\\ & \qquad - F_2\bigg(t,v^+_{\Omega^o}[\mu^o](\phi(t)) + V_{\partial\Omega^i[\phi]}[\mu^i\circ \phi^{({-}1)}](\phi(t)) +\rho^o ,\nonumber\\ & \qquad \qquad V_{\partial\Omega^i[\phi]}[\eta^i\circ \phi^{({-}1)}](\phi(t)) +\rho^i \bigg ) \quad\forall t \in \partial\Omega^i \end{align}
for all 
$(\phi ,\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in \mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$. We incidentally observe that by the definition of 
$\Omega ^i[\phi ]$ we have that 
$\partial \Omega ^i[\phi ]=\phi (\partial \Omega ^i)$. Then, by the definition of 
$M$, we can deduce the following result.
Proposition 5.1 Let A be as in (4.8). Let assumptions (4.16) and (4.21) hold. Let
\[ (\phi,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i) \in \mathcal{A}^{\Omega^o}_{\partial\Omega^i} \times C^{0,\alpha}(\partial\Omega^o)_0 \times C^{0,\alpha}(\partial\Omega^i) \times C^{0,\alpha}(\partial\Omega^i)_0 \times \mathbb{R}^2. \]
Then the pair of functions
\[ (U^o_{\Omega^i[\phi]}[\mu^o,\mu^i\circ \phi^{({-}1)},\eta^i\circ \phi^{({-}1)},\rho^o,\rho^i],U^i_{\Omega^i[\phi]}[\mu^o,\mu^i\circ \phi^{({-}1)},\eta^i\circ \phi^{({-}1)},\rho^o,\rho^i]) \]
defined by (4.1) is a solution of problem (1.6) if and only if
\begin{equation} M[\phi,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] = (0,0,0). \end{equation}In particular, equation
\begin{equation} M[\phi_0,\mu^o,\mu^i,\eta^i,\rho^o,\rho^i] = (0,0,0) \end{equation}
is equivalent to the system (4.13) and has a solution 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ (recall that 
$\phi _0\equiv \text {id}_{\partial \Omega ^i}$).
Proof. We first observe that, by the regularity of 
$\phi \in \mathcal {A}^{\Omega ^o}_{\partial \Omega ^i}$, if 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$, then
\[ (\mu^o,\mu^i,\eta^i,\rho^o,\rho^i) \in C^{0,\alpha}(\partial\Omega^o)_0 \times C^{0,\alpha}(\partial\Omega^i) \times C^{0,\alpha}(\partial\Omega^i)_0 \times \mathbb{R}^2. \]
Moreover, since 
$\overline {\Omega ^i[\phi ]} \subset \Omega ^o$, we can apply lemma 4.1 with 
$\Omega =\Omega ^i[\phi ]$. Then by the jump relations for the single-layer potential (cf. theorem 3.4(iii)), by a change of variable on 
$\phi (\partial \Omega ^i)$ and by the definition of 
$M$ (cf. (5.1)), we obtain that the pair of functions
\begin{align*} & U^o_{\Omega^i[\phi]}[\mu^o,\mu^i\circ \phi^{({-}1)},\eta^i\circ \phi^{({-}1)},\rho^o,\rho^i] = (v^+_{\Omega^o} [\mu^o] + v^-_{\Omega^i[\phi]}[\mu^i\circ \phi^{({-}1)}] + \rho^o)_{| \overline{\Omega^o} {\setminus} \Omega^i[\phi]}, \\ & U^i_{\Omega^i[\phi]}[\mu^o,\mu^i\circ \phi^{({-}1)},\eta^i\circ \phi^{({-}1)},\rho^o,\rho^i] = v^+_{\Omega^i[\phi]}[\eta^i\circ \phi^{({-}1)}] + \rho^i \end{align*}
is a solution of problem (1.6) if and only if (5.2) is satisfied. Finally, since 
$\phi _0\equiv \text {id}_{\partial \Omega ^i}$ (cf. (1.5)) and by the definition of 
$J_A$ (cf. (4.11)), we obtain that, for all 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$, equation (5.3) is equivalent to the system (4.13). Then the existence of a solution 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0)$ of (5.3) follows by proposition 4.7.
By proposition 5.1, the study of problem (1.6) is reduced to that of equation (5.2). We now wish to apply the Implicit Function Theorem for real analytic maps in Banach spaces (cf. Deimling [Reference Deimling9, theorem 15.3]) to equation (5.2) around the value 
$\phi _0$. As a first step we have to analyse the regularity of the map 
$M$.
In what follows we will assume the following:
\begin{align} \bullet & \,\mbox{The superposition operators } \mathcal{N}_{F_1} \mbox{ and } \mathcal{N}_{F_2} \mbox{ are real analytic from } \nonumber\\ &\quad (C^{0,\alpha}(\partial\Omega^i))^2 \mbox{ into } C^{0,\alpha}(\partial\Omega^i). \end{align} For conditions on 
$F_1$ and 
$F_2$ which imply the validity of assumption (5.4), we refer to Valent [Reference Valent32, chapter II]. We now show that 
$M$ is real analytic.
Proposition 5.2 Let assumption (5.4) holds. Then the map 
$M$ is real analytic from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$.
Proof. We only prove the analyticity of 
$M_2$. The analyticity of 
$M_1$ and of 
$M_3$ can be proved similarly and it is left to the reader. Therefore, we now analyse 
$M_2$. The map from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^i)$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$(\phi ,\mu ^i)$ to the function of the variable 
$t\in \partial \Omega ^i$ defined by
\begin{align*} &\left( \frac{1}{2} I + W^\ast_{\partial\Omega^i[\phi]} \right) [\mu^i \circ \phi^{({-}1)}] (\phi(t)) = \frac{1}{2}\mu^i(t) + W^\ast_{\partial\Omega^i[\phi]} [ \mu^i\circ \phi^{({-}1)}] (\phi(t)) \\ &\quad= \frac{1}{2}\mu^i(t) + \int_{\partial\Omega^i} (\nu_{\Omega^i[\phi]}(\phi(t))) \cdot \nabla S_n(\phi(t)-\phi(s))) \,\mu^i(s)\, \tilde{\sigma}_n[\phi](s) \,{\rm d}\sigma_s \end{align*}
is real analytic by the real analyticity result for the dependence of layer potentials upon perturbation of the support and of the density of Lanza de Cristoforis and Rossi [Reference Lanza de Cristoforis and Rossi21, theorem 3.12] and Lanza de Cristoforis [Reference Lanza de Cristoforis18, proposition 7] (see also lemma 2.1). The map from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$(\phi ,\mu ^o)$ to the function of the variable 
$t \in \partial \Omega ^i$ defined by
\[ \nu_{\Omega^i[\phi]}(\phi(t)) \cdot \nabla v^+_{\Omega^o}[\mu^o](\phi(t)) = \int_{\partial\Omega^o} (\nu_{\Omega^i[\phi]}(\phi(t)) \cdot \nabla S_n(\phi(t)-y)) \, \mu^o(y) \,{\rm d}\sigma_y \]
can be proven to be real analytic by the properties of integral operators with real analytic kernels and no singularities (see Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]). For the third term of 
$M_2$ we proceed in this way. The map from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$(\phi ,\mu ^o)$ to the function of the variable 
$t\in \partial \Omega ^i$ defined by
\[ v^+_{\Omega^o}[\mu^o](\phi(t)) = \int_{\partial\Omega^o} S_n(\phi(t)-y) \, \mu^o(y) \,{\rm d}\sigma_y \]
can be proven to be real analytic by the properties of integral operators with real analytic kernels and no singularities (see Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]). The map from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^i)$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes 
$(\phi ,\mu ^o)$ to the function of the variable 
$t\in \partial \Omega ^i$ defined by
\[ V_{\partial\Omega^i[\phi]}[\mu^i\circ \phi^{({-}1)}](\phi(t)) = \int_{\phi(\partial\Omega^i)} S_n(\phi(t)-y) \, \mu^i\circ \phi^{({-}1)}(y) \,{\rm d}\sigma_y \]
is real analytic by a result of real analytic dependence for the single-layer potential upon perturbation of the support and of the density (see Lanza de Cristoforis and Rossi [Reference Lanza de Cristoforis and Rossi21, theorem 3.12], Lanza de Cristoforis [Reference Lanza de Cristoforis18, proposition 7]). Similarly we can treat 
$V_{\partial \Omega ^i[\phi ]}[\eta ^i\circ \phi ^{(-1)}](\phi (\cdot ))$. Hence, by the real analyticity of the composition of real analytic maps and by (5.4), we conclude that the map from 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i} \times C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^i)$ that takes a sextuple 
$(\phi ,\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ to the function
\begin{align*} & \mathcal{N}_{F_1}\bigg( v^+_{\Omega^o}[\mu^o](\phi({\cdot}))_{|\partial\Omega^i} + V_{\partial\Omega^i[\phi]}[\mu^i\circ \phi^{({-}1)}](\phi({\cdot})) +\rho^o ,\\ &\quad V_{\partial\Omega^i[\phi]}[\eta^i\circ \phi^{({-}1)}](\phi({\cdot})) +\rho^i \bigg) \end{align*}
is real analytic. As a consequence 
$M_2$ is real analytic.
It will be convenient to consider 
$F_1$, 
$F_2$ as two components of a vector field on 
$\partial \Omega ^i\times \mathbb {R}^2$. We denote by 
$F$ the function from 
$\partial \Omega ^i\times \mathbb {R}^2$ to 
$\mathbb {R}^2$ defined by
\[ F(t,\zeta_1,\zeta_2) = (F_1(t,\zeta_1,\zeta_2),F_2(t,\zeta_1,\zeta_2)) \quad \forall (t,\zeta_1,\zeta_2) \in \partial\Omega^i\times{\mathbb{R}}^2. \]
Clearly, we can extend the definition of the superposition operator (cf. § 2) in a natural way, i.e., by setting
\[ \mathcal{N}_F : (C^{0,\alpha}(\partial\Omega^i))^2 \to (C^{0,\alpha}(\partial\Omega^i))^2, \quad \mathcal{N}_F \equiv (\mathcal{N}_{F_1},\mathcal{N}_{F_2}). \]
Now let 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ be as in proposition 4.8. By standard calculus in Banach space, we have the following formula regarding the first-order differential of 
$\mathcal {N}_F$:
\[ {\rm d}\mathcal{N}_F(v^+_{\Omega^o}[\mu^o_0]_{|\partial\Omega^i} + V_{\Omega^i}[\mu^i_0] + \rho^o_0 , V_{\Omega^i}[\eta^i_0] +\rho^i_0) .(h_1,h_2) = A_{\mathcal{N}_F,0} \begin{pmatrix} h_1 \\ h_2 \end{pmatrix} \]
for all 
$(h_1,h_2) \in (C^{0,\alpha }(\partial \Omega ^i))^2$, where
\begin{equation} A_{\mathcal{N}_F,0} \equiv \begin{pmatrix} \mathcal{N}_{\partial_{\zeta_1}F_1}(\alpha^1_0 , \alpha^2_0) & \mathcal{N}_{\partial_{\zeta_2}F_1}(\alpha^1_0 , \alpha^2_0) \\ \mathcal{N}_{\partial_{\zeta_1}F_2}(\alpha^1_0 , \alpha^2_0) & \mathcal{N}_{\partial_{\zeta_2}F_2}(\alpha^1_0 , \alpha^2_0) \end{pmatrix} \end{equation}
and 
$\alpha ^1_0$ and 
$\alpha ^2_0$ are the functions from 
$\partial \Omega ^i$ to 
$\mathbb {R}$ defined by
\begin{equation} \alpha^1_0 \equiv v^+_{\Omega^o}[\mu^o_0]_{|\partial\Omega^i} + V_{\Omega^i}[\mu^i_0] + \rho^o_0, \quad \alpha^2_0 \equiv V_{\Omega^i}[\eta^i_0] +\rho^i_0. \end{equation}
We will require that the matrix 
$A_{\mathcal {N}_F,0}$ given by (5.5)–(5.6) satisfies assumption (4.8). In particular, we notice that assumption (5.4) implies the validity of the first of the three conditions of (4.8) for the matrix 
$A_{\mathcal {N}_F,0}$. In order to apply the Implicit Function Theorem (cf. Deimling [Reference Deimling9, theorem 15.3]) to equation (5.2) we need to prove the invertibility of the partial differential of 
$M$.
Proposition 5.3 Let assumptions (4.16) and (5.4) hold. Let 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ be as in proposition 4.8. Let 
$A_{\mathcal {N}_F,0}$ be as in (5.5)–(5.6) and assume that satisfies assumption (4.8). Then the partial differential of 
$M$ with respect to 
$(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)$ evaluated at the point 
$(\phi _0,\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0)$, which we denote by
\begin{equation} \partial_{(\mu^o, \mu^i, \eta^i, \rho^o,\rho^i)} M[\phi_0,\mu^o_0, \mu^i_0 ,\eta^i_0,\rho^o_0,\rho^i_0], \end{equation}
is an isomorphism from 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ to 
$C^{0,\alpha }(\partial \Omega ^o) \times (C^{0,\alpha }(\partial \Omega ^i))^2$.
Proof. By standard calculus in Banach spaces, we can verify that the partial differential (5.7) is the linear and continuous operator defined by
\begin{align*} & \partial_{(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)} M_1[\phi_0,\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0]. (\tilde{\mu}^o,\tilde{\mu}^i,\tilde{\eta}^i,\tilde{\rho}^o,\tilde{\rho}^i)(x) \\ &\quad = \left( -\frac{1}{2} I + W^\ast_{\Omega^o} \right) [\tilde{\mu}^o] (x) + \nu_{\Omega^o}(x) \cdot \nabla v^-_{\Omega^i}[\tilde{\mu}^i] (x) \quad \forall x \in \partial\Omega^o \\ & \partial_{(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)} M_2[\phi_0,\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0]. (\tilde{\mu}^o,\tilde{\mu}^i,\tilde{\eta}^i,\tilde{\rho}^o,\tilde{\rho}^i)(t) \\ &\quad = \left( \frac{1}{2} I + W^\ast_{\Omega^i} \right) [\tilde{\mu}^i] (t) + \nu_{\Omega^i}(t) \cdot \nabla v^+_{\Omega^o}[\mu^o](t) \\ & \qquad -\partial_{\zeta_1}F_1\left(t,v^+_{\Omega^o}[\mu^o_0](t) + V_{\Omega^i}[\mu^i_0](t) + \rho^o_0 , V_{\Omega^i}[\eta^i_0](t) +\rho^i_0 \right) \\ & \qquad \times \left(v^+_{\Omega^o}[\tilde{\mu}^o](t) + V_{\Omega^i}[\tilde{\mu}^i](t) + \tilde{\rho}^o\right) \\ & \qquad -\partial_{\zeta_2}F_1\left(t,v^+_{\Omega^o}[\mu^o_0](t) + V_{\Omega^i}[\mu^i_0](t) + \rho^o_0 , V_{\Omega^i}[\eta^i_0](t) +\rho^i_0 \right) \\ & \qquad \times \left( V_{\Omega^i}[\tilde{\eta}^i](t) + \tilde{\rho}^i\right) \quad \forall t \in \partial\Omega^i \\ & \partial_{(\mu^o,\mu^i,\eta^i,\rho^o,\rho^i)} M_3[\phi_0,\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0]. (\tilde{\mu}^o,\tilde{\mu}^i,\tilde{\eta}^i,\tilde{\rho}^o,\tilde{\rho}^i)(t) \\ &\quad = \left( -\frac{1}{2} I + W^\ast_{\Omega^i} \right) [\tilde{\eta}^i] (t) \\ & \qquad -\partial_{\zeta_1}F_2\left(t,v^+_{\Omega^o}[\mu^o_0](t) + V_{\Omega^i}[\mu^i_0](t) + \rho^o_0 , V_{\Omega^i}[\eta^i_0](t) +\rho^i_0 \right) \, \\ & \qquad \times \left(v^+_{\Omega^o}[\tilde{\mu}^o](t) + V_{\Omega^i}[\tilde{\mu}^i](t) + \tilde{\rho}^o\right) \\ & \qquad -\partial_{\zeta_2}F_2\left(t,v^+_{\Omega^o}[\mu^o_0](t) + V_{\Omega^i}[\mu^i_0](t) + \rho^o_0 , V_{\Omega^i}[\eta^i_0](t) +\rho^i_0 \right) \, \\ & \qquad \times \left( V_{\Omega^i}[\tilde{\eta}^i](t) + \tilde{\rho}^i\right) \quad \forall t \in \partial\Omega^i \end{align*}
for all 
$(\tilde {\mu }^o,\tilde {\mu }^i,\tilde {\eta }^i,\tilde {\rho }^o,\tilde {\rho }^i) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$. Then, by proposition 4.3 with 
$A=A_{\mathcal {N}_F,0}$ (cf. (4.11)) and since 
$A_{\mathcal {N}_F,0}$ satisfies (4.8), 
$\partial _{(\mu ^o,\mu ^i,\eta ^i,\rho ^o,\rho ^i)} M[\phi _0,\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0]$ is an isomorphism of Banach spaces.
By propositions 5.1–5.3, and by applying the Implicit Function Theorem for real analytic functions in Banach spaces (cf. Deimling [Reference Deimling9, theorem 15.3]) to equation (5.2), we deduce the following real analyticity result for the dependence of the densities in the integral representation fomula for the solutions of problem (1.6) upon the perturbation of the shape of the inclusion 
$\Omega ^i$.
Theorem 5.4 Let assumptions (4.16) and (5.4) hold. Let 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0) \in C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$ be as in proposition 4.8. Let 
$A_{\mathcal {N}_F,0}$ be as in (5.5)–(5.6) and assume that satisfies assumption (4.8). Then, there exist two open neighbourhoods 
$Q_0$ of 
$\phi _0$ in 
$\mathcal {A}^{\Omega ^o}_{\partial \Omega ^i}$ and 
$U_0$ of 
$(\mu ^o_0,\mu ^i_0,\eta ^i_0,\rho ^o_0,\rho ^i_0)$ in 
$C^{0,\alpha }(\partial \Omega ^o)_0 \times C^{0,\alpha }(\partial \Omega ^i) \times C^{0,\alpha }(\partial \Omega ^i)_0 \times \mathbb {R}^2$, and a real analytic map 
$\Lambda \equiv (M^o,M^i,N^i,R^o,R^i): Q_0 \to U_0$ such that the set of zeros of 
$M$ in 
$Q_0 \times U_0$ coincides with the graph of the function 
$\Lambda$. In particular,
\[ \Lambda[\phi_0]=(M^o[\phi_0],M^i[\phi_0],N^i[\phi_0],R^o[\phi_0],R^i[\phi_0])= (\mu^o_0,\mu^i_0,\eta^i_0,\rho^o_0,\rho^i_0). \]
We are now ready to exhibit a family of solutions of problem (1.6).
Definition 5.5 Let assumptions (4.16) and (5.4) hold. Let 
$A_{\mathcal {N}_F,0}$ be as in (5.5)–(5.6) and assume that satisfies assumption (4.8). Let 
$Q_0$ and 
$\Lambda \equiv (M^o,M^i,N^i,R^o,R^i)$ be as in theorem 5.4. Then, for each 
$\phi \in Q_0$ we set
\begin{align*}
u^o_\phi(x) &=
U^o_{\Omega^i[\phi]}[M^o[\phi],M^i[\phi]\circ
\phi^{({-}1)},N^i[\phi]\circ
\phi^{({-}1)},R^o[\phi],\\
&\quad \times R^i[\phi]](x) \quad\forall x \in
\overline{\Omega^o} {\setminus} \Omega^i[\phi], \\
u^i_\phi(x) &=
U^i_{\Omega^i[\phi]}[M^o[\phi],M^i[\phi]\circ
\phi^{({-}1)},N^i[\phi]\circ
\phi^{({-}1)},R^o[\phi],R^i[\phi]](x) \quad \forall x \in
\overline{\Omega^i[\phi]}, \end{align*}
where the pair 
$(U^o_{\Omega ^i[\phi ]}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ],U^i_{\Omega ^i[\phi ]}[\cdot ,\cdot ,\cdot ,\cdot ,\cdot ])$ is defined by (4.1).
By propositions 4.8, 5.1, and theorem 5.4, we deduce the following.
Theorem 5.6 Let assumptions (4.16) and (5.4) hold. Let 
$A_{\mathcal {N}_F,0}$ be as in (5.5)–(5.6) and assume that satisfies assumption (4.8). Let 
$Q_0$ be as in theorem 5.4 and let 
$(u^o_\phi ,u^i_\phi )$ be as in definition 5.5. Then, for all 
$\phi \in Q_0$, 
$(u^o_\phi ,u^i_\phi ) \in C^{1,\alpha }(\overline {\Omega ^o} {\setminus} \Omega ^i[\phi ]) \times C^{1,\alpha }(\overline {\Omega ^i[\phi ]})$ is a solution of problem (1.6). In particular 
$(u^o_{\phi _0},u^i_{\phi _0})= (u^o_0,u^i_0)$ is a solution of problem (1.3).
We are now ready to prove our main result, where we show that suitable restrictions of the functions 
$u^o_\phi$ and 
$u^i_\phi$ depend real analytically on the parameter 
$\phi$ which determines the domain perturbation.
Theorem 5.7 Let assumptions (4.16) and (5.4) hold. Let 
$A_{\mathcal {N}_F,0}$ be as in (5.5)–(5.6) and assume that satisfies assumption (4.8). Let 
$Q_0$ be as in theorem 5.4 and let 
$(u^o_\phi ,u^i_\phi )$ be as in definition 5.5. Then, the following statements hold.
(i) Let
$\Omega _\mathtt {int}$ be a bounded open subset of $\Omega ^o$. Let $Q_\mathtt {int} \subseteq Q_0$ be an open neighbourhood of $\phi _0$ such that
Then the map from\[ \overline{\Omega_\mathtt{int}} \subset {\Omega^i[\phi]} \quad \forall \phi \in Q_\mathtt{int}. \]
$Q_\mathtt {int}$ to $C^{1,\alpha }(\overline {\Omega _\mathtt {int}})$ that takes $\phi$ to $u^i_{\phi | \overline {\Omega _\mathtt {int}}}$ is real analytic.(ii) Let
$\Omega _\mathtt {ext}$ be a bounded open subset of $\Omega ^o$. Let $Q_\mathtt {ext} \subseteq Q_0$ be an open neighbourhood of $\phi _0$ such that
Then the map from\[ \overline{\Omega_\mathtt{ext}} \subset \Omega^o {\setminus} \overline{\Omega^i[\phi]} \quad \forall \phi \in Q_\mathtt{ext}. \]
$Q_\mathtt {ext}$ to $C^{1,\alpha }(\overline {\Omega _\mathtt {ext}})$ that takes $\phi$ to $u^o_{\phi | \overline {\Omega _\mathtt {ext}}}$ is real analytic.
Proof. We prove (i). By definition 5.5, by (4.1) and by lemma 2.1, we have
\begin{align*} u^i_\phi(x) &= U^i_{\Omega^i[\phi]}[M^o[\phi],M^i[\phi]\circ \phi^{({-}1)},N^i[\phi]\circ \phi^{({-}1)},R^o[\phi],R^i[\phi]](x)\\ & = \int_{\partial\Omega^i} S_n(x-\phi(s)) \, N^i[\phi](s) \, \tilde{\sigma}_n[\phi](s) \,{\rm d}\sigma_s + R^i[\phi] \quad \forall x \in \overline{\Omega^i[\phi]} \end{align*}
and for all 
$\phi \in Q_0$. By the assumption 
$Q_\mathtt {int} \subseteq Q_0$ and theorem 5.4, we know that the map from 
$Q_\mathtt {int}$ to 
$\mathbb {R}$ that takes 
$\phi$ to 
$R^i[\phi ]$ is real analytic. Moreover, by the real analyticity of 
$N^i[\cdot ]$ (cf. theorem 5.4) and by the properties of integral operators with real analytic kernels and no singularities (see Lanza de Cristoforis and Musolino [Reference Lanza de Cristoforis and Musolino20, proposition 4.1]), we can prove that the map from 
$Q_\mathtt {int}$ to 
$C^{1,\alpha }(\overline {\Omega _\mathtt {int}})$ that takes 
$\phi$ to the function 
$\int _{\partial \Omega ^i} S_n(x-\phi (s)) \, N^i[\phi ](s) \, \tilde {\sigma }_n[\phi ](s) \,\textrm {d}\sigma _s$ of the variable 
$x \in \overline {\Omega _\mathtt {int}}$ is real analytic (see also lemma 2.1). Hence, we deduce the validity of (i). The proof of (ii) is similar and it is left to the reader.
Acknowledgments
The authors are members of the ‘Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni’ (GNAMPA) of the ‘Istituto Nazionale di Alta Matematica’ (INdAM). R.M. acknowledges the support of the Project ‘Variational methods for stationary and evolution problems with singularities and interfaces’ (PRIN 2017) funded by the Italian Ministry of Education, University, and Research. P.M. acknowledges the support of the Project BIRD191739/19 ‘Sensitivity analysis of partial differential equations in the mathematical theory of electromagnetism’ (University of Padova), of the ‘INdAM GNAMPA Project 2020 – Analisi e ottimizzazione asintotica per autovalori in domini con piccoli buchi’, and of the grant ‘Challenges in Asymptotic and Shape Analysis – CASA’ (Ca’ Foscari University of Venice).
Open access
