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

On the Fundamental Eigenvalue Gap of Sturm-Liouville Operators

Mohammed Ahrami, Zakaria El Allali, and Evans M. Harrell II
Abstract.

We use methods of direct optimization as in [9] to find the minimizers of the fundamental gap of Sturm-Liouville operators on an interval, under the constraint that the potential is of single-well form and that the weight function is of single-barrier form, and under similar constraints expressed in terms of convexity.

Key words and phrases:
Fundamental spectral gap, eigenvalue estimates, Sturm-Liouville operators, single-well potential, Dirichlet boundary conditions
2010 Mathematics Subject Classification:
34L15, 34B27, 35J60, 35B05.

1. Introduction

The goal of this article is to find optimal estimates, under constraints on the form of the coefficient functions in (1), for the fundamental eigenvalue gap Γ:=λ2λ1assignΓsubscript𝜆2subscript𝜆1\Gamma:=\lambda_{2}-\lambda_{1}roman_Γ := italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT of the Sturm-Liouville equation

(1) H(p,q)u:=ddx(p(x)dudx)+V(x)u=λw(x)uassign𝐻𝑝𝑞𝑢𝑑𝑑𝑥𝑝𝑥𝑑𝑢𝑑𝑥𝑉𝑥𝑢𝜆𝑤𝑥𝑢H(p,q)u:=-\frac{d}{dx}\left(p(x)\frac{du}{dx}\right)+V(x)u=\lambda w(x)uitalic_H ( italic_p , italic_q ) italic_u := - divide start_ARG italic_d end_ARG start_ARG italic_d italic_x end_ARG ( italic_p ( italic_x ) divide start_ARG italic_d italic_u end_ARG start_ARG italic_d italic_x end_ARG ) + italic_V ( italic_x ) italic_u = italic_λ italic_w ( italic_x ) italic_u

on a finite interval, with self-adjoint boundary conditions, According to [21], §8.4, one may impose any separated homogeneous boundary conditions of the form

u(0)cosα(pu)(0)sinα𝑢0𝛼𝑝superscript𝑢0𝛼\displaystyle u(0)\cos\alpha-(pu^{\prime})(0)\sin\alphaitalic_u ( 0 ) roman_cos italic_α - ( italic_p italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ( 0 ) roman_sin italic_α =0absent0\displaystyle=0= 0
u(π)cosβ(pu)(π)sinβ𝑢𝜋𝛽𝑝superscript𝑢𝜋𝛽\displaystyle u(\pi)\cos\beta-(pu^{\prime})(\pi)\sin\betaitalic_u ( italic_π ) roman_cos italic_β - ( italic_p italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ( italic_π ) roman_sin italic_β =0,absent0\displaystyle=0,= 0 ,

where 0α,β<πformulae-sequence0𝛼𝛽𝜋0\leq\alpha,\beta<\pi0 ≤ italic_α , italic_β < italic_π, to make H𝐻Hitalic_H self-adjoint; the interval has been standardized as [0,π]0𝜋[0,\pi][ 0 , italic_π ] without loss of generality. To keep the exposition simple we restrict to Dirichlet boundary conditions, i.e., u(0)=u(π)=0𝑢0𝑢𝜋0u(0)=u(\pi)=0italic_u ( 0 ) = italic_u ( italic_π ) = 0, although our techniques also work with other self-adjoint boundary conditions, with suitable changes. Some further simplifications will be imposed below.

The quantity ΓΓ\Gammaroman_Γ is of interest as the ionization energy in quantum theory, and sharp bounds for natural categories of potential energies V(x)𝑉𝑥V(x)italic_V ( italic_x ) that prevent its collapse, especially single-well and convex V(x)𝑉𝑥V(x)italic_V ( italic_x ), have been studied since the 1980s; cf. [1, 2, 3, 4, 5, 6, 10, 13, 14, 18] and references therein. The problem of maximizing ΓΓ\Gammaroman_Γ not only casts light on the physical problem of ionization, but is interesting as a mathematical problem in its own right.

Most prior work has assumed or emphasized the case where both p(x)𝑝𝑥p(x)italic_p ( italic_x ) and the weight w(x)𝑤𝑥w(x)italic_w ( italic_x ) are held constant. Liouville transformations allow one to convert (1) into equivalent equations with different p(x)𝑝𝑥p(x)italic_p ( italic_x ), V(x)𝑉𝑥V(x)italic_V ( italic_x ), and w(x)𝑤𝑥w(x)italic_w ( italic_x ), but in general only one of these three functions can be eliminated. In the first three sections of this article we standardize with p(x)1𝑝𝑥1p(x)\equiv 1italic_p ( italic_x ) ≡ 1, and recall how the Liouville transformation works in an appendix. In summary, these remarks allow us to concentrate on the problem

(2) {u′′+V(x)u=λw(x)u,x[0,π]u(0)=u(π)=0casessuperscript𝑢′′𝑉𝑥𝑢𝜆𝑤𝑥𝑢𝑥0𝜋𝑢0𝑢𝜋0missing-subexpression\left\{\begin{array}[]{lr}-u^{\prime\prime}+V(x)u=\lambda w(x)u,&x\in[0,\pi]\\ u(0)=u(\pi)=0&\end{array}\right.{ start_ARRAY start_ROW start_CELL - italic_u start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT + italic_V ( italic_x ) italic_u = italic_λ italic_w ( italic_x ) italic_u , end_CELL start_CELL italic_x ∈ [ 0 , italic_π ] end_CELL end_ROW start_ROW start_CELL italic_u ( 0 ) = italic_u ( italic_π ) = 0 end_CELL start_CELL end_CELL end_ROW end_ARRAY

in the following two sections.

In [2], for the problem (2) with Dirichlet conditions, Ashbaugh and Benguria proved that the optimal lower bound for ΓΓ\Gammaroman_Γ for symmetric single-well potentials is achieved if and only if V𝑉Vitalic_V is constant on (0,π)0𝜋(0,\pi)( 0 , italic_π ). In [18] Lavine considered the class of convex potentials on [0,π]0𝜋[0,\pi][ 0 , italic_π ] and proved, with either Dirichlet or Neumann boundary conditions, that the constant potential function minimizes ΓΓ\Gammaroman_Γ. Later Horváth [13] returned with Lavine’s methods to the problem of single-well potentials, but without symmetry assumptions, and again showed that the constant potential was optimal with some restrictions on the transition point, and in 2015 Yu and Yang [20] extended Horváth’s result by allowing other transition points and both Dirichlet and Neumann conditions. More recently El Allali and Harrell [9] used direct optimization methods to prove sharp lower bounds for ΓΓ\Gammaroman_Γ with general single-well potential V(x)𝑉𝑥V(x)italic_V ( italic_x ), without any restriction on the transition point a[0,π]𝑎0𝜋a\in[0,\pi]italic_a ∈ [ 0 , italic_π ], and obtained similar results in the case where the potential is convex. El Allali and Harrell were furthermore able to analyze the case where V=V0+V1𝑉subscript𝑉0subscript𝑉1V=V_{0}+V_{1}italic_V = italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, where V0subscript𝑉0V_{0}italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is a fixed background potential energy and V1subscript𝑉1V_{1}italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is assumed either single-well or convex. In contrast to the earlier studies of single-well potentials, which restrict the transition point in one way or other, the minimizing single-well potentials they found are in general step functions and not necessarily constant, unless extra conditions are imposed. In the classic case where p=1𝑝1p=1italic_p = 1 they recovered with different arguments the result of Lavine that ΓΓ\Gammaroman_Γ is uniquely minimized among convex V𝑉Vitalic_V by the constant, and in the case of single-well potentials, with no restrictions on the position of the minimum, they obtained a new, sharp bound, that Γ>2.04575Γ2.04575\Gamma>2.04575\dotsroman_Γ > 2.04575 …. Some further related articles are Huang’s discussion of the eigenvalue gap [15] and eigenvalue ratio [16] for the vibrating string with symmetric densities, i.e., allowing variable p(x)𝑝𝑥p(x)italic_p ( italic_x ), and the works of Ashbaugh and Benguria [4], Huang and Law [17], and Horvàth and Kiss [14], which include other expressions related to the low-lying eigenvalues such as eigenvalue ratios like λnλ1subscript𝜆𝑛subscript𝜆1\frac{\lambda_{n}}{\lambda_{1}}divide start_ARG italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG.

2. Simple properties of the fundamental gap ΓΓ\Gammaroman_Γ

We shall use expressions such as λk(V,w)subscript𝜆𝑘𝑉𝑤\lambda_{k}(V,w)italic_λ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_V , italic_w ) and Γ(V,w)=λ2(V,w)λ1(V,w)Γ𝑉𝑤subscript𝜆2𝑉𝑤subscript𝜆1𝑉𝑤\Gamma(V,w)=\lambda_{2}(V,w)-\lambda_{1}(V,w)roman_Γ ( italic_V , italic_w ) = italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V , italic_w ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V , italic_w ) to indicate the dependence on the fundamental gap on coefficients in (2) with respect to which we wish to optimize. In this section we review and slightly extend some useful observations about ΓΓ\Gammaroman_Γ that are familiar from previous sources such as [2, 18, 9]. Most importantly, there is an explicit formula for the first derivative with respect to perturbations of V𝑉Vitalic_V and w𝑤witalic_w:

Lemma 2.1.

Suppose that V(.,t)V(.,t)italic_V ( . , italic_t ) and w(.,t)w(.,t)italic_w ( . , italic_t ) are one-parameter families of real-valued, locally L1superscript𝐿1L^{1}italic_L start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT functions, with infV(x,κ)>infimum𝑉𝑥𝜅\inf V(x,\kappa)>-\inftyroman_inf italic_V ( italic_x , italic_κ ) > - ∞, Cw(x,κ)1C𝐶𝑤𝑥𝜅1𝐶C\geq w(x,\kappa)\geq\frac{1}{C}italic_C ≥ italic_w ( italic_x , italic_κ ) ≥ divide start_ARG 1 end_ARG start_ARG italic_C end_ARG for some C>0𝐶0C>0italic_C > 0, and Vκ(x,κ)𝑉𝜅𝑥𝜅\frac{\partial V}{\partial\kappa}(x,\kappa)divide start_ARG ∂ italic_V end_ARG start_ARG ∂ italic_κ end_ARG ( italic_x , italic_κ ) and wκ(x,κ)L1(0,π)𝑤𝜅𝑥𝜅superscript𝐿10𝜋\frac{\partial w}{\partial\kappa}(x,\kappa)\in L^{1}(0,\pi)divide start_ARG ∂ italic_w end_ARG start_ARG ∂ italic_κ end_ARG ( italic_x , italic_κ ) ∈ italic_L start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( 0 , italic_π ). Then

dλn(κ)dκ=λn0πwκ(x,κ)un2(x,κ)𝑑x+0πVκ(x,κ)un2(x,κ)𝑑x.𝑑subscript𝜆𝑛𝜅𝑑𝜅subscript𝜆𝑛superscriptsubscript0𝜋𝑤𝜅𝑥𝜅superscriptsubscript𝑢𝑛2𝑥𝜅differential-d𝑥superscriptsubscript0𝜋𝑉𝜅𝑥𝜅superscriptsubscript𝑢𝑛2𝑥𝜅differential-d𝑥\frac{d\lambda_{n}(\kappa)}{d\kappa}=-\lambda_{n}\int_{0}^{\pi}\frac{\partial w% }{\partial\kappa}(x,\kappa)u_{n}^{2}(x,\kappa)dx+\int_{0}^{\pi}\frac{\partial V% }{\partial\kappa}(x,\kappa)u_{n}^{2}(x,\kappa)dx.divide start_ARG italic_d italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_κ ) end_ARG start_ARG italic_d italic_κ end_ARG = - italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT divide start_ARG ∂ italic_w end_ARG start_ARG ∂ italic_κ end_ARG ( italic_x , italic_κ ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , italic_κ ) italic_d italic_x + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT divide start_ARG ∂ italic_V end_ARG start_ARG ∂ italic_κ end_ARG ( italic_x , italic_κ ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , italic_κ ) italic_d italic_x .
Proof.

Because Vκ𝑉𝜅\frac{\partial V}{\partial\kappa}divide start_ARG ∂ italic_V end_ARG start_ARG ∂ italic_κ end_ARG and wκ𝑤𝜅\frac{\partial w}{\partial\kappa}divide start_ARG ∂ italic_w end_ARG start_ARG ∂ italic_κ end_ARG are relatively bounded perturbations, Kato’s theory of analytic perturbations applies, and since the eigenvalues with separated homogeneous boundary conditions are simple, this justifies the use of a formal expansion to calculate the effect of the perturbation, à la Feynman-Hellmann: Denoting u˙=uκ˙𝑢𝑢𝜅\dot{u}=\frac{\partial u}{\partial\kappa}over˙ start_ARG italic_u end_ARG = divide start_ARG ∂ italic_u end_ARG start_ARG ∂ italic_κ end_ARG, differentiation of (2) with respect to κ𝜅\kappaitalic_κ gives

u˙n′′+(λ˙nw+λnw˙V˙)un+(λnwV)u˙n=0.subscriptsuperscript˙𝑢′′𝑛subscript˙𝜆𝑛𝑤subscript𝜆𝑛˙𝑤˙𝑉subscript𝑢𝑛subscript𝜆𝑛𝑤𝑉subscript˙𝑢𝑛0\dot{u}^{\prime\prime}_{n}+(\dot{\lambda}_{n}w+\lambda_{n}\dot{w}-\dot{V})u_{n% }+(\lambda_{n}w-V)\dot{u}_{n}=0.over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + ( over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w + italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over˙ start_ARG italic_w end_ARG - over˙ start_ARG italic_V end_ARG ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + ( italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w - italic_V ) over˙ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 0 .

We multiply the equation above by un(x,κ)subscript𝑢𝑛𝑥𝜅u_{n}(x,\kappa)italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_x , italic_κ ) and integrate with respect to x𝑥xitalic_x from 00 to π𝜋\piitalic_π. This yields

0πu˙n′′un𝑑x+0π(λnwV)u˙nun𝑑x=0π(λ˙nw+λnw˙V˙)un2𝑑x.superscriptsubscript0𝜋subscriptsuperscript˙𝑢′′𝑛subscript𝑢𝑛differential-d𝑥superscriptsubscript0𝜋subscript𝜆𝑛𝑤𝑉subscript˙𝑢𝑛subscript𝑢𝑛differential-d𝑥superscriptsubscript0𝜋subscript˙𝜆𝑛𝑤subscript𝜆𝑛˙𝑤˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥\int_{0}^{\pi}\dot{u}^{\prime\prime}_{n}u_{n}dx+\int_{0}^{\pi}(\lambda_{n}w-V)% \dot{u}_{n}u_{n}dx=-\int_{0}^{\pi}(\dot{\lambda}_{n}w+\lambda_{n}\dot{w}-\dot{% V})u_{n}^{2}dx.∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w - italic_V ) over˙ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x = - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT ( over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w + italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over˙ start_ARG italic_w end_ARG - over˙ start_ARG italic_V end_ARG ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x .

Observing that un(λnwV)=un′′subscript𝑢𝑛subscript𝜆𝑛𝑤𝑉superscriptsubscript𝑢𝑛′′u_{n}(\lambda_{n}w-V)=-u_{n}^{\prime\prime}italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w - italic_V ) = - italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT,

0πu˙n′′un𝑑x0πu˙nun′′𝑑x=0π(λ˙nw+λnw˙V˙)un2𝑑x.superscriptsubscript0𝜋subscriptsuperscript˙𝑢′′𝑛subscript𝑢𝑛differential-d𝑥superscriptsubscript0𝜋subscript˙𝑢𝑛subscriptsuperscript𝑢′′𝑛differential-d𝑥superscriptsubscript0𝜋subscript˙𝜆𝑛𝑤subscript𝜆𝑛˙𝑤˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥\int_{0}^{\pi}\dot{u}^{\prime\prime}_{n}u_{n}dx-\int_{0}^{\pi}\dot{u}_{n}u^{% \prime\prime}_{n}dx=-\int_{0}^{\pi}(\dot{\lambda}_{n}w+\lambda_{n}\dot{w}-\dot% {V})u_{n}^{2}dx.∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x = - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT ( over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w + italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over˙ start_ARG italic_w end_ARG - over˙ start_ARG italic_V end_ARG ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x .

Integrating by parts twice yields

[u˙nun]0π0πu˙nun𝑑x[u˙nun]0π+0πu˙nun𝑑x=0π(λ˙nw+λnw˙V˙)un2𝑑x,superscriptsubscriptdelimited-[]subscriptsuperscript˙𝑢𝑛subscript𝑢𝑛0𝜋superscriptsubscript0𝜋subscriptsuperscript˙𝑢𝑛subscriptsuperscript𝑢𝑛differential-d𝑥superscriptsubscriptdelimited-[]subscriptsuperscript˙𝑢𝑛subscriptsuperscript𝑢𝑛0𝜋superscriptsubscript0𝜋subscriptsuperscript˙𝑢𝑛subscriptsuperscript𝑢𝑛differential-d𝑥superscriptsubscript0𝜋subscript˙𝜆𝑛𝑤subscript𝜆𝑛˙𝑤˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥[\dot{u}^{\prime}_{n}u_{n}]_{0}^{\pi}-\int_{0}^{\pi}\dot{u}^{\prime}_{n}u^{% \prime}_{n}dx-[\dot{u}^{\prime}_{n}u^{\prime}_{n}]_{0}^{\pi}+\int_{0}^{\pi}% \dot{u}^{\prime}_{n}u^{\prime}_{n}dx=-\int_{0}^{\pi}(\dot{\lambda}_{n}w+% \lambda_{n}\dot{w}-\dot{V})u_{n}^{2}dx,[ over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x - [ over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_d italic_x = - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT ( over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w + italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over˙ start_ARG italic_w end_ARG - over˙ start_ARG italic_V end_ARG ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x ,

with boundary conditions un(0,t)=un(π,t)=0subscript𝑢𝑛0𝑡subscript𝑢𝑛𝜋𝑡0u_{n}(0,t)=u_{n}(\pi,t)=0italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( 0 , italic_t ) = italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_π , italic_t ) = 0. This gives

0π(λ˙nw+λnw˙V˙)un2𝑑x=0,superscriptsubscript0𝜋subscript˙𝜆𝑛𝑤subscript𝜆𝑛˙𝑤˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥0-\int_{0}^{\pi}(\dot{\lambda}_{n}w+\lambda_{n}\dot{w}-\dot{V})u_{n}^{2}dx=0,- ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT ( over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w + italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over˙ start_ARG italic_w end_ARG - over˙ start_ARG italic_V end_ARG ) italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x = 0 ,

so that

λ˙n0πun2w𝑑x=λn0πw˙un2𝑑x+0πV˙un2𝑑x.subscript˙𝜆𝑛superscriptsubscript0𝜋superscriptsubscript𝑢𝑛2𝑤differential-d𝑥subscript𝜆𝑛superscriptsubscript0𝜋˙𝑤superscriptsubscript𝑢𝑛2differential-d𝑥superscriptsubscript0𝜋˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥\dot{\lambda}_{n}\int_{0}^{\pi}u_{n}^{2}wdx=-\lambda_{n}\int_{0}^{\pi}\dot{w}u% _{n}^{2}dx+\int_{0}^{\pi}\dot{V}u_{n}^{2}dx.over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_w italic_d italic_x = - italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_w end_ARG italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_V end_ARG italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x .

Noting that 0πwun2𝑑x=1,superscriptsubscript0𝜋𝑤superscriptsubscript𝑢𝑛2differential-d𝑥1\displaystyle\int_{0}^{\pi}wu_{n}^{2}dx=1,∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT italic_w italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x = 1 ,

λ˙n=λn0πw˙un2𝑑x+0πV˙un2𝑑x.subscript˙𝜆𝑛subscript𝜆𝑛superscriptsubscript0𝜋˙𝑤superscriptsubscript𝑢𝑛2differential-d𝑥superscriptsubscript0𝜋˙𝑉superscriptsubscript𝑢𝑛2differential-d𝑥\dot{\lambda}_{n}=-\lambda_{n}\int_{0}^{\pi}\dot{w}u_{n}^{2}dx+\int_{0}^{\pi}% \dot{V}u_{n}^{2}dx.over˙ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = - italic_λ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_w end_ARG italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x + ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT over˙ start_ARG italic_V end_ARG italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x .

We next adapt the monotonicity argument of [2, 9] to incorporate the weight:

Lemma 2.2.

Consider the problem (2) with the same assumptions on V𝑉Vitalic_V and w𝑤witalic_w as in Lemma 2.1. Without loss of generality we standardize the first two normalized eigenfunctions so that u1,2(x)>0subscript𝑢12𝑥0u_{1,2}(x)>0italic_u start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT ( italic_x ) > 0 for 0<x<ϵ0𝑥italic-ϵ0<x<\epsilon0 < italic_x < italic_ϵ for some ϵitalic-ϵ\epsilonitalic_ϵ. Then:

  1. (1)

    u2u1subscript𝑢2subscript𝑢1\dfrac{u_{2}}{u_{1}}divide start_ARG italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG is decreasing on (0,π)0𝜋(0,\pi)( 0 , italic_π ).

  2. (2)

    The equation |u1(x)|=|u2(x)|subscript𝑢1𝑥subscript𝑢2𝑥|u_{1}(x)|=|u_{2}(x)|| italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) | = | italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) | has either one or two solutions on (0,π)0𝜋(0,\pi)( 0 , italic_π ).

  3. (3)

    There exist two points xsubscript𝑥x_{-}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT and x+subscript𝑥x_{+}italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT 0x<x+π0subscript𝑥subscript𝑥𝜋0\leq x_{-}<x_{+}\leq\pi0 ≤ italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT < italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ≤ italic_π, at least one of which is interior to (0,π)0𝜋(0,\pi)( 0 , italic_π ), such that u12(x)>u22(x)superscriptsubscript𝑢12𝑥superscriptsubscript𝑢22𝑥u_{1}^{2}(x)>u_{2}^{2}(x)italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on (x,x+)subscript𝑥subscript𝑥(x_{-},x_{+})( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) and u12(x)u22(x)superscriptsubscript𝑢12𝑥superscriptsubscript𝑢22𝑥u_{1}^{2}(x)\leq u_{2}^{2}(x)italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ≤ italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on (x,x+)csuperscriptsubscript𝑥subscript𝑥𝑐(x_{-},x_{+})^{c}( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT.

  4. (4)

    The equation λ1|u12(x)|=λ2|u22(x)|subscript𝜆1superscriptsubscript𝑢12𝑥subscript𝜆2superscriptsubscript𝑢22𝑥\lambda_{1}|u_{1}^{2}(x)|=\lambda_{2}|u_{2}^{2}(x)|italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) | = italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) | has either one or two solutions on (0,π)0𝜋(0,\pi)( 0 , italic_π ).

  5. (5)

    There exist two points x^subscript^𝑥\widehat{x}_{-}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT and x^+subscript^𝑥\widehat{x}_{+}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT 0x^<x^+π0subscript^𝑥subscript^𝑥𝜋0\leq\widehat{x}_{-}<\widehat{x}_{+}\leq\pi0 ≤ over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ≤ italic_π, at least one of which is interior to (0,π)0𝜋(0,\pi)( 0 , italic_π ), such that λ1u12(x)>λ2u22(x)subscript𝜆1superscriptsubscript𝑢12𝑥subscript𝜆2superscriptsubscript𝑢22𝑥\lambda_{1}u_{1}^{2}(x)>\lambda_{2}u_{2}^{2}(x)italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on (x^,x^+)subscript^𝑥subscript^𝑥(\widehat{x}_{-},\widehat{x}_{+})( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) and λ1u12(x)λ2u22(x)subscript𝜆1superscriptsubscript𝑢12𝑥subscript𝜆2superscriptsubscript𝑢22𝑥\lambda_{1}u_{1}^{2}(x)\leq\lambda_{2}u_{2}^{2}(x)italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ≤ italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on (x^,x^+)csuperscriptsubscript^𝑥subscript^𝑥𝑐(\widehat{x}_{-},\widehat{x}_{+})^{c}( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT.

Proof.

We first show that (u2u1)<0superscriptsubscript𝑢2subscript𝑢10\left(\dfrac{u_{2}}{u_{1}}\right)^{\prime}<0( divide start_ARG italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT < 0 for 0<x<x00𝑥subscript𝑥00<x<x_{0}0 < italic_x < italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, where u2(x0)=0subscript𝑢2subscript𝑥00u_{2}(x_{0})=0italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = 0 and hence that there can be at most one value x(0,x0)subscript𝑥0subscript𝑥0x_{-}\in(0,x_{0})italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ∈ ( 0 , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) for which u1(x)=u2(x)subscript𝑢1subscript𝑥subscript𝑢2subscript𝑥u_{1}(x_{-})=u_{2}(x_{-})italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) = italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ).

The Wronskian is by definition

W(x)=u1(x)u2(x)u2(x)u1(x).𝑊𝑥subscript𝑢1𝑥superscriptsubscript𝑢2𝑥subscript𝑢2𝑥superscriptsubscript𝑢1𝑥W(x)=u_{1}(x)u_{2}^{\prime}(x)-u_{2}(x)u_{1}^{\prime}(x).italic_W ( italic_x ) = italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) - italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) .

Thus with the weight in (2),

W(x)=(λ1λ2)w(x)u1(x)u2(x),superscript𝑊𝑥subscript𝜆1subscript𝜆2𝑤𝑥subscript𝑢1𝑥subscript𝑢2𝑥W^{\prime}(x)=(\lambda_{1}-\lambda_{2})w(x)u_{1}(x)u_{2}(x),italic_W start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) = ( italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_w ( italic_x ) italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) ,

and for the the quotient v(x):=u2(x)u1(x)assign𝑣𝑥subscript𝑢2𝑥subscript𝑢1𝑥v(x):=\dfrac{u_{2}(x)}{u_{1}(x)}italic_v ( italic_x ) := divide start_ARG italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) end_ARG,

v(x)=u1(x)u2(x)u2(x)u1(x)u12(x)=W(x)u12(x).superscript𝑣𝑥subscript𝑢1𝑥superscriptsubscript𝑢2𝑥subscript𝑢2𝑥superscriptsubscript𝑢1𝑥superscriptsubscript𝑢12𝑥𝑊𝑥superscriptsubscript𝑢12𝑥v^{\prime}(x)=\frac{u_{1}(x)u_{2}^{\prime}(x)-u_{2}(x)u_{1}^{\prime}(x)}{u_{1}% ^{2}(x)}=\frac{W(x)}{u_{1}^{2}(x)}.italic_v start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) = divide start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) - italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) end_ARG = divide start_ARG italic_W ( italic_x ) end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) end_ARG .

Hence

(3) v(x)=1u12(x)0x(λ1λ2)w(t)u1(t)u2(t)𝑑t<0,superscript𝑣𝑥1superscriptsubscript𝑢12𝑥superscriptsubscript0𝑥subscript𝜆1subscript𝜆2𝑤𝑡subscript𝑢1𝑡subscript𝑢2𝑡differential-d𝑡0v^{\prime}(x)=\frac{1}{u_{1}^{2}(x)}\int_{0}^{x}(\lambda_{1}-\lambda_{2})w(t)u% _{1}(t)u_{2}(t)dt<0,italic_v start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_x ) = divide start_ARG 1 end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) end_ARG ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT ( italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_w ( italic_t ) italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_t ) italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_t ) italic_d italic_t < 0 ,

since λ1<λ2subscript𝜆1subscript𝜆2\lambda_{1}<\lambda_{2}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and u1,u2>0subscript𝑢1subscript𝑢20u_{1},u_{2}>0italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT > 0 on (0,x0)0subscript𝑥0(0,x_{0})( 0 , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ).

Suppose that there exist distinct α1,2(0,x0)subscript𝛼120subscript𝑥0\alpha_{1,2}\in(0,x_{0})italic_α start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT ∈ ( 0 , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) such that

u2(αi)=u1(αi),i=1,2.formulae-sequencesubscript𝑢2subscript𝛼𝑖subscript𝑢1subscript𝛼𝑖𝑖12u_{2}(\alpha_{i})=u_{1}(\alpha_{i}),\leavevmode\nobreak\ \leavevmode\nobreak\ % \leavevmode\nobreak\ \leavevmode\nobreak\ i=1,2\leavevmode\nobreak\ .italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_i = 1 , 2 .

Then v(α1)=v(α2)𝑣subscript𝛼1𝑣subscript𝛼2v(\alpha_{1})=v(\alpha_{2})italic_v ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = italic_v ( italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ). By Rolle’s theorem there exists ξ(α1,α2)(0,x0)𝜉subscript𝛼1subscript𝛼20subscript𝑥0\xi\in(\alpha_{1},\alpha_{2})\subset(0,x_{0})italic_ξ ∈ ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⊂ ( 0 , italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) such that:::: v(ξ)=0superscript𝑣𝜉0v^{\prime}(\xi)=0italic_v start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_ξ ) = 0, but this contradicts (3).

Since u2subscript𝑢2u_{2}italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT vanishes at a unique point x0subscript𝑥0x_{0}italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and the same argument can be carried out after the change of variables xπx𝑥𝜋𝑥x\to\pi-xitalic_x → italic_π - italic_x and an adjustment of the sign of u2subscript𝑢2u_{2}italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, it follows that (u2u1)subscript𝑢2subscript𝑢1\left(\dfrac{u_{2}}{u_{1}}\right)( divide start_ARG italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ) is strictly monotonic on (x0,π)subscript𝑥0𝜋(x_{0},\pi)( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_π ) and that there is at most one value x+(x0,π)subscript𝑥subscript𝑥0𝜋x_{+}\in(x_{0},\pi)italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ∈ ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_π ) for which u1(x+)=u2(x+)subscript𝑢1subscript𝑥subscript𝑢2subscript𝑥u_{1}(x_{+})=u_{2}(x_{+})italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) = italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ).

At least one of the points x±(0,π)subscript𝑥plus-or-minus0𝜋x_{\pm}\in(0,\pi)italic_x start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ∈ ( 0 , italic_π ), because if x=0subscript𝑥0x_{-}=0italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT = 0 and x+=πsubscript𝑥𝜋x_{+}=\piitalic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT = italic_π then u1(x)>|u2(x)|subscript𝑢1𝑥subscript𝑢2𝑥u_{1}(x)>|u_{2}(x)|italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) > | italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) | for all x(0,π)𝑥0𝜋x\in(0,\pi)italic_x ∈ ( 0 , italic_π ), which would contradict u12=u22subscriptnormsubscript𝑢12subscriptnormsubscript𝑢22\|u_{1}\|_{2}=\|u_{2}\|_{2}∥ italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∥ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = ∥ italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ∥ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT. ∎

3. Characterization of optimizers

In this section, we determine the explicit form of the gap-minimizing potential and density function of problem (2), closely following the strategy of [9].

3.1. The class of single-well potentials and single-barrier densities

Definition 3.1.

Let 1<M1𝑀1<M\leq\infty1 < italic_M ≤ ∞. The function V𝑉Vitalic_V is called a single-well function if V𝑉Vitalic_V is non-increasing on [0,a]0𝑎[0,a][ 0 , italic_a ] and non-decreasing on [a,π]𝑎𝜋[a,\pi][ italic_a , italic_π ], for some a[0,π]𝑎0𝜋a\in[0,\pi]italic_a ∈ [ 0 , italic_π ]. The point a𝑎aitalic_a is called a transition point (with no assumption of uniqueness). The notation below will be used through this article:

SW[0,π],M={V(x):0V(x)M,where V is a single-well function on [0,π]}.𝑆subscript𝑊0𝜋𝑀conditional-set𝑉𝑥formulae-sequence0𝑉𝑥𝑀where 𝑉 is a single-well function on 0𝜋SW_{[0,\pi],M}=\{V(x):0\leq V(x)\leq M,\textrm{where }V\textrm{ is a single-% well function on }[0,\pi]\}.italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT = { italic_V ( italic_x ) : 0 ≤ italic_V ( italic_x ) ≤ italic_M , where italic_V is a single-well function on [ 0 , italic_π ] } .
Definition 3.2.

Let 0<N<N><0subscript𝑁subscript𝑁0<N_{<}\leq N_{>}<\infty0 < italic_N start_POSTSUBSCRIPT < end_POSTSUBSCRIPT ≤ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT < ∞. The function w𝑤witalic_w is called a single-barrier density if w𝑤witalic_w is non-decreasing on [0,b]0𝑏[0,b][ 0 , italic_b ] and non-increasing on [b,π]𝑏𝜋[b,\pi][ italic_b , italic_π ] for some b[0,π]𝑏0𝜋b\in[0,\pi]italic_b ∈ [ 0 , italic_π ]. The following notation will be used through this article:

SB[0,π],N<,N>={w(x):N<w(x)N>,where w is a single-barrier density on [0,π]}.𝑆subscript𝐵0𝜋𝑁subscript𝑁conditional-set𝑤𝑥formulae-sequencesubscript𝑁𝑤𝑥subscript𝑁where 𝑤 is a single-barrier density on 0𝜋SB_{[0,\pi],N<,N_{>}}=\left\{w(x):N_{<}\leq w(x)\leq N_{>},\textrm{where }w% \textrm{ is a single-barrier density on }[0,\pi]\right\}.italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT = { italic_w ( italic_x ) : italic_N start_POSTSUBSCRIPT < end_POSTSUBSCRIPT ≤ italic_w ( italic_x ) ≤ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT , where italic_w is a single-barrier density on [ 0 , italic_π ] } .
Definition 3.3.

Consider the Sturm Liouville problem

(4) u′′+(V0(x)+V(x))u=λw(x)usuperscript𝑢′′subscript𝑉0𝑥𝑉𝑥𝑢𝜆𝑤𝑥𝑢\displaystyle{-u^{\prime\prime}+(V_{0}(x)+V(x))u=\lambda w(x)u}- italic_u start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT + ( italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + italic_V ( italic_x ) ) italic_u = italic_λ italic_w ( italic_x ) italic_u
u(0)=u(π)=0𝑢0𝑢𝜋0\displaystyle u(0)=u(\pi)=0italic_u ( 0 ) = italic_u ( italic_π ) = 0 ,

where V𝑉Vitalic_V is a single-well function and w𝑤witalic_w is a single-barrier density. The background potential V0subscript𝑉0V_{0}italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is assumed bounded and measurable. If there exist VW[0,π],Msubscript𝑉subscript𝑊0𝜋𝑀V_{*}\in W_{[0,\pi],M}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and wSB[0,π],N<,N>subscript𝑤𝑆subscript𝐵0𝜋subscript𝑁subscript𝑁w_{*}\in SB_{[0,\pi],N_{<},N_{>}}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N start_POSTSUBSCRIPT < end_POSTSUBSCRIPT , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT such that

Γ(V,w)=inf(Γ(V,w),VSW[0,π],M,wSB[0,π],N<,N>),\Gamma(V_{*},w_{*})=\inf{(\Gamma(V,w),V\in SW_{[0,\pi],M},\leavevmode\nobreak% \ w\in SB_{[0,\pi],N<,N_{>}})},roman_Γ ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) = roman_inf ( roman_Γ ( italic_V , italic_w ) , italic_V ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT , italic_w ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ,

then we call the function Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT an optimal potential and the function wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT an optimal density for problem (4).

As in [9], we will use compactness of the sets SW𝑆𝑊SWitalic_S italic_W and SB𝑆𝐵SBitalic_S italic_B, following from a theorem of Helly:

Proposition 3.1.

For any sequence fnΛsubscript𝑓𝑛Λf_{n}\in\Lambdaitalic_f start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∈ roman_Λ, (Λ=SWΛ𝑆𝑊\Lambda=SWroman_Λ = italic_S italic_W or SB𝑆𝐵SBitalic_S italic_B with any fixed positive M,N<,>𝑀subscript𝑁M,N_{<,>}italic_M , italic_N start_POSTSUBSCRIPT < , > end_POSTSUBSCRIPT), there exist a subsequence fnksubscript𝑓subscript𝑛𝑘f_{n_{k}}italic_f start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT and a function fsubscript𝑓f_{\star}italic_f start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT such that fnk(x)f(x)Λsubscript𝑓subscript𝑛𝑘𝑥subscript𝑓𝑥Λf_{n_{k}}(x)\longrightarrow f_{\star}(x)\in\Lambdaitalic_f start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_x ) ⟶ italic_f start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ( italic_x ) ∈ roman_Λ for a.e. x𝑥xitalic_x.

For the proof, see [9], Proposition 2.12.12.12.1.

Corollary 3.1.

There exist a potential VSW[0,π],Msubscript𝑉𝑆subscript𝑊0𝜋𝑀V_{*}\in SW_{[0,\pi],M}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and a density wSB[0,π],N<,N>subscript𝑤𝑆subscript𝐵0𝜋𝑁subscript𝑁w_{*}\in SB_{[0,\pi],N<,N_{>}}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT that minimize Γ[V,w]Γ𝑉𝑤\Gamma[V,w]roman_Γ [ italic_V , italic_w ].

Proof.

According to Proposition 3.1, since the gap Γ[V,w]Γ𝑉𝑤\Gamma[V,w]roman_Γ [ italic_V , italic_w ] is positive, there exist minimizing sequences (Vn)nΛsubscriptsubscript𝑉𝑛𝑛Λ(V_{n})_{n\in\mathbb{N}}\subset\Lambda( italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_n ∈ blackboard_N end_POSTSUBSCRIPT ⊂ roman_Λ and (wn)nΛsubscriptsubscript𝑤𝑛𝑛Λ(w_{n})_{n\in\mathbb{N}}\subset\Lambda( italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_n ∈ blackboard_N end_POSTSUBSCRIPT ⊂ roman_Λ such that

limn+(λ2(Vn,wn)λ1(Vn,wn))=inf{λ2(V,w)λ1(V,w);V,wΛ}.subscript𝑛subscript𝜆2subscript𝑉𝑛subscript𝑤𝑛subscript𝜆1subscript𝑉𝑛subscript𝑤𝑛infimumsubscript𝜆2𝑉𝑤subscript𝜆1𝑉𝑤𝑉𝑤Λ\lim_{n\rightarrow+\infty}(\lambda_{2}(V_{n},w_{n})-\lambda_{1}(V_{n},w_{n}))=% \inf\{\lambda_{2}(V,w)-\lambda_{1}(V,w);\leavevmode\nobreak\ \leavevmode% \nobreak\ V,w\in\Lambda\}.roman_lim start_POSTSUBSCRIPT italic_n → + ∞ end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) ) = roman_inf { italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V , italic_w ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V , italic_w ) ; italic_V , italic_w ∈ roman_Λ } .

By Proposition 3.1 we may pass to subsequences in ΛΛ\Lambdaroman_Λ that converge pointwise a.e. to limits VΛsubscript𝑉ΛV_{*}\in\Lambdaitalic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ roman_Λ and wΛsubscript𝑤Λw_{*}\in\Lambdaitalic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ roman_Λ:

limn+Vn=Vandlimn+wn=wformulae-sequencesubscript𝑛subscript𝑉𝑛subscript𝑉andsubscript𝑛subscript𝑤𝑛subscript𝑤\lim_{n\rightarrow+\infty}V_{n}=V_{*}\quad\textsf{and}\quad\lim_{n\rightarrow+% \infty}w_{n}=w_{*}roman_lim start_POSTSUBSCRIPT italic_n → + ∞ end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT and roman_lim start_POSTSUBSCRIPT italic_n → + ∞ end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT

By the dominated convergence theorem these sequences also converge in L1(0,π)superscript𝐿10𝜋L^{1}(0,\pi)italic_L start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( 0 , italic_π ). Hence by continuity of ΓΓ\Gammaroman_Γ with respect to relatively bounded perturbations,

limn+(λ2(Vn,wn)λ1(Vn,wn))=inf{λ2(V,w)λ1(V,w);V,wΛ}=λ2(V,w)λ1(V,w).subscript𝑛subscript𝜆2subscript𝑉𝑛subscript𝑤𝑛subscript𝜆1subscript𝑉𝑛subscript𝑤𝑛infimumsubscript𝜆2𝑉𝑤subscript𝜆1𝑉𝑤𝑉𝑤Λsubscript𝜆2subscript𝑉subscript𝑤subscript𝜆1subscript𝑉subscript𝑤\lim_{n\rightarrow+\infty}(\lambda_{2}(V_{n},w_{n})-\lambda_{1}(V_{n},w_{n}))=% \inf\{\lambda_{2}(V,w)-\lambda_{1}(V,w);\leavevmode\nobreak\ \leavevmode% \nobreak\ V,w\in\Lambda\}=\lambda_{2}(V_{*},w_{*})-\lambda_{1}(V_{*},w_{*}).roman_lim start_POSTSUBSCRIPT italic_n → + ∞ end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) ) = roman_inf { italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V , italic_w ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V , italic_w ) ; italic_V , italic_w ∈ roman_Λ } = italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) .

Theorem 3.1.

For any piecewise continuous, strictly positive weight function w𝑤witalic_w, the optimal potential Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is a step function. For any piecewise continuous potential function V𝑉Vitalic_V, the optimal weight wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is a step function. The same is true for jointly optimal Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT and wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT. In each case the optimizers have the following characterization:

  1. (1)

    V(x)=0subscript𝑉𝑥0V_{*}(x)=0italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = 0 a.e. on a connected component of {x:u22(x)>u12(x)}conditional-set𝑥superscriptsubscript𝑢22𝑥superscriptsubscript𝑢12𝑥\{x:u_{2}^{2}(x)>u_{1}^{2}(x)\}{ italic_x : italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } and on the complement of that interval, V(x)=max(V)subscript𝑉𝑥subscript𝑉V_{*}(x)=\max(V_{*})italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) a.e.

  2. (2)

    w(x)=N>subscript𝑤𝑥subscript𝑁w_{*}(x)=N_{>}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT a.e. on a connected component of {x:λ2u22(x)>λ1u12(x)}conditional-set𝑥subscript𝜆2superscriptsubscript𝑢22𝑥subscript𝜆1superscriptsubscript𝑢12𝑥\{x:\lambda_{2}u_{2}^{2}(x)>\lambda_{1}u_{1}^{2}(x)\}{ italic_x : italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } and on the complement of that interval, w(x)=min(w)subscript𝑤𝑥subscript𝑤w_{*}(x)=\min(w_{*})italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = roman_min ( italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) a.e.

Proof.

Let VSW[0,π],Msubscript𝑉𝑆subscript𝑊0𝜋𝑀V_{*}\in SW_{[0,\pi],M}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and wSB[0,π],N<,N>subscript𝑤𝑆subscript𝐵0𝜋𝑁subscript𝑁w_{*}\in SB_{[0,\pi],N<,N_{>}}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT be the minimizing potential and density guaranteed by the lemma.

We can characterize Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT with the same argument as was used in [9], since Lemma 2.2 allows for the possibility of variable weight. We repeat it here to make the proof self-contained, and to point out a key difference. By Lemma 2.2 there exist x±subscript𝑥plus-or-minusx_{\pm}italic_x start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT: 0x<x+π0subscript𝑥subscript𝑥𝜋0\leqslant x_{-}<x_{+}\leqslant\pi0 ⩽ italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT < italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ⩽ italic_π, for which

u22(x)>u12(x) on (0,x)(x+,π)(one of these intervals may be vacuous)superscriptsubscript𝑢22𝑥superscriptsubscript𝑢12𝑥 on 0subscript𝑥subscript𝑥𝜋(one of these intervals may be vacuous)\displaystyle u_{2}^{2}(x)>u_{1}^{2}(x)\textrm{ on }(0,x_{-})\cup(x_{+},\pi)% \quad\textrm{(one of these intervals may be vacuous)}italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ∪ ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_π ) (one of these intervals may be vacuous)
u12(x)>u22(x) on (x,x+).superscriptsubscript𝑢12𝑥superscriptsubscript𝑢22𝑥 on subscript𝑥subscript𝑥\displaystyle u_{1}^{2}(x)>u_{2}^{2}(x)\textrm{ on }(x_{-},x_{+}).italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) .

We define a family of single-well potentials by perturbing Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT so that

V(x,κ)=κV1(x)+(1κ)V(x)t[0,1].formulae-sequence𝑉𝑥𝜅𝜅subscript𝑉1𝑥1𝜅subscript𝑉𝑥𝑡01V(x,\kappa)=\kappa V_{1}(x)+(1-\kappa)V_{*}(x)\quad t\in[0,1].italic_V ( italic_x , italic_κ ) = italic_κ italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) + ( 1 - italic_κ ) italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) italic_t ∈ [ 0 , 1 ] .

Next, we work out the explicit form of Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT in two cases, beginning with

(i) xa<x+subscript𝑥𝑎subscript𝑥x_{-}\leq a<x_{+}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ≤ italic_a < italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. For definiteness we arrange by reflecting if necessary that 0<xax+0subscript𝑥𝑎subscript𝑥0<x_{-}\leq a\leq x_{+}0 < italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ≤ italic_a ≤ italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, V(x+)V(x)subscript𝑉subscript𝑥subscript𝑉subscript𝑥V_{\star}(x_{+})\geq V_{*}(x_{-})italic_V start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) ≥ italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ), and Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is nondecreasing for xx+𝑥subscript𝑥x\geq x_{+}italic_x ≥ italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. In this case we proceed in two stages. First, let

V1(x)={V(x) on (0,a)V(x+) on (a,π).subscript𝑉1𝑥casessubscript𝑉subscript𝑥 on 0𝑎subscript𝑉subscript𝑥 on 𝑎𝜋V_{1}(x)=\left\{\begin{array}[]{c}V_{*}(x_{-})\textrm{ on }(0,a)\\ V_{*}(x_{+})\textrm{ on }(a,\pi).\end{array}\right.italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) on ( 0 , italic_a ) end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( italic_a , italic_π ) . end_CELL end_ROW end_ARRAY

We observe that V1SW[0,π],Msubscript𝑉1𝑆subscript𝑊0𝜋𝑀V_{1}\in SW_{[0,\pi],M}italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT has been chosen so that V(x,κ)SW[0,π],M𝑉𝑥𝜅𝑆subscript𝑊0𝜋𝑀V(x,\kappa)\in SW_{[0,\pi],M}italic_V ( italic_x , italic_κ ) ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and V1(x)subscript𝑉1𝑥V_{1}(x)italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) has the opposite sign to u22(x,0)u12(x,0)superscriptsubscript𝑢22𝑥0superscriptsubscript𝑢12𝑥0u_{2}^{2}(x,0)-u_{1}^{2}(x,0)italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) a.e.

By the Feynman-Hellmann formula, at κ=0𝜅0\kappa=0italic_κ = 0,

(5) d(λ2(0)λ1(0))dκ=0π[V1(x)V(x)][u22(x,0)u12(x,0)]𝑑x.𝑑subscript𝜆20subscript𝜆10𝑑𝜅superscriptsubscript0𝜋delimited-[]subscript𝑉1𝑥subscript𝑉𝑥delimited-[]superscriptsubscript𝑢22𝑥0superscriptsubscript𝑢12𝑥0differential-d𝑥\frac{d(\lambda_{2}(0)-\lambda_{1}(0))}{d\kappa}=\int_{0}^{\pi}[V_{1}(x)-V_{*}% (x)][u_{2}^{2}(x,0)-u_{1}^{2}(x,0)]dx.divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 0 ) ) end_ARG start_ARG italic_d italic_κ end_ARG = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT [ italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) ] [ italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) ] italic_d italic_x .

Since Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is a minimizer and

[V1(x)V(x)][u22(x,0)u12(x,0)]dx0,delimited-[]subscript𝑉1𝑥subscript𝑉𝑥delimited-[]superscriptsubscript𝑢22𝑥0superscriptsubscript𝑢12𝑥0𝑑𝑥0[V_{1}(x)-V_{*}(x)][u_{2}^{2}(x,0)-u_{1}^{2}(x,0)]dx\leqslant 0,[ italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) ] [ italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) ] italic_d italic_x ⩽ 0 ,
0d(λ2(0)λ1(0))dt0,0𝑑subscript𝜆20subscript𝜆10𝑑𝑡00\leqslant\frac{d(\lambda_{2}(0)-\lambda_{1}(0))}{dt}\leqslant 0,0 ⩽ divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 0 ) ) end_ARG start_ARG italic_d italic_t end_ARG ⩽ 0 ,

which implies that the integrand in (5) equals 0 a.e. Thus V1(x)=V(x)subscript𝑉1𝑥subscript𝑉𝑥V_{1}(x)=V_{*}(x)italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) a.e. on [0,π]0𝜋[0,\pi][ 0 , italic_π ]. I.e., V(x)=V(x+)χ(a,π)subscript𝑉𝑥subscript𝑉subscript𝑥subscript𝜒𝑎𝜋V_{*}(x)=V_{*}(x_{+})\chi_{(a,\pi)}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) italic_χ start_POSTSUBSCRIPT ( italic_a , italic_π ) end_POSTSUBSCRIPT a.e. But if a>x𝑎subscript𝑥a>x_{-}italic_a > italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT the alternative choice

V1(x)={0 on (0,x)V(x+) on (x,π)subscript𝑉1𝑥cases0 on 0subscript𝑥subscript𝑉subscript𝑥 on subscript𝑥𝜋V_{1}(x)=\left\{\begin{array}[]{c}0\textrm{ on }(0,x_{-})\\ V_{*}(x_{+})\textrm{ on }(x_{-},\pi)\end{array}\right.italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL 0 on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) end_CELL end_ROW end_ARRAY

is also valid, ensuring that V(x,κ)SW[0,π],M𝑉𝑥𝜅𝑆subscript𝑊0𝜋𝑀V(x,\kappa)\in SW_{[0,\pi],M}italic_V ( italic_x , italic_κ ) ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and that V1(x)subscript𝑉1𝑥V_{1}(x)italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) has the opposite sign to u22(x,0)u12(x,0)superscriptsubscript𝑢22𝑥0superscriptsubscript𝑢12𝑥0u_{2}^{2}(x,0)-u_{1}^{2}(x,0)italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) a.e. We are thus led to the conclusion that V(x)=V(x+)χ(x,π)subscript𝑉𝑥subscript𝑉subscript𝑥subscript𝜒subscript𝑥,𝜋V_{*}(x)=V_{*}(x_{+})\chi_{(x_{,}\pi)}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) italic_χ start_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT , end_POSTSUBSCRIPT italic_π ) end_POSTSUBSCRIPT a.e. In particular, for M>0𝑀0M>0italic_M > 0 the unique transition point is a=x𝑎subscript𝑥a=x_{-}italic_a = italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT.

(ii) Secondly, suppose that a<x𝑎subscript𝑥a<x_{-}italic_a < italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT; the case when x+<asubscript𝑥𝑎x_{+}<aitalic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT < italic_a, is similar. Let

V1(x)={V(a) on (0,x)V(x+) on (x,π).subscript𝑉1𝑥casessubscript𝑉𝑎 on 0subscript𝑥subscript𝑉subscript𝑥 on subscript𝑥𝜋V_{1}(x)=\left\{\begin{array}[]{c}V_{*}(a)\textrm{ on }(0,x_{-})\\ \,\,\,\,V_{*}(x_{+})\textrm{ on }(x_{-},\pi).\end{array}\right.italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_a ) on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) . end_CELL end_ROW end_ARRAY

Then V1SW[0,π],Msubscript𝑉1𝑆subscript𝑊0𝜋𝑀V_{1}\in SW_{[0,\pi],M}italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT with

V1(x)V(x)subscript𝑉1𝑥subscript𝑉𝑥\displaystyle V_{1}(x)-V_{*}(x)italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) 0 on (0,x)(x+,π)absent0 on 0subscript𝑥subscript𝑥𝜋\displaystyle\geq 0\textrm{ on }(0,x_{-})\cup(x_{+},\pi)≥ 0 on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ∪ ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_π )
0 on (x,x+).absent0 on subscript𝑥subscript𝑥\displaystyle\leq 0\textrm{ on }(x_{-},x_{+}).≤ 0 on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) .

We note that V(x,t)SW[0,π],M𝑉𝑥𝑡𝑆subscript𝑊0𝜋𝑀V(x,t)\in SW_{[0,\pi],M}italic_V ( italic_x , italic_t ) ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT, and by the optimality of Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT,

0d(λ2(0)λ1(0))dt0.0𝑑subscript𝜆20subscript𝜆10𝑑𝑡00\leqslant\frac{d(\lambda_{2}(0)-\lambda_{1}(0))}{dt}\leqslant 0.0 ⩽ divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 0 ) ) end_ARG start_ARG italic_d italic_t end_ARG ⩽ 0 .

Hence, as before we conclude that V1(x)=V(x)subscript𝑉1𝑥subscript𝑉𝑥V_{1}(x)=V_{*}(x)italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) on [0,π]0𝜋[0,\pi][ 0 , italic_π ]. Indeed, since we have concluded that V(a)=V(x)subscript𝑉𝑎subscript𝑉subscript𝑥V_{*}(a)=V_{*}(x_{-})italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_a ) = italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ), we may as well redefine a𝑎aitalic_a as xsubscript𝑥x_{-}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT, which reduces case (ii) to case (i). In conclusion all optimal Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT must be step functions with a unique jump coinciding with either xsubscript𝑥x_{-}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT or x+subscript𝑥x_{+}italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT.

In [9], where the weight was constant, it was possible to conclude that Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT was of the form MχI(x)𝑀subscript𝜒𝐼𝑥M\chi_{I}(x)italic_M italic_χ start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ( italic_x ) for an interval I𝐼Iitalic_I by relying on the fact that adding a constant to the potential function does not change the gap ΓΓ\Gammaroman_Γ. With a variable weight we are only able to conclude that in general V=CχI(x)subscript𝑉𝐶subscript𝜒𝐼𝑥V_{*}=C\chi_{I}(x)italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT = italic_C italic_χ start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT ( italic_x ) for an undetermined constant CM𝐶𝑀C\leq Mitalic_C ≤ italic_M.

It remains to characterize the optimal weight wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT by a similar argument applied to w𝑤witalic_w for fixed V𝑉Vitalic_V. Just as before, we can use Lemma 2.2 to conclude that there exist x^±subscript^𝑥plus-or-minus\widehat{x}_{\pm}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT: 0x^<x^+π0subscript^𝑥subscript^𝑥𝜋0\leqslant\widehat{x}_{-}<\widehat{x}_{+}\leqslant\pi0 ⩽ over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ⩽ italic_π, satisfying

λ2u22(x)>λ1u12(x) on (0,x^)(x^+,π)subscript𝜆2superscriptsubscript𝑢22𝑥subscript𝜆1superscriptsubscript𝑢12𝑥 on 0subscript^𝑥subscript^𝑥𝜋\displaystyle\lambda_{2}u_{2}^{2}(x)>\lambda_{1}u_{1}^{2}(x)\textrm{ on }(0,% \widehat{x}_{-})\cup(\widehat{x}_{+},\pi)italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on ( 0 , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ∪ ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_π )
λ1u12(x)>λ2u22(x) on (x^,x^+).subscript𝜆1superscriptsubscript𝑢12𝑥subscript𝜆2superscriptsubscript𝑢22𝑥 on subscript^𝑥subscript^𝑥\displaystyle\lambda_{1}u_{1}^{2}(x)>\lambda_{2}u_{2}^{2}(x)\textrm{ on }(% \widehat{x}_{-},\widehat{x}_{+}).italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) > italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) .

Given an optimal wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT, we extend it to a family of weights by

w(x,κ)=tw1(x)+(1κ)w(x)t[0,1].formulae-sequence𝑤𝑥𝜅𝑡subscript𝑤1𝑥1𝜅subscript𝑤𝑥𝑡01w(x,\kappa)=tw_{1}(x)+(1-\kappa)w_{*}(x)\quad t\in[0,1].italic_w ( italic_x , italic_κ ) = italic_t italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) + ( 1 - italic_κ ) italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) italic_t ∈ [ 0 , 1 ] .

There are again two cases, beginning with

(i) x^b<x^+subscript^𝑥𝑏subscript^𝑥\widehat{x}_{-}\leq b<\widehat{x}_{+}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ≤ italic_b < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. As before we can arrange a convenient orientation by reflecting if necessary. We thus posit that 0<x^bx^+0subscript^𝑥𝑏subscript^𝑥0<\widehat{x}_{-}\leq b\leq\widehat{x}_{+}0 < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ≤ italic_b ≤ over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, w(x^+)w(x^)subscript𝑤subscript^𝑥subscript𝑤subscript^𝑥w_{\star}(\widehat{x}_{+})\leq w_{*}(\widehat{x}_{-})italic_w start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) ≤ italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ), and wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is nonincreasing for xx^+𝑥subscript^𝑥x\geq\widehat{x}_{+}italic_x ≥ over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. If we let

w1(x)={w(x^) on (0,b)w(x^+) on (b,π),subscript𝑤1𝑥casessubscript𝑤subscript^𝑥 on 0𝑏subscript𝑤subscript^𝑥 on 𝑏𝜋w_{1}(x)=\left\{\begin{array}[]{c}w_{*}(\widehat{x}_{-})\textrm{ on }(0,b)\\ w_{*}(\widehat{x}_{+})\textrm{ on }(b,\pi),\end{array}\right.italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) on ( 0 , italic_b ) end_CELL end_ROW start_ROW start_CELL italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( italic_b , italic_π ) , end_CELL end_ROW end_ARRAY

then w(x,κ)SB[0,π],N<,N>𝑤𝑥𝜅𝑆subscript𝐵0𝜋𝑁subscript𝑁w(x,\kappa)\in SB_{[0,\pi],N<,N_{>}}italic_w ( italic_x , italic_κ ) ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT and w1(x)subscript𝑤1𝑥w_{1}(x)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) has the opposite sign to λ1u12(x,0)λ2u22(x,0)subscript𝜆1superscriptsubscript𝑢12𝑥0subscript𝜆2superscriptsubscript𝑢22𝑥0\lambda_{1}u_{1}^{2}(x,0)-\lambda_{2}u_{2}^{2}(x,0)italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) a.e.

Fixing V𝑉Vitalic_V and differentiating the eigenvalues at κ=0𝜅0\kappa=0italic_κ = 0 by the Feynman-Hellman formula,

d(λ2λ1)dκ=0π[w1(x)w(x)][λ1u12(x,0)λ2u22(x,0)]𝑑x.𝑑subscript𝜆2subscript𝜆1𝑑𝜅superscriptsubscript0𝜋delimited-[]subscript𝑤1𝑥subscript𝑤𝑥delimited-[]subscript𝜆1superscriptsubscript𝑢12𝑥0subscript𝜆2superscriptsubscript𝑢22𝑥0differential-d𝑥\frac{d(\lambda_{2}-\lambda_{1})}{d\kappa}=\int_{0}^{\pi}[w_{1}(x)-w_{*}(x)][% \lambda_{1}u_{1}^{2}(x,0)-\lambda_{2}u_{2}^{2}(x,0)]dx.divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_ARG start_ARG italic_d italic_κ end_ARG = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT [ italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) ] [ italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) ] italic_d italic_x .

By the optimality of wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT and since

[w1(x)w(x)][λ1u12(x,0)λ2u22(x,0)]dx0,delimited-[]subscript𝑤1𝑥subscript𝑤𝑥delimited-[]subscript𝜆1superscriptsubscript𝑢12𝑥0subscript𝜆2superscriptsubscript𝑢22𝑥0𝑑𝑥0[w_{1}(x)-w_{*}(x)][\lambda_{1}u_{1}^{2}(x,0)-\lambda_{2}u_{2}^{2}(x,0)]dx% \leqslant 0,[ italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) ] [ italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) ] italic_d italic_x ⩽ 0 ,
0d(λ2(0)λ1(0))dκ00𝑑subscript𝜆20subscript𝜆10𝑑𝜅00\leqslant\frac{d(\lambda_{2}(0)-\lambda_{1}(0))}{d\kappa}\leqslant 00 ⩽ divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 0 ) ) end_ARG start_ARG italic_d italic_κ end_ARG ⩽ 0

implies that w1(x)=w(x)subscript𝑤1𝑥subscript𝑤𝑥w_{1}(x)=w_{*}(x)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) a.e. on [0,π]0𝜋[0,\pi][ 0 , italic_π ].
In conclusion the optimal wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT must be a step function with at most one jump, located at b𝑏bitalic_b. If b>x^𝑏subscript^𝑥b>\widehat{x}_{-}italic_b > over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT, then the alternative choice

w1(x)={N> on (0,x^)w(x^+) on (x^,π)subscript𝑤1𝑥casessubscript𝑁 on 0subscript^𝑥subscript𝑤subscript^𝑥 on subscript^𝑥𝜋w_{1}(x)=\left\{\begin{array}[]{c}N_{>}\textrm{ on }(0,\widehat{x}_{-})\\ w_{*}(\widehat{x}_{+})\textrm{ on }(\widehat{x}_{-},\pi)\end{array}\right.italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT on ( 0 , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) end_CELL end_ROW end_ARRAY

still ensures that w(x,κ)SW[0,π],M𝑤𝑥𝜅𝑆subscript𝑊0𝜋𝑀w(x,\kappa)\in SW_{[0,\pi],M}italic_w ( italic_x , italic_κ ) ∈ italic_S italic_W start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_M end_POSTSUBSCRIPT and w1(x)subscript𝑤1𝑥w_{1}(x)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) and has the opposite sign to λ1u12(x,0)λ2u22(x,0)subscript𝜆1superscriptsubscript𝑢12𝑥0subscript𝜆2superscriptsubscript𝑢22𝑥0\lambda_{1}u_{1}^{2}(x,0)-\lambda_{2}u_{2}^{2}(x,0)italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) - italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x , 0 ) a.e. We are thus led to the conclusion that w(x)=N>subscript𝑤𝑥subscript𝑁w_{*}(x)=N_{>}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT a.e. for x<x^𝑥subscript^𝑥x<\widehat{x}_{-}italic_x < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT and w(x)=w(x+)subscript𝑤𝑥subscript𝑤subscript𝑥w_{*}(x)=w_{*}(x_{+})italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) a.e. for x>x^𝑥subscript^𝑥x>\widehat{x}_{-}italic_x > over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT. In particular, for N>>N<subscript𝑁subscript𝑁N_{>}>N_{<}italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT > italic_N start_POSTSUBSCRIPT < end_POSTSUBSCRIPT the unique transition point is b=x^𝑏subscript^𝑥b=\widehat{x}_{-}italic_b = over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT.

(ii) Secondly, suppose that b<x^𝑏subscript^𝑥b<\widehat{x}_{-}italic_b < over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT; the case when x^+<bsubscript^𝑥𝑏\widehat{x}_{+}<bover^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT < italic_b is similar. Let

w1(x)={N> on (0,x^)w(x^+) on (x^,π).subscript𝑤1𝑥casessubscript𝑁 on 0subscript^𝑥subscript𝑤subscript^𝑥 on subscript^𝑥𝜋w_{1}(x)=\left\{\begin{array}[]{c}N_{>}\textrm{ on }(0,\widehat{x}_{-})\\ \,\,\,\,w_{*}(\widehat{x}_{+})\textrm{ on }(\widehat{x}_{-},\pi).\end{array}\right.italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT on ( 0 , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) . end_CELL end_ROW end_ARRAY

We have w1SB[0,π],N<,N>subscript𝑤1𝑆subscript𝐵0𝜋𝑁subscript𝑁w_{1}\in SB_{[0,\pi],N<,N_{>}}italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT with

w1(x)w(x)subscript𝑤1𝑥subscript𝑤𝑥\displaystyle w_{1}(x)-w_{*}(x)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) - italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) 0 on (0,x^)(x^+,π)absent0 on 0subscript^𝑥subscript^𝑥𝜋\displaystyle\leq 0\textrm{ on }(0,\widehat{x}_{-})\cup(\widehat{x}_{+},\pi)≤ 0 on ( 0 , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ∪ ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_π )
0 on (x^,x^+).absent0 on subscript^𝑥subscript^𝑥\displaystyle\geq 0\textrm{ on }(\widehat{x}_{-},\widehat{x}_{+}).≥ 0 on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) .

Noting that w(x,t)SB[0,π],N<,N>𝑤𝑥𝑡𝑆subscript𝐵0𝜋𝑁subscript𝑁w(x,t)\in SB_{[0,\pi],N<,N_{>}}italic_w ( italic_x , italic_t ) ∈ italic_S italic_B start_POSTSUBSCRIPT [ 0 , italic_π ] , italic_N < , italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_POSTSUBSCRIPT and using the optimality of wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT,

0d(λ2(0)λ1(0))dκ0.0𝑑subscript𝜆20subscript𝜆10𝑑𝜅00\leqslant\frac{d(\lambda_{2}(0)-\lambda_{1}(0))}{d\kappa}\leqslant 0.0 ⩽ divide start_ARG italic_d ( italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 ) - italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( 0 ) ) end_ARG start_ARG italic_d italic_κ end_ARG ⩽ 0 .

We conclude that w1(x)=w(x)subscript𝑤1𝑥subscript𝑤𝑥w_{1}(x)=w_{*}(x)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) = italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) a.e. on [0,π]0𝜋[0,\pi][ 0 , italic_π ]. As with the characterization of V𝑉Vitalic_V, this means that we can now suppose that b=x^𝑏subscript^𝑥b=\widehat{x}_{-}italic_b = over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT is a transition point, reducing this case to case (i).

Theorem 3.2.

The eigenvalues of the Sturm-Liouville problem (2) correspond to the real roots of the transcendental equation

ηtan(z(πx^))=ztan[η(x^x)+arctan(ηttan(tx))]ifλ>max(V)min(w),formulae-sequence𝜂𝑧𝜋subscript^𝑥𝑧𝜂subscript^𝑥subscript𝑥𝜂𝑡𝑡subscript𝑥𝑖𝑓𝜆subscript𝑉subscript𝑤\eta\tan(z(\pi-\widehat{x}_{-}))=-z\tan\left[\eta(\widehat{x}_{-}-x_{-})+% \arctan\left(\frac{\eta}{t}\tan(tx_{-})\right)\right]\leavevmode\nobreak\ % \leavevmode\nobreak\ \leavevmode\nobreak\ if\leavevmode\nobreak\ \leavevmode% \nobreak\ \lambda>\frac{\max(V_{\star})}{\min(w_{\star})},italic_η roman_tan ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) = - italic_z roman_tan [ italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + roman_arctan ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ] italic_i italic_f italic_λ > divide start_ARG roman_max ( italic_V start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ) end_ARG start_ARG roman_min ( italic_w start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ) end_ARG ,

where η:=λN>max(V)assign𝜂𝜆subscript𝑁subscript𝑉\eta:=\sqrt{\lambda N_{>}-\max(V_{\ast})}italic_η := square-root start_ARG italic_λ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT - roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) end_ARG, z:=λmin(w)max(V)assign𝑧𝜆subscript𝑤subscript𝑉z:=\sqrt{\lambda\min(w_{\ast})-\max(V_{\ast})}italic_z := square-root start_ARG italic_λ roman_min ( italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) - roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) end_ARG, and t:=λN>assign𝑡𝜆subscript𝑁t:=\sqrt{\lambda N_{>}}italic_t := square-root start_ARG italic_λ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_ARG.

Proof.

By Theorem 3.1, the optimal potential Vsubscript𝑉V_{*}italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT must be of the form

V(x)={0 on (0,x)max(V) on (x,π),subscript𝑉𝑥cases0 on 0subscript𝑥subscript𝑉 on subscript𝑥𝜋V_{*}(x)=\left\{\begin{array}[]{c}0\textrm{ on }(0,x_{-})\\ \max(V_{\ast})\textrm{ on }(x_{-},\pi),\end{array}\right.italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL 0 on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) , end_CELL end_ROW end_ARRAY

and the optimal density wsubscript𝑤w_{*}italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT must be of the form

w(x)={N> on (0,x^)min(w) on (x^,π).subscript𝑤𝑥casessubscript𝑁 on 0subscript^𝑥subscript𝑤 on subscript^𝑥𝜋w_{*}(x)=\left\{\begin{array}[]{c}N_{>}\textrm{ on }(0,\widehat{x}_{-})\\ \min(w_{\ast})\textrm{ on }(\widehat{x}_{-},\pi).\end{array}\right.italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT on ( 0 , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL roman_min ( italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) . end_CELL end_ROW end_ARRAY

The eigenfunctions are given by

u(x)={α1sin(tx) on (0,x)β1sin(η(xx))+β2cos(η(xx)) on (x,x^)α2sin(z(πx)) on (x^,π).𝑢𝑥casessubscript𝛼1𝑡𝑥 on 0subscript𝑥subscript𝛽1𝜂𝑥subscript𝑥subscript𝛽2𝜂𝑥subscript𝑥 on subscript𝑥subscript^𝑥subscript𝛼2𝑧𝜋𝑥 on subscript^𝑥𝜋u(x)=\left\{\begin{array}[]{c}\alpha_{1}\sin(tx)\textrm{ on }(0,x_{-})\\ \beta_{1}\sin(\eta(x-x_{-}))+\beta_{2}\cos(\eta(x-x_{-}))\textrm{ on }(x_{-},% \widehat{x}_{-})\\ \alpha_{2}\sin(z(\pi-x))\textrm{ on }(\widehat{x}_{-},\pi).\par\end{array}\right.italic_u ( italic_x ) = { start_ARRAY start_ROW start_CELL italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_sin ( italic_t italic_x ) on ( 0 , italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_sin ( italic_η ( italic_x - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) + italic_β start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_cos ( italic_η ( italic_x - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) on ( italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin ( italic_z ( italic_π - italic_x ) ) on ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_π ) . end_CELL end_ROW end_ARRAY

Where η=λN>max(V)𝜂𝜆subscript𝑁subscript𝑉\eta=\sqrt{\lambda N_{>}-\max(V_{\ast})}italic_η = square-root start_ARG italic_λ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT - roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) end_ARG , z=λmin(w)max(V)𝑧𝜆subscript𝑤subscript𝑉z=\sqrt{\lambda\min(w_{\ast})-\max(V_{\ast})}italic_z = square-root start_ARG italic_λ roman_min ( italic_w start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) - roman_max ( italic_V start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) end_ARG , t=λN>𝑡𝜆subscript𝑁t=\sqrt{\lambda N_{>}}italic_t = square-root start_ARG italic_λ italic_N start_POSTSUBSCRIPT > end_POSTSUBSCRIPT end_ARG and α1,α2,β1,β2subscript𝛼1subscript𝛼2subscript𝛽1subscript𝛽2\alpha_{1},\alpha_{2},\beta_{1},\beta_{2}italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_β start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are real constants, with β10subscript𝛽10\beta_{1}\neq 0italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≠ 0.

Continuity of u𝑢uitalic_u at xsubscript𝑥x_{-}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT gives

β2=α1sin(tx)subscript𝛽2subscript𝛼1𝑡subscript𝑥\beta_{2}=\alpha_{1}\sin(tx_{-})italic_β start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_sin ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT )

and continuity of usuperscript𝑢u^{\prime}italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT at xsubscript𝑥x_{-}italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT gives

β1=tα1ηcos(tx).subscript𝛽1𝑡subscript𝛼1𝜂𝑡subscript𝑥\beta_{1}=\frac{t\alpha_{1}}{\eta}\cos(tx_{-}).italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG italic_t italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG roman_cos ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) .

Then

β2β1=ηttan(tx).subscript𝛽2subscript𝛽1𝜂𝑡𝑡subscript𝑥\frac{\beta_{2}}{\beta_{1}}=\frac{\eta}{t}\tan(tx_{-}).divide start_ARG italic_β start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG = divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) .

Continuity of u𝑢uitalic_u at x^subscript^𝑥\widehat{x}_{-}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT gives

α1(tηcos(tx)sin(η(x^x))+sin(tx)cos(η(x^x)))=α2sin(z(πx^)).subscript𝛼1𝑡𝜂𝑡subscript𝑥𝜂subscript^𝑥subscript𝑥𝑡subscript𝑥𝜂subscript^𝑥subscript𝑥subscript𝛼2𝑧𝜋subscript^𝑥\alpha_{1}\left(\frac{t}{\eta}\cos(tx_{-})\sin(\eta(\widehat{x}_{-}-x_{-}))+% \sin(tx_{-})\cos(\eta(\widehat{x}_{-}-x_{-}))\right)=\alpha_{2}\sin(z(\pi-% \widehat{x}_{-})).italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( divide start_ARG italic_t end_ARG start_ARG italic_η end_ARG roman_cos ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) roman_sin ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) + roman_sin ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) roman_cos ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ) = italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) .

It follows that

tα1ηcos(tx)1+(ηttan(tx))2sin(η(x^x)+arctan(ηttan(tx)))=α2sin(z(πx^)),𝑡subscript𝛼1𝜂𝑡subscript𝑥1superscript𝜂𝑡𝑡subscript𝑥2𝜂subscript^𝑥subscript𝑥𝜂𝑡𝑡subscript𝑥subscript𝛼2𝑧𝜋subscript^𝑥\frac{t\alpha_{1}}{\eta}\cos(tx_{-})\sqrt{1+\left(\frac{\eta}{t}\tan(tx_{-})% \right)^{2}}\sin\left(\eta(\widehat{x}_{-}-x_{-})+\arctan\left(\frac{\eta}{t}% \tan(tx_{-})\right)\right)=\alpha_{2}\sin(z(\pi-\widehat{x}_{-})),divide start_ARG italic_t italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG roman_cos ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) square-root start_ARG 1 + ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG roman_sin ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + roman_arctan ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ) = italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ,

and by the continuity of usuperscript𝑢u^{\prime}italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT at x^subscript^𝑥\widehat{x}_{-}over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ,

α1η(tηcos(tx)cos(η(x^x))sin(tx)sin(η(x^x)))=α2zcos(z(πx^)).subscript𝛼1𝜂𝑡𝜂𝑡subscript𝑥𝜂subscript^𝑥subscript𝑥𝑡subscript𝑥𝜂subscript^𝑥subscript𝑥subscript𝛼2𝑧𝑧𝜋subscript^𝑥\alpha_{1}\eta\left(\frac{t}{\eta}\cos(tx_{-})\cos(\eta(\widehat{x}_{-}-x_{-})% )-\sin(tx_{-})\sin(\eta(\widehat{x}_{-}-x_{-}))\right)=-\alpha_{2}z\cos(z(\pi-% \widehat{x}_{-})).italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_η ( divide start_ARG italic_t end_ARG start_ARG italic_η end_ARG roman_cos ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) roman_cos ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) - roman_sin ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) roman_sin ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ) = - italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_z roman_cos ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) .

Therefore,

α1tcos(tx)1+(ηttan(tx))2cos(η(x^x)+arctan(ηttan(tx)))=α2zcos(z(πx^)).subscript𝛼1𝑡𝑡subscript𝑥1superscript𝜂𝑡𝑡subscript𝑥2𝜂subscript^𝑥subscript𝑥𝜂𝑡𝑡subscript𝑥subscript𝛼2𝑧𝑧𝜋subscript^𝑥\alpha_{1}t\cos(tx_{-})\sqrt{1+\left(\frac{\eta}{t}\tan(tx_{-})\right)^{2}}% \cos\left(\eta(\widehat{x}_{-}-x_{-})+\arctan\left(\frac{\eta}{t}\tan(tx_{-})% \right)\right)=-\alpha_{2}z\cos(z(\pi-\widehat{x}_{-})).italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_t roman_cos ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) square-root start_ARG 1 + ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG roman_cos ( italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + roman_arctan ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ) = - italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_z roman_cos ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) .

Then

ηtan(z(πx^))=ztan[η(x^x)+arctan(ηttan(tx))].𝜂𝑧𝜋subscript^𝑥𝑧𝜂subscript^𝑥subscript𝑥𝜂𝑡𝑡subscript𝑥\eta\tan(z(\pi-\widehat{x}_{-}))=-z\tan\left[\eta(\widehat{x}_{-}-x_{-})+% \arctan\left(\frac{\eta}{t}\tan(tx_{-})\right)\right].italic_η roman_tan ( italic_z ( italic_π - over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) = - italic_z roman_tan [ italic_η ( over^ start_ARG italic_x end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) + roman_arctan ( divide start_ARG italic_η end_ARG start_ARG italic_t end_ARG roman_tan ( italic_t italic_x start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ) ] .

This ends the proof. ∎

4. Liouville transform of Sturm-Liouville operators

In this section we apply Lavine’s estimate on the fundamental gap to the Sturm-Liouville equation (2),

u′′+V(x)u=λw(x)u,x[0,π].superscript𝑢′′𝑉𝑥𝑢𝜆𝑤𝑥𝑢𝑥0𝜋\displaystyle\begin{array}[]{lr}-u^{\prime\prime}+V(x)u=\lambda w(x)u,&x\in[0,% \pi].\\ \end{array}start_ARRAY start_ROW start_CELL - italic_u start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT + italic_V ( italic_x ) italic_u = italic_λ italic_w ( italic_x ) italic_u , end_CELL start_CELL italic_x ∈ [ 0 , italic_π ] . end_CELL end_ROW end_ARRAY

The eigenvalues of (2) coincide with these of the corresponding eigenvalue problem in Liouville normal form,

(6) d2ηdξ2+(λψ(ξ))η=0on[0,L],superscript𝑑2𝜂𝑑superscript𝜉2𝜆𝜓𝜉𝜂0𝑜𝑛0𝐿\displaystyle\frac{d^{2}\eta}{d\xi^{2}}+(\lambda-\psi(\xi))\eta=0\leavevmode% \nobreak\ \leavevmode\nobreak\ on\leavevmode\nobreak\ \leavevmode\nobreak\ [0,% L],divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_η end_ARG start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + ( italic_λ - italic_ψ ( italic_ξ ) ) italic_η = 0 italic_o italic_n [ 0 , italic_L ] ,

where ψ(ξ)𝜓𝜉\psi(\xi)italic_ψ ( italic_ξ ) is the Liouville potential defined by

(7) ψ(ξ)=w′′4w25(w)216w3+Vw,𝜓𝜉superscript𝑤′′4superscript𝑤25superscriptsuperscript𝑤216superscript𝑤3𝑉𝑤\displaystyle\psi(\xi)=\frac{w^{\prime\prime}}{4w^{2}}-\frac{5(w^{\prime})^{2}% }{16w^{3}}+\frac{V}{w},italic_ψ ( italic_ξ ) = divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 5 ( italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_V end_ARG start_ARG italic_w end_ARG ,

with L=0πw(t)𝑑t𝐿superscriptsubscript0𝜋𝑤𝑡differential-d𝑡L=\int_{0}^{\pi}\sqrt{w(t)}dtitalic_L = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t. (For background on the Liouville transform we refer to [22]). In particular we have::::

Γ[V,w]=Γ[ψ].Γ𝑉𝑤Γdelimited-[]𝜓\Gamma[V,w]=\Gamma[\psi].roman_Γ [ italic_V , italic_w ] = roman_Γ [ italic_ψ ] .
Proposition 4.1.

If the Liouville potential ψ𝜓\psiitalic_ψ given by (7) of the Sturm-Liouville problem (6) is convex, then

Γ[V,w]3π2(0πw(t)𝑑t)2,Γ𝑉𝑤3superscript𝜋2superscriptsuperscriptsubscript0𝜋𝑤𝑡differential-d𝑡2\Gamma[V,w]\geqslant\frac{3\pi^{2}}{(\int_{0}^{\pi}\sqrt{w(t)}dt)^{2}},roman_Γ [ italic_V , italic_w ] ⩾ divide start_ARG 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,

and equality is obtained if and only if ψ𝜓\psiitalic_ψ is constant.

Proof.

If the Liouville potential ψ𝜓\psiitalic_ψ of the Sturm-Liouville problem (6) is convex, then by [18], we have

Γ[V,w]3π2L2.Γ𝑉𝑤3superscript𝜋2superscript𝐿2\Gamma[V,w]\geqslant\frac{3\pi^{2}}{L^{2}}.roman_Γ [ italic_V , italic_w ] ⩾ divide start_ARG 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

Because the interval [0,π]0𝜋[0,\pi][ 0 , italic_π ] is tranformed to [0,L]0𝐿[0,L][ 0 , italic_L ] with L=0πw(t)𝑑t𝐿superscriptsubscript0𝜋𝑤𝑡differential-d𝑡L=\int_{0}^{\pi}\sqrt{w(t)}dtitalic_L = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t. Thus

Γ[V,w]3π2(0πw(t)𝑑t)2.Γ𝑉𝑤3superscript𝜋2superscriptsuperscriptsubscript0𝜋𝑤𝑡differential-d𝑡2\Gamma[V,w]\geqslant\frac{3\pi^{2}}{(\int_{0}^{\pi}\sqrt{w(t)}dt)^{2}}.roman_Γ [ italic_V , italic_w ] ⩾ divide start_ARG 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

Remark 4.1.

To see that convexity of V𝑉Vitalic_V is not required for Proposition 4.1, consider the example V(x)=x2𝑉𝑥superscript𝑥2V(x)=-x^{2}italic_V ( italic_x ) = - italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Let w(x)=x2const𝑤𝑥superscript𝑥2constw(x)=x^{2}\neq\textrm{const}italic_w ( italic_x ) = italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≠ const, then ψ′′(ξ)=3x225x60superscript𝜓′′𝜉3superscript𝑥225superscript𝑥60\psi^{{}^{\prime\prime}}(\xi)=\frac{3x^{2}-25}{x^{6}}\geq 0italic_ψ start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ( italic_ξ ) = divide start_ARG 3 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 25 end_ARG start_ARG italic_x start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG ≥ 0 on [5,6]56[5,6][ 5 , 6 ]. On this interval the Liouville potential ψ𝜓\psiitalic_ψ is convex and hence Proposition 4.1 is applicable, so

Γ[V,w]3π2(56w(t)𝑑t)2,Γ𝑉𝑤3superscript𝜋2superscriptsuperscriptsubscript56𝑤𝑡differential-d𝑡2\Gamma[V,w]\geqslant\frac{3\pi^{2}}{(\int_{5}^{6}\sqrt{w(t)}dt)^{2}},roman_Γ [ italic_V , italic_w ] ⩾ divide start_ARG 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( ∫ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,

i.e.

Γ[V,w]0.978803.Γ𝑉𝑤0.978803\Gamma[V,w]\geqslant 0.978803\dots.roman_Γ [ italic_V , italic_w ] ⩾ 0.978803 … .
Proposition 4.2.

Consider the Sturm-Liouville problem (2) with positive density function wC2(0,π)𝑤superscript𝐶20𝜋w\in C^{2}(0,\pi)italic_w ∈ italic_C start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 0 , italic_π ) and continuous convex potential V𝑉Vitalic_V on [0,π]0𝜋[0,\pi][ 0 , italic_π ]. If the fundamental gap satisfies

Γ[V,w]=3π2(0πw(t)𝑑t)2Γ𝑉𝑤3superscript𝜋2superscriptsuperscriptsubscript0𝜋𝑤𝑡differential-d𝑡2\Gamma[V,w]=3\pi^{2}\left(\int_{0}^{\pi}\sqrt{w(t)}dt\right)^{-2}roman_Γ [ italic_V , italic_w ] = 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t ) start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT

Then

w(3)4w2w′′w2w310w(3)16w3+15w316w4+VwVww2=0.superscript𝑤34superscript𝑤2superscript𝑤′′superscript𝑤2superscript𝑤310superscript𝑤316superscript𝑤315superscript𝑤316superscript𝑤4superscript𝑉𝑤𝑉superscript𝑤superscript𝑤20\frac{w^{(3)}}{4w^{2}}-\frac{w^{\prime\prime}w^{\prime}}{2w^{3}}-\frac{10w^{(3% )}}{16w^{3}}+\frac{15w^{\prime 3}}{16w^{4}}+\frac{V^{\prime}}{w}-\frac{Vw^{% \prime}}{w^{2}}=0.divide start_ARG italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 10 italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 15 italic_w start_POSTSUPERSCRIPT ′ 3 end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w end_ARG - divide start_ARG italic_V italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = 0 .

for all x[0,π]𝑥0𝜋x\in[0,\pi]italic_x ∈ [ 0 , italic_π ].

Proof.

If the fundamental gap satisfies

Γ[V,w]=3π2(0πw(t)𝑑t)2Γ𝑉𝑤3superscript𝜋2superscriptsuperscriptsubscript0𝜋𝑤𝑡differential-d𝑡2\Gamma[V,w]=3\pi^{2}\left(\int_{0}^{\pi}\sqrt{w(t)}dt\right)^{-2}roman_Γ [ italic_V , italic_w ] = 3 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT square-root start_ARG italic_w ( italic_t ) end_ARG italic_d italic_t ) start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT

then the Liouville potential ψ𝜓\psiitalic_ψ is constant, in which case ψ=0superscript𝜓0\psi^{\prime}=0italic_ψ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 0. Therefore

w(3)4w2w′′w2w310w(3)16w3+15w316w4+VwVww2=0.superscript𝑤34superscript𝑤2superscript𝑤′′superscript𝑤2superscript𝑤310superscript𝑤316superscript𝑤315superscript𝑤316superscript𝑤4superscript𝑉𝑤𝑉superscript𝑤superscript𝑤20\frac{w^{(3)}}{4w^{2}}-\frac{w^{\prime\prime}w^{\prime}}{2w^{3}}-\frac{10w^{(3% )}}{16w^{3}}+\frac{15w^{\prime 3}}{16w^{4}}+\frac{V^{\prime}}{w}-\frac{Vw^{% \prime}}{w^{2}}=0.divide start_ARG italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 10 italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 15 italic_w start_POSTSUPERSCRIPT ′ 3 end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w end_ARG - divide start_ARG italic_V italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = 0 .

Proposition 4.3.

The Liouville potential ψ𝜓\psiitalic_ψ corresponding to (2) is convex if

w(4)4w2superscript𝑤44superscript𝑤2\displaystyle\frac{w^{(4)}}{4w^{2}}divide start_ARG italic_w start_POSTSUPERSCRIPT ( 4 ) end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG w(3)ww3w′′22w3+3w′′w22w41016[w(4)w33w(3)ww4]superscript𝑤3superscript𝑤superscript𝑤3superscript𝑤′′22superscript𝑤33superscript𝑤′′superscript𝑤22superscript𝑤41016delimited-[]superscript𝑤4superscript𝑤33superscript𝑤3superscript𝑤superscript𝑤4\displaystyle-\frac{w^{(3)}w^{\prime}}{w^{3}}-\frac{w^{\prime\prime 2}}{2w^{3}% }+\frac{3w^{\prime\prime}w^{\prime 2}}{2w^{4}}-\frac{10}{16}\left[\frac{w^{(4)% }}{w^{3}}-\frac{3w^{(3)}w^{\prime}}{w^{4}}\right]- divide start_ARG italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 10 end_ARG start_ARG 16 end_ARG [ divide start_ARG italic_w start_POSTSUPERSCRIPT ( 4 ) end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ]
+1516[3w2w′′w44w4w5]+V′′wVww2Vw+Vw′′w2+2w2Vw301516delimited-[]3superscript𝑤2superscript𝑤′′superscript𝑤44superscript𝑤4superscript𝑤5superscript𝑉′′𝑤superscript𝑉superscript𝑤superscript𝑤2superscript𝑉superscript𝑤𝑉superscript𝑤′′superscript𝑤22superscript𝑤2𝑉superscript𝑤30\displaystyle+\frac{15}{16}\left[\frac{3w^{\prime 2}w^{\prime\prime}}{w^{4}}-% \frac{4w^{\prime 4}}{w^{5}}\right]+\frac{V^{\prime\prime}}{w}-\frac{V^{\prime}% w^{\prime}}{w^{2}}-\frac{V^{\prime}w^{\prime}+Vw^{\prime\prime}}{w^{2}}+\frac{% 2w^{\prime 2}V}{w^{3}}\geq 0+ divide start_ARG 15 end_ARG start_ARG 16 end_ARG [ divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 4 italic_w start_POSTSUPERSCRIPT ′ 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG ] + divide start_ARG italic_V start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w end_ARG - divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_V italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 2 italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT italic_V end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ≥ 0

on [0,π]0𝜋[0,\pi][ 0 , italic_π ].

Proof

To analyze the convexity of ψ𝜓\psiitalic_ψ, we will apply the nonnegativity criterion for the second derivative. To differentiate the Liouville potential we make use of the chain rule. In particular,

dψdξ=gdψdx.𝑑𝜓𝑑𝜉𝑔𝑑𝜓𝑑𝑥\frac{d\psi}{d\xi}=\sqrt{g}\frac{d\psi}{dx}.divide start_ARG italic_d italic_ψ end_ARG start_ARG italic_d italic_ξ end_ARG = square-root start_ARG italic_g end_ARG divide start_ARG italic_d italic_ψ end_ARG start_ARG italic_d italic_x end_ARG .

If g=1𝑔1g=1italic_g = 1 then ψ𝜓\psiitalic_ψ is convex, so d2ψdx20superscript𝑑2𝜓𝑑superscript𝑥20\displaystyle\frac{d^{2}\psi}{dx^{2}}\geqslant 0divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ end_ARG start_ARG italic_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ⩾ 0.
As a consequence

dψdξ=w(3)4w2w′′w2w310w(3)16w3+15w316w4+VwVww2.𝑑𝜓𝑑𝜉superscript𝑤34superscript𝑤2superscript𝑤′′superscript𝑤2superscript𝑤310superscript𝑤316superscript𝑤315superscript𝑤316superscript𝑤4superscript𝑉𝑤𝑉superscript𝑤superscript𝑤2\frac{d\psi}{d\xi}=\frac{w^{(3)}}{4w^{2}}-\frac{w^{\prime\prime}w^{\prime}}{2w% ^{3}}-\frac{10w^{(3)}}{16w^{3}}+\frac{15w^{\prime 3}}{16w^{4}}+\frac{V^{\prime% }}{w}-\frac{Vw^{\prime}}{w^{2}}.divide start_ARG italic_d italic_ψ end_ARG start_ARG italic_d italic_ξ end_ARG = divide start_ARG italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 10 italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 15 italic_w start_POSTSUPERSCRIPT ′ 3 end_POSTSUPERSCRIPT end_ARG start_ARG 16 italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w end_ARG - divide start_ARG italic_V italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

This yields that

d2ψdξ2superscript𝑑2𝜓𝑑superscript𝜉2\displaystyle\frac{d^{2}\psi}{d\xi^{2}}divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ end_ARG start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG =w(4)4w2w(3)ww3w′′22w3+3w′′w22w41016[w(4)w33w(3)ww4]absentsuperscript𝑤44superscript𝑤2superscript𝑤3superscript𝑤superscript𝑤3superscript𝑤′′22superscript𝑤33superscript𝑤′′superscript𝑤22superscript𝑤41016delimited-[]superscript𝑤4superscript𝑤33superscript𝑤3superscript𝑤superscript𝑤4\displaystyle=\frac{w^{(4)}}{4w^{2}}-\frac{w^{(3)}w^{\prime}}{w^{3}}-\frac{w^{% \prime\prime 2}}{2w^{3}}+\frac{3w^{\prime\prime}w^{\prime 2}}{2w^{4}}-\frac{10% }{16}\left[\frac{w^{(4)}}{w^{3}}-\frac{3w^{(3)}w^{\prime}}{w^{4}}\right]= divide start_ARG italic_w start_POSTSUPERSCRIPT ( 4 ) end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_w start_POSTSUPERSCRIPT ′ ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 10 end_ARG start_ARG 16 end_ARG [ divide start_ARG italic_w start_POSTSUPERSCRIPT ( 4 ) end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ]
+1516[3w2w′′w44w4w5]+V′′wVww2Vw+Vw′′w2+2w2Vw3.1516delimited-[]3superscript𝑤2superscript𝑤′′superscript𝑤44superscript𝑤4superscript𝑤5superscript𝑉′′𝑤superscript𝑉superscript𝑤superscript𝑤2superscript𝑉superscript𝑤𝑉superscript𝑤′′superscript𝑤22superscript𝑤2𝑉superscript𝑤3\displaystyle\quad+\frac{15}{16}\left[\frac{3w^{\prime 2}w^{\prime\prime}}{w^{% 4}}-\frac{4w^{\prime 4}}{w^{5}}\right]+\frac{V^{\prime\prime}}{w}-\frac{V^{% \prime}w^{\prime}}{w^{2}}-\frac{V^{\prime}w^{\prime}+Vw^{\prime\prime}}{w^{2}}% +\frac{2w^{\prime 2}V}{w^{3}}.+ divide start_ARG 15 end_ARG start_ARG 16 end_ARG [ divide start_ARG 3 italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 4 italic_w start_POSTSUPERSCRIPT ′ 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG ] + divide start_ARG italic_V start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w end_ARG - divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_V start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_w start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_V italic_w start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 2 italic_w start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT italic_V end_ARG start_ARG italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG .

References

  • [1] M. Ahrami and Z. El Allali, Lower bounds on the fundamental spectral gap with Robin boundary conditions, 2021 UNC Greensboro PDE Conference. Electron. J. Diff. Eqns. Conf. 26 (2022), pp 1-11.
  • [2] M. Ashbaugh and R. Benguria, Optimal lower bound for the gap between the first two eigenvalues of one-dimensional Schrödinger operators with symmetric single-well potentials, Proc. Amer. Math. Soc. 105 (1989), 419–424.
  • [3] M. S. Ashbaugh and R. Svirsky, Periodic potentials with minimal energy bands, Proc. Amer. Math. Soc., 114 (1992) 69–77.
  • [4] M. Ashbaugh and R. Benguria, Optimal bounds for ratios of eigenvalues of one-dimensional Schrödinger operators with Dirichlet boundary conditions and positive potentials, Comm. Math. Phys., 124 (1989), 403–415.
  • [5] M. S. Ashbaugh, E. M. Harrell II, and R. Svirsky, On minimal and maximal eigenvalues gaps and their causes, Pac. J. Math. 147 (1991) 1–24.
  • [6] Y. H. Cheng, S. Y. Kung, C. K. Law and W. C. Lian, The dual eigenvalue problems for the Sturm-Liouville system, Computers and Mathematics with Applications, 60 (2010) 2556–2563.
  • [7] M. S. Ashbaugh and E. M. Harrell II, Perturbation theory for shape resonances and large barrier potentials, Commun. Math. Phys. 83 (1982) 151–170.
  • [8] J. L. Doob, Measure theory, Graduate Texts in Mathematics (143), 1994 Springer.
  • [9] Z. El Allali and E. M. Harrell II , Optimal bounds on the fundamental spectral gap with single-well potentials, Proc. Amer. Math. Soc., 150,(2022), 57–587.
  • [10] G. Bognàr and O. Dosly, The ratio of eigenvalues of the Dirichlet eigenvalue problem for equations with one-dimensional p-Laplacian, Abstract and Applied Analyis, 2010 (2010).
  • [11] T. Kato, Perturbation theory for linear operator, 1980 Springer-Verlag.
  • [12] M. J. Huang, On the eigenvalue ratio for vibrating strings, Proc. Amer. Math. Soc., 127 (1999) 1805–1813.
  • [13] M. Horváth, On the first two eigenvalues of Sturm-Liouville operators, Proc. Amer. Math. Soc. 131 4 (2002) 1215–1224.
  • [14] M. Horváth and M. Kiss, A bound for ratios of eigenvalues of Schrödinger operators with single-well potentials, Proc. Amer. Math. Soc., 134 (2006) 1425–1434.
  • [15] M. J. Huang, The eigenvalue gap for vibrating strings with symmetric densities, Acta Math. Hungar., 117 (2007) 341–348.
  • [16] M. J. Huang, A note on the eigenvalues ratio of vibrating strings, Acta Math. Hungar., 123 (2009) 265–271.
  • [17] M. J. Huang and C. K. Law, Eigenvalue ratios for the regular Sturm-Liouville system, Proc. Amer. Math. Soc., 124 (1996) 1427–1436.
  • [18] R. Lavine, The eigenvalue gap for one-dimensional convex potentials, Proceedings of the American Mathematical Society 1994; 121, 815–821.
  • [19] R. Svirsky, Maximally resonant potentials subject to p-Norm constraints, Pac. J. Math. 129, 357–374 (1987).
  • [20] X. J. Yu and C. F. Yang, The gap between the first two eigenvalues of Schrödinger operators with single-well potential, Appl. Math. Comp. 268 (2015) 275–283.
  • [21] J. Weidmann, Linear operators in Hilbert spaces, volume 68 of Graduate Texts in Mathematics. Springer-Verlag, New York-Berlin, 1980. Translated from the German by Joseph Szücs.
  • [22] Garrett Birkhoff and Gian-Carlo Rota, Ordinary Differential Equations, 4th Edition. New York: Wiley, 1989. See Chapter 10, section 9.
  • [23] A. Zettl Strum-Liouville theory, American Mathematical Society, 2005.