Next Article in Journal
Multivariate Structural Equation Modeling Techniques for Estimating Reliability, Measurement Error, and Subscale Viability When Using Both Composite and Subscale Scores in Practice
Next Article in Special Issue
On Aspects of Continuous Approximation of Diatomic Lattice
Previous Article in Journal
Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation
Previous Article in Special Issue
Utilization of a Genetic Algorithm to Identify Optimal Geometric Shapes for a Seismic Protective Barrier
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs

by
Kirill D. Cherednichenko
1,
Yulia Yu. Ershova
2 and
Alexander V. Kiselev
1,*
1
Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
2
Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(8), 1161; https://doi.org/10.3390/math12081161
Submission received: 18 March 2024 / Revised: 4 April 2024 / Accepted: 6 April 2024 / Published: 12 April 2024
(This article belongs to the Special Issue Multiscale Mathematical Modeling)

Abstract

:
Norm-resolvent convergence with an order-sharp error estimate is established for Neumann Laplacians on thin domains in R d ,   d 2 , converging to metric graphs in the limit of vanishing thickness parameter in the “resonant” case. The vertex matching conditions of the limiting quantum graph are revealed as being closely related to those of the δ type.
MSC:
35C20; 35J25; 47A10; 47A20; 81Q35; 74K99

1. Introduction

In the book [1] (see also references therein, in particular [2,3,4]), Neumann Laplacians A ε on thin manifolds, converging to metric graphs as ε 0 , were studied; see, e.g., Figure 1. Here, ε represents the “thickness” of the manifolds in those parts where they converge to the graph edges. The named works attacked the question of spectral (and, in the case of [1], norm-resolvent) convergence of such partial differential operators to a graph Laplacian with certain matching conditions at the graph vertices. This latter ordinary differential operator is introduced as follows. Denoting by E the set of edges e of the limiting graph G , each e E can be identified with the interval [ 0 , l e ] , where l e is the length of e. We use notation L 2 ( e ) for the associated Hilbert space L 2 ( 0 , l e ) . Similarly, we denote L 2 ( G ) : = e L 2 ( e ) . The graph Laplacian A G is generated by the differential expression u , u H 2 ( e ) on each edge e E separately (see [5] for details), subject to the vertex conditions discussed below.
It was proved in [2,3,4] that, within any compact K C , the spectra of A ε converge in the Hausdorff sense to the spectrum of a graph Laplacian A G . In the book [1], the claimed convergence was enhanced to the norm-resolvent type, with an explicit control of the error as O ( ε γ ) , where γ > 0 depends on whether the ambient space is two-dimensional. The matching conditions at the vertices of the limiting graph turn out to be any one of the following:
(i)
Kirchhoff (i.e., standard) if the vertex volumes are decaying, as ε 0 , faster than the edge volumes;
(ii)
“Resonant”, which can be equivalently described in terms of δ -type matching conditions with coupling constants proportional to the spectral parameter z , see [2,3], if the vertex and edge volumes are of the same order;
(iii)
“Dirichlet-decoupled” (i.e., the graph Laplacian becomes completely decoupled) if the vertex volumes vanish slower than the edge ones.
We also refer the reader to the frequently overlooked papers [6,7], where an alternative approach to the asymptotic analysis of thin networks, aimed at capturing the resonant and scattering features, is developed; see also the references therein. In the present paper we consider the case of scalar PDEs, although our approach can be immediately generalised to the case of PDE systems, in particular, in the context of thin elastic structures with applications to, e.g., pentamodes (see [8,9]) by utilising the results of [10].
In the present paper, we are primarily interested in the most non-trivial resonant case (ii). We provide a straightforward, alternative to [1], proof of the fact that the Neumann Laplacians A ε in this case converge in a norm-resolvent sense to a linear operator acting in the Hilbert space L 2 ( G ) C N , where N is the number of vertices. The mentioned limiting operator is in fact the one first pointed out in [4] as a self-adjoint operator whose spectrum coincides with the Hausdorff limit of spectra for the family A ε .
On the technical side, our approach can be seen as a modification of the one developed by us in [10,11]; see also [12,13,14]. We specifically point out that the framework originally developed for high-contrast homogenisation admits a natural generalisation to setups where the “contrast” is achieved by purely geometric means, including (but not limited to) thin networks, as in the present paper. This seems to widen significantly the range of dimension-reduction-type models that are amenable to this kind of analysis and thus to establish a transparent connection between previously unrelated physical contexts, providing for a possibility to develop new types of media in materials science.
We obtain a better error bound than [1] (Section 6.7) (and skip a visit to the “plumber’s shop” of [1,15]). Our estimate in the planar case is O ( ε | log ε | ) and in the case of R d ,   d 3 , it is O ( ε ) . Unlike that of [1], our method does not allow us to study the full set of asymptotic regimes in (iii). This is due to the fact that the argument of our paper [11] is based on the Dirichlet-to-Neumann (DN) machinery. There is a possibility to modify the approach by invoking Neumann-to-Dirichlet maps instead, which would have two advantages: one could consider all rates of vertex volume decay in (iii), and certain geometric smoothness requirements could be somewhat relaxed. Nevertheless, in this paper we stick with the DN version of the approach in order to align the exposition with that of [11]. The alternative strategy will be followed up elsewhere, both in the present context and in the setting of [11].
The above results of course imply the Hausdorff spectral convergence, at the same time yielding a sharp estimate on its rate. Moreover, in contrast to [1], our approach allows one to consider “high-frequency” regimes, i.e., setups where the spectral parameter (which in the wave propagation context may represent the square of the frequency) is no longer constrained to a compact set but is still constrained by some negative power of the small parameter ε . (In non-dimensional terms this corresponds to the wavelength being of the order of some positive power of ε ). This leads to a sequence of “effective”, dimensionally reduced, models of the thin structure, which are sequentially applicable for a set of (asymptotic) frequency intervals. The complexity of the dimension reduction process for these models increases along the sequence. While an initial result regarding the high-frequency situation is presented below (which suffices to reveal a metamaterial in a periodic thin network of [16]), we postpone the full analysis to a future publication.
Alongside the high-frequency analysis, yet another sequence of models will be revealed by a version of the same argument. This corresponds to the transition from the resonant regime to the Dirichet-decoupled one and will allow us to reconcile the asymptotic analysis of [1] with that of [7] by introducing “transitional” models of increasing complexity. In these transitional regimes, from the point of view of [1], one gets arbitrarily close to the Dirichlet-decoupled situation, see (iii) above, whereas from the point of view of [7] (where the vertex volumes do not decay at all), one faces a highly non-trivial picture of resonant scattering.
Aiming at better clarity, in the present paper we restrict ourselves to the case where (without much loss of generality) the edge subdomains are assumed straight and uniformly thin, whereas the vertex subdomains are smooth with, possibly, the exception of the points where they meet the edges (see Section 2 for further details). These assumptions appear very natural in view of the possible examples shown in Figure 1 and Figure 2.
The ultimate section of the paper is dedicated to the spectral analysis of the effective graph Laplacian. We show that effective matching conditions at graph vertices, having been commonly treated as impedance type (namely, linear in the spectral parameter), are in fact of δ type after an application of a unitary gauge. This latter, however, introduces constant magnetic potentials on the edges of the limit graph. It is notable that the presence of such a magnetic field leads to a phase transition to a medium exhibiting double-negative metamaterial properties; see [16]. This observation will permit us to construct thin periodic networks with negative group velocity without the need for an external magnetic field, by relying instead on the resonant geometric properties of the network.

2. Problem Setup and Preliminaries

In what follows, we consider a prototypical setup only, which already presents all the challenges appearing in the general case. We also refer the reader to [1], where the most general setup is meticulously introduced.
For the limiting graph, the following notation will be used: the metric graph G will be identified with the set of edges E, so each individual edge is denoted by e E and is associated with an interval [ 0 , l e ] . We denote by V the set of graph vertices and treat each v V as the set of edge endpoints meeting at v. The graph G is assumed to be oriented throughout.
Proceeding to the setup for the Neumann Laplacian on a thin-graph-convergent structure, let a connected (open) domain Q be the union of the “vertex” part Q V , the “edge” part Q E , and the interface boundary Γ v ε between the two, where Q E will be assumed to be a finite collection of ε -thin cylinders (for d = 2 , rectangular boxes), Q E = e Q e of lengths l e . For each e, the domain Q e is assumed to be, up to a linear spatial transform (i.e., shift and rotation), defined by (for d = 2 )
Q e = x R 2 : x 1 ( 0 , l e ) , x 2 ( 0 , ε ) = ( 0 , l e ) × ( 0 , ε ) .
In dimensions d 3 , the set Q e is defined as the direct product of the interval ( 0 , l e ) and a smooth cross-section Q e c of volume ε d 1 .
It is further assumed that Q V = v Q v , where each of the disjoint domains Q v is connected and has piecewise smooth boundary Q v , which can be decomposed as Q v = Γ ˜ v ε ¯ Γ v ε ¯ , where Γ v ε is the interface between the vertex and edge domains and Γ ˜ v ε is the remaining part of Q v . Henceforth, we often drop the dependence on ε in the notation for those geometric elements where it can be removed by a rescaling or is otherwise irrelevant in the analysis to follow.
The contact part Γ v ε is further decomposed into a union of flat plates (for d = 2 , straight segments): Γ v ε = Γ e v ε . Here, the union is taken over all edge domains Q e connected to Q v so that Q ¯ e Q ¯ v = Γ e v ε ¯ . In what follows, we will refer to the segments Γ e v ε as contact plates. Since operators of Zaremba (or mixed) boundary value problem [17] will be used below, we further require that the contact plates Γ e v ε meet Γ ˜ v ε at angles strictly less than π ; see [18,19] for further details. We will further assume that the curves Γ ˜ v ε are smooth.
Furthermore, we assume that for all vertices v the domains Q v have volumes of order ε d 1 and are obtained by suitable transformations of an ε -independent domain Q v 0 containing the origin. More precisely, we assume that for all v and ε > 0 one has
Q v = ε d 1 d T v ε ( Q v 0 ) + b v ε ,
where b v ε R d , and the domain T v ε ( Q v 0 ) is obtained from Q v 0 by a suitable homeomorphism T v ε that (a) preserves volume and (b) maintains the angles between the contact plates and the adjacent parts of the boundary to be uniformly less than π , cf. [17,18].
For example, if Q v is star-shaped with respect to the origin (recall that a domain Q v is said to be star-shaped with respect to x v Q v if for all x Q v one has { x v + t ( x x v ) , t [ 0 , 1 ] } Q v ), the deformation T v ε can be constructed by “cut-and-glue surgery” applied to Q v 0 followed by a suitable “radial” scaling, as follows. The domain Q v 0 is first split into several conical (sectorial for d = 2 ) subdomains by making cuts alongside the cone “generatrices” (“radii” for d = 2 ) from the origin to the boundary points (endpoints for d = 2 ) of all subsets of Q v 0 pertaining to the contact plates of Q v . Each cone is then transformed by a suitable change of variable so that (A) the images of the cone bases pertaining to the contact plates of Q v have linear size ε 1 / d and (B) the image of the complementary part of Q v 0 is reattached to the lateral boundaries of the cones in A. Note that this procedure in general changes the domain volume by a small quantity. In order to restore the volume, a suitable “radial” scaling (centred at the origin) is applied, ensuring that it does not affect any of the mentioned conical subsets.
We remark that the representation (1) guarantees that we are in the “resonant” (or “borderline”) case of [2,3], i.e., that the volumes of the contact plates are proportional to the volumes of vertex domains.
The requirement that the homeomorphism T v ε in (1) be volume preserving is not essential and can be removed, as long as the ratios of the volumes of the vertex domains to the volumes of the contact plates are bounded below. In this case, the corresponding final convergence estimates, similar to those we derive below (see Theorem 3), will in general depend on the said ratios, as is evident from the Formula (14).
On the domain Q, we consider a family of self-adjoint operators A ε defined by their sesquilinear forms and corresponding to the differential expression Δ subject to Neumann boundary conditions. Other types of boundary conditions can be considered as well, including those of Robin and Dirichlet [6,20], in which case the lower edge of the spectrum is ε -dependent, which in the context of homogenisation corresponds to the so-called high-frequency regime or, using the terminology of M. S. Birman, to “the neighbourhood of the edge of a gap” [21,22].
Consider σ > 0 and an arbitrary non-negative function R = R ( ε ) ,   ε ( 0 , 1 ) such that R ( ε ) ε 1 / 2 | log ε | γ 0 as ε 0 , where γ = 1 for d = 2 and γ = 0 otherwise. In what follows, we will deal with the family of resolvents ( A ε z ) 1 . We shall always assume that z C is separated from the spectrum of the original operator family. In particular cases, we will assume that either z is constrained to the ε -growing compact set
K σ ε : = B R ( ε ) ( 0 ) { z C | dist ( z , R ) σ } ,
where B R ( ε ) ( 0 ) denotes the ball of radius R ( ε ) centered at the origin, or
z K σ : = z C | z K a compact set in C , dist ( z , R ) σ .
In particular, for all ε small enough, one has K σ K σ ε . After we have established the operator-norm asymptotics of ( A ε z ) 1 for z K σ , the result is extended by analyticity to a compact set K σ ext whose distance to the spectrum of the leading order of the asymptotics is bounded below by σ , and the same statement holds in the case of K σ ε .
Proceeding similarly to, e.g., [23] in the related area of critical-contrast homogenisation and facilitated by the abstract framework of [24], we consider A ε as operators of transmission problems (see [25] and references therein) relative to the internal boundary  Γ ε : = e , v Γ e v ε . The transmission problem is formulated as, given a function f L 2 ( Q ) , finding the weak solution u L 2 ( Q ) of the boundary value problem
Δ u ( x ) z u ( x ) = f ( x ) , x Q V or x Q E , u v ( x ) = u e ( x ) , u v n u e n = 0 on Γ e v ε , u n = 0 on Q .
Here, u v : = u | Q v , u e : = u | Q e for all admissible e and v (i.e., when Γ e v ε ), and n represents the exterior normal on Q and the “edge-inward" normal (i.e., directed from Q v to Q e ) on any of contact plates Γ e v ε . By a classical argument, the solution of the above problem is shown to be equal to ( A ε z ) 1 f . Note that the boundary value problem for the Neumann Laplacian on Q is given by the first and third lines in (2). While the second line is, strictly speaking, redundant, it proves important in order to view the problem as a transmission one relative to Γ ε .
It remains to be seen that the linear operator of the transmission problem (2) defined via the technique in [24], which we briefly recall below, is the same operator A ε ; the proof of this fact follows easily by combining [24] and the main estimate of [25].
Following the approach of [24] (cf. [26,27] and references therein for alternative approaches), which is based on the ideas of the classical Birman–Kreĭn–Višik theory (see [28,29,30]), the linear operator of the transmission boundary value problem is introduced as follows. Let H : = L 2 ( Γ ε ) = e , v L 2 ( Γ e v ε ) , and consider the harmonic lift operators Π V and Π E defined on ϕ H via
Π V ϕ : = u ϕ , where Δ u ϕ = 0 , u ϕ L 2 ( Q V ) , u ϕ | Γ ε = ϕ , Π E ϕ : = u ϕ , where Δ u ϕ = 0 , u ϕ L 2 ( Q E ) , u ϕ | Γ ε = ϕ ,
subject to Neumann boundary conditions on Q . These operators are first defined on ϕ C 0 2 ( Γ ε ) , in which case the corresponding solutions u ϕ can be seen as classical [17]. The results of [18] allow one to extend both harmonic lifts to bounded (in fact, compact) operators on H , in which case u ϕ are to be treated as distributional solutions of the respective boundary value problems. The solution operator Π : H L 2 ( Q ) = L 2 ( Q V ) L 2 ( Q E ) is defined as follows:
Π ϕ : = Π V ϕ Π E ϕ .
Consider the self-adjoint operator family A 0 (choosing not to reflect the ε -dependence in the notation) to be the Dirichlet decoupling of the operator family A ε , i.e., the operator of the boundary value problem on both Q V and Q E , where the Dirichlet boundary conditions are imposed on Γ ε together with Neumann boundary conditions on Q . The operator A 0 is generated by the same differential expression as A ε . Clearly, one has A 0 = A 0 V A 0 E relative to the orthogonal decomposition L 2 ( Q ) = L 2 ( Q V ) L 2 ( Q E ) ; all three operators A 0 ,   A 0 V and A 0 E are self-adjoint and positive-definite. Moreover, by [18,31] there exists a bounded inverse A 0 1 . Note that dom A 0 ran Π = ; see [24].
Furthermore, denoting by Γ ˜ 0 V and Γ ˜ 0 E the left inverse of Π V and Π E , respectively, one introduces the trace operator Γ 0 V (respectively, Γ 0 E ) as the null extension of Γ ˜ 0 V (respectively, Γ ˜ 0 E ) to the domain dom A 0 V ran Π V (respectively, dom A 0 E ran Π E ). In the same way, we introduce the operator Γ ˜ 0 and its null extension Γ 0 to the domain dom A 0 ran Π .
The solution operators S z V , S z E of the boundary value problems
Δ u ϕ z u ϕ = 0 , u ϕ dom A 0 V ran Π V , Γ 0 V u ϕ = ϕ , Δ u ϕ z u ϕ = 0 , u ϕ dom A 0 E ran Π E , Γ 0 E u ϕ = ϕ
are defined as linear mappings from ϕ to u ϕ , respectively. These operators are bounded from L 2 ( Γ ε ) to L 2 ( Q V ) and L 2 ( Q E ) , respectively, and admit the following representations:
S z E = ( 1 z ( A 0 E ) 1 ) 1 Π E , S z V = ( 1 z ( A 0 V ) 1 ) 1 Π V .
The solution operator S z from L 2 ( Γ ε ) to L 2 ( Q V ) L 2 ( Q E ) is now defined as S z = S z V S z E ; it admits the representation S z = ( 1 z ( A 0 ) 1 ) 1 Π and is bounded.
Having introduced orthogonal projections P V and P E from L 2 ( Q ) onto L 2 ( Q V ) and L 2 ( Q E ) , respectively, one has the obvious identities
S z V = P V S z , S z E = P E S z , Π V = P V Π , Π E = P E Π .
Fix self-adjoint (and, in general, unbounded) operators Λ E , Λ V defined on domains dom Λ E ,   dom Λ V L 2 ( Γ ε ) (in what follows these operators will be chosen as DN maps of Zaremba problems on Q E and Q V , respectively, and well-defined on H 1 ( Γ ε ) , where H 1 ( Γ ε ) is the standard Sobolev space pertaining to the internal boundary Γ ε ). Still following [24], we define the “second boundary operators” Γ 1 E and Γ 1 V to be linear operators on the domains
dom Γ 1 E : = dom A 0 E Π E dom Λ E , dom Γ 1 V : = dom A 0 V Π V dom Λ V .
The action of Γ 1 E ( V ) is set by:
Γ 1 E : ( A 0 E ) 1 f Π E ϕ Π E * f + Λ E ϕ , Γ 1 V : ( A 0 V ) 1 f Π V ϕ Π V * f + Λ V ϕ
for all f L 2 ( Q E ) , ϕ dom Λ E and f L 2 ( Q V ) , ϕ dom Λ V , respectively.
Alongside Γ 1 E ( V ) , introduce a self-adjoint operator Λ on dom Λ H and then the following “boundary” operator Γ 1 :
dom Γ 1 : = dom A 0 1 Π dom Λ ,
Γ 1 : A 0 1 f Π ϕ Π * f + Λ ϕ f L 2 ( Q ) , ϕ dom Λ .
We remark that the operators Γ 1 , Γ 1 E ( V ) thus defined are assumed to be neither closed nor indeed closable.
In our setup, we make the following concrete choice of the operators Λ E ( V ) : in what follows, they are the DN maps pertaining to the components Q E and Q V , respectively. More precisely, for the problem
Δ u ϕ = 0 , u ϕ L 2 ( Q E ) , u ϕ | Γ ε = ϕ , n u ϕ | Q = 0 ,
the operator Λ E maps the boundary values ϕ of u ϕ to the negative traces of its normal derivative
n u ϕ | Γ ε : = u · n | Γ ε ,
where n = n E is as above the “edge-inward” normal. This operator is well defined by its sesquilinear form as a self-adjoint operator on L 2 ( Γ ε ) (see, e.g., [32,33]), and dom Λ E H 1 ( Γ ε ) by [18]. The above definition is inspired by [24]. Note that the operator thus defined is negative the classical D N map of, e.g., [23].
On the vertex part Q V , we consider the problem
Δ u ϕ = 0 , u ϕ L 2 ( Q V ) , u ϕ | Γ ε = ϕ , n u ϕ | Q = 0 ,
and define Λ V as the operator mapping the boundary values ϕ of u ϕ to the negative traces of its normal derivative n u ϕ | Γ ε , where n = n V is again the “edge-inward” normal. The self-adjointness of Λ V on dom Λ V H 1 ( Γ ε ) follows by an unchanged argument.
Finally we introduce the operator Λ which on ϕ H 1 ( Γ ε ) is the sum Λ ϕ = Λ V ϕ + Λ E ϕ . It is also a self-adjoint operator on dom Λ H 1 ( Γ ε ) . This can be ascertained either by the argument of [25], in which case it is defined as the inverse of a compact self-adjoint operator on the orthogonal complement to constants L 2 ( Γ ε ) { c 1 } , extended to { c 1 } by zero or, alternatively, from its definition by a closed sesquilinear form.
The choice of Λ E ( V ) made above allows us to consider Γ 1 on the domain dom A 0 Π dom Λ . One then writes [24] the second Green identity in the following form:
A u , v L 2 ( Q ) u , A v L 2 ( Q ) = Γ 1 u , Γ 0 v L 2 ( Γ ε ) Γ 0 u , Γ 1 v L 2 ( Γ ε )
for all u , v dom Γ 1 = dom A 0 Π dom Λ , where the operator A is the null extension (see [34]) of the operator A 0 onto dom Γ 1 . Thus the triple ( H , Γ 0 , Γ 1 ) is closely related to a boundary quasi-triple of [26] (see also [35]) for the transmission problem considered; cf. [27] for an alternative approach.
The calculation of Π * in [24] shows that Π * = Γ 1 A 0 1 and therefore Γ 1 as introduced above acts as follows:
Γ 1 : u = P E u + P V u n P e u | Γ ε n P V u | Γ ε ,
where P E and P V are the orthogonal projections of L 2 ( Q ) onto L 2 ( Q E ) , L 2 ( Q V ) , respectively. Therefore, transmission problem at hand (at least, formally so far) corresponds to the interface (or “matching”) condition Γ 1 u = 0 .
Definition 1 
([24]). The operator-valued function M ( z ) defined on the domain dom Λ for z ρ ( A 0 ) (and in particular, for z K σ ε ) by the formula
M ( z ) ϕ = Γ 1 S z ϕ = Γ 1 ( 1 z A 0 1 ) 1 Π ϕ
is referred to as the M-function of the problem (2).
The following result of [24] summarises the properties of the M -function which we will need in what follows.
Proposition 1 
([24], Theorem 3.3). 1. One has the following representation:
M ( z ) = Λ + z Π * ( 1 z A 0 1 ) 1 Π , z ρ ( A 0 ) .
2. M ( z ) is an analytic operator-function with values in the set of closed operators in L 2 ( Γ ε ) densely defined on the z-independent domain dom Λ .
3. For z , ζ ρ ( A 0 ) the operator M ( z ) M ( ζ ) is bounded and
M ( z ) M ( ζ ) = ( z ζ ) S z ¯ * S ζ
In particular, M ( z ) = ( z ) S z ¯ * S z ¯ and ( M ( z ) ) * = M ( z ¯ ) .
4. For u z ker ( A z I ) { dom A 0 Π dom Λ } , the following formula holds:
M ( z ) Γ 0 u z = Γ 1 u z .
Alongside M ( z ) , we define M V ( z ) and M E ( z ) , which pertain to the vertex Q V and edge Q E parts of the domain Q , respectively, by the formulae
M V ( z ) ϕ = Γ 1 V S z V ϕ = Γ 1 V 1 z ( A 0 V ) 1 1 Π V ϕ , M E ( z ) ϕ = Γ 1 E S z E ϕ = Γ 1 E 1 z ( A 0 E ) 1 1 Π E ϕ .
As before, in the notation we suppress the dependence on the parameter ε for brevity.
The value of the fact that for ϕ H 1 ( Γ ε ) one has M ( z ) ϕ = M E ( z ) ϕ + M V ( z ) ϕ is clear: in contrast to A ε , which cannot be additively decomposed into “independent” terms pertaining to the vertex and edge parts of the medium Q owing to the transmission interface conditions on Γ ε , the M-function is additive (see, e.g., [11,36], where this property was observed and exploited in the related settings of homogenisation and scattering, respectively). In what follows, we will observe that the resolvent ( A ε z ) 1 can be expressed in terms of M ( z ) via a version of the celebrated Kreĭn formula, thus reducing the asymptotic analysis of the resolvent to that of the corresponding M-function (see, e.g., [26,37] for alternative approaches to derivation of the Kreĭn formula in our setting).
Alongside the transmission problem (2), the boundary conditions of which can be now (so far, formally) represented as u dom A 0 Π L 2 ( Γ ε ) , Γ 1 u = 0 , in what follows we will require a wider class of problems of this type. This class is formally given by the transmission conditions
u dom A dom A 0 Π L 2 ( Γ ε ) , β 0 Γ 0 u + β 1 Γ 1 u = 0 ,
where β 1 is a bounded operator on L 2 ( Γ ε ) and β 0 is a linear operator defined on the domain dom β 0 dom Λ .
In general, the operator β 0 Γ 0 + β 1 Γ 1 is not defined on the domain dom A . This problem is being taken care of by the following assumption, which will be satisfied throughout:
β 0 + β 1 Λ defined on dom Λ is closable in H .
We remark that by Proposition 1 the operators β 0 + β 1 M ( z ) are then closable for all z ρ ( A 0 ) , and the domains of their closures coincide with dom β 0 + β 1 Λ ¯ .
For any f H ,   ϕ dom Λ , the equality
( β 0 Γ 0 + β 1 Γ 1 ) ( A 0 1 f + Π ϕ ) = β 1 Π * f + ( β 0 + β 1 Λ ) ϕ
shows that the operator β 0 Γ 0 + β 1 Γ 1 is correctly defined on A 0 1 H Π dom Λ dom A . Denoting B : = β 0 + β 1 Λ ¯ with the domain dom B dom Λ , one checks that H B : = A 0 1 H Π dom B is a Hilbert space with respect to the norm
u B 2 : = f H 2 + ϕ H 2 + B ϕ H 2 , u = A 0 1 f + Π ϕ .
It is then proved [24] (Lemma 4.1) that β 0 Γ 0 + β 1 Γ 1 extends to a bounded operator from H B to H . For the sake of convenience, same notation β 0 Γ 0 + β 1 Γ 1 is preserved for this extension.
We will make use of the following version of the celebrated Kreĭn formula.
Proposition 2 
([24], Theorem 5.1). Let z ρ ( A 0 ) be such that the operator β 0 + β 1 M ( z ) ¯ defined on dom B is boundedly invertible. Then
R β 0 , β 1 ( z ) : = ( A 0 z ) 1 + S z Q β 0 , β 1 ( z ) S z ¯ * , w h e r e Q β 0 , β 1 : = ( β 0 + β 1 M ( z ) ¯ ) 1 β 1 ,
is the resolvent of a closed densely defined operator A β 0 , β 1 with the domain
dom A β 0 , β 1 = { u H B | ( β 0 Γ 0 + β 1 Γ 1 ) u = 0 } = ker ( β 0 Γ 0 + β 1 Γ 1 ) .
In particular, the (self-adjoint) operator of the transmission problem (2), which corresponds to the choice β 0 = 0 , β 1 = I , admits the following characterisation in terms of its resolvent:
R 0 , I ( z ) = ( A 0 z ) 1 S z M 1 ( z ) S z ¯ * .
In this case, one clearly has H B = A 0 1 H Π dom Λ and dom A 0 , I = { u H B | Γ 1 u = 0 } , which, together with the discussion at the beginning of this section, yields A 0 , I = A ε .
We remark that the operators β 0 and β 1 above can be assumed z-dependent, as this change does not impact the corresponding proofs of [24]. In this case, however, the corresponding operator-function R β 0 , β 1 ( z ) is shown to be the resolvent of a z-dependent operator family. Within the self-adjoint setup of the present paper, R β 0 , β 1 ( z ) is guaranteed to represent a generalised resolvent in the sense of [38,39,40].

3. Auxiliary Estimates

In this section, we collect a number of auxiliary statements required in our proof of the main result.
We start with the analysis of the operators Π V and S z V introduced in Section 2. First, we note that each of these operators admits a decomposition into an orthogonal sum over N vertex domains { Q v } of Q. It therefore suffices to consider a single vertex domain Q v (we recall for readers’ convenience that the volume of this domain is assumed to be decaying with ε 0 ). Its boundary Q v contains a disjoint set of straight segments belonging to the internal boundary Γ ε , which are, in line with what has been said above, denoted as Γ e v ε ; the union of the latter is Γ v ε .
The decoupled operator A 0 has L 2 ( Q v ) as its invariant subspace. We will denote by A 0 ( v ) its self-adjoint restriction, A 0 ( v ) : = A 0 | L 2 ( Q v ) . By construction, the operator A 0 ( v ) is the Laplacian with the so-called Zaremba, or mixed Neumann–Dirichlet, boundary condition [17,41]. More precisely, it is subject to the Dirichlet boundary condition on Γ v ε and to Neumann boundary condition on its complement Γ ˜ v ε . Clearly, this operator is boundedly invertible; moreover, the following statement holds.
Proposition 3 
(see [31,42]). There exists a constant C > 0 such that for all ε one has
( A 0 ( v ) ) 1 C | Q v | | log ε | , d = 2 , ε 2 d , d 3 ε | log ε | γ ,
where, as before, γ = 1 for d = 2 and γ = 0 otherwise.
Remark 1. 
The above proposition holds under more general conditions than those we impose. Namely, the domain Q v is only required to be Lipschitz and no conditions whatsoever are imposed on the geometry of the set Γ v ε .
Next, we turn our attention to the solution operator S z ( v ) : = S z V | Γ v ε and the corresponding harmonic lift Π ( v ) : = Π V | Γ v ε . The two are clearly related by the formula
S z ( v ) = 1 z ( A 0 ( v ) ) 1 1 Π ( v ) .
In order to bound the norm of Π ( v ) , we can follow, e.g., the following approach. First, consider the corresponding Zaremba problem on Q v 0 . We proceed by relating the norm of the corresponding Poisson operator to the least Steklov eigenvalue of the bi-Laplacian, following the blueprint of [43], based in turn on Fichera’s duality principle; see [44]. Since the boundary of Q v 0 is non-smooth, in doing so we follow the generalisations developed in [45,46], with obvious modifications required when passing from the Dirichlet to Zaremba setup. The estimate for the said Steklov eigenvalue is then taken from the norm of the compact embedding of H 2 ( Q v 0 ) to the traces of normal derivatives on the contact plates; see, e.g., [47]. Rescaling back to Q v , we obtain the following auxiliary result.
Lemma 1. 
There exists C > 0 such that Π ( v )     C for all ε .
By Proposition 3, the above lemma yields the following estimate for the solution operator S z ( v ) .
Lemma 2. 
For ε ( 0 , 1 ) , uniformly in z K σ ε (and, in particular, for z K σ ) one has
S z ( v ) = 1 + O | z | ε | log ε | γ Π ( v ) = Π ( v ) + O ( | z | ε | log ε | γ ) ,
where the error bounds are understood in the uniform operator norm topology.
Our next step is the analysis of the “part” of the DN map Λ V pertaining to the vertex domain Q v . We will denote by Λ v V its self-adjoint restriction Λ V | L 2 ( Γ v ε ) = : Λ v V .
First, we note that the spectrum of Λ v V (which can be termed as the Steklov spectrum of the sloshing problem pertaining to A 0 ( v ) ; see [48]) is discrete and accumulates to negative infinity. The point λ 1 = 0 is the least (by absolute value) Steklov eigenvalue with ψ v = | Γ v ε | 1 / 2 1 | Γ v ε being the corresponding eigenvector. For the second eigenvalue λ 2 , one has the following estimate; see, e.g., [32] and references therein.
Lemma 3. 
There exists C > 0 such that | λ 2 |     C ε 1 .
Introduce the N-dimensional orthogonal projection
P : = v · , ψ v ψ v ,
define P : = 1 P and consider the operator P M ( z ) P , which is well defined since P dom Λ dom Λ . By a straightforward estimate for sesquilinear forms, see [11] (Section 3.2), and taking into account (4) applied to M V ( z ) and combined with Proposition 3 as well as Lemmata 1, 3, one has the following statement.
Lemma 4. 
There exists C > 0 such that for all ε ( 0 , 1 ) one has
( P M ( z ) P ) 1 C ε , z K σ ε ,
where the operator P M ( z ) P is considered as a linear (unbounded) operator in P H .
We conclude this section by noting that, alternatively, one can derive the bound on the Poisson operator Π ( v ) claimed in Lemma 1 by employing the scaling property of the Dirichlet-to-Neumann map Λ v V , similar to the argument of [32] referenced above, combined with a standard estimate on the solutions to the classical Neumann problem.

4. Norm-Resolvent Asymptotics

We will make use of the Kreĭn Formula (7) to obtain a norm-resolvent asymptotics of the family A ε . In doing so, we will compute the asymptotics of M 1 ( z ) based on a Schur–Frobenius-type inversion formula, having first rewritten M ( z ) as a 2 × 2 operator matrix relative to the orthogonal decomposition of the Hilbert space H = P H P H . In the study of operator matrices, we rely upon the material of [49]; see also references therein.
The operator M ( z ) admits the block matrix representation
M ( z ) = A B E D with A , B , E bounded .
For the inversion of M ( z ) , we then use the Schur–Frobenius inversion formula [49] (Theorem 2.3.3)
A B E D ¯ 1 = A 1 + A 1 B ¯ S ¯ 1 E A 1 A 1 B ¯ S ¯ 1 S ¯ 1 E A 1 S ¯ 1 with S : = D E A 1 B .
Note that by Proposition 1, one has M ( z ) = ( z ) S z ¯ * S z ¯ . Moreover, since S z = ( 1 z A 0 1 ) 1 Π , one has
S z ¯ * S z ¯ = Π * ( 1 z A 0 1 ) 1 ( 1 z ¯ A 0 1 ) 1 Π ,
and therefore, for some constants c 1 , c 2 > 0 ,
S z ¯ * S z ¯ P ϕ , P ϕ H = ( 1 z ¯ A 0 1 ) 1 Π P ϕ 2 c 1 Π P ϕ 2 c 1 Π V P ϕ 2 c 2 P ϕ H 2
for all ϕ H , z K σ ε , where we have used the fact that the operator A 0 is bounded below by a positive constant. It follows that A 1 = ( P M ( z ) P ) 1 is bounded.
Proceeding exactly as in [11] based on the estimate provided by Lemma 4 which now reads
D 1 C ε ,
we use S 1 = ( I D 1 E A 1 B ) 1 D 1 to obtain S 1 = O ( ε ) .
Returning to (8), one obtains
M ( z ) 1 = A B E D 1 = A 1 0 0 0 + O ( ε )
with a uniform estimate for the remainder term. Comparing our result with (6) of Proposition 2 with β 0 : = P and β 1 : = P , one arrives at the following.
Theorem 1. 
There exists C > 0 such that for all ε ( 0 , 1 ) and z K σ ε (in partucular, for all z K σ ) one has the estimate
( A ε z ) 1 ( A β 0 , β 1 z ) 1 C | z | 2 ε
for a universal constant C and β 0 = P , β 1 = P , where the operator A β 0 , β 1 is defined in Proposition 2.
Proof. 
The proof is identical to that of [11] (Theorem 3.1); we include it here for the sake of completeness. For the resolvent ( A ε z ) 1 , the Formula (7) is applicable, in which for M ( z ) 1 we use (9). As for the resolvent A P , P z 1 , Proposition 2 with β 0 = P ,   β 1 = P is clearly applicable. Moreover, for this choice of β 0 ,   β 1 , the operator
Q P , P ( z ) = P + P M ( z ) ¯ 1 P
in (6) is easily computable (e.g., by the Schur–Frobenius inversion formula of [49], see (8)). The mentioned computation is facilitated by the fact that P + P M ( z ) is triangular ( A = P M ( z ) P ,   B = P M ( z ) P ,   E = 0 ,   D = I in (8)) with respect to the decomposition H = P H P H , yielding
Q P , P ( z ) = P P M ( z ) P 1 P ,
and the claim follows. □
Already the estimate of Theorem 1 establishes norm-resolvent convergence of the family A ε to an operator which by (10) is a relative (i.e., with respect to the difference of resolvents) finite-rank perturbation of the decoupled operator A 0 . However, in the case z K σ , which we will assume henceforth, it is possible to obtain a further simplification of this answer, relating the leading-order asymptotic term to a self-adjoint operator on the limiting metric graph. This procedure follows the blueprint of our paper [11]. We next briefly outline the related argument. For the case of z not constrained to a compact, a similar argument yields a sequence of dimensionally reduced models, as mentioned in the Introduction.
Note first that ( A 0 z ) 1 = ( A 0 V z ) 1 ( A 0 E z ) 1 is easily analysed. Indeed, by Proposition 3 one has
( A 0 V z ) 1 = O ε | log ε | γ .
Furthermore, the operator ( A 0 E z ) 1 , by separation of variables, is O ( ε 2 ) -close to the Dirichlet Laplacian on the space
H G : = v L 2 ( [ 0 , l e ] × 1 ε ) ,
where 1 e : = ε ( d 1 ) / 2 1 is the normalised constant function in the variable transverse to the edge e.
We remark that this is the only place where we use the assumption about the geometric shape of the edge parts Q e of the thin structure. This can be generalised to the setup of [1], allowing for curvature and non-uniform thickness, leading to Laplace–Beltrami operators on the edges of the limiting graph.
The operator ( A 0 z ) 1 is therefore close, uniformly in z K σ , to an operator that is unitary equivalent to the resolvent of A 0 G , where A 0 G is the Dirichlet-decoupled graph Laplacian pertaining to the graph G. The related error estimates are the same as in (11). The finite-dimensional second term on the right-hand side of (10) is therefore expected to encode the matching conditions at the vertices of the limiting graph G. In order to see this, one passes over to the generalised resolvent R ε ( z ) : = P E ( A ε z ) 1 P E , which is shown to admit the following asymptotics.
Theorem 2. 
The operator family R ε ( z ) admits the following asymptotics in the operator-norm topology for z K σ :
R ε ( z ) R eff ( z ) = O ( ε ) ,
where R eff ( z ) is the solution operator for the following spectral BVP on the edge domain Q E :
Δ u z u = f , f L 2 ( Q E ) , β 0 ( z ) Γ 0 E u + β 1 Γ 1 E u = 0 ,
with β 0 ( z ) = P P B ( z ) P , B ( z ) : = M V ( z ) and β 1 = P .
The boundary condition in (12) can be written in the more conventional form
P u | Γ ε = 0 , P n u = P B ( z ) P u | Γ ε .
Equivalently,
R ε ( z ) A P P B ( z ) P , P E z 1 = O ( ε ) ,
where A P P B ( z ) P , P E , for any fixed z , is the operator in L 2 ( Q E ) defined by Proposition 2 relative to the triple ( H , Π E , Λ E ) , where the term “triple" is understood in the sense of [24]. This operator is maximal anti-dissipative for z C + and maximal dissipative for z C ; see [40].
The proof of the theorem follows immediately from Theorem 1; see [11] (Theorem 3.6), together with the observation that
P M ( z ) P = P M E ( z ) P + P M V ( z ) P .
The next step of our argument is to introduce the truncated (reduced) boundary space H ˘ in order to make all the ingredients finite-dimensional. In view of clarity, in what follows we consistently supply the (finite-dimensional) “truncated” spaces and operators pertaining to them by the breve overscript.
We put H ˘ : = P H (noting that in our setup H ˘ is N-dimensional, where N is the number of vertices; see Section 2). Introduce the truncated Poisson operator on H ˘ by Π ˘ E : = Π E | H ˘ and the truncated DN map Λ ˘ E : = P Λ E | H ˘ . Then, the following statement holds.
Proposition 4 
([11], Theorem 3.7). 1. The formula
R eff ( z ) = A 0 E z 1 S ˘ z E M ˘ E ( z ) P B ( z ) P 1 ( S ˘ z ¯ E ) *
holds, where S ˘ z E is the solution operator of the problem
Δ u ϕ z u ϕ = 0 , u ϕ dom A 0 E ran Π ˘ E , Γ 0 E u ϕ = ϕ , ϕ H ˘ ,
and M ˘ E is the M-operator defined in accordance with (3), (5) relative to the triple ( H ˘ , Π ˘ E , Λ ˘ E ) .
2. The “effective” generalised resolvent R eff ( z ) is represented as the generalised resolvent of the problem
Δ u z u = f , f L 2 ( Q E ) , u dom A 0 E ran Π ˘ E , P n u | Γ ε = P B ( z ) P u | Γ ε .
3. The triple ( H ˘ , Γ ε ˘ 0 E , Γ ε ˘ 1 E ) is the classical boundary triple [50,51] for the operator A max defined by the differential expression Δ on the domain dom A max = dom A 0 E ran Π ˘ E . Here, Γ ε ˘ 0 E and Γ ε ˘ 1 E are defined on dom A max as the operator of the boundary trace on Γ and P n u , respectively.
We now consider the operator P B ( z ) P in (13); since B = M V by definition, we invoke the estimates derived in Section 3 to obtain
P B P = P Λ V P z P Π V * Π V P + O ( ε | log ε | γ ) = z Π ˘ V * Π ˘ V + O ( ε | log ε | γ ) ,
with a uniform estimate for the remainder term. Here, the truncated Poisson operator Π ˘ V is introduced as Π ˘ V : = Π V | H ˘ relative to the same truncated boundary space as above, H ˘ = P H . As a result, we obtain
R eff ( z ) R hom ( z ) = O ε | log ε | γ ,
with
R hom ( z ) : = ( A 0 E z ) 1 S ˘ z E M ˘ E ( z ) + z Π ˘ V * Π ˘ V 1 ( S ˘ z ¯ E ) * .
By a classical result of [40] (see also [38,39]), the operator R eff ( z ) is a generalised resolvent, so it defines a z-dependent family of closed densely defined operators in L 2 ( Q E ) , which are maximal anti-dissipative for z C + and maximal dissipative for z C . Writing the resolvent ( A ε z ) 1 in the matrix form relative to the orthogonal decomposition L 2 ( Q ) = P E L 2 ( Q ) P V L 2 ( Q ) = L 2 ( Q E ) L 2 ( Q V ) then yields the following result.
Theorem 3. 
The resolvent ( A ε z ) 1 admits the following asymptotics in the uniform operator-norm topology:
A ε z 1 = R eff ( z ) + O ε | log ε | γ ,
where the operator R eff ( z ) has the following representation relative to the decomposition L 2 ( Q E ) L 2 ( Q V ) :
R eff ( z ) = R eff ( z ) K z ¯ R eff ( z ¯ ) ( A 0 E z ¯ ) 1 * Π ˘ V * Π ˘ V K z R eff ( z ) ( A 0 E z ) 1 Π ˘ V K z K z ¯ R eff ( z ¯ ) ( A 0 E z ¯ ) 1 * Π ˘ V * .
Here, K z : = Γ 0 E | N z with N z : = ran S z E P , z C ± , and the generalised resolvent R eff ( z ) is defined by (14).
The above theorem provides us with the simplest possible leading-order term R eff ( z ) of the asymptotic expansion for ( A ε z ) 1 . However, it is not yet obvious whether it is the resolvent of some self-adjoint operator in the space L 2 ( Q E ) Π ˘ V H ˘ L 2 ( Q ) . It turns out that this is indeed so, which is seen via the following explicit construction.
Put L 2 ( G ) : = e L 2 ( 0 , l e ) , H 2 ( G ) : = e H 2 ( 0 , l e ) . For all u H 2 ( G ) , denote by u e v the limit of u e ( x ) : = u | e ( x ) at the vertex v. Let H eff : = L 2 ( G ) C N , and set
dom A hom = { ( u , β ) H eff : u H 2 ( G ) , u e v = u e v = : u v for any v and e , e incident to v , and β = κ u V } ,
where u V is the N-dimensional vector of { u v } v V and κ is the diagonal matrix
κ : = diag | Q v 0 | 1 / 2 .
The action of the operator is set by
A hom u β = u κ 1 n u | V , u β dom A hom ,
where n u | V is the N-dimensional vector { e v n u e | v } v V , i.e., the vector for which each element is represented by the sum of edge-inward normal derivatives of the function u over all the edges incident to the vertex v. We write e v if and only if the edge e is incident to the vertex v.
The main result of the present work, which is obtained by computing explicitly the resolvent of (16)–(18) and comparing it with (15) (see details of a similar computation in [13]), is formulated next.
Theorem 4. 
The resolvent ( A ε z ) 1 admits the following estimate in the uniform operator norm topology, uniform in z K σ :
A ε z 1 Θ A hom z 1 Θ * = O ε | log ε | γ ,
where Θ is a partial isometry from H eff onto L 2 ( Q ) , acting as follows:
  • For every edge e G , e = [ 0 , l e ] , it embeds u H 2 ( e ) into L 2 ( Q e ) as u ( x ) × ε ( d 1 ) / 2 1 ( y ) , where y is the variable in the direction transverse to that of x;
  • For every vertex v G , it embeds the value u v , i.e., the common value of u H 2 ( G ) at the vertex v, into L 2 ( Q v ) as ε ( d 1 ) / 2 u v 1 .

5. Analysis of Vertex Matching Conditions

In the present section, we continue our study of the operator (16)–(18) associated with an arbitrary metric graph G, with a view to analyse its spectral structure. We will show that the matching conditions at graph vertices associated with the spectral problem for the mentioned operator, albeit closely resembling δ -type conditions with coupling constants linear in the spectral parameter z , can be in fact represented (up to a unitary gauge) by δ -type matching conditions for all z 0 (while for z = 0 they coincide with the classical Kirchhoff condition). We will assume throughout that this graph contains no loops, in line with the assumptions imposed on the thin network studied above. We will further assume without loss of generality that the graph G is connected and that the matrix κ is invertible.
Since the operator A hom can be viewed as a self-adjoint out-of-space extension (see, e.g., [52] and references therein) of a symmetric differential operator on the metric graph G, it is amenable to the classical boundary triples theory; see [53,54]. We recall that for a closed and densely defined symmetric operator A on a separable Hilbert space H with domain dom A , a boundary triple is defined as follows.
Definition 2 
([55]). A triple ( K , Γ 0 , Γ 1 ) consisting of an auxiliary Hilbert space K and linear mappings Γ 0 , Γ 1 defined everywhere on dom A * is called a boundary triple for A * if the following conditions are satisfied:
(1) 
The abstract Green’s formula is valid
( A * u , v ) H ( u , A * v ) H = ( Γ 1 u , Γ 0 v ) K ( Γ 0 u , Γ 1 v ) K , u , v dom A *
(2) 
For any Y 0 , Y 1 K there exist u dom A * , such that Γ 0 u = Y 0 , Γ 1 u = Y 1 . In other words, the mapping u Γ 0 u Γ 1 u from dom A * to K K is surjective.
It can be shown (see [55]) that a boundary triple for A * exists, although it is not unique.
Definition 3. 
Let T = ( K , Γ 0 , Γ 1 ) be a boundary triple of A * . The Weyl function of A * corresponding to T and denoted by M ( z ) , z C R , is an analytic operator-function with a positive imaginary part for z C + (i.e., an operator R-function) with values in the algebra of bounded operators on K such that
M ( z ) Γ 0 u = Γ 1 u u ker ( A * z I ) .
For z C R one has ( M ( z ) ) * = ( M ( z ¯ ) ) and ( z ) ( M ( z ) ) > 0 .
A comparison with the assertion 4 of Proposition 1 shows that the Weyl function M ( z ) in the context of the boundary triples theory is intimately related to the object introduced in Definition 1. The overall setup leading to its construction is, however, different and is based on the explicit choice of the boundary operators Γ 0 and Γ 1 .
Definition 4. 
An extension A of a closed densely defined symmetric operator A is called almost solvable and is denoted by A = A B , if there exist a boundary triple ( K , Γ 0 , Γ 1 ) for A * and a bounded operator B : K K defined everywhere in K such that
u dom A B Γ 1 u = B Γ 0 u .
This definition implies that dom A B dom A * and A B is a restriction of A * to the linear set dom A B : = { u dom A * : Γ 1 u = B Γ 0 u } . In this context, the operator B plays the rôle of a parameter for the family of extensions { A B B : K K } . It can be shown (see [36] for references) that the resolvent set of A B is non-empty (i.e., A B is maximal), both A B and ( A B ) * = A B * are restrictions of A * to their domains, and A B and B are self-adjont or dissipative simultaneously.
Under the additional assumption that A is simple (or, in other words, completely non-self-adjoint), that is, it has no reducing self-adjoint “parts”, the spectrum of A B coincides, counting multiplicities, with the set of points z 0 C into which ( M ( z 0 ) B ) 1 does not admit analytic continuation. In the general case, however, the spectrum is a union of the “zeroes” of the operator-valued function M ( z ) B introduced above and the spectrum of the self-adjoint “part” of the symmetric operator A in its Wold decomposition [56].
Our immediate aim is to construct a convenient boundary triple for the operator A eff . In doing so, we rely upon the framework (in a particular case of a loop-graph with exactly one vertex) of the paper [57].
We define the symmetric operator A as follows (cf. (16)):
dom A = { ( u , β ) H eff : u H 2 ( G ) , u e v = u e v = : u v for any v and e , e incident to v , n u | V = 0 , and β = κ u V } ,
where u V is, as above, the N-dimensional vector of { u v } v V , κ is the diagonal matrix (17), and n u | V is the N-dimensional vector of { e v n u e | v } v V . The action of the operator A is set by (18). The operator thus defined is clearly symmetric in L 2 ( G ) C N ; its adjoint A * is defined by the same expression (18) on the domain
dom A * = { ( u , β ) H eff : u H 2 ( G ) , u e v = u e v = : u v for any v and e , e incident to v } .
We have the following lemma.
Lemma 5. 
Let K = C N , Γ 0 ( u , β ) : = n u | V , and Γ 1 ( u , β ) : = κ 1 β u V . The triple ( K , Γ 0 , Γ 1 ) is a boundary triple for the operator A * . The operator A eff is a self-adjoint almost solvable extension of A , corresponding to the matrix B = 0 with respect to this boundary triple.
The proof of the above lemma is obtained via integration by parts; for details, see [57].
We shall further require two ordinary differential operators on the metric graph G together with their boundary triples. Consider A max δ to be the operator generated by the negative Laplacian Δ on L 2 ( G ) , defined on the domain
dom A max δ = { u H 2 ( G ) : u v : = u e v = u e v for any v and e , e incident to v } .
In the paper [58], it is shown that this is a natural choice of a maximal operator if one seeks to consider the so-called δ -type matching conditions at the graph vertices, i.e., matching conditions of the type
u continuous at every vertex v and n u | V = η 2 u V ,
where η is a diagonal matrix. The conditions (20) reduce to Kirchhoff, or standard, matching conditions under the choice η = 0 . The natural boundary triple for A max δ is ( K , Γ 0 δ , Γ 1 δ ) , where K = C N , Γ 0 δ u = u V , and Γ 1 δ u = n u | V . The corresponding Weyl function (“M-matrix”) admits the form [58 M δ ( z ) : = { m v v δ } v , v V , where
m v v δ ( k ) = k e v cot k l e , v = v , k e v , e v ( sin k l e ) 1 , v v .
Here, k : = z is such that k 0 .
Next, we consider the magnetic Laplacian on the graph G subject to δ -type matching at the vertices. Namely, we assume that for all edges e the action of the operator is described as
d d x + i τ e 2 ,
where τ e is an edgewise-constant magnetic potential. In order to introduce δ -type matching, consider the co-normal derivatives
n τ u e | v = u e | v + i τ e u e | v i f v is the left endpoint of e , ( u e | v + i τ e u e | v ) otherwise .
We say that the magnetic Laplacian on G is subject to δ -type matching at the vertices if
n τ u | v : = n τ u e | v = n τ u e | v for any e , e v , and Σ u | V = η 2 n τ u | V v V ,
where n τ u | V and Σ u | V denote the vectors { n τ u | v } v V and { e v u e | v } v V , respectively.
By an argument similar to that of [58], see also [13], one easily checks that ( K , Γ ^ 0 , Γ ^ 1 ) is a boundary triple for the magnetic Laplacian A max δ defined on
dom A max δ = { u H 2 ( G ) : n τ u | v : = n τ u e | v = n τ u e | v for any v and e , e incident to v } ,
if K = C N , Γ ^ 0 u : = { n τ u | v } v V , and, finally, Γ ^ 1 u : = { e v u e | v } v V . The corresponding M-matrix M ^ ( z ) : = { m ^ v v } v , v V admits the form [58]
m ^ v v ( k ) = k 1 e v cot k l e , v = v , k 1 e v , e v exp ( i σ e ( v , v ) τ e l e ) ( sin k l e ) 1 , v v .
Here k = z such that k 0 and σ e ( v , v ) = 1 if e is directed from v to v , σ e ( v , v ) = 1 otherwise.
We will now fix the values of the magnetic potential as follows: τ e : = π / l e . The operator A max δ corresponding to this choice will be henceforth denoted by A ^ max . Its M-matrix M ^ relative to the triple ( K , Γ ^ 0 , Γ ^ 1 ) admits the form
m ^ v v ( k ) = k 1 e v cot k l e , v = v , k 1 e v , e v ( sin k l e ) 1 , v v ,
which coincides with (21) up to the factor z 1 . We remark that the operator A ^ max and any of its self-adjoint restrictions can be unitary transformed into a regular (non-magnetic) Laplacian on the same graph G by a standard gauge transform; this will, however, be reflected in the corresponding change of matching conditions at the vertices. Motivated by applications to electromagnetic wave propagation [16], here we prefer to proceed with the magnetic setup.
Returning to the analysis of the operator A eff , we now have the following lemma.
Lemma 6. 
Relative to the boundary triple of Lemma 5, the Weyl M-matrix of the operator A eff admits the form
M ( z ) = M δ 1 ( z ) 1 z κ 2 = 1 z ( M ^ 1 ( z ) + κ 2 ) .
Proof. 
In view of Lemma 5, consider the vector ( u , β ) dom A * such that
A * u β = z u β .
By the definition of A * , one equivalently has u = z u on G and κ 1 n u | V = z β . Abbreviating Γ 0 ( u , β ) = : α , one therefore has, in view of the definition of Γ 0 (see Lemma 5):
β = 1 z κ 1 α .
Furthermore, taking now into account the definition of the second boundary operator Γ 1 , one has
Γ 1 ( u , β ) = κ 1 β u V = 1 z κ 2 α u V .
Note that the function u, by the definition of A * , must be continuous at every v V and therefore belongs to the domain of the operator A max δ . Therefore, u ker ( A max δ z ) and thus one has Γ 1 δ u = M δ ( z ) Γ 0 δ u . By construction, one has Γ 0 δ u = u | V and Γ 1 δ u = n u | V , whence u V = M δ 1 ( z ) α . Ultimately,
M ( z ) α = Γ 1 ( u , β ) = 1 z κ 2 α M δ 1 ( z ) α ,
as claimed. □
Let E ( · ) denote the orthogonal operator spectral measure of the self-adjoint operator A eff and E ( · ) the orthogonal spectral measure of the self-adjoint operator A ^ κ 2 , where the latter is defined as the almost solvable extension of ( A max δ ) * corresponding to the parameterising operator B = κ 2 . In other words, it is the magnetic Laplacian on the graph G subject to the condition τ e = π / l e for all e E and δ -type matching conditions (23) at the graph vertices with η = κ . We have the following statement.
Theorem 5. 
For any c 0 > 0 , the operators A eff E ( c 0 , ) and A ^ κ 2 E ( c 0 , ) are unitary equivalent.
Proof. 
In the generic case when A is simple (in particular when all l e are rationally independent), the claim follows immediately from Lemma 6. Indeed, in this case no positive z can be an eigenvalue of the operator A and the same applies to the operator ( A max δ ) * . Therefore, any reducing self-adjoint “part” of either symmetric operator can only be zero. In the general case, if z 0 > 0 is an eigenvalue of A , then the corresponding eigenfunction solves u = z 0 u on G subject to κ 1 n u | V = z 0 β ,   u V = κ 1 β ,   n u | V = 0 and therefore u V = n u | V = 0 . Thus, any eigenfunction must be of the form c e sin z 0 x on each e; moreover, one must also have sin z 0 l e = 0 for all e such that c e 0 and n u | V = 0 . For each of these, the function that is edgewise transformed as c e sin z 0 x c e exp ( i τ e x ) cos z 0 x is shown to be an eigenfunction of ( A max δ ) * corresponding to the eigenvalue z 0 . The same argument applied in the opposite direction completes the proof. □
Remark 2. 
The analysis of the unitary equivalence linking A eff E ( c 0 , ) and A ^ κ 2 E ( c 0 , ) is an exciting possible development from the point of view of classical functional analysis, since it appears to be a natural graph-based generalisation of the classical Hilbert transform. This can be seen, in particular, from the explicit calculation in the case of an infinite chain graph [14].
The above theorem shows that the spectral analysis of the operator A eff in relation to its non-zero spectrum reduces to that of the magnetic graph Laplacian with (non-trivial) δ -type matching condition at the graph vertices, which could come as a surprise given that the eigenvalue problem for A eff yields, for the first component of the eigenvector ( u , β ) , the equation u = z u on G subject to n u | V = z κ 2 u | V , which on the face of it is a δ -type matching condition, albeit with coupling constants proportional to the spectral parameter z. It would seem natural, therefore, for the self-adjoint operator A eff to be a H 1 singular perturbation of the graph Laplacian with standard boundary conditions. Instead, our result shows that it is in fact a more singular H 2 perturbation. The underlying reasons of this peculiar behaviour are discussed in detail in [59], where the relationship of operators of the class considered with those in the area of zero-range potentials with internal structure, as introduced by B. S. Pavlov [60,61], is explained. We also point out that the assertion of Theorem 5 has been observed in a particular case of the cycle graph with one vertex (the loop) in [14] in the context of a high-contrast homogenisation problem on the real line.
We conclude the spectral analysis of the operator A eff by a brief discussion of its kernel and the comparison of the latter with that of the operator A ^ κ 2 . It turns out that, unlike what happens with its non-zero spectrum, the kernel ker A eff of the operator A eff is that of the Kirchhoff graph Laplacian.
Note first that ker A eff necessarily belongs to the non-simple (i.e., self-adjoint) part of A . Indeed, for ( u , β ) ker A eff it has to satisfy
u = 0 on G , n u | V = 0 , u V = κ 1 β ,
where the boundary conditions are equivalent to Γ 0 ( u , β ) = 0 , Γ 1 ( u , β ) = 0 . On the other hand, this is precisely the condition for u to be in the kernel of a graph Laplacian with Kirchhoff matching conditions. One therefore infers from [5] that dim ker A eff = 1 (in our case of connected graphs), and the elements of ker A eff are constants on G.
The kernel of  A ^ κ 2 is spanned by functions u such that
u = 0 on G and Σ u | V = κ 2 n τ u | V .
It is easily checked that a non-trivial solution to this problem could exist only if n τ u | V = 0 ,   Σ u | V = 0 or, in other words, if Γ ^ 0 u = Γ ^ 1 u = 0 . This means that, precisely as in the case of A eff , the kernel of A ^ κ 2 necessarily belongs to the non-simple (i.e., self-adjoint) part of the symmetric operator ( A max δ ) * . The question of its existence and dimension admits a simple answer in terms of the graph topology. It is clear that it is trivial in the case when G is a tree; in general, its dimension is shown to be equal to the cyclomatic number χ of the graph G. In particular, this yields unitary equivalence of A eff and A ^ κ 2 in the case where G contains exactly one cycle.
Remark 3. 
The result of the present section seems to have been overlooked in a number of now-classical papers dealing with Sturm–Liouville problems on an interval with boundary conditions depending on a spectral parameter; see, e.g., [62,63,64,65,66,67].
Conjecture 1. 
The above discussion raises the question of which definition of δ -type interaction on a graph is motivated physically, i.e., whether it is the one emerging from the analysis of thin networks as the operator A hom , see (16)–(18), or the traditional (see [5]) definition (22)–(23). At first sight, the difference between the two operators is insignificant: it is only in their kernels. However, it can happen to be of paramount importance if, e.g., one considers an ε -periodic graph with δ -type matching conditions, in which case the homogenisation procedure [11,12] will lead to drastically different outcomes for the two related setups, as it relies upon a “threshold effect” [68] in the behaviour of the least eigenvalue of the operator on the fundamental cell for small quasimomenta.

Author Contributions

Writing—original draft, K.D.C., Y.Y.E. and A.V.K.; Writing—review & editing, K.D.C., Y.Y.E. and A.V.K. All authors have read and agreed to the published version of the manuscript.

Funding

K.D.C. and Y.Y.E. are grateful for the financial support of EPSRC Grant EP/L018802/2. KDC and AVK are grateful for the financial support of EPSRC Grant EP/V013025/1. Y.Y.E. and A.V.K. are grateful to IIMAS–UNAM for the hospitality and financial support during the research visit when part of this work was carried out.

Data Availability Statement

No new data were generated or analysed during this study.

Acknowledgments

The authors acknowledge valuable comments by the referees, which have helped eliminate a number of inaccuracies.

Conflicts of Interest

There are no conflicts of interest to be declared.

References

  1. Post, O. Spectral Analysis on Graph-Like Spaces; Lecture Notes in Mathematics 2039; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  2. Exner, P.; Post, O. Convergence of spectra of graph-like thin manifolds. J. Geom. Phys. 2005, 54, 77–115. [Google Scholar] [CrossRef]
  3. Kuchment, P.; Zeng, H. Convergence of spectra of mesoscopic systems collapsing onto a graph. J. Math. Anal. Appl. 2001, 258, 671–700. [Google Scholar] [CrossRef]
  4. Kuchment, P.; Zeng, H. Asymptotics of spectra of Neumann Laplacians in thin domains. Contemp. Math. 2003, 327, 199–213. [Google Scholar]
  5. Berkolaiko, G.; Kuchment, P. Introduction to Quantum Graphs. In Mathematical Surveys and Monographs; American Mathematical Society: Providence, RI, USA, 2012. [Google Scholar]
  6. Mikhailova, A.B.; Pavlov, B.S.; Prokhorov, L.V. Intermediate Hamiltonian via Glazman’s splitting and analytic perturbation for meromorphic matrix-functions. Math. Nachr. 2007, 280, 1376–1416. [Google Scholar] [CrossRef]
  7. Pavlov, B.S. Krein formula with compensated singularities for the ND-mapping and the generalized Kirchhoff condition at the Neumann Schrödinger junction. Russ. J. Math. Phys. 2008, 15, 364–388. [Google Scholar] [CrossRef]
  8. Kadic, M.; Bückmann, T.; Stenger, N.; Thiel, M.; Wegener, M. On the practicability of pentamode mechanical metamaterials. Appl. Phys. Lett. 2012, 100, 191901. [Google Scholar] [CrossRef]
  9. Milton, G.W.; Cherkaev, A.V. Which Elasticity Tensors are Realizable? J. Eng. Mater. Technol. 1995, 117, 483. [Google Scholar] [CrossRef]
  10. Cherednichenko, K.; Kiselev, A.V.; Velčić, I.; Žubrinić, J. Effective behaviour of critical-contrast PDEs: Micro-resonances, frequency conversion, and time dispersive properties. II. arXiv 2023, arXiv:2307.01125. [Google Scholar] [CrossRef]
  11. Cherednichenko, K.D.; Ershova, Y.Y.; Kiselev, A.V. Effective behaviour of critical-contrast PDEs: Micro-resonances, frequency conversion, and time dispersive properties. I. Commun. Math. Phys. 2020, 375, 1833–1884. [Google Scholar] [CrossRef]
  12. Cherednichenko, K.D.; Ershova, Y.Y.; Kiselev, A.V. Time-dispersive behaviour as a feature of critical contrast media. SIAM J. Appl. Math. 2019, 79, 690–715. [Google Scholar] [CrossRef]
  13. Cherednichenko, K.; Ershova, Y.; Kiselev, A.; Naboko, S. Unified approach to critical-contrast homogenisation with explicit links to time-dispersive media. Trans. Mosc. Math. Soc. 2019, 80, 295–342. [Google Scholar] [CrossRef]
  14. Cherednichenko, K.D.; Kiselev, A.V. Norm-resolvent convergence of one-dimensional high-contrast periodic problems to a Kronig-Penney dipole-type model. Commun. Math. Phys. 2017, 349, 441–480. [Google Scholar] [CrossRef]
  15. Rubinstein, J.; Schatzman, M. Variational Problems on multiply connected thin strips I: Basic estimates and convergence of the Laplacian spectrum. Arch. Ration. Mech. Anal. 2001, 160, 271–308. [Google Scholar] [CrossRef]
  16. Kiselev, A.V.; Ryadovkin, K.S. Phase transition in a periodic tubular structure. arXiv 2023, arXiv:2303.00872; to appear in: SIAM J. Appl. Math. (2024). [Google Scholar]
  17. Zaremba, S. Sur un probléme mixte relatif à l’équation de Laplace (French). In Bulletin International de l’Académie des Sciences de Cracovie. Classe des Sciences Mathématiques et Naturelles, Serie A: Sciences Mathématiques; Uniwersytet Jagielloński: Krakow, Poland, 1910; pp. 313–344. [Google Scholar]
  18. Brown, R.M. The mixed problem for Laplace’s equation in a class of Lipschitz domains. Commun. Partial Diff. Eqns. 1994, 19, 1217–1233. [Google Scholar] [CrossRef]
  19. Taylor, J.L.; Ott, K.A.; Brown, R.M. The mixed problem in Lipschitz domains with general decompositions of the boundary. Trans. Am. Math. Soc. 2013, 365, 2895–2930. [Google Scholar] [CrossRef]
  20. Grieser, D. Spectra of graph neighborhoods and scattering. Proc. Lond. Math. Soc. 2008, 97, 718–752. [Google Scholar] [CrossRef]
  21. Birman, M.S. On the averaging procedure for periodic operators in a neighborhood of an edge of an internal gap. St. Petersburg Math. J. 2004, 15, 507–513. [Google Scholar] [CrossRef]
  22. Craster, R.V.; Kaplunov, J.; Pichugin, A.V. High-frequency homogenization for periodic media. Proc. R. Soc. Lond. Ser. A 2010, 466, 2341–2362. [Google Scholar] [CrossRef]
  23. Friedlander, L. On the density of states of periodic media in the large coupling limit. Commun. Part. Diff. Eq. 2002, 27, 355–380. [Google Scholar] [CrossRef]
  24. Ryzhov, V. Linear operators and operator functions associated with spectral boundary value problems. Oper. Adv. Appl. 2020, 276, 576–626. [Google Scholar]
  25. Schechter, M. A generalization of the problem of transmission. Ann. Scuola Norm. Sup. Pisa 1960, 14, 207–236. [Google Scholar]
  26. Behrndt, J.; Langer, M. Boundary value problems for elliptic partial differential operators on bounded domains. J. Func. Anal. 2007, 243, 536–565. [Google Scholar] [CrossRef]
  27. Brown, M.; Marletta, M.; Naboko, S.; Wood, I. Boundary triples and M-functions for non-selfadjoint operators with applications to elliptic PDEs and block operator matrices. J. Lond. Math. Soc. 2008, 77, 700–718. [Google Scholar] [CrossRef]
  28. Birman, M.S. On the self-adjoint extensions of positive definite operators. Math. Sb. 1956, 38, 431–450. [Google Scholar]
  29. Kreĭn, M.G. Theory of self-adjoint extensions of semibounded hermitian operators and applications II. Mat. Sb. 1947, 21, 365–404. [Google Scholar]
  30. Višik, M.I. On general boundary problems for elliptic differential equations (Russian). Trudy Moskov. Mat. ObŠc. 1952, 1, 187–246. [Google Scholar]
  31. Denzler, J. Bounds for the heat diffusion through windows of given area. J. Math. Anal. Appl. 1998, 217, 405–422. [Google Scholar] [CrossRef]
  32. Agranovich, M.S. On a mixed Poincaré-Steklov type spectral problem in a Lipschitz domain. Russ. J. Math. Phys. 2006, 13, 239–244. [Google Scholar] [CrossRef]
  33. Grubb, G. The mixed boundary value problem, Krein resolvent formulas and spectral asymptotic estimates. J. Math. Anal. Appl. 2011, 382, 339–363. [Google Scholar] [CrossRef]
  34. Ryzhov, V. Weyl-Titchmarsh function of an abstract boundary value problem, operator colligations, and linear systems with boundary control. Complex Anal. Oper. Theory 2009, 3, 289–322. [Google Scholar] [CrossRef]
  35. Behrndt, J.; Rohleder, J. Spectral analysis of selfadjoint elliptic differential operators, Dirichlet-to-Neumann maps, and abstract Weyl functions. Adv. Math. 2015, 285, 1301–1338. [Google Scholar] [CrossRef]
  36. Cherednichenko, K.D.; Kiselev, A.V.; Silva, L.O. Scattering theory for non-selfadjoint extensions of symmetric operators. Oper. Theory Adv. Appl. 2020, 276, 194–230. [Google Scholar]
  37. Abels, H.; Grubb, G.; Wood, I.G. Extension theory and Kreĭn-type resolvent formulas for non-smooth boundary value problems. J. Funct. Anal. 2014, 266, 4037–4100. [Google Scholar] [CrossRef]
  38. Neumark, M. Spectral functions of a symmetric operator (Russian). Bull. Acad. Sci. URSS. Sér. Math. 1940, 4, 277–318. [Google Scholar]
  39. Neumark, M. Positive definite operator functions on a commutative group (Russian). Bull. Acad. Sci. URSS Sér. Math. 1943, 7, 237–244. [Google Scholar]
  40. Strauss, A.V. Generalised resolvents of symmetric operators (Russian). Izv. Akad. Nauk SSSR Ser. Mat. 1954, 18, 51–86. [Google Scholar]
  41. Girouard, A.; Polterovich, I. Spectral geometry of the Steklov problem (survey article). J. Spec. Theory 2017, 7, 321–359. [Google Scholar] [CrossRef] [PubMed]
  42. Denzler, J. Windows of given area with minimal heat diffusion. Trans. Am. Math. Soc. 1999, 351, 569–580. [Google Scholar] [CrossRef]
  43. Kuttler, J.R. Remarks on a Stekloff eigenvalue problem. SIAM J. Numer. Anal. 1972, 9, 1–5. [Google Scholar] [CrossRef]
  44. Fichera, G. Su un principio di dualità per talune formole di maggiorazione relative alle equazioni differenziali (Italian). Atti Accad. Naz. Lincei Rend. Cl. Sci. Fis. Mat. Nat. 1956, 19, 411–418. [Google Scholar]
  45. Auchmuty, G. The S.V.D. of the Poisson kernel. J. Fourier Anal. Appl. 2017, 23, 1517–1536. [Google Scholar] [CrossRef]
  46. Bucur, D.; Ferrero, A.; Gazzola, F. On the first eigenvalue of a fourth order Steklov problem. Calc. Var. 2009, 35, 103–131. [Google Scholar] [CrossRef]
  47. Ferrero, A.; Gazzola, F.; Weth, T. On a fourth order Steklov eigenvalue problem. Analysis 2005, 25, 315–332. [Google Scholar] [CrossRef]
  48. Levitin, M.; Parnovski, L.; Polterovich, I.; Scher, D. Sloshing, Steklov and corners: Asymptotics of sloshing eigenvalues. J. Anal. Math. 2022, 146, 65–125. [Google Scholar] [CrossRef]
  49. Tretter, C. Spectral Theory of Block Operator Matrices and Applications; Imperial College Press: London, UK, 2008. [Google Scholar]
  50. Derkach, V.A.; Malamud, M.M. Generalised resolvents and the boundary value problems for Hermitian operators with gaps. J. Funct. Anal. 1991, 95, 1–95. [Google Scholar] [CrossRef]
  51. Gorbachuk, V.I.; Gorbachuk, M.L. Boundary Value Problems for Operator Differential Equations; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1991. [Google Scholar]
  52. Štraus, A.V. Functional models and generalized spectral functions of symmetric operators. St. Petersburg Math. J. 1999, 10, 733–784. [Google Scholar]
  53. Behrndt, J.; Hassi, S.; de Snoo, H. Boundary Value Problems, Weyl Functions, and Differential Operators; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  54. Schmüdgen, K. Unbounded Self-Adjoint Operators on Hilbert Space; Graduate Texts in Mathematics 265; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  55. Kočubeĭ, A.N. Extensions of symmetric operators and of symmetric binary relations. Math. Notes 1975, 17, 25–28. [Google Scholar] [CrossRef]
  56. Gorbachuk, M.I.; Gorbachuk, V.I.M.G. Krein’s Lectures on Entire Operators; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  57. Cherednichenko, K.; Ershova, Y.; Naboko, S. Functional model for generalised resolvents and its application to time-dispersive media. arXiv 2021, arXiv:2111.05387. [Google Scholar]
  58. Ershova, Y.; Karpenko, I.I.; Kiselev, A.V. Isospectrality for graph Laplacians under the change of coupling at graph vertices. J. Spectr. Theory 2016, 6, 43–66. [Google Scholar] [CrossRef]
  59. Cherednichenko, K.D.; Ershova, Y.Y.; Kiselev, A.V.; Ryzhov, V.A.; Silva, L.O. Asymptotic analysis of operator families and applications to resonant media. Oper. Theory Adv. Appl. 2023, 291, 239–311. [Google Scholar]
  60. Pavlov, B.S. A model of zero-radius potential with internal structure (Russian). Teoret. Mat. Fiz. 1984, 59, 345–353. [Google Scholar]
  61. Pavlov, B.S. The theory of extensions, and explicitly soluble models. Russ. Math. Surv. 1987, 42, 127–168. [Google Scholar] [CrossRef]
  62. Evans, W.D. A non-selfadjoint differential operator in Lσ2[a,b). Quart. J. Math. 1970, 21, 371–383. [Google Scholar] [CrossRef]
  63. Fulton, C.T. Two-point boundary value problems with eigenvalue parameter contained in the boundary conditions. Proc. R. Soc. Edinb A. 1977, 77, 293–308. [Google Scholar] [CrossRef]
  64. Hinton, D.B. An expansion theorem for an eigenvalue problem in the boundary conditions. Quart. J. Math. 1977, 30, 33–42. [Google Scholar] [CrossRef]
  65. Schneider, A. A note on eigenvalue problems with eigenvalue parameter in the boundary conditions. Math. Z. 1974, 136, 163–167. [Google Scholar] [CrossRef]
  66. Shkalikov, A.A. Boundary problems for ordinary differential equations with parameter in the boundary conditions. J. Soviet. Math. 1983, 33, 1311–1342. [Google Scholar] [CrossRef]
  67. Walter, J. Regular eigenvalue problems with eigenvalue parameter in the boundary condition. Math. Z. 1973, 133, 301–312. [Google Scholar] [CrossRef]
  68. Birman, M.S.; Suslina, T.A. Second order periodic differential operators. Threshold properties and homogenisation. St. Petersburg Math. J. 2004, 15, 639–714. [Google Scholar] [CrossRef]
Figure 1. An example of a thin network in the case of a resonant scaling of the edge thickness and vertex diameter. The (non-dimensional) parameter ε is the ratio of the actual size to the wavelength. As ε 0 , the structure converges to a graph; the related notation is introduced at the start of Section 2.
Figure 1. An example of a thin network in the case of a resonant scaling of the edge thickness and vertex diameter. The (non-dimensional) parameter ε is the ratio of the actual size to the wavelength. As ε 0 , the structure converges to a graph; the related notation is introduced at the start of Section 2.
Mathematics 12 01161 g001
Figure 2. An example of a vertex subdomain Q v and an edge subdomain Q e separated by a contact plate Γ e v ε of linear size ε . The volumes of Q v and Q e are assumed to be of the same order ε d 1 , which is also the area of the contact plate. As in Figure 1, the (non-dimensional) parameter ε is the ratio of the actual size to the wavelength.
Figure 2. An example of a vertex subdomain Q v and an edge subdomain Q e separated by a contact plate Γ e v ε of linear size ε . The volumes of Q v and Q e are assumed to be of the same order ε d 1 , which is also the area of the contact plate. As in Figure 1, the (non-dimensional) parameter ε is the ratio of the actual size to the wavelength.
Mathematics 12 01161 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cherednichenko, K.D.; Ershova, Y.Y.; Kiselev, A.V. Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs. Mathematics 2024, 12, 1161. https://doi.org/10.3390/math12081161

AMA Style

Cherednichenko KD, Ershova YY, Kiselev AV. Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs. Mathematics. 2024; 12(8):1161. https://doi.org/10.3390/math12081161

Chicago/Turabian Style

Cherednichenko, Kirill D., Yulia Yu. Ershova, and Alexander V. Kiselev. 2024. "Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs" Mathematics 12, no. 8: 1161. https://doi.org/10.3390/math12081161

APA Style

Cherednichenko, K. D., Ershova, Y. Y., & Kiselev, A. V. (2024). Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs. Mathematics, 12(8), 1161. https://doi.org/10.3390/math12081161

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop