Tunneling in soft waveguides: closing a book

We investigate the spectrum of a soft quantum waveguide in two dimensions of the generalized ‘bookcover’ shape, that is, Schrödinger operator with the potential in the form of a ditch consisting of a finite curved part and straight asymptotes which are parallel or almost parallel pointing in the same direction. We show how the eigenvalues accumulate when the angle between the asymptotes tends to zero. In case of parallel asymptotes the existence of a discrete spectrum depends on the ditch profile. We prove that it is absent in the weak-coupling case, on the other hand, it exists provided the transverse potential is strong enough. We also present a numerical example in which the critical strength can be assessed.


Introduction
Quantum waveguides have been a subject of intense investigation in the previous three or four decades; for a survey (except the last few years) and an extensive bibliography we refer to [EK15].They serve as models of numerous physical systems, and at the arXiv:2307.01536v1[math.SP] 4 Jul 2023 same time, they represent a source of interesting mathematical problems.The nature of the interaction responsible for the guiding may differ.Often the walls of the waveguide are hard, meaning mathematically that the Hamiltonian is Dirichlet Laplacian in an appropriate spatial region.Alternatively singular Schrödinger operators are used with an attractive interaction supported by curves, surfaces, or more complicated sets of lower dimensionality in the configuration space.
Recently an alternative, more realistic model of 'soft' waveguides attracted attention, cf.[Ex20] and further developments in [KKK21, EL21, EKL22, EV23], where the interaction term is a regular potential 'ditch' the axis of which is a fixed curve.In such a description the tunneling between different parts of the guide is not suppressed while its transverse width need not be zero.The mentioned work dealt with the twodimensional situation; there are extensions to three dimensions describing soft quantum tubes [Ex22] and layers [EKP20, KK23], or even to situations where the ditch is replaced by an array of potential wells [Ex23].
The most remarkable feature of these systems are relations between their geometry and the energy spectrum, in particular, the fact that bending gives rise to an effective attractive interaction able to produce localized states in an infinitely extended ditch.This property is shared with the other waveguide types mentioned above, however, in contrast to them, the existence results obtained so far lack the universal character of their 'hard' or 'singular' counterparts.For two-dimensional Schrödinger operator with a δ interaction on a non-straight curve, for instance, Birman-Schwinger method yields the discrete spectrum existence under rather weak regularity and asymptotic straightness requirements [EI01], while the sufficient condition obtained by the same method in the soft case [Ex20] is much weaker.
As an alternative, one may apply the variational method to the operator directly; the difficulty with this approach is that in general it is not easy to construct a suitable trial function.In [KKK21] a simple case, usually dubbed bookcover in the literature [SM90], was treated; the existence of bound states was proved for any bending angle θ ∈ (0, π) between the straight parts of the guide and any channel profile such that the transverse part of the operator has a discrete eigenvalue.Two questions have been left open, both related to the tunneling between the guide arms, namely the accumulation rate of the eigenvalues as the 'book' closes, θ → π−, and the existence of the discrete spectrum in the limiting case, θ = π.
A considerably more general existence result was obtained by variational method in [EV23] allowing the finite non-straight part to be built over a C 3 smooth curve, while the semi-infinite asymptotic parts were supposed to be straight as before.The above two questions remained open again, and goal of this paper is to address them in the less restrictive geometrical setting of [EV23] where the 'spine' of the book may have an arbitrary shape provided it is finite and sufficiently regular.As in [EL21] we characterize the ditch profile by a measure type potential which allows us to treat on the same footing regular and singular interactions; recall that for the latter the said two questions are also open.
In the next section we collect the necessary notions and hypotheses, we also find the essential spectrum which differs in the parallel and non-parallel case due to the different character of the tunneling in these two situations.Section 3 is devoted to 'closing' the book.Using the angle β = π − θ as a more natural parameter in this case, we show the curvature-induced eigenvalues fill the difference between the two essential spectra in the limit β → 0, their number in any subinterval behaving as ∼ β −1 .
After that we turn to the situation when the asymptotes are parallel, β = 0.In Section 4 we consider the weak-coupling case in which the channel profile is shallow; we show that the discrete spectrum is then empty.The opposite extreme is the topic of Section 5.The strong-coupling asymptotics offer more than one regime and we discuss three of them, a singular δ potential, a rectangular well of a non-flat bottom, and a sufficiently deep and narrow profile; in all of them we prove that strong coupling gives rise to discrete eigenvalues.As a byproduct, we generalize in Sec.5.2 the well-known result about the discrete spectrum existence [EK15, Thm.1.1] to a class of Schrödinger operators in curved Dirichlet strips.As finding the critical interaction strength for U-shaped soft waveguides and their generalizations is out of reach analytically, we present in Section 6 a numerical example of an U-shaped channel with a specific transversal profile to illustrate bound states.Furthermore, we present computations regarding the dependence of the critical interaction strength on the geometry of the waveguide.

Problem setting
The object of our interest are two-dimensional Schrödinger operators with a potential in the shape of a 'ditch' of a finite width and fixed profile.Let us describe the setting in more precise terms starting from the channel axis Γ.We suppose it is an infinite and smooth planar curve without self-intersections, naturally parametrized by its arc length s, that is, the graph of a C 1 function Γ : R → R 2 such that Γ(s) = 1, where the dot conventionally denotes the derivative with respect to s.With an abuse of notation we employ the same symbol for the map Γ and for its range.If Γ is C 2 , its signed curvature is γ : γ(s) = ( Γ2 Γ1 − Γ1 Γ2 )(s) being naturally independent of the Cartesian coordinates used.The curve Γ is supposed to satisfy the following assumption: (a) Γ is piecewise C 2 -smooth and its curvature is of compact support.assa Without loss of generality we may suppose that the straight parts of Γ are the halflines Γ ± := x, ±(ρ + x tan( 1 2 β) : x ≥ 0 for a positive ρ and β ∈ [0, 1 2 π), symmetric with respect to the x axis, parametrized as for some s 0 > 2ρ.The middle part of the curve referring to s ∈ (−s 0 , s 0 ) will be denoted as Γ int ; without loss of generality again, we may suppose that it is contained in the open left halfplane, x < 0. Recall that θ(s 2 , s 1 ) := s 2 s 1 γ(s) ds is the angle between the tangents at the points Γ(s j ), j = 1, 2, hence (2.2) ?angle?
we will be interested in the situations when β = 0 and in the asymptotic regime β → 0+.
Let us turn to the potential.We suppose that it is restricted transversally being a subset of the strip of a halfwidth a > 0 built over Γ.We suppose that a is small enough so that we can parametrize Ω a using the parallel (often called Fermi) coordinates, given by the arc length s and the distance u along the normal N (s) = (− Γ2 (s), Γ1 (s)) to Γ at the point s, in such a way that (b) the map R × I a → Ω a defined by (2.4) is a diffeomorphism, assb where I a := (−a, a); note that it can be satisfied, in particular, only if a < ρ.
The channel profile will be determined by the one-dimensional Schrödinger operator The potential v is always supposed to be nontrivial ; when the interaction strength will be of importance we will replace h v by h λv with λ > 0. Mostly we will not strive for the maximum generality focusing primarily on two situations; in the first assc needless to say one then has D(h v ) = H 2 (R).Moreover, we have σ ess (h v ) = [0, ∞) and v := inf σ(h v ) < 0 is a simple eigenvalue.Using then the map (2.4) we define The second case of interest is the leaky curve when instead of (2.6b) we consider the operator assd Alternatively one can say that the self-adjoint operator H Γ,α acts as the negative Laplacian on functions ψ ∈ H 1 (R 2 ) ∩ H 2 loc (R 2 \ Γ) which have at points of Γ the jump of the normal derivative equal to −αψ(x); formally one can write H Γ,α = −∆ − αδ(x − Γ).It may be regarded as a singular version of the previous case, namely the family H Γ,Vε corresponding to the scaled transverse potentials v ε BEHL17].The analogue of the operator (2.5), formally − d 2 dx 2 − αδ(x), will be denoted as h α ; its ground state is in this case known explicitly, α := inf σ(h α ) = − 1 4 α 2 .While the two situations are of primary interest to us, they are particular cases of a more general channel profile.In the spirit of [EL21] one can consider the operator (e) H Γ,µ associated with the quadratic form , where dµ(x) = (1 − uγ(s)) dsdµ ⊥ (u) corresponds to a finite positive Borel measure µ ⊥ supported on the interval (−a, a); asse the assumptions (c) and (d) correspond respectively to dµ ⊥ (u) = v(u)du and dµ ⊥ (u) = αδ, where δ is the Dirac measure supported at the origin.The analogue of (2.5), the selfadjoint operator associated with the form In addition to (2.5) we need also its double-well counterpart and its analogues h α,ρ and h µ ⊥ ,ρ corresponding to assumptions (d) and (e), respectively.The essential spectrum of these operators is, of course, R + and their ground-state eigenvalues will be denoted as v,ρ and analogously in the other two situations.
rop:spectess Proposition 2.1.Assume (a), (b) and (e), then Proof.Under the assumptions (c) or (d) the essential spectrum was determined in [Ex20, Prop.3.1] and [EI01, Prop.5.1] for β > 0, the case β = 0 in the bookcover situation was dealt with in [KKK21, Prop.2].The more general waveguide satisfying the present assumptions can be treated similarly; the core of the argument is the existence of an increasing family of rectangles on which the potential channel is straight (or the parallel potential channels are straight for β = 0) so that one can separate the variables and construct a suitable Weyl sequence combining transversally the groundstate eigenfunction of h µ ⊥ or h µ ⊥ ,ρ , restricted to intervals of a growing length, with a mollifier in the longitudinal direction.The interaction term in h µ ⊥ ,ρ is mirror-symmetric, hence its ground state coincides with that of the single well and the Neumann boundary condition at the distance ρ from its 'center'.By the bracketing argument [RS78, Sec.XIII.15] the function ρ → µ ⊥ ,ρ is non-decreasing.In fact, it is increasing with the indicated limit; one way to see that is to use the formula , where ϕ µ ⊥ ,ρ is the real-valued eigenfunction associated with µ ⊥ ,ρ , obtained by an easy modification of the argument used in [DH93] (an alternative is to employ an Agmon-type estimate which shows that the convergence is even exponentially fast).

Closing the book
s:closing Let us now investigate what happens if we are 'closing the book', in other words, how the spectrum of H Γ,µ behaves asymptotically in the limit β → 0; our aim is to show that the eigenvalues will then fill the gap between the corresponding essential spectrum thresholds, µ ⊥ ,ρ and µ ⊥ .Specifically, we are going to show that for any ν ∈ ( µ ⊥ ,ρ , µ ⊥ ) there are eigenvalues of H Γ,µ in the interval ( µ ⊥ ,ρ , ν) for β small enough and to estimate their accumulation rate.
Proof.It is sufficient to find the appropriate number of linearly independent test functions for which the quadratic form associated with H Γ,µ satisfies q[φ] < ν φ 2 .We will construct such functions with the support in the part of the plane where Ω a is straight, choosing them even with the respect to the axis of the angle between the two asymptotes, φ(x, y) = φ(x, −y).In the upper halfplane we choose a pair of points of Γ with the coordinates s 1 , which can without loss of generality refer to the place where the straight parts begins, and s 2 > s 1 .The support of φ for y ≥ 0 will be a semiinfinite strip Σ the boundaries of which will be lines normal to Γ at the points s 1 , s 2 ; we divide it into the rectagular part Σ r and the wedge-shaped Σ w as sketched in Fig. 1.
The regions used in the proof of Theorem 3.1 It is obvious that in Σ one can use the parallel coordinates; those referring to Cartesian (x, y) will be denoted as s(x, y) and u(x, y).With this licence we use the following test function Ansatz, where ϕ ρ β is the eigenfunction associated with the eigenvalue µ ⊥ ,ρ β of operator (2.7) (for simplicity we drop the subscript µ ⊥ ), ρ β = ρ sec β 2 , and (3.2) eq:DiscSpectrumDef the functions (3.1) are obviously admissible test functions belonging to H 1 (R 2 ).
Let us now express the quadratic form starting from its kinetic term, ∇φ 2 .Using the symmetry of φ and the fact that the parallel coordinates in Σ are obtained by rotation of the Cartesian ones, we get where L := s 2 − s 1 , the norms refer to L 2 (s 1 , s 2 ) and L 2 (R), respectively, and when integrating over u we used the fact that ϕ ρ β is an even function.For the interaction part of the quadratic form we get similarly and the norm of the test function is Putting now (3.3)-(3.5)together, we can express the shifted quadratic form as where we have used the fact that ϕ ρ β is by assumption the ground-state eigenfunction of h µ ⊥ ,ρ β .The last expression is negative provided where we set η ρ β := for the sake of brevity; this quantity is bounded as a function of the angle β and lim β→0 η ρ β = η ρ > 0. The right-hand side of inequality (3.6), which we denote as R ν (L, β), is continuous with respect to the parameters and for small enough β it is positive tending to ν − µ ⊥ ,ρ as β → 0. The left-hand side is in view of (3.2) nothing else than the normalized quadratic form of the Dirichlet Laplacian on an interval of length L, hence the maximum number n ν of mutually orthogonal trial functions making the form q[•] − ν • 2 negative has to satisfy the inequality πnν L 2 < R ν (L, β).The indicated properties of the right-hand side imply, in particular, that there are positive a ν , b ν such that R ν (L, β) > a ν holds whenever L tan β 2 < b ν , and consequently, we get for n ν the bound which is the result we have set out to prove.

Parallel asymptotes, weak coupling
s:weak Let us now turn to the situation when the 'book is closed', β = 0. We know from Proposition 2.1 how the essential spectrum looks in this case; what we want to know is whether the discrete spectrum is nonempty.It appear that the answer depends on the strength of the transversal coupling.To this aim we introduce a coupling constant λ > 0 into the picture; modifying assumption (e) we consider operator H Γ,λµ associated with the quadratic form ).First we focus on the situation when the coupling is weak.hm:weakcoupl Theorem 4.1.Assume (a), (b) and (e).Then σ disc (H Γ,λµ ) is empty for all λ small enough.
We will need the following claim extending the well-known result of [Si76]: eak1Dmeasure Lemma 4.1.Let h λη be the Schrödinger operator in L 2 (R) associated with the quadratic , where η is a compactly supported finite positive Borel measure.For all sufficiently small positive λ, the operator has exactly one eigenvalue (λ) which behaves asymptotically as (4.1) ?weak1D?
Proof.As in the regular case, on can use the Birman-Schwinger principle: by [BEK Š94, Lemma 2.3] we have to check that the operator λR ηη (κ), the trace of the resolvent − d 2 dx 2 + κ 2 −1 at the space L 2 (R, dη), has for small λ a simple eigenvalue one.The trace acts on f ∈ L 2 (R, dη) as R ηη (κ)f = G κ f η-almost everywhere with G κ = 1 2κ e −κ|•| defines the kernel of the one-dimensional Laplacian resolvent.As in [Si76] we split the operator into two parts, R(κ) = L κ + M κ , where L κ is the integral operator with the kernel 1 2κ whose range is one-dimensional consisting of constant functions, and M κ with the kernel 1 2κ e −κ|x−y| − 1 .The last opertor is obviously bounded which means that I − λM κ is invertible for λ small enough; using then the identity we see that the task is reduced to finding Ker I − (I − λM κ ) −1 λL κ ηη , and since we have noted that the operator is one-dimensional, the explicit form of L κ yields the equation The rest of the proof proceeds as in the regular case; to complete the proof one has to check that the right-hand side is a real analytic function of λ up to the point λ = 0 and to expand it there in the leading order.
Proof of Theorem 4.1.We use again a bracketing argument and impose additional Neumann boundary conditions on the y-axis obtaining thus the inequality where the first operator in the direct sum refers to the curved part in the left halfplane and the second to the parallel channels in the right one.It is thus sufficient to check that none of the two operators has spectral points below λµ ⊥ ,ρ .The spectrum of H a Γ,λµ is purely absolutely continuous covering the interval [ λµ ⊥ ,ρ , ∞), and choosing for η in Lemma 4.1 the 'double-well measure', dη(u 3) eq:right_est Consider next the operator H c Γ,λµ and compare it with H Γ,λµ where Γ a closed curve being a union of Γ int and its mirror image with respect to the y-axis.By construction Γ is C 1 -smooth, and since it is a finite loop, the negative spectrum of H Γ,λµ is discrete.In view of the symmetry, the operator allows for a parity decomposition in the x direction, H Γ,λµ = H sym Γ,λµ ⊕ H asym Γ,λµ and it holds H sym Γ,λµ ≤ H asym Γ,λµ , in particular, inf σ(H Γ,λµ ) = inf σ(H sym Γ,λµ ) = inf σ(H c Γ,λµ ).By assumption, the measure µ ⊥ is positive, and the same is, of course, true for µ.This allows us to use the result of [KL14] by which H Γ,λµ has for all λ small enough a single eigenvalues which behaves asymptotically as follows, as λ → 0, (4.4) eq:loop_est where C µ is a positive constant depending on the measure µ.A comparison of (4.3) and (4.4), taking into account that µ ⊥ (R) > 0, implies in view of (4.2) that for sufficiently small λ the spectrum of H Γ,λµ below λµ ⊥ ,ρ is empty; this concludes the proof.
sign_change ?Remark 4.1.The assumption (e) includes positivity of the measure.This may not be necessary, for instance, a regular potential of assumption (c) might be signchanging.The claim of Theorem 4.1 remains nevertheless valid in such a situation even if R v(x) dx = 0. Indeed, we know from [Si76] that the weakly coupled state then still exists and it satisfies − (λ) = 1 4 λ 2 R 2 v(x)|x−y|v(y) dx dy+O(λ 3 ), hence the essential spectrum threshold behaves as O(λ 4 ).We cannot apply the result of [KL14] directly to get (4.4) because the asymptotics was obtained under the positivity assumption -see, however, Remark 3.5 of the paper -it is easy to bypass this limitation.If µ ⊥ is signchanging, we can write it as a difference of two positive measures with disjoint supports, ⊥ , and the same decomposition applies to µ.Since H Γ,λµ ≥ H Γ,λµ (+) , the eigenvalue 0 (λ) is bounded from below by the right-hand side of (4.4) with µ(R 2 ) replaced by µ (+) (R 2 ) tending exponentially fast to zero, so for small enough λ it is again the essential spectrum threshold which dominates.

Parallel asymptotes, strong coupling
s:strong In contrast to the weak coupling, the strong one offers a wider variety of asymptotic regimes referring to particular subclasses of potentials.In this section we focus on three of them.
Proof.The claim follows from the strong-coupling asymptotics of σ disc (H Γ,α ) derived in [EY01], see also [EK15,Cor. 10.3.1].It was stated there for curves the asymptotes of which were not parallel, however, this assumption was not used in its proof; the vital assumptions are satisfied since Γ ∈ C 4 (R) and its curvature is compactly supported.The jth eigenvalues of H Γ,α then behaves asymptotically as (5.1) leakyasympt where j < 0 is the jth eigenvalue of the comparison operator S Γ .The curvature-induced potential of S Γ is nonzero, attractive, bounded and compactly supported, hence σ disc (S Γ ) is nonempty and finite.

An interlude: Dirichlet guides
ss:dirichlet To consider another strong-coupling asymptotic, we shall extend a known result about the discrete spectrum of curved Dirichlet strips [EK15, Thm.1.1].Comparing the assumption concerning the Dirichlet strips in that book with those about the potential support used here, we see that (i) of [EK15, Sec.1.1] guaranteeing the existence of parallel coordinates coincides with (b), while (a) is weaker than (ii) because it requires the curvature only to be piecewise continuous.For the needs of the present paper it would be sufficient to have the curvature compactly supported which would imply assumptions (iii)-(v) but we will not limit ourselves here to this case and prove the result for strips Ω a which are only asymptotically straight.
The object to consider is a Dirichlet waveguide with a generally non-flat bottom, that is, the operator H D,v := −∆ Ω a D − V (x), where V (x) is the potential (2.6a) corresponding to v satisfying assumption (c).Let h D,v be the operator on L 2 (I a ) with the domain H 2 (I a ) ∩ H 1 0 (I a ) acting as (2.5), in other words, the transverse Schrödinger operator on I a with Dirichlet boundary conditions at ±a.The spectrum of this operator is simple and purely discrete; we denote D,v := inf σ(h D,v ) and χ 0 will be the eigenfunction corresponding to this eigenvalue.ichletGuides Theorem 5.1.Let Ω a satisfy assumptions (a) and (b), and v satisfy assumption (c).
Proof.The argument follows the proof of the mentioned Theorem 1.1 in [EK15] except that we weaken the regularity requirements on Γ in the spirit of [K Š12].The main tool is the quadratic form of the operator H D,v − D,v I which can be, using the parallel coordinates, written as (5.2) eq:HQWExtendedQuad for any ψ ∈ H 1 (R × I a ); we recall that g(s, u) = (1 − uγ(s)).To find the essential spectrum threshold, we use Neumann bracketing, adding Neumann condition at the segments perpendicular to Γ at s = ±s 0 .This allows us to estimates H D,v from below by the direct sum The middle part refers to a compact region being thus irrelevant from the point of the essential spectrum.The shifted quadratic forms of the other two parts can be estimated from below as follows, , where the norms refer to L 2 ((s 0 , ∞) × I a ) and L 2 ((−∞, −s 0 ) × I a ), respectively, and γ ± are the restrictions of γ to the appropriate intervals, and we have used the inequality The parts in the square bracket are nothing but quadratic forms of the operators on the respective parts of the straight channel and as such they they bounded from below by D,v ψ 2 ± , so we have ± , and since γ(s) tends to zero as |s| → ∞ by assumption, the left-hand side can be made arbitrarily close to zero by choosing s 0 large enough.
The prove the second claim, one has to find a trial function from the form domain of H D,v which makes the left-hand side of (5.2) negative.We use the idea proposed in this context in [GJ92] and choose the function in the form ψ := φ ς χ 0 + εf , where f will be chose later and φ ς is a suitable mollifier, say where g −1/2 (s) = a −a g −1/2 (s, u) du is the transverse average of the inverse Jacobian.By Fubini theorem the order of integration in the last three terms makes no difference.Eigenfunctions of h D,v may be chosen real and by assumption we have an easy integration by parts with (1 − uγ(s))χ 0 (u)du using the fact, that χ 0 (±a) = 0, shows that the sum of the three terms is zero.The remaining term is easily estimated; we get the bound (5.5) ?eq:HQWExtendedFir showing that by the choice ς the positive contribution from the trial function tails can be made arbitrarily small.To complete the construction, we have to choose the function f .We pick it from C ∞ 0 ((−s 1 , s 1 ) × I a ) in which case ψ ∈ H 1 (R × I a ) and (5.6) GJ where q(•, •) is the corresponding sesquilinear form.Using the fact that the supports of φ ς χ 0 and f are disjoint, we evaluate easily the linear term coefficient, Since f is picked from an infinitely dimensional space, we can obviously choose it in such a way that the coefficient is nonzero and, say, negative.For small ε > 0 the linear term dominates over the quadratic one; then we can choose ε so that the sum of the last two terms on the right-hand side of (5.6) is negative, and subsequently to choose ς in order not to spoil the negativity of the whole estimating expression.
Note that we have identified only the threshold of σ ess (H D,v ).Under stronger regularity assumptions one can prove that the essential spectrum covers the whole interval [ D,v , ∞) but we will not need it here.

Making the ditch deeper
? ss:deep ?Let us return to our problem in the situation with parallel asymptotes, β = 0.A natural way to achieve strong coupling for a regular potential channel consists of modifying its depth without changing the profile of its bottom.For the sake of brevity we denote by χ a the characteristic function of the set Ω a ⊂ R 2 , with an abuse of notation we will use the same symbol the characteristic function of the interval I a ⊂ R and for the measure χ a (x)dx on R; then we consider the behavior of the operator H Γ,V +λχa .
Proof.Consider the operator family {H Γ,V +λχa +λI : λ ≥ 0}.The corresponding family of quadratic forms is monotonously increasing, then it follows from [RS80, Thm.S.14] that these operators converge to H D,V as λ → ∞ in the generalized strong resolvent sense [W00, Sec.9.3].By the same argument, h Γ,v+λχa + λI → h D,v as λ → ∞ in the generalized strong resolvent sense.We will use the last result to prove (5.7).Its first relation follows form Proposition 2.1.To get the second one, we note that by the standard double-well estimate [Si84, Thm.1.5] we have where the constant c is given by the associated Agmon metric.
The second claim follows from Theorem 5.1.As the curvature γ is by assumption nonzero and compactly supported, the assumptions are satisfied and the limiting operator H D,V has at least one eigenvalue below D,v , hence H Γ,V +λχa + λI has an eigenvalue below µ ⊥ +λχa,ρ , and consequently, below −λ + D,v for λ large enough.
Then there is a function g 0 with the indicated properties such that σ disc (H Γ,V g(λ) ) is nonempty for any g ≥ g 0 and all λ large enough.
Proof.Instead of (5.8), consider the two-parameter family of operators H Γ,V ξ,λ with the profile potential v ξ,λ (x) := ξλv(λx).As we have already mentioned, for a fixed ξ > 0 these operators converge in the norm-resolvent sense to H Γ,α ξ with α ξ := ξ v 1 as λ → ∞, cf.[EI01,BEHL17].The difficulty to deal with is that the essential spectrum threshold depend on ξ.Given a below bounded self-adjoint operator A, we denote by µ k (A) the corresponding numbers obtained from the minimax principle ‡ [RS78, Sec.XIII.1]and by µ ∞ (A) the limiting value of this sequence (for operators with a finite discrete spectrum we consider here, the sequence is constant from some index on).Then we introduce the operators HΓ,V ξ,λ := H Γ,V ξ,λ − µ ∞ (H Γ,V ξ,λ ) and HΓ,α ξ := H Γ,α ξ − µ ∞ (H Γ,α ξ ) ‡ This is a traditional notation and we are sure there is no danger of confusion with the symbol used to describe measure-type potentials.
Note that one could get a better idea of the function g 0 looking into the error terms of the tow limits involved, but we will not follow this route here.

A numerical example
s:num While we have been able to establish the weak/strong coupling dichotomy, finer properties are out of reach for an analytical treatment, for instance, it is not easy to find the critical strength needed to produce the discrete spectrum.They can be dealt with numerically as we are going to illustrate here.As an example, we consider the U-shaped channel with a polynomial profile, a multiple of where Ω := (−ρ − a, −ρ + a) (ρ − a, ρ + a) and α ≥ 2 is an even integer.Such a profile obviously tends to a rectangular well as α → ∞ but we work with finite values as it makes the numerics easier.Naturally, we have to solve the equation in a finite region.We do it for both the Dirichlet and Neumann boundaries which allows us to control the precision using the gap between the two.In computations, we made use of the Spectral Method [T00, Chap.9].Looking for the critical coupling constant of the potential λv α , we plot in Figure 2 the dependence of two dimensionless quantities, − |λ|A, where A := 1 π a −a v α (x) dx, and a/ρ, the relative width of the channel, for two values of α.The former is, of course, an upper estimate, asymptotically exact, of the number of bound states minus one in the one-dimensional potential well of that profile.Predictably, the curves are monotonously decreasing since the tunnelling between the parallel channels plays a more prominent role as a increases and a stronger coupling is needed to compete with it.Also the dependence on α makes sense, as the dimensionless 'volume' increases as α grows.We  also see that for larger α the curve becomes steeper as we approach the point a = ρ as the 'residual' barrier is thinner then.
For illustration we also plot in Figure 3 the eigenfunction in the situation when the coupling is stronger than critical; the thin black line in the picture indicates the potential support.

Concluding remarks
While the shape of the curved part in our result is quite general, modulo the regularity requirement, the guide is supposed to be straight outside a compact.A question arising naturally is what happens if Γ approached parallel lines only asymptotically; one expects that the general picture would remain similar as here.Other type of asymptotic behaviors, on the other hand, may lead to different spectral properties.One can ask, for instance, whether there can be geometrically induced eigenfunctions supported away from the 'bend' of Γ; one expects this to happen, say, when the double guide is locally bent or the gap between two channels is locally diminished at an appropriate place.

Figure 2 .
Figure 2. The critical interaction strength vs. the relative channel width iticalLambda