Asymptotic formulas of the eigenvalues for the linearization of the scalar field equation

Abstract

We establish asymptotic formulas for all the eigenvalues of the linearization problem of the Neumann problem for the scalar field equation in a finite interval

\[ \begin{cases} \varepsilon^2u_{xx}-u+u^3=0, & 0< x<1,\\ u_x(0)=u_x(1)=0. \end{cases} \]

In the previous paper of the third author [T. Wakasa and S. Yotsutani, J. Differ. Equ. 258 (2015), 3960–4006] asymptotic formulas for the Allen–Cahn case $\varepsilon ^2u_{xx}+u-u^3=0$ were established. In this paper, we apply the method developed in the previous paper to our case. We show that all the eigenvalues can be classified into three categories, i.e., near $-3$ eigenvalues, near $0$ eigenvalues and the other eigenvalues. We see that the number of the near $-3$ eigenvalues (resp. the near $0$ eigenvalues) is equal to the number of the interior and boundary peaks (resp. the interior peaks) of a solution for the nonlinear problem. The main technical tools are various asymptotic formulas for complete elliptic integrals.

1. Introduction and main results

We are concerned with the Neumann problem of a reaction–diffusion equation

(1.1)$$ \begin{cases} \varepsilon^2u_{xx}(x)+f(u(x))=0, & 0< x<1,\\ u_x(0)=u_x(1)=0 \end{cases} $$

and its linearized eigenvalue problem associated with a positive solution $u(x)$

(1.2)$$ \begin{cases} \varepsilon^2\varphi_{xx}(x)+f_u(u(x))\varphi(x)={-}\lambda\varphi(x), & 0< x<1,\\ \varphi_x(0)=\varphi_x(1)=0, \end{cases} $$

where $\varepsilon >0$ is a parameter and $f$ is smooth enough. We denote $F(u)$ by

\[ F(u)=\int_0^uf(s)\,{\rm d}s. \]

Precise information of eigenvalues is crucial not only in a stability analysis of a stationary solution of the associated parabolic problem but also in a study of dynamics. The goal of the present paper is to obtain a rather explicit expression of all the eigenvalues of (1.2) in the case

\[ f(u)={-}u+u^3. \]

The equation $\varepsilon ^2u_{xx}-u+u^p=0$, $p>1$, is called the scalar field equation which appears in the study of a standing wave of a nonlinear Schrödinger equation [2], of elastic curves in differential geometry [6] and of the Gierer–Meinhardt model in mathematical biology [5]. This equation has attracted much attention for these three decades. A complete bifurcation diagram of the positive solutions for (1.1) with $f(u)=-u+u^p$, $p>1$, was obtained in [7, 16]. See Fig. 1. Let $p>1$ and $\varepsilon _n:=\sqrt {p-1}/n\pi$. For each $n\ge 1$, (1.1) has exactly two $n$-mode solutions $u^{\pm }_{n,\varepsilon }$ if $0<\varepsilon <\varepsilon _n$. Here,

\[ {u^+_{n,\varepsilon}} (0)=\min_{0\le x\le 1}{u^+_{n,\varepsilon}} (x),\quad {u^-_{n,\varepsilon}} (0)=\max_{0\le x\le 1}{u^-_{n,\varepsilon}} (x). \]

In particular, if $n=2k+1$, then both $u^+_{n,\varepsilon }$ and $u^-_{u,\varepsilon }$ have $k$ interior peaks and $1$ boundary peak. On the other hand, if $n=2k$, then $u^+_{n,\varepsilon }$ (resp. $u^-_{n,\varepsilon }$) has $k$ interior peaks and no boundary peak (resp. $k-1$ interior peaks and $2$ boundary peaks). Each peak becomes sharp as $\varepsilon \to 0$.

Figure 1. The complete bifurcation diagram for (1.1) with $f(u)=-u+u^3$.

Hereafter, we consider the case $f(u)=-u+u^3$. Then

\[ \varepsilon_n=\frac{\sqrt{2}}{n\pi}. \]

All the solutions of (1.1) can be written in terms of elliptic functions:

Proposition 1.1 Let $n\ge 1$ be fixed and let $0<\varepsilon <\varepsilon _n$. Let $k_{\varepsilon }=k_{n,\varepsilon }\in (0,1)$ be a solution of

(1.3)$$ \sqrt{2-k^2}K(k)=\frac{1}{n\varepsilon}. $$

Note that $k_{\varepsilon }$ is uniquely determined. Then ${u^+_{n,\varepsilon }}$ and ${u^-_{n,\varepsilon }}$ can be written as follows:

(1.4)$$ u^{{\pm}}_{n,\varepsilon}(x)=\sqrt{\frac{2}{2-k_{\varepsilon}^2}}{\rm DN}_n^{{\pm}}(x,k_{\varepsilon}), $$

where

\[ {\rm DN}_n^+(x,k_{\varepsilon}):={\rm dn}(K(k_{\varepsilon})(1+nx),k_{\varepsilon}),\quad {\rm DN}_n^-(x,k_{\varepsilon}):={\rm dn}(nK(k_{\varepsilon})x,k_{\varepsilon}). \]

The proof of proposition 1.1 is given in §2.

In this paper, we frequently use Jacobi's elliptic functions ${\rm sn}(x,k)$ and ${\rm dn}(x,k)$ and complete elliptic integrals $K(k)$, $E(k)$ and $\Pi (\nu,k)$. The definitions and various properties are summarized in Appendix A of the present paper.

The linearized eigenvalue problem (1.2) associated with $u^{\pm }_{n,\varepsilon }$ is as follows:

(LP_±)$$ \begin{cases} {\varepsilon }^2\varphi _{xx}(x) + f_u(u_{n,\varepsilon }^ \pm (x))\varphi (x) = -\lambda \varphi (x), & {0 < x < 1,} \\ \varphi _x(0) = \varphi _x(1) = 0. \end{cases} $$

Hereafter, $\lambda ^{\pm }_{j,\varepsilon }$, $j\ge 0$ denotes the $j+1$-th eigenvalue of (LP_±) and $\varphi ^{\pm }_{j,\varepsilon }$ denotes an associated eigenfunction. It is well known that each eigenvalue $\lambda ^{\pm }_{j,\varepsilon }$ is simple.

The following three exact eigenvalues of (LP_±) were obtained in the previous paper of the third author [11]:

Proposition 1.2 Let $n\ge 1$,

$$\begin{align*} & {\rm SN}_n^+(x,k):={\rm sn}(K(k)(1+nx),k),\quad {\rm SN}_n^-(x,k):={\rm sn}(nK(k)x,k),\\ & {\rm CN}_n^+(x,k):={\rm cn}(K(k)(1+nx),k),\quad {\rm CN}_n^-(x,k):={\rm cn}(nK(k)x,k). \end{align*}$$

The following hold$:$

1. (i) The problem (${\rm LP}_+$) has the following two pairs of eigenvalues and eigenfunctions$:$

   $$\begin{align*} & \lambda^+_{0,\varepsilon}={-}1-\frac{2\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}}{2-k_{\varepsilon}^2},\quad \varphi^+_{0,\varepsilon}(x)=1\\& \quad -\left(1+k_{\varepsilon}^2-\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}\right){\rm SN}_n^+(x,k_{\varepsilon})^2,\\ & \lambda^+_{n,\varepsilon}={-}1+\frac{2\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}}{2-k_{\varepsilon}^2},\quad \varphi^+_{n,\varepsilon}(x)={-}1\\& \quad +\left(1+k_{\varepsilon}^2+\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}\right){\rm SN}_n^+(x,k_{\varepsilon})^2. \end{align*}$$
   Moreover, if $n$ is even, then (LP$_+$) has one more pair$:$
   \[ \lambda^+_{n/2,\varepsilon}={-}\frac{3(1-k_{\varepsilon}^2)}{2-k_{\varepsilon}^2},\quad \varphi^+_{n/2,\varepsilon}(x)={\rm SN}_n^+(x,k_{\varepsilon}){\rm DN}_n^+(x,k_{\varepsilon}). \]

2. (ii) The problem (${\rm LP}_-$) has the following pairs of eigenvalues and eigenfunctions$:$

   $$\begin{align*} & \lambda^-_{0,\varepsilon}={-}1-\frac{2\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}}{2-k_{\varepsilon}^2},\quad \varphi^-_{0,\varepsilon}(x)=1\\& \quad-\left(1+k_{\varepsilon}^2-\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}\right){\rm SN}_n^-(x,k_{\varepsilon})^2,\\ & \lambda^-_{n,\varepsilon}={-}1+\frac{2\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}}{2-k_{\varepsilon}^2},\quad \varphi^-_{n,\varepsilon}(x)={-}1\\& \quad +\left(1+k_{\varepsilon}^2+\sqrt{1-k_{\varepsilon}^2+k_{\varepsilon}^4}\right){\rm SN}_n^-(x,k_{\varepsilon})^2. \end{align*}$$
   Moreover, if $n$ is even, then (LP$_-$) has one more pair$:$
   \[ \lambda^-_{n/2,\varepsilon}={-}\frac{3}{2-k_{\varepsilon}^2},\quad \varphi^-_{n/2,\varepsilon}(x)={\rm CN}_n^-(x,k_{\varepsilon}){\rm DN}_n^-(x,k_{\varepsilon}). \]

Singularly perturbed problems have attracted great attention, because of a wide variety of applications. It follows from (1.3) that a singularly perturbed problem corresponds to the case where $k$ is close to $1$. The main result of the paper is the following asymptotic formula of the other eigenvalues as $k\to 1:$

Theorem 1.3 Let $n\ge 1$ and $j>0$ be fixed. Let $k_{\varepsilon }$ be as in proposition 1.1. The following hold for (${\rm LP}_{\pm }$)$:$

1. (i) For $0< j< n/2$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}={-}3+3(1-k_{\varepsilon}^2)+\frac{3}{4}\left(\sin^2\frac{j\pi}{n}-5\right)(1-k_{\varepsilon}^2)^2\!+\!o\left((1-k_{\varepsilon}^2)^2\right)\ \textrm{as}\ k_{\varepsilon}\!\to\! 1. \]

2. (ii) For $n/2< j< n$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}={-}3\left(\sin^2\frac{j\pi}{n}\right)(1-k_{\varepsilon}^2)+o(1-k_{\varepsilon}^2)\ \textrm{as}\ k_{\varepsilon}\to 1. \]

3. (iii) For $j>n$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}=1+\frac{(j-n)^2\pi^2}{n^2}\frac{1}{K(k_{\varepsilon})^2}+o\left(\frac{1}{K(k_{\varepsilon})^2}\right)\ \textrm{as}\ k_{\varepsilon}\to 1, \]
   where $K(k_{\varepsilon })\to \infty$ as $k_{\varepsilon }\to 1$.

Combining proposition 1.2 and theorem 1.3, we can relate locations of eigenvalues with the number of the peaks of $u^{\pm }_{n,\varepsilon }$ as follows: All the eigenvalues can be classified into three categories, i.e.,

(1.5)$$\begin{cases} \textrm{(a) near ${-}3$ eigenvalues which converge to ${-}3$ as $k_{\varepsilon}\to 1$,}\\ \textrm{(b) near $0$ eigenvalues which converge to $0$ as $k_{\varepsilon}\to 1$,}\\ \textrm{(c) the other eigenvalues which converge to $1$ as $k_{\varepsilon}\to 1$.} \end{cases} $$

The number of the near $-3$ eigenvalues is equal to the number of the interior and boundary peaks of $u^{\pm }_{n,\varepsilon }$. The number of the near $0$ eigenvalues is equal to the number of the interior peaks of $u^{\pm }_{n,\varepsilon }$.

Example 1.4

1. (i) The solution $u^+_{10,\varepsilon }$ has $5$ interior peaks and no boundary peak. Hence, (${\rm LP}_+$) has $5$ eigenvalues near $-3$ and $5$ eigenvalues near $0$ when $\varepsilon >0$ is small. On the other hand, $u^-_{10,\varepsilon }$ has $4$ interior peaks and 2 boundary peaks. Hence (${\rm LP}_-$) has 6 eigenvalues near $-3$ and 4 eigenvalues near $0$ when $\varepsilon >0$ is small.

2. (ii) Each of $u^+_{11,\varepsilon }$ and $u^-_{11,\varepsilon }$ has $5$ interior peaks and 1 boundary peak. Hence, (${\rm LP}_{\pm }$) has 6 eigenvalues near $-3$ and 5 eigenvalues near $0$ when $\varepsilon >0$ is small.

   In figures 2–4 profiles of eigenfunctions on (${\rm LP}_\pm$) with $u^\pm _{10,\varepsilon }$ and profiles of eigenfunctions on (${\rm LP}_+$) with $u^+_{11,\varepsilon }$ are shown. In each figure (a-0), (a-1) and (a-2) belong to category (a) in (LP_±), and so on.

Figure 2. Profiles of eigenfunctions for (${\rm LP}_+$) with $u_{10,\varepsilon }^+(x)$; (a-0) $\varphi ^+_0$, (a-1) $\varphi ^+_1$, (a-2) $\varphi ^+_2$, (b-0) $\varphi ^+_5$, (b-1) $\varphi ^+_6$, (b-2) $\varphi ^+_7$, (c-0) $\varphi ^+_{10}$, (c-1) $\varphi ^+_{11}$ and (c-2) $\varphi ^+_{12}$ .

Figure 3. Profiles of eigenfunctions for (LP$_-$) with $u_{10,\varepsilon }^-(x)$; (a-0) $\varphi ^-_0$, (a-1) $\varphi ^-_1$, (a-2) $\varphi ^-_2$, (a-5) $\varphi ^-_5$, (b-1) $\varphi ^-_6$, (b-2) $\varphi ^-_7$, (c-0) $\varphi ^-_{10}$, (c-1) $\varphi ^-_{11}$ and (c-2) $\varphi ^-_{12}$.

Figure 4. Profiles of eigenfunctions for (${\rm LP}_+$) with $u_{11,\varepsilon }^+(x)$; (a-0) $\varphi ^+_0$, (a-1) $\varphi ^+_1$, (a-2) $\varphi ^+_2$, (a-5) $\varphi ^+_5$, (b-1) $\varphi ^+_6$, (b-2) $\varphi ^+_7$, (c-0) $\varphi ^+_{11}$, (c-1) $\varphi ^+_{12}$ and (c-2) $\varphi ^+_{13}$.

Let us consider the stretched problem

\[ \begin{cases} \tilde{u}_{yy}-\tilde{u}+\tilde{u}^3=0, & -\frac{x_{\ell,\varepsilon}}{\varepsilon}< y<\frac{1-x_{\ell,\varepsilon}}{\varepsilon},\\ \tilde{u}_y\left(-\frac{x_{\ell,\varepsilon}}{\varepsilon}\right)=\tilde{u}_y(\frac{1-x_{\ell,\varepsilon}}{\varepsilon})=0 \end{cases} \]

and the limit problem

(1.6)$$ \begin{cases} \bar{u}_{yy}-\bar{u}+\bar{u}^3=0, & -\infty< y<\infty,\\ \bar{u}\in L^2(\mathbb{R}). \end{cases}$$

Here $\tilde {u}(y):=u(x)$, $y:=(x-x_{\ell,\varepsilon })/\varepsilon$ and $x_{\ell,\varepsilon }\in (0,1)$ is a position of an $\ell$-th interior peak. Then $\bar {u}(y)=\frac {\sqrt {2}}{\cosh (y)}$ is a solution of (1.6) which has a one-peak at $y=0$. The associated eigenvalue problem is as follows:

(1.7)$$ \begin{cases} \bar{\varphi}_{yy}+({-}1+3\bar{u}^2)\bar{\varphi}={-}\lambda\bar{\varphi}, & -\infty< y<\infty,\\ \bar{\varphi}\in H^1(\mathbb{R}). \end{cases} $$

It is known that the spectral set consists of two eigenvalues $-3$ and $0$ and the continuous spectrum $[1,\infty )$. The near $-3$ eigenvalues of (LP$_{\pm }$) come from the $-3$ eigenvalue of (1.7). Figures 2–4 indicate that an eigenfunction in category (a) can be approximated by a linear combination of translations of a compressed first eigenfunction of (1.7). The near $0$ eigenvalues of (LP$_{\pm }$) come from a translation invariance of the one-peak solution $\bar {u}(y)$, which is indicated by a shape of an eigenfunction of category (b), e.g., (b-0), (b-1) and (b-2) of figure 2. Specifically, an eigenfunction is close to $cu'(x)$, $c\in \mathbb {R}$, near each peak of $u(x)$. The other eigenvalues converge to $1$ which is an end point of $[1,\infty )$. The spectrum $1$ of (1.7) is not an eigenvalue, since a corresponding eigenfunction is not in $H^1(\mathbb {R})$. Indeed, a graph of an eigenfunction of category (c) in figures 2–4 is close to $\cos (m\pi x)$ for most $x\in (0,1)$, and hence its stretched function is not in $H^1(\mathbb {R})$ or $0$.

Remark 1.5 If $n$ is odd, then the eigenvalue sets of (${\rm LP}_+$) and (${\rm LP}_-$) are the same, i.e., $\lambda ^+_{j,\varepsilon }=\lambda ^-_{j,\varepsilon }$ for all $j\ge 0$, since $u^-_{n,\varepsilon }(x)=u^+_{n,\varepsilon }(1-x)$. However, if $n$ is even, then $\lambda ^+_{j,\varepsilon }=\lambda ^-_{j,\varepsilon }$ for $j\neq n/2$, but $\lambda ^+_{n/2,\varepsilon }\neq \lambda ^-_{n/2,\varepsilon }$. In the even case the number of all the peaks of $u_{n,\varepsilon }^+$ is different from that of $u_{n,\varepsilon }^-$.

Corollaries 1.6 and 1.7 are asymptotic formulas with respect to ${\rm e}^{-2/n\varepsilon }$.

Corollary 1.6 Let $n\ge 1$ be fixed. The following hold for (LP_±):

1. (i) $\displaystyle \lambda ^{\pm }_{0,\varepsilon }=-3+48\,{\rm e}^{-\frac {2}{n\varepsilon }}+o({\rm e}^{-\frac {2}{n\varepsilon }})$ as $\varepsilon \to 0$.

2. (ii) If $n$ is even, then as $\varepsilon \to 0$,

   \[ \lambda^+_{n/2,\varepsilon}={-}48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}),\quad \lambda^-_{n/2,\varepsilon}={-}3+48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}). \]

3. (iii) $\displaystyle \lambda ^{\pm }_{n,\varepsilon }=1-48\,{\rm e}^{-\frac {2}{n\varepsilon }}+o({\rm e}^{-\frac {2}{n\varepsilon }})$ as $\varepsilon \to 0$.

Corollary 1.7 Let $n\ge 1$ and $j\ge 0$ be fixed. The following hold for (${\rm LP}_{\pm }$):

1. (i) For $0< j< n/2$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}={-}3+\left\{1+4\left(\sin^2\frac{j\pi}{n}\right){\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}})\right\}(\lambda^{{\pm}}_{0,\varepsilon}+3) \ \textrm{as}\ \varepsilon\to 0. \]

2. (ii) For $n/2< j< n$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}={-}48\left(\sin^2\frac{j\pi}{n}\right){\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}})\ \textrm{as}\ \ \varepsilon\to 0. \]

3. (iii) For $j>n$,

   \[ \lambda_{j,\varepsilon}^{{\pm}}=1+(j-n)^2\pi^2\varepsilon^2+o(\varepsilon^2) \ \textrm{as}\ \varepsilon\to 0. \]

It follows from corollary 1.7 (i) that $\lambda ^{\pm }_{j,\varepsilon }$, $0< j< n/2$, can be also written as follows:

\[ \lambda^{{\pm}}_{j,\varepsilon}={-}3+48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}) \ \textrm{as}\ \ \varepsilon\to 0. \]

However a $j$-dependence does not appear in the first two terms.

Let us recall known results. It is difficult to obtain exact expressions of eigenvalues of an elliptic differential operator even if a domain is a finite interval. Only few examples are known for elliptic differential operators with variable coefficients. A model case is the linearization problem for the Allen–Cahn equation (1.1) with $f(u)=u-u^3$. The problem (1.1) with $f(u)=u-u^3$ has $n$-mode solutions $u^{\pm }_{n,\varepsilon }$ for small $\varepsilon >0$ which has $n$ transition layers in the interval. In [4] the so-called very slow dynamics of a transition layer solution of an associated parabolic problem were studied. Then the first $n$ eigenvalues of (LP_±), $\{\lambda ^{\pm }_{j,\varepsilon }\}_{j=0}^{n-1}$, play a crucial role. It was shown in [4, Corollary 4.2] that for $0\le j< n$,

(1.8)$$ \lambda^{{\pm}}_{j,\varepsilon}=O({\rm e}^{{-}d/\varepsilon})\ {\rm with\ some}\ d>0. $$

These near $0$ eigenvalues come from translation invariance of each transition layer. Later exact expressions of three special eigenvalues $\lambda ^{\pm }_{j,\varepsilon }$, $j=0,n,2n$, were obtained in [9]. Using these three exact eigenvalues, one can obtain

$$\begin{align*} & \lambda^{{\pm}}_{0,\varepsilon}={-}96\,{\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}}+o({\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0,\\ & \lambda^{{\pm}}_{n,\varepsilon}= \frac{3}{2}-12\,{\rm e}^{-\frac{1}{\sqrt{2}n\varepsilon}}+o({\rm e}^{-\frac{1}{\sqrt{2}n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0,\\ & \lambda^{{\pm}}_{2n,\varepsilon}= 2+96\,{\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}}+o({\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0. \end{align*}$$

Note that these three eigenvalues correspond to $\lambda ^{\pm }_{j,\varepsilon }$, $j=0,n/2,n$, in our case $f(u)=-u+u^3$. Then an exact representation formula of eigenvalues for general $f$, which is lemma 2.2 in our case, was obtained in [12]. The authors of [12] applied it to the case $f(u)=\sin u$ and obtained an asymptotic formula for every eigenvalue. After that in [14] asymptotic formulas of all the eigenvalues for the Allen–Cahn case $f(u)=u-u^3$ were established. Specifically,

$$\begin{align*} & \textrm{for}\ 0< j< n,\ \lambda^{{\pm}}_{j,\varepsilon}={-}96\left(\cos^2\frac{j\pi}{2n}\right){\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}}+o({\rm e}^{-\frac{\sqrt{2}}{n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0,\\ & \textrm{for}\ n< j<2n,\ \lambda^{{\pm}}_{j,\varepsilon}=\frac{3}{2}-12\left(\cos\frac{(j-n)\pi}{n}\right){\rm e}^{-\frac{1}{\sqrt{2}n\varepsilon}}+o({\rm e}^{-\frac{1}{\sqrt{2}n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0,\\ & \textrm{for}\ j>2n, \lambda^{{\pm}}_{j,\varepsilon}=2+(j-2n)^2\pi^2\varepsilon^2+o(\varepsilon^2)\ \textrm{as}\ \varepsilon\to 0. \end{align*}$$

These formulas show that (1.8) holds for $0\le j< n$.

Detailed studies of all the eigenfunctions were made for the case $f(u)=\sin u$ in [13] and for the case $f(u)=u-u^3$ in [15]. The authors of [8] obtained exact eigenvalues and eigenfunctions of a one-dimensional Gel'fand problem $f(u)=e^u$, using the method developed in [12]. In this paper, we apply the method developed in [14] to the scalar field equation (1.1) with $f(u)=-u+u^3$ to establish asymptotic formulas. This method seems to work for other nonlinearities and other boundary conditions. They may be future works.

The paper consists of five sections. In §2 we prove proposition 1.1. We obtain asymptotic formulas for three special eigenvalues $\lambda ^{\pm }_{0,\varepsilon }$, $\lambda ^{\pm }_{n/2,\varepsilon }$ and $\lambda ^{\pm }_{n,\varepsilon }$ (corollary 1.6). We recall an exact representation formulas for the case $f(u)=-u+u^3$ in lemma 2.2. In particular, the other eigenvalues are given by a unique solution of

\[ \mathcal{A}_0(k,\lambda_j)=\frac{j\pi}{n}\ \textrm{for}\ j\neq 0,n/2,n, \]

where $\mathcal {A}_0$ is given in lemma 2.2. Let $\mu :=(2-k^2)\lambda$. We use $\mu$ instead of $\lambda$, since various formulas becomes simple. We define $\mathcal {A}(k,\mu ):=A_0(k,\lambda )$. We see in lemma 2.3 that the characteristic function $\mathcal {A}(k,\mu )$ can be written in terms of the complete elliptic integral of the third kind. In § 3 we study the shape of the graph of $\mathcal {A}(k,\mu )$, since the graph of $\mathcal {A}(k,\mu )$ is directly related to the $j+1$-th eigenvalue by the equation

(1.9)$$ \mathcal{A}(k,\mu)=\frac{j\pi}{n}\ \textrm{for}\ j\neq 0,n/2,n. $$

Specifically, we show that $\mathcal {A}(k,\mu )$ is increasing in $\mu$ and that for each $j$, (1.9) has a unique solution. See figure 5 in § 2 for a graph of $\mathcal {A}(k,\mu )$. In § 4 we prove theorem 1.3 and corollary 1.7, using results proved in § 3. In Appendix A we recall definitions and basic properties of Jacobi's elliptic functions ${\rm sn}(x,k)$, ${\rm cn}(x,k)$ and ${\rm dn}(x,k)$. Then we recall basic properties of the complete elliptic integrals $K(k)$, $E(k)$ and $\Pi (\nu,k)$. We also recall various asymptotic formulas for $K(k)$ and $\Pi (\nu,k)$ which are used in proofs of lemmas in § 3 and theorem 1.3.

Figure 5. A graph of $\mathcal {A}(k,\mu )$ with $k=3/4$. $\mathcal {A}(k,\mu )$ is defined on $\Sigma$ and monotone increasing in $\mu$.

2. Preliminaries

2.1 Exact solutions of the nonlinear problem

Proof of proposition 1.1 Let $y={\rm dn}(x,k)$. Then, $y(x)$ satisfies $y''-(2-k^2)y+2y^3=0$. We can find a solution, assuming $u(x)=c_0y(c_1x)$ and fixing $c_0$ and $c_1$. However, we take a different approach here.

We apply a change of variables $\tau =\sqrt {1-k^2s^2}$ to (A.1). We have

(2.1)$$ x=\int_0^{{\rm sn}(x,k)}\frac{{\rm d}s}{\sqrt{(1-s^2)(1-k^2s^2)}} =\int_{{\rm dn}(x,k)}^1\frac{{\rm d}\tau}{\sqrt{(1-\tau^2)(k^2-1+\tau^2)}}, $$

where (A.2) is used. Let ${\rm dn}^{-1}(y,k)$ denote the inverse function of $y={\rm dn}(x,k)$. Substituting $x={\rm dn}^{-1}(y,k)$ into (2.1), we obtain

(2.2)$$ {\rm dn}^{{-}1}(y,k)=\int_y^1\frac{{\rm d}\tau}{\sqrt{(1-\tau^2)(k^2-1+\tau^2)}}. $$

We consider only the case $u^-_{n,\varepsilon }(x)$, since the case $u^+_{n,\varepsilon }(x)$ is similar. For simplicity, $u$ stands for $u^-_{n,\varepsilon }$. By (1.1) we have

\[ \frac{\varepsilon^2u'^2}{2}-\frac{u^2}{2}+\frac{u^4}{4}={-}\frac{\alpha^2}{2}+\frac{\alpha^4}{4}, \]

where

\[ \alpha:=u(0)=\max_{0\le x\le 1}u(x) ({>}1). \]

Since $2\varepsilon ^2u'^2=(\alpha ^2-u^2)(\alpha ^2-2+u^2)$, we have

\[ \frac{x}{\sqrt{2}\varepsilon}=\int_0^x\frac{-u'\,{\rm d}x}{\sqrt{(\alpha^2-u^2)(\alpha^2-2+u^2)}}, \]

where we use the fact that $u'(x)<0$ for $0< x<1/n$. Let $k$ be defined by the relation

\[ \alpha=\sqrt{\frac{2}{2-k^2}}. \]

Using a change of variables $\tau :=u/\alpha$, we have

$$\begin{align*} \frac{x}{\sqrt{2}\varepsilon}& =\int_{u(x)}^{\alpha}\frac{{\rm d}u}{\sqrt{(\alpha^2-u^2)(\alpha^2-2+u^2)}} =\frac{1}{\alpha}\int_{\frac{u(x)}{\alpha}}^1\frac{{\rm d}\tau}{\sqrt{(1-\tau^2)\left(1-\frac{2}{\alpha^2}+\tau^2\right)}}\\ & = \sqrt{\frac{2-k^2}{2}}\int_{\sqrt{\frac{2-k^2}{2}}u(x)}^1\frac{{\rm d}\tau}{\sqrt{(1-\tau^2)(k^2-1+\tau^2)}} \\ & \quad =\sqrt{\frac{2-k^2}{2}}{\rm dn}^{{-}1}\left(\sqrt{\frac{2-k^2}{2}}u(x),k\right), \end{align*}$$

where (2.2) is used in the last equality. Thus,

(2.3)$$ u(x)=\sqrt{\frac{2}{2-k^2}}{\rm dn}\left(\frac{x}{\varepsilon\sqrt{2-k^2}},k\right). $$

Since $K(k)$ is a half-period of ${\rm dn}(x,k)$ and $1/n$ is a half-period of $u(x)$, we have

\[ \frac{1}{n\varepsilon\sqrt{2-k^2}}=K(k), \]

and hence (1.3) is obtained. Substituting $\varepsilon =1/\sqrt {2-k^2}nK(k)$ into (2.3), we have (1.4).

2.2 Three special eigenvalues

Proposition 1.2 says that there are three (resp. two) exact eigenvalues of (LP_±) if $n$ is even (resp. odd). These three or two eigenvalues are special and play a key role in this paper. Proposition 1.2 can be proved by direct calculation. We define $\lambda _{\pm }(k)$ by

(2.4)$$ \lambda_+:={-}1+\frac{2\sqrt{1-k^2+k^4}}{2-k^2},\quad \lambda_-:={-}1-\frac{2\sqrt{1-k^2+k^4}}{2-k^2}. $$

Then it is obvious that $\lambda ^{\pm }_{0,\varepsilon }=\lambda _-$ and $\lambda ^{\pm }_{n,\varepsilon }=\lambda _+$.

Lemma 2.1 Let $k_{n,\varepsilon }$ be the unique solution of (1.3). Then,

\[ 1-k^2_{n,\varepsilon}=16\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}) \ \textrm{as}\ \ \varepsilon\to 0. \]

Proof. By lemma A.3 we have

\[ (1-k^2)K(k)+\frac{1}{2}(1-k^2)\log(1-k^2)-2(1-k^2)\log 2=o(1) \ \textrm{as}\ \ k\to 1, \]

and hence $(1-k^2)K(k)\to 0$ as $k\to 1$. By (1.3) we have

\[ -\frac{2}{\sqrt{2-k^2}n\varepsilon}={-}\frac{2}{n\varepsilon}+\frac{2(1-k^2)K(k)}{\sqrt{2-k^2}+1}={-}\frac{2}{n\varepsilon}+o(1)\ \textrm{as}\ \ \varepsilon\to 0. \]

Using $K(k)=1/\sqrt {2-k^2}n\varepsilon$, we have

$$\begin{align*} 1-k^2=\exp\left({-}2K(k)+4\log 2+o(1)\right) & =16\exp\left(\frac{-2}{\sqrt{2-k^2}n\varepsilon}+o(1)\right)\\ & = 16\exp\left(\frac{-2}{n\varepsilon}+o(1)\right) \\ & =16\exp\left({-\frac{2}{n\varepsilon}}\right)\left(1+o(1)\right) \ \textrm{as}\ \varepsilon\to 0. \end{align*}$$

We obtain the desired result.

Using lemma 2.1, we obtain asymptotic expansions of three special eigenvalues.

Proof of corollary 1.6 First we prove (i) and (iii). By lemma 2.1 we have

\[ \frac{2\sqrt{1-k^2+k^4}}{2-k^2}=2-3(1-k^2)+o(1-k^2) =2-48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}})\ \textrm{as}\ \varepsilon\to 0. \]

Then

$$\begin{align*} & \lambda^{{\pm}}_{0,\varepsilon}=\lambda_-{=}-1-\frac{2\sqrt{1-k^2+k^4}}{2-k^2}={-}3+48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o\left({\rm e}^{-\frac{2}{n\varepsilon}}\right)\ \textrm{as}\ \varepsilon\to 0,\\ & \lambda^{{\pm}}_{n,\varepsilon}=\lambda_+{=}-1+\frac{2\sqrt{1-k^2+k^4}}{2-k^2}=1-48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o\left({\rm e}^{-\frac{2}{n\varepsilon}}\right)\ \textrm{as}\ \varepsilon\to 0. \end{align*}$$

Next we prove (ii). If $n$ is even, then

$$\begin{align*} & \lambda^+_{n/2,\varepsilon}={-}\frac{3(1-k^2)}{2-k^2}={-}3(1-k^2)+o(1-k^2)={-}48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}) \ \textrm{as}\ \varepsilon\to 0,\\ & \lambda^-_{n/2,\varepsilon}={-}\frac{3}{2-k^2}={-}3+3(1-k^2)\\ & \quad +o(1-k^2)={-}3+48\,{\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}) \ \textrm{as}\ \varepsilon\to 0. \end{align*}$$

2.3 Exact representation of the other eigenvalues

Before we consider the scalar field equation, let us briefly explain a theory for a general nonlinear term developed in Wakasa–Yotsutani [12].

Let $u(x)$ be a solution of (1.1), and let $\alpha$ denote the maximum value of $u(x)$. Substituting $\varphi (x)=\sqrt {\psi (x)}$ into the equation in (1.2), we see that $\psi (x)$ satisfies

\[ \varepsilon^2(2\psi\psi''-\psi'^2)+4(f'(u)+\lambda)\psi^2=0. \]

Let us consider the following function:

(2.5)$$ \Psi(x):=\varepsilon^2(2\psi\psi''-\psi'^2)+4(f'(u)+\lambda)\psi^2. $$

We have

\[ \Psi'(x)=2\psi\left\{\varepsilon^2\psi'''+4(f'(u)+\lambda)\psi'+2f''(u)u'\psi\right\}. \]

We look for a solution of the equation

(2.6)$$ \varepsilon^2\psi'''+4(f'(u)+\lambda)\psi'+2f''(u)u'\psi=0 $$

of the form $\psi (x)=h(u(x))$, where $h(\,\cdot \,)$ is an unknown positive function. Substituting $\psi (x)=h(u(x))$ into (2.6), we obtain a key equation

(2.7)$$ 2(F(\alpha)-F(u))h'''(u)-3f(u)h''(u)+(3f'(u)+4\lambda)h'(u)+2f''(u)h(u)=0, $$

where we use the following relations:

\[ \varepsilon^2u'^2=2(F(\alpha)-F(u))\ \textrm{and}\ \varepsilon^2u''={-}f(u). \]

We assume that an exact expression $h(u)$ of a positive solution of (2.7) can be obtained. Then, $h(u)$ is also a solution of (2.6). Thus, $\Psi '(x)\equiv 0$ and $\Psi$ is constant. Substituting $h(u)$ into (2.5), we see that there exists $\rho \in \mathbb {R}$ such that

(2.8)$$ \Psi(x)=(F(\alpha)-F(u))(4hh''-2h'^2)-2f(u)hh'+4(f'(u)+\lambda)h^2=:4\rho. $$

Now, we construct an eigenfunction of the form

(2.9)$$ \varphi(x)=\sqrt{h(u(x))}W(\theta(x)), $$

using a solution $h(u)$ of (2.7). Here, $W(\,\cdot \,)$ and $\theta (x)$ are defined later. Using (2.8), we have

(2.10)$$\begin{align} & \varepsilon^2\frac{{\rm d}^2\sqrt{h(u(x))}}{{\rm d}x^2}+f'(u)\sqrt{h(u(x))}+\lambda\sqrt{h(u(x))}\nonumber\\ & = \frac{1}{4h^{3/2}}\left\{ (F(\alpha)-F(u))(4hh''-2h'^2)-2f(u)hh'+4(f'(u)+\lambda)h^2\right\} =\frac{\rho}{h^{3/2}}. \end{align}$$

Substituting (2.9) into (1.2), we have

(2.11)$$\begin{align} 0 & =\varepsilon^2\left\{\frac{{\rm d}^2\sqrt{h}}{{\rm d}x^2}W+2\frac{{\rm d}\sqrt{h}}{{\rm d}x}W'\theta'+\sqrt{h}\left(W''\theta'^2+W'\theta''\right)\right\} +f'(u)\sqrt{h}W+\lambda\sqrt{h}W\nonumber\\ & =\varepsilon^2\sqrt{h}\theta'^2\left( W''+\frac{2\frac{{\rm d}\sqrt{h}}{{\rm d}x}\theta'+\sqrt{h}\theta''}{\sqrt{h}\theta'^2}W'\right) +\left(\varepsilon^2\frac{{\rm d}^2\sqrt{h}}{{\rm d}x^2}+f'(u)\sqrt{h}+\lambda\sqrt{h}\right)W\nonumber\\ & =\varepsilon^2\sqrt{h}\theta'^2\left( W''+\frac{1}{h\theta'^2}\frac{{\rm d}}{{\rm d}x}\left(h\theta'\right)W'+\frac{\rho}{\varepsilon^2h^2\theta'^2}W\right), \end{align}$$

where (2.10) is used in the last equality. We assume that $\rho >0$. We define $\theta (x)$ by a solution of

(2.12)$$ h(u(x))\theta'(x)=\frac{\sqrt{\rho}}{\varepsilon}. $$

Specifically, the following function is a solution of (2.12):

\[ \theta(x)=\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho}{\rm d}\xi}{h(u(\xi))}+\theta_0. \]

Moreover, it follows from (2.11) that $W''(\theta )+W(\theta )=0$, and hence $W(\theta )=C \cos (\theta +\theta _1)$. By (2.9) we obtain

(2.13)$$ \varphi(x)=C\sqrt{h(u(x))}\cos\left(\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho_0}\,{\rm d}\xi}{h(u(\xi))}+x_0\right), $$

where $x_0:=\theta _0+\theta _1$. In other words, a general solution of $\varepsilon ^2\varphi ''+f'(u)\varphi +\lambda \varphi =0$ can be written in terms of a solution $u(x)$ of the nonlinear problem (1.1). Readers can find more details of (2.7) and (2.13) in [8, 12, 14].

In our case the maximum value of $u^{\pm }_{n,\varepsilon }(x)$, which is denoted by $\alpha$, is given by

(2.14)$$ \alpha:=\sqrt{\frac{2}{2-k^2}} $$

and (2.7) becomes

(2.15)$$ 2\left(-\frac{\alpha^2}{2}+\frac{\alpha^4}{4}+\frac{u^2}{2}-\frac{u^4}{4}\right)h''' +3(u-u^3)h''+\left\{ 3({-}1+3u^2)+4\lambda\right\}h'+12uh=0. $$

There seems to be no solution formula of (2.7) for general nonlinearities. In the case (2.15) every coefficient is a polynomial of $u$. We substitute $h(u)=u^p$ into (2.15). Then it does not vanish, and the top term becomes

(2.16)$$ -\frac{1}{2}(p-4)(p+1)(p+6)u^{p+1}. $$

This suggests that a bi-quadratic polynomial is a candidate of a solution, since (2.16) vanishes for $p=4$. Actually, $h$ defined by (2.19) is a solution of (2.15). We evaluate (2.8) at a maximum point of $u$. Then,

\[ 2\rho_0={-}f(\alpha)h(\alpha)h'(\alpha)+2(f'(\alpha)+\lambda)h(\alpha)^2, \]

which yields (2.18). Using $h(u)$ and $\rho _0$, we can construct a general solution of $\varepsilon ^2\varphi ''+f(u)\varphi +\lambda \varphi =0$. We take the Neumann boundary condition into account. Let $x_0=0$ in (2.13). Since $1/n$ is a half-period of $u^{\pm }_{n,\varepsilon }(x)$, the following holds:

\[ \frac{n}{\varepsilon}\int_0^{1/n}\frac{\sqrt{\rho_0}\,{\rm d}\xi}{h(u^{{\pm}}_{n,\varepsilon}(\xi))}=j\pi\ \textrm{for}\ j\in\{1,2,\ldots\}. \]

Eigenvalues can be found by solving the above equation.

Finding a solution of (2.7) is key in this theory. A solution for (2.7) was found for

$$\begin{align*} f(u)& =\sin u\ \textrm{in}\ [12],\quad f(u)=u-u^3\ \textrm{in}\ [14],\quad f(u)={\rm e}^{{\pm} u}\ \textrm{in}\ [8]\ \textrm{and}\ f(u)\\ & =\sinh u\ \textrm{in}\ [1]. \end{align*}$$

We summarize these results in the following lemma:

Lemma 2.2 Let $j\neq 0,n$ and $j\neq n/2$ if $n$ is even. Let $k_{\varepsilon }$ be the unique solution of (1.3). We define $\mathcal {A}_0$ by

(2.17)$$\mathcal{A}_0(k,\lambda):=\frac{1}{\varepsilon}\int_0^{1/n}\frac{\sqrt{\rho_0(\lambda,k)}{\rm d}\xi}{|h(u^{{\pm}}_{n,\varepsilon}(\xi),\lambda,k)|} $$

for $\lambda \in \{\lambda \in \mathbb {R};\ \rho _0(\lambda,k)>0\}$, where

(2.18)$$ \rho_0(\lambda,k):=\frac{16}{81}\lambda\left(\lambda+\frac{3}{2-k^2}\right)\left(\lambda+\frac{3-3k^2}{2-k^2}\right)\left(\lambda-\lambda_+\right)\left(\lambda-\lambda_-\right), $$

(2.19)$$ h(u,\lambda,k)=u^4-2\left(\frac{\lambda}{3}+1\right)u^2+\frac{4}{9}\lambda^2+\frac{4}{3}\lambda +\frac{4(1-k^2)}{(2-k^2)^2}. $$

1. (i) Let $\alpha$ be defined by (2.14). If $\rho _0(\lambda,k)>0$, then either $h(u,\lambda,k)>0$ or $h(u,\lambda,k)<0$ for $u\in \left [\sqrt {2-\alpha ^2},\alpha \right ]$.

   (ii) If the equation for $\lambda$

   (2.20)$$ \mathcal{A}_0(k_{\varepsilon},\lambda)=\frac{j\pi}{n} $$
   has a solution $\lambda \in \{\lambda \in \mathbb {R}|\ \rho _0(\lambda,k)>0\}$, then $\lambda$ is a $j+1$-th eigenvalue of $({\rm LP}_{\pm })$ and
   \[ \varphi^{{\pm}}_{j,\varepsilon}(x)=\sqrt{|h(u^{{\pm}}_{n,\varepsilon}(x),\lambda,k_{\varepsilon})|}\cos\left(\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho_0(\lambda,k_{\varepsilon})}\,{\rm d}\xi}{|h(u^{{\pm}}_{n,\varepsilon}(\xi),\lambda,k_{\varepsilon})|}\right) \]
   is an associated eigenfunction.

Proof. Here, we directly prove the lemma.

1. (i) By direct calculation we can check that $h(u)$ satisfies

   (2.21)$$ (F(\alpha)-F(u))(2hh''-h'^2)-f(u)hh'+2(f'(u)+\lambda)h^2=2\rho_0. $$
   Suppose the contrary, i.e., there exists $\bar {u}\in [\sqrt {2-\alpha ^2},\alpha ]$ such that $h(\bar {u})=0$. Since $F(\alpha )-F(\bar {u})\ge 0$, it follows from (2.21) that
   \[ 0\ge -(F(\alpha)-F(\bar{u}))h'^2=2\rho_0>0, \]
   which is a contradiction. This contradiction leads to the assertion (i).

2. (ii) Without loss of generality, we can assume that $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. Let

   \[ \varphi(x):=\sqrt{h(u^{{\pm}}_{n,\varepsilon}(x))}\cos\left(\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho_0}\,{\rm d}\xi}{h(u^{{\pm}}_{n,\varepsilon}(\xi))}\right). \]
   For simplicity, we write $u^{\pm }_{n,\varepsilon }$ and $h(u^{\pm }_{n,\varepsilon })$ as $u$ and $h$, respectively. Then,
   $$\begin{align*} & \varepsilon^2\varphi_{xx}+f'(u)\varphi+\lambda\varphi\\ & =\frac{1}{2h^{3/2}}\left\{\frac{\varepsilon^2u'^2}{2}(2hh''-h'^2)+\varepsilon^2u''hh'-2\rho_0+2f'(u)h^2+2\lambda h^2 \right\}\\& \quad \cos\left(\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho_0}\,{\rm d}\xi}{h(u(\xi))}\right)\\ & =\frac{1}{2h^{3/2}}\left\{ (F(\alpha)-F(u))(2hh''-h'^2)-f(u)hh'-2\rho_0+2f'(u)h^2+2\lambda h^2 \right\}\\& \quad\cos\left(\frac{1}{\varepsilon}\int_0^x\frac{\sqrt{\rho_0}\,{\rm d}\xi}{h(u(\xi))}\right)=0. \end{align*}$$
   By (2.20) we see that $\varphi (x)$ has exactly $j$ zero(s) in $0\le x\le 1$. It follows from Sturm–Liouville theory that $\varphi (x)$ is a $j+1$-th eigenfunction and $\lambda$ is the associated eigenvalue.

Since $\sqrt {2-\alpha ^2}\le u^{\pm }_{n,\varepsilon }(x)\le \alpha$ for $0\le\! x\le\! 1$ (proposition 1.1), either $h(u^{\pm }_{n,\varepsilon }(x))\!>\!0$ or $h(u^{\pm }_{n,\varepsilon }(x))<0$ holds in the case $\rho _0(\lambda,k)>0$. Hence, $\mathcal {A}_0(k,\lambda )$ is well defined. It follows from lemma 2.2 that the $j+1$-th eigenvalue is determined by (2.20) and that the other eigenvalues are on $\{\rho _0>0\}$. Hence, it is important to study the function $\mathcal {A}_0(k,\lambda )$.

Let

\[ \mu:=(2-k^2)\lambda. \]

We sometimes use $\mu$ instead of $\lambda$, since various formulas become simple. We define

(2.22)$$\begin{align} & \mu_+:=(2-k^2)\lambda_+{=}-2+k^2+2\sqrt{1-k^2+k^4},\quad\mu_-:=(2-k^2)\lambda_-{=}-2\nonumber\\ & \quad +k^2-2\sqrt{1-k^2+k^4}. \end{align}$$

We will show that the other eigenvalues are on $\{\lambda \in \mathbb {R}|\ \rho _0(\lambda,k)>0\}$. The set $\{(k,\lambda )|\ k\in (0,1), \rho _0(\lambda,k)>0\}$ corresponds to

(2.23)$$ \Sigma:=\left\{(k,\mu)|\ k\in (0,1),\ \rho_0\left(\frac{\mu}{2-k^2},k\right)>0\right\}. $$

The set $\Sigma$ is split into three components, i.e., $\Sigma =\Sigma _0\cup \Sigma _1\cup \Sigma _2$, where

(2.24)$$ \begin{aligned} \Sigma_0 & :=\{(k,\mu)|\ k\in (0,1),\ \mu\in (\mu_-,-3)\},\\ \Sigma_1 & :=\{(k,\mu)|\ k\in (0,1),\ \mu\in ({-}3+3k^2,0)\},\\ \Sigma_2 & :=\{(k,\mu)|\ k\in (0,1),\ \mu\in (\mu_+,\infty)\}. \end{aligned}$$

Let

(2.25)$$ \mathcal{A}(k,\mu):=\mathcal{A}_0(k,\lambda), $$

where $\mathcal {A}_0$ is defined by (2.17). Then the characteristic equation (2.20) becomes

\[ \mathcal{A}(k,\mu)=\frac{j\pi}{n}. \]

Figure 5 shows the graph of $\mathcal {A}(k,\,\cdot \,)$. In § 3 we rigorously study the graph of $\mathcal {A}(k,\,\cdot \,)$. Specifically, end points of $\mathcal {A}(k,\,\cdot \,)$ are obtained in lemma 3.1 and the monotone increase of $\mathcal {A}(k,\mu )$ in $\mu$ is proved in lemma 3.4.

We would like to obtain a more simple expression of $\mathcal {A}(k,\mu )$. The following lemma indicates that $\mathcal {A}(k,\mu )$ can be decomposed into two complete elliptic integrals of the third kind.

Lemma 2.3 Suppose that the assumptions of lemma 2.2 hold. Then the following hold$:$

1. (i) If $(k,\mu )\in \Sigma _0\cup \Sigma _1$, then

   (2.26)$$ \mathcal{A}(k,\mu)=|{\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)-{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k)|, $$
   where ${\rm sgn}(\,\cdot \,)$ denotes the sign function,
   (2.27)$$\begin{align} \mathcal{M}(\nu,k)& :=\sqrt{\frac{(1+\nu)(k^2+\nu)}{\nu}}\Pi(\nu,k),\quad \Pi(\nu,k)\nonumber\\ & :=\int_0^1\frac{{\rm d}s}{(1+\nu s^2)\sqrt{(1-s^2)(1-k^2s^2)}}, \end{align}$$
   (2.28)$$ \nu_{{\pm}}(k,\mu):=\frac{3k^2}{2}\frac{\mu-3k^2\pm\sqrt{-3\mu^2+6(k^2-2)\mu+9k^4}}{\mu(\mu+3-3k^2)}. $$

2. (ii) If $(k,\mu )\in \Sigma _2$, then

   \[ \mathcal{A}(k,\mu)=\frac{1}{9k^4}\sqrt{\mu(\mu+3)(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}\tilde{\Pi}(a,b,k), \]
   where
   (2.29)$$\begin{align} \tilde{\Pi}(a,b,k):=\int_0^1\frac{{\rm d}s}{\sqrt{(1-s^2)(1-k^2s^2)}\{a+(b-s^2)^2\}}, \end{align}$$
   (2.30)$$\begin{align} a(k,\mu):=\frac{1}{12k^2}(\mu-\mu_+)(\mu-\mu_-),\quad b(k,\mu):=\frac{3k^2-\mu}{6k^2}. \end{align}$$

Proof. (i) We consider the case $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. In both cases ${u^+_{n,\varepsilon }}$ and ${u^-_{n,\varepsilon }}$ by lemma 2.2 we see that

\[ \mathcal{A}_0(k,\lambda)=\frac{\sqrt{\rho_0}}{\varepsilon}\int_0^{1/n}\frac{{\rm d}\xi}{h({u^-_{n,\varepsilon}}(\xi))}, \]

because

\[ \int_0^{1/n}\frac{{\rm d}\xi}{h({u^+_{n,\varepsilon}}(\xi))}=\int_0^{1/n}\frac{{\rm d}\xi}{h({u^-_{n,\varepsilon}}(\xi))}. \]

Using the change of variables $w={u^-_{n,\varepsilon }}(\xi )$, we have

\[ \mathcal{A}_0(k,\lambda)=\frac{\sqrt{\rho_0}}{\varepsilon}\int_{\sqrt{2-\alpha^2}}^{\alpha}\frac{{\rm d}w}{\left|\frac{{\rm d}{u^-_{n,\varepsilon}}}{{\rm d}x}(\xi)\right|h({u^-_{n,\varepsilon}}(\xi))} =\sqrt{\frac{\rho_0}{2}}\int_{\sqrt{2-\alpha^2}}^{\alpha}\frac{{\rm d}w}{\sqrt{F(\alpha)-F(w)}h(w)}, \]

where we use

\[ \varepsilon^2\left(\frac{{\rm d}{u^-_{n,\varepsilon}}}{{\rm d}x}\right)^2=2(F(\alpha)-F({u^-_{n,\varepsilon}})). \]

In order to bring the integral closer to the expression of the complete elliptic integral of the third kind, we let $s=w/\alpha$. Since

\[ 2(F(\alpha)-F(\alpha s))=\frac{\alpha^4}{2}(1-s^2)(s^2-1+k^2), \]

we have

(2.31)$$\begin{align} \mathcal{A}_0(k,\lambda)=\frac{\sqrt{2\rho_0}}{\alpha}\int_{\sqrt{1-k^2}}^1\frac{{\rm d}s}{\sqrt{(1-s^2)(s^2-1+k^2)}h(\alpha s)}. \end{align}$$

We define $\sigma _{\pm }$ by

(2.32)$$\begin{align} \sigma_{{\pm}}:=\frac{\alpha^2}{\frac{\lambda}{3}+1\pm\sqrt{\frac{3k^4-(2-k^2)^2(\lambda^2+2\lambda)}{3(2-k^2)^2}}}. \end{align}$$

Then

(2.33)$$\begin{align} \frac{1}{h(\alpha s)}=\frac{\sqrt{3}}{2\alpha^2\sqrt{-(\lambda-\lambda_+)(\lambda-\lambda_-)}} \left(\frac{\sigma_-}{1-\sigma_- s^2}-\frac{\sigma_+}{1-\sigma_+ s^2}\right). \end{align}$$

Note that $-(\lambda -\lambda _+)(\lambda -\lambda _-)>0$ for $(k,\mu )\in \Sigma _0\cup \Sigma _1$. Using (2.33), we have

$$\begin{align*} \mathcal{A}_0(k,\lambda)& =\frac{\sqrt{6\rho_0(k,\lambda)}}{2\alpha^3\sqrt{-(\lambda-\lambda_+)(\lambda-\lambda_-)}}\int_{\sqrt{1-k^2}}^1\frac{1}{\sqrt{(1-s^2)(s^2-1+k^2)}}\\ & \quad \left(\frac{\sigma_-}{1-\sigma_- s^2}-\frac{\sigma_+}{1-\sigma_+ s^2} \right){\rm d}s. \end{align*}$$

Using (2.18), which is the definition of $\rho _0$, and

\[ \left(\lambda+\frac{3}{2-k^2}\right)\left(\lambda+\frac{3-3k^2}{2-k^2}\right)=\frac{4\alpha^4}{9\sigma_+\sigma_-}, \]

we have

\[ \mathcal{A}_0(k,\lambda)=\sqrt{\frac{(2-k^2)(-\lambda)}{3\sigma_+\sigma_-}}\left(\sigma_-\Pi_0(\sigma_-,k)-\sigma_+\Pi_0(\sigma_+,k)\right), \]

where

\[ \Pi_0(\sigma,k):=\int_{\sqrt{1-k^2}}^1\frac{{\rm d}s}{(1-\sigma s^2)\sqrt{(1-s^2)(s^2-1+k^2)}}. \]

Using the change of variables $s:=\sqrt {1-k^2\tau ^2}$, we have

\[ \Pi_0(\sigma,k)=\int_0^1\frac{{\rm d}\tau}{(1-\sigma+\sigma k^2\tau^2)\sqrt{(1-\tau^2)(1-k^2\tau^2)}}=\frac{1}{1-\sigma}\Pi\left(\frac{k^2\sigma}{1-\sigma},k\right). \]

Therefore,

\[ \mathcal{A}_0(k,\lambda)=\sqrt{\frac{(2-k^2)(-\lambda)}{3\sigma_+\sigma_-}}\left(\frac{\sigma_-}{1-\sigma_-}\Pi\left(\frac{k^2\sigma_-}{1-\sigma_-},k\right) -\frac{\sigma_+}{1-\sigma_+}\Pi\left(\frac{k^2\sigma_+}{1-\sigma_+},k\right)\right). \]

Since $\lambda =\mu /(2-k^2)$, we can check that

(2.34)$$\begin{align} \nu_+{=}\frac{k^2\sigma_-}{1-\sigma_-}\ \ \textrm{and}\ \ \nu_-{=}\frac{k^2\sigma_+}{1-\sigma_+}. \end{align}$$

Here $\sigma _{\pm }$ are given by (2.32) and $\nu _{\pm }$ are given by (2.28). Using

\[ \sigma_+{=}\frac{\nu_-}{\nu_-{+}k^2}\ \ \textrm{and}\ \sigma_-{=}\frac{\nu_+}{\nu_+{+}k^2}, \]

by lemma 2.4 below we have

(2.35)$$\begin{align} \mathcal{A}(k,\mu) & =\sqrt{\frac{-\mu}{3}}\sqrt{\frac{(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_+\nu_-}} \left(\frac{\nu_+}{k^2}\Pi(\nu_+,k)-\frac{\nu_-}{k^2}\Pi(\nu_-,k)\right)\nonumber\\ & =\frac{1}{k^2}\sqrt{\frac{-\mu}{3}}\left( {\sqrt{\left|\frac{\nu_+(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_-}\right|}}{\rm sgn}(\nu_+)\Pi(\nu_+,k)\right.\nonumber\\ & \quad\left.-{\sqrt{\left|\frac{\nu_-(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_+}\right|}}{\rm sgn}(\nu_-)\Pi(\nu_-,k) \right) \nonumber\\ & ={\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)-{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k). \end{align}$$

Here $\mathcal {M}$ and $\Pi$ are defined by (2.27). When $h({u^-_{n,\varepsilon }}(x))<0$ for $x\in [0,1]$, we have

\[ \mathcal{A}_0(k,\lambda)=\frac{\sqrt{\rho_0}}{\varepsilon}\int_0^{1/n}\frac{{\rm d}\xi}{-h({u^-_{n,\varepsilon}}(\xi))} ={-}\left({\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)\!-\!{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k)\right)\!>\!0. \]

We obtain (2.26).

(ii) We can check that if $(k,\mu )\in \Sigma _2$, then $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. By the same way as in (i) we have (2.31). Using the change of variables $s:=\sqrt {1-k^2\tau ^2}$, we have

\[ \mathcal{A}(k,\mu)=\frac{\sqrt{2\rho_0}}{\alpha^5k^4}\int_0^1\frac{{\rm d}\tau}{\sqrt{(1-\tau^2)(1-k^2\tau^2)}\left\{a+(b-\tau^2)^2\right\}}, \]

where

$$\begin{align*} & a=\frac{1}{\alpha^4k^4}\left(\frac{\lambda^2+2\lambda}{3}-\frac{k^4}{(2-k^2)^2}\right) =\frac{1}{12k^4}(\mu-\mu_+)(\mu-\mu_-),\\ & b=\frac{1}{\alpha^2k^2}\left(\alpha^2-1-\frac{\lambda}{3}\right) =\frac{3k^2-\mu}{6k^2}. \end{align*}$$

In this case the polynomial $a+(b-\tau ^2)^2$ cannot be factored into two real quadratic polynomials. We have

\[ \mathcal{A}(k,\mu)=\frac{1}{9k^4}\sqrt{\mu(\mu+3)(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}\tilde{\Pi}(a,b,k). \]

The proof of (ii) is complete.

The following equality, which is somewhat nontrivial, is left in the proof of lemma 2.3.

Lemma 2.4 Let $\sigma _{\pm }$ and $\nu _{\pm }$ be defined by (2.32) and (2.28), respectively. If $(k,\mu )\in \Sigma _0\cup \Sigma _1$, then

(2.36)$$\begin{align} \frac{1}{k^2}\sqrt{\frac{-\mu}{3}}\sqrt{\left|\frac{\nu_{{\pm}}(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_{{\mp}}}\right|} =\sqrt{\frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}}}. \end{align}$$

Proof. By (2.28) we can check that

(2.37)$$\begin{align} \nu_+{+}\nu_-{=}\frac{3k^2(\mu-3k^2)}{\mu(\mu+3-2k^2)},\quad \nu_+\nu_-{=}\frac{9k^4}{\mu(\mu+3-3k^2)}. \end{align}$$

Using (2.37), we have

(2.38)$$\begin{align} & -\frac{\mu}{3k^4}\frac{(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_{{\pm}}\nu_{{\mp}}} ={-}\frac{\mu}{3k^4}\frac{\nu_+\nu_-{+}k^2(\nu_+{+}\nu_-)+k^4}{\nu_+\nu_-}\nonumber\\ & \quad ={-}\frac{1}{27k^4}\mu(\mu+3)(\mu+3-3k^2). \end{align}$$

On the other hand, we consider the right-hand side of (2.36). Let $\gamma _{\pm }:=\alpha ^2/\sigma _{\pm }$. Then $\gamma _{\pm }$ satisfies

(2.39)$$\begin{align} \gamma_{{\pm}}^2-2\left(\frac{\lambda}{3}+1\right)\gamma_{{\pm}} +\frac{4}{9}\lambda^2+\frac{4}{3}\lambda+\frac{4(1-k^2)}{(2-k^2)^2}=0. \end{align}$$

Since

\[ \frac{1}{\nu_{{\pm}}}=\frac{1}{k^2}\left(\frac{1}{\sigma_{{\mp}}}-1\right)=\frac{1}{k^2}\left(\frac{2-k^2}{2}\gamma_{{\mp}}-1\right), \]

we have

$$\begin{align*} \frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}^3} & =\frac{1}{\nu_{{\pm}}}\left(1+\frac{1}{\nu_{{\pm}}}\right)\left(1+\frac{k^2}{\nu_{{\pm}}}\right)\\ & =\frac{1}{k^2}\left(\frac{2-k^2}{2}\gamma_{{\mp}}-1\right)\left(1+\frac{2-k^2}{2k^2}\gamma_{{\mp}}-\frac{1}{k^2}\right)\\ & \quad \left(1+\frac{2-k^2}{2}\gamma_{{\mp}}-1\right)\\ & =\frac{(2-k^2)^3}{8k^4}\left(\gamma^2_{{\mp}}-2\gamma_{{\mp}}+\frac{4(1-k^2)}{(2-k^2)^2} \right)\gamma_{{\mp}}. \end{align*}$$

By (2.39) we have

(2.40)$$\begin{align} \frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}^3} & =\frac{(2-k^2)^3}{8k^4}\left(\frac{2}{3}\lambda\gamma_{{\mp}}-\frac{4}{9}\lambda^2-\frac{4}{3}\lambda\right)\gamma_{{\mp}} \nonumber\\ & =\frac{(2-k^2)^3}{12k^4}\lambda\left\{\gamma_{{\mp}}^2-2\left(\frac{\lambda}{3}+1\right)\gamma_{{\mp}}\right\}\nonumber\\ & ={-}\frac{(2-k^2)^3}{12k^4}\lambda\left(\frac{4}{9}\lambda^2+\frac{4}{3}\lambda+\frac{4(1-k^2)}{(2-k^2)^2}\right) \nonumber\\ & ={-}\frac{1}{27k^4}\mu(\mu+3)(\mu+3-3k^2). \end{align}$$

It follows from (2.40) and (2.38) that

\[ -\frac{\mu}{3k^4}\frac{(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_{{\pm}}\nu_{{\mp}}} =\frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}^3}. \]

Multiplying both sides by $\nu _{\pm }^2$, we have

\[ -\frac{\mu}{3k^4}\frac{\nu_{{\pm}}(\nu_+{+}k^2)(\nu_-{+}k^2)}{\nu_{{\mp}}} =\frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}}. \]

The assertion holds.

Now the proof of lemma 2.3 is complete.

3. Fundamental properties of $\mathcal {A}$

Let $\mathcal {A}(k,\mu )$ be defined by (2.25), and let $\Sigma$ be defined by (2.23).

Lemma 3.1 Let $k\in (0,1)$ be fixed and $(k,\mu )\in \Sigma$, and let $\mu _{\pm }$ be defined by (2.22). Then the following hold$:$

1. (i) $\mathcal {A}(k,\mu )\to 0$ as $\mu \to \mu _-$.

2. (ii) $\mathcal {A}(k,\mu )\to {\pi }/{2}$ as $\mu \to -3$.

3. (iii) $\mathcal {A}(k,\mu )\to {\pi }/{2}$ as $\mu \to -3+3k^2$.

4. (iv) $\mathcal {A}(k,\mu )\to \pi$ as $\mu \to 0$.

5. (v) $\mathcal {A}(k,\mu )\to \pi$ as $\mu \to \mu _+$.

6. (vi) $\mathcal {A}(k,\mu )\to \infty$ as $\mu \to \infty$.

Proof.

1. (i) In this case we can check that $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. Let $\nu _{\pm }$ be defined by (2.28). We see that there exists $\nu _*$ such that

   \[ \lim_{\mu\to\mu_-}\nu_+{=}\nu_*,\quad \lim_{\mu\to\mu_-}\nu_-{=}\nu_*. \]
   By direct calculation we can check that $-1<\nu _*<0$. Let $\mathcal {M}$ be defined by (2.27). Since
   \[ \lim_{\nu_+\to\nu_*}\mathcal{M}(\nu_+,k)=\mathcal{M}(\nu_*,k),\quad \lim_{\nu_-\to\nu_*}\mathcal{M}(\nu_-,k)=\mathcal{M}(\nu_*,k), \]
   we have
   \[ \lim_{\mu\to\mu_-}\mathcal{A}(k,\mu)=0. \]

2. (ii) In this case we can check that $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. We see that $\nu _+\to -k^2$ ($\mu \to -3$) and $\nu _-\to -1$ ($\mu \to -3$). By lemma A.9 we have

   \[ \lim_{\nu_+\to -k^2}{\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)=0,\quad \lim_{\nu_-\to-1}{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k)={-}\frac{\pi}{2}. \]
   Then,
   \[ \lim_{\mu\to -3}\mathcal{A}(k,\mu)=0-\left(-\frac{\pi}{2}\right)=\frac{\pi}{2}. \]

3. (iii) In this case we can check that $h(u^{\pm }_{n,\varepsilon }(x))<0$ for $x\in [0,1]$. We see that $\nu _+\to -k^2$ $(\mu \to -3+3k^2)$ and $\nu _-\to \infty$ ($\mu \to -3+3k^2$). By lemma A.9 we have

   \[ \lim_{\nu_+\to -k^2}{\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)=0,\quad \lim_{\nu_-\to\infty}{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k)=\frac{\pi}{2}. \]
   Then,
   \[ \lim_{\mu\to -3+3k^2}\mathcal{A}(k,\mu)={-}\left(0-\frac{\pi}{2}\right)=\frac{\pi}{2}. \]

4. (iv) In this case we can check that $h(u^{\pm }_{n,\varepsilon }(x))<0$ for $x\in [0,1]$. We see that $\nu _+\to -1$ ($\mu \to 0$) and $\nu _-\to \infty$ ($\mu \to 0$). By lemma A.9 we have

   \[ \lim_{\nu_+\to -1}{\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k)={-}\frac{\pi}{2},\quad \lim_{\nu_-\to\infty}{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k)=\frac{\pi}{2}. \]
   Then,
   \[ \lim_{\mu\to 0}\mathcal{A}(k,\mu)={-}\left(-\frac{\pi}{2}-\frac{\pi}{2}\right)=\pi. \]

5. (v) In this case we can check that $h(u^{\pm }_{n,\varepsilon }(x))>0$ for $x\in [0,1]$. Let $a,b$ be defined by (2.30). By lemmas 2.3 and A.7 we have

   $$\begin{align*} \mathcal{A}(k,\mu) & =\frac{1}{9k^4}\sqrt{\mu(\mu+3)(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}\tilde{\Pi}(a,b,k)\\ & =\frac{2}{3\sqrt{3}k^2}\sqrt{\mu(\mu+3)(\mu+3-3k^2)}\sqrt{a}\tilde{\Pi}(a,b,k)\\ & \to \frac{2}{3\sqrt{3}k^2}\sqrt{\mu_+(\mu_+{+}3)(\mu_+{+}3-3k^2)}\\ & \qquad \frac{\pi}{2\sqrt{b_0(1-b_0)(1-k^2b_0)}}\ \ \textrm{as}\ \ \mu\to\mu_+. \end{align*}$$
   Here, $b_0:=(3k^2-\mu _+)/6k^2$. Since
   \[ \frac{2}{3\sqrt{3}k^2}\sqrt{\mu_+(\mu_+{+}3)(\mu_+{+}3-3k^2)} =\frac{2}{3\sqrt{3}k^2}\sqrt{(1-k^2+k^4)\mu_+{+}3k^4(2-k^2)}, \]
   \[ \frac{\pi}{2\sqrt{b_0(1-b_0)(1-k^2b_0)}} =\frac{3\sqrt{3}\pi k^2}{2\sqrt{(1-k^2+k^4)\mu_+{+}3k^4(2-k^2)}}, \]
   we have
   \[ \lim_{\mu\to\mu_+}\mathcal{A}(k,\mu)=\pi. \]

6. (vi) Let $\lambda _{\pm }$ be defined by (2.4). We see that $\lambda _+:(0,1)\to \mathbb {R}$ is increasing in $k$,

   \[ \lim_{k\to 0}\lambda_+(k)=0,\ \ \lim_{k\to 1}\lambda_+(k)=1. \]
   Hence, $0<\lambda _+<1$ for $k\in (0,1)$.

   Let $a(k,\mu )$, $b(k,\mu )$ be defined by (2.30). Since $\lambda >\lambda _+>0$, we have see $\mu >\mu _+>0$ and

   \[ (b(k,\mu)-1)^2-b(k,\mu)^2=\frac{\mu}{3k^2}>0. \]
   Hence,
   \[ 0< a(k,\mu)+b(k,\mu)^2< a(k,\mu)+(b(k,\mu)-1)^2. \]
   Using the inequality
   \[ \tilde{\Pi}(a,b,k)\ge\min\left\{\frac{1}{a+b^2},\frac{1}{a+(b-1)^2}\right\}K(k) \]
   for any $a>0$, $b\in (0,1)$, $k\in (0,1)$, we have
   $$\begin{align*} \mathcal{A}(k,\mu) & =\frac{1}{9k^4}\sqrt{\mu(\mu+3)(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}\tilde{\Pi}(a(k,\mu),b(k,\mu),k)\\ & \ge\frac{1}{9k^4}\sqrt{\mu(\mu+3)(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}\\ & \qquad\frac{K(k)}{a(k,\mu)+(b(k,\mu)-1)^2}\\ & =\sqrt{\frac{(\mu+3-3k^2)(\mu-\mu_+)(\mu-\mu_-)}{\mu(\mu+3)}}K(k)\\ & \ge c\sqrt{\mu-\mu_+}K(k), \end{align*}$$
   because there exists $c>0$ such that $\sqrt {\frac {(\mu +3-3k^2)(\mu -\mu _-)}{\mu (\mu +3)}}>c$ for $\mu >\mu _+$. Since
   \[ \mathcal{A}(k,\mu)\ge c\sqrt{\mu-\mu_+}K(k)\to \infty\ \ \textrm{as}\ \ \mu\to\infty, \]
   the assertion (vi) holds.

Let $\nu _{\pm }$ be defined by (2.28). We need two lemmas to prove a monotonicity of $\mathcal {A}(k,\mu )$ in $\mu$.

Lemma 3.2 Let $\nu _{\pm }$ be defined by (2.28). Then the following hold$:$

$$\begin{align*} & \frac{\partial}{\partial\mu}\left[\left\{\frac{k^2}{2}\left(\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)+1\right\} \left(-\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)\right] =\frac{-\{\mu^2+3(2-k^2)\mu+6(1-k^2)\}}{3\sqrt{3}k^2\sqrt{-\mu^2+(2k^2-4)\mu+3k^4}},\\ & \frac{\partial}{\partial\mu}\left( -\frac{1}{\nu_+}+\frac{1}{\nu_-}\right) =\frac{-(\mu-k^2+2)}{\sqrt{3}k^2\sqrt{-\mu^2+(2k^2-4)\mu+3k^4}}. \end{align*}$$

Proof. By (2.28) we can check that

\[ \frac{1}{\nu_+}+\frac{1}{\nu_-}=\frac{\mu-3k^2}{3k^2},\quad - \frac{1}{\nu_+}+\frac{1}{\nu_-}=\frac{\sqrt{-\mu^2+(2k^2-4)\mu+3k^4}}{\sqrt{3}k^2}. \]

Then by direct calculation we obtain the conclusion. We omit details.

We define

(3.1)$$\begin{align} \mathcal{B}(k,\mu):=\left\{\mu^2+3(2-k^2)\mu+6(1-k^2)\right\}K(k)-3(\mu-k^2+2)E(k), \end{align}$$

since a sign of $\mathcal {B}(k,\mu )$ is important in the study of ${\partial \mathcal {A}}/{\partial \mu }$.

Lemma 3.3 The following holds$:$

\[ \mathcal{B}(k,0)<0\ \ \textrm{for}\ \ k\in (0,1). \]

Proof. It is obvious that $\mathcal {B}(0,0)=0$. By lemma A.1 (i) and (ii) we have

$$\begin{align*} \frac{\partial\mathcal{B}}{\partial k}(k,0) & ={-}12kK(k)+6(1-k^2)\frac{E(k)-(1-k^2)K(k)}{k(1-k^2)}\\ & \quad +6kE(k)+3(k^2-2)\frac{E(k)-K(k)}{k}\\ & =9k(E(k)-K(k))<0 \end{align*}$$

for $k\in (0,1)$. Thus, the assertion holds.

Lemma 3.4 Let $\Sigma _i$, $i=0,1,2$, be defined by (2.24). For each $i=0,1,2$,

(3.2)$$\begin{align} \frac{\partial\mathcal{A}}{\partial\mu}(k,\mu)>0\ \ \textrm{for}\ \ (k,\mu)\in\Sigma_i. \end{align}$$

Proof. We consider the case $(k,\mu )\in \Sigma _0$. Let $\nu _{\pm }$ be defined by (2.28). By (2.37) and (2.40) we can see that

(3.3)$$\begin{align} -1<\nu_-{<}\nu_+{<}-k^2,\ \frac{(1+\nu_-)(k^2+\nu_-)}{\nu_-^2}<0,\ \frac{(1+\nu_+)(k^2+\nu_+)}{\nu_+^2}<0. \end{align}$$

Then, $\mathcal {A}(k,\mu )=-\mathcal {M}(\nu _+,k)+\mathcal {M}(\nu _-,k)$. By lemma A.1(iv) we have

$$\begin{align*} \frac{\partial\mathcal{A}}{\partial\mu} & ={-}\frac{\partial\mathcal{M}}{\partial\nu_+}\frac{\partial\nu_+}{\partial\mu} +\frac{\partial\mathcal{M}}{\partial\nu_-}\frac{\partial\nu_-}{\partial\mu} \\ & ={-}\frac{1}{2}\sqrt{\frac{(1+\nu_+)(k^2+\nu_+)}{\nu_+}}\left(\frac{-K(k)}{\nu_+(1+\nu_+)}+\frac{E(k)}{(1+\nu_+)(k^2+\nu_+)}\right)\frac{\partial\nu_+}{\partial\mu}\\ & \quad +\frac{1}{2}\sqrt{\frac{(1+\nu_-)(k^2+\nu_-)}{\nu_-}}\left(\frac{-K(k)}{\nu_-(1+\nu_-)}+\frac{E(k)}{(1+\nu_-)(k^2+\nu_-)}\right)\frac{\partial\nu_-}{\partial\mu}. \end{align*}$$

By (3.3) we have

$$\begin{align*} \frac{\partial\mathcal{A}}{\partial\mu} & ={-}\frac{1}{2}\sqrt{\frac{\nu_+^3}{(1+\nu_+)(k^2+\nu_+)}}\left(\frac{k^2+\nu_+}{\nu_+^3}K(k)-\frac{1}{\nu_+^2}E(k) \right)\frac{\partial\nu_+}{\partial\mu}\\ & \quad +\frac{1}{2}\sqrt{\frac{\nu_-^3}{(1+\nu_-)(k^2+\nu_-)}}\left(\frac{k^2+\nu_-}{\nu_-^3}K(k)-\frac{1}{\nu_-^2}E(k)\right)\frac{\partial\nu_-}{\partial\mu}. \end{align*}$$

By (2.40) we have

\[ \frac{(1+\nu_{{\pm}})(k^2+\nu_{{\pm}})}{\nu_{{\pm}}^3}={-}\frac{1}{27k^4}\mu(\mu+3)(\mu+3-3k^2)=:\mathcal{R}. \]

Then,

$$\begin{align*} \frac{\partial\mathcal{A}}{\partial\mu} & =\frac{1}{2\sqrt{\mathcal{R}}}\left\{\left( -\frac{k^2+\nu_+}{\nu_+^3}\frac{\partial\nu_+}{\partial\mu} +\frac{k^2+\nu_-}{\nu_-^3}\frac{\partial\nu_-}{\partial\mu}\right)K(k)\right.\\ & \left.\quad +\,\left(\frac{1}{\nu_+^2}\frac{\partial\nu_+}{\partial\mu} -\frac{1}{\nu_-^2}\frac{\partial\nu_-}{\partial\mu} \right)E(k)\right\}\\ & =\frac{1}{2\sqrt{\mathcal{R}}}\left[\frac{\partial}{\partial\mu}\left\{\left( -\frac{k^2}{2}\left(\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)-1\right)\left(-\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)\right\}K(k) \right.\\ & \quad \left. +\,\frac{\partial}{\partial\mu}\left(-\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)E(k)\right]. \end{align*}$$

By lemma 3.2 we have

\[ \frac{\partial\mathcal{A}}{\partial\mu} =\frac{\mathcal{B}(k,\mu)}{6\sqrt{3}k^2\sqrt{\mathcal{R}}\sqrt{-\mu^2+(2k^2-4)\mu+3k^4}}, \]

where $\mathcal {B}(k,\mu )$ is defined by (3.1). Since $(k,\mu )\in \Sigma _0$, we see that $-3(\mu -k^2+2)>0$. By lemma A.2 we have

$$\begin{align*} \mathcal{B}(k,\mu) & =\left\{\mu^2+3(2-k^2)\mu+6(1-k^2) \right\}K(k) -3(\mu-k^2+2)E(k)\\ & >\left\{\mu^2+3(2-k^2)\mu+6(1-k^2)-3(\mu-k^2+2)(1-k^2)\right\}K(k)\\ & =\left\{\mu^2+3\mu+3k^2(1-k^2)\right\}K(k)>0. \end{align*}$$

Thus, (3.2) holds for $(k,\mu )\in \Sigma _0$.

We consider the case $(k,\mu )\in \Sigma _1$. By (2.37) and (2.40) we can see that

\[ \nu_+{<}0<\nu_-,\ \frac{(1+\nu_-)(k^2+\nu_-)}{\nu_-^2}>0,\ \frac{(1+\nu_+)(k^2+\nu_+)}{\nu_+^2}<0. \]

Then, $\mathcal {A}(k,\mu )=\mathcal {M}(\nu _-,k)+\mathcal {M}(\nu _+,k)$. By a similar calculation as in (i) we have

$$\begin{align*} \frac{\partial\mathcal{A}}{\partial\mu} & =\frac{1}{2\sqrt{\mathcal{R}}}\left[\frac{\partial}{\partial\mu}\left\{\left( \frac{k^2}{2}\left(\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)+1\right)\left(-\frac{1}{\nu_+}+\frac{1}{\nu_-}\right)\right\}K(k) \right.\\ & \quad \left.+\,\frac{\partial}{\partial\mu} \left(\frac{1}{\nu_+}-\frac{1}{\nu_-}\right)E(k) \right]\\ & =\frac{-\mathcal{B}(k,\mu)}{6\sqrt{3}k^2\sqrt{\mathcal{R}}\sqrt{ -\mu^2+(2k^2-4)\mu+3k^4}}. \end{align*}$$

We show that $\mathcal {B}(k,\mu )<0$ for $(k,\mu )\in \Sigma _1$. Since $\mathcal {B}(k,\mu )$ is convex in $\mu$, it is enough to show that $\mathcal {B}(k,-3+3k^2)<0$ and $\mathcal {B}(k,0)<0$. By lemma A.2 we have

$$\begin{align*} \mathcal{B}(k,-3+3k^2) & ={-}3(1-k^2)K(k)+3(1-2k^2)E(k)\\ & \le -3\frac{1-k^2}{1-\frac{k^2}{2}}E(k)+3(1-2k^2)E(k) ={-}\frac{3k^2(3-2k^2)}{2-k^2}E(k)<0 \end{align*}$$

for $k\in (0,1)$. By lemma 3.3 we see that $\mathcal {B}(k,0)<0$ for $k\in (0,1)$. Thus, (3.2) holds for $(k,\mu )\in \Sigma _1$.

We consider the case $(k,\mu )\in \Sigma _2$. In this case by (2.35) we see that $\mathcal {A}(k,\mu )=\mathcal {M}(\nu _+,k)-\mathcal {M}(\nu _-,k)$, where we consider $\mathcal {M}(\nu _+,k)$ and $\mathcal {M}(\nu _-,k)$ as complex-valued functions. Then by a similar calculation as in (i) we obtain

\[ \frac{\partial\mathcal{A}}{\partial\mu} =\frac{-\mathcal{B}(k,\mu)}{6\sqrt{3}k^2\sqrt{\mathcal{R}}\sqrt{-\mu^2+(2k^2-4)\mu+3k^4}}, \]

Since $\mathcal {R}<0$ and $-\mu ^2+(2k^2-4)\mu +3k^4<0$, we have

\[ \frac{\partial\mathcal{A}}{\partial\mu} =\frac{\mathcal{B}(k,\mu)}{6\sqrt{3}k^2\sqrt{|\mathcal{R}|}\sqrt{|-\mu^2+(2k^2-4)\mu+3k^4|}}. \]

Thus, we show that $\mathcal {B}(k,\mu )>0$ for $\mu >\mu _+$. Since $\mu -k^2+2>0$ for $\mu >\mu _+$, by lemma A.2 we have

$$\begin{align*} \mathcal{B}(k,\mu) & \ge\left\{\mu^2+3(2-k^2)\mu+6(1-k^2)\right\}K(k) -3(\mu-k^2+2)\left(1-\frac{1}{2}k^2\right)K(k)\\ & =\left\{\mu^2+3\left(1-\frac{1}{2}k^2\right)\mu-\frac{3}{2}k^4\right\}K(k)\\ & >\left\{\left( k^2-2+2\sqrt{1-k^2+k^4}\right)^2 \right.\\ & \quad \left. +\,3\left(1-\frac{1}{2}k^2\right)\left( k^2-2+2\sqrt{1-k^2+k^4}\right)-\frac{3}{2}k^4\right\}K(k)\\ & =2\sqrt{1-k^2+k^4}\left\{\sqrt{1-k^2+k^4}-\left(1-\frac{1}{2}k^2\right)\right\}K(k) >0. \end{align*}$$

Thus, (3.2) holds for $(k,\mu )\in \Sigma _2$.

4. Asymptotic formulas

In this section let $\lambda _j(k)$, $j\ge 0$, denote the $j+1$-th eigenvalue of (LP$_{\pm }$) and let $\mu _j(k):=(2-k^2)\lambda _j(k)$.

We know the following:

• • $\mathcal {A}(k,\mu )$ is defined on $\Sigma _0\cup \Sigma _1\cup \Sigma _2$ (lemma 2.2),

• • $\mathcal {A}(k,\mu )$ is increasing in $\mu$ (lemma 3.4),

• • the range of $\mathcal {A}(k,\,\cdot \,)$ is $\mathbb {R}_{+}\setminus \{\pi /2,\pi \}$, where $\mathbb {R}_{+}:=\{x|\ x>0\}$ (lemma 3.1).

The equation

\[ \mathcal{A}(k,\mu_j)=\frac{j\pi}{n},\ j\neq 0,n/2,n \]

has a unique solution $\mu _j$ in $\Sigma _0$ for $0< j< n/2$, in $\Sigma _1$ for $n/2< j< n$ and in $\Sigma _2$ for $j>n$. Then, $\lambda _j=\mu _j/(2-k^2)$ is the $j+1$-th eigenvalue for $j\neq 0,n/2,n$. In this section we obtain an asymptotic expansion of $\lambda _j$ as $k\to 1$ for $j\neq 0,n/2,n$.

Proof of theorem 1.3(i) Let $\lambda _{\pm }$ be defined by (2.4). We consider the case $\lambda _-<\lambda _j<-3/(2-k^2)$, which corresponds to the case $0< j< n/2$. Because of the definition of $\lambda _-$ we have

\[ \frac{2}{(2-k^2)(2-k^2+\sqrt{1-k^2+k^4})} =\frac{\frac{\lambda_-}{3}+1}{1-k^2}\le\frac{\frac{\lambda_j}{3}+1}{1-k^2}\le\frac{\frac{-3}{3(2-k^2)}+1}{1-k^2} =\frac{1}{2-k^2}, \]

and hence

(4.1)$$\begin{align} \lim_{k\to 1}\frac{\frac{\lambda_j}{3}+1}{1-k^2}=1. \end{align}$$

We define $r_j(k)$ by

\[ r_j:=\frac{\lambda_j+3}{\lambda_-{+}3}-1. \]

It follows from corollary 1.6(i) that $\lambda _-+3>0$ for $k$ close to $1$. Since $r_j={(\lambda _j-\lambda _-)}/{(\lambda _-+3)}>0$, We have

$$\begin{align*} & 0\le\frac{\lambda_j+3}{\lambda_-{+}3}-1 \le\frac{\frac{-3}{2-k^2}+3}{-1-\frac{2\sqrt{1-k^2+k^4}}{2-k^2}+3}-1\\ & = \frac{1-k^2}{2(k^2+\sqrt{1-k^2+k^4})}=\left(\frac{1}{4}+o(1)\right)(1-k^2) \ \ \textrm{as}\ \ k\to 1. \end{align*}$$

Therefore,

\[ 0\le\frac{r_j}{1-k^2}\le\frac{1}{4}+o(1) \ \ \textrm{as}\ \ k\to 1. \]

Then there exists $r^*_j\in [0,{1}/{4}]$ such that

\[ \lim_{m\to\infty}\frac{r_j(k_m)}{1-k_m^2}=r^*_j,\ \textrm i.e., \ r_j=(r^*_j+o(1))(1-k_m^2)\ \ \textrm{as}\ \ m\to\infty \]

with a suitable monotonically increasing sequence $\{k_m\}_{m=1}^{\infty }$ satisfying $k_m\to 1$ as $m\to \infty$.

Let $\sigma _{\pm }(k_m)$ be defined by (2.32) with $k_m$ and $\lambda _j(k_m)$. Let $\nu _{\pm }(k_m)$ be defined by (2.28) with $k_m$ and $\mu _j(k_m)$. Next we calculate the limit $\lim _{m\to \infty }{(1+\nu _{\pm })}/{(1-k_m^2)}$. By (2.4) we see that

(4.2)$$\begin{align} \lambda_-{+}3=3(1-k^2)-\frac{15}{4}(1-k^2)^2+o((1-k^2)^2)\ \ \textrm{as}\ \ k\to 1. \end{align}$$

Using (4.2), we have

(4.3)$$\begin{align} & \lambda_j-\lambda_-{=}r_j(\lambda_-{+}3) =(r^*_j+o(1))(1-k_m^2)\cdot (3+o(1))(1-k_m^2) \nonumber\\ & \quad =(3r^*_j+o(1))(1-k^2_m)^2. \end{align}$$

Since

\[ 3k^4-(2-k^2)^2(\lambda^2+2\lambda)={-}(2-k^2)^2(\lambda-\lambda_+)(\lambda-\lambda_-), \]

by (4.1) and (4.3) we have

(4.4)$$\begin{align} \frac{1}{(1-k_m^2)\sigma_-} & =\frac{2-k_m^2}{2}\left(\frac{\frac{\lambda_j}{3}+1}{1-k_m^2}-\frac{1}{1-k_m^2}\sqrt{\frac{-(2-k_m^2)^2(\lambda_j-\lambda_+)(\lambda_j-\lambda_-)}{3(2-k_m^2)^2}} \right)\nonumber\\ & =\frac{2-k_m^2}{2}\left(\frac{\frac{\lambda_j}{3}+1}{1-k_m^2}-\frac{1}{1-k_m^2}\right.\nonumber\\ & \left.\quad \sqrt{\frac{-(1+o(1))^2({-}4+o(1))3(r^*_j+o(1))(1-k_m^2)^2}{3}} \right)\nonumber\\ & =\frac{1-2\sqrt{r^*_j}}{2}+o(1) \ \ \textrm{as}\ \ m\to\infty. \end{align}$$

Moreover,

(4.5)$$\begin{align} \frac{1}{\sigma_-}=\left(\frac{1-2\sqrt{r^*_j}}{2}+o(1) \right)(1-k_m^2)\to 0\ \ \textrm{as}\ \ m\to\infty. \end{align}$$

Thus, by (4.5), (4.4) and (2.34) we have

(4.6)$$\begin{align} \nu^*_+:=\lim_{m\to\infty}\frac{1+\nu_+}{1-k_m^2} =\lim_{m\to\infty}\frac{\frac{1}{(1-k_m^2)\sigma_-}-1}{\frac{1}{\sigma_-}-1} =\frac{1+2\sqrt{r^*_j}}{2}. \end{align}$$

Similarly we have

(4.7)$$\begin{align} \nu^*_-:=\lim_{m\to\infty}\frac{1+\nu_-}{1-k_m^2} =\lim_{m\to\infty}\frac{\frac{1}{(1-k_m^2)\sigma_+}-1}{\frac{1}{\sigma_+}-1} =\frac{1-2\sqrt{r^*_j}}{2}. \end{align}$$

By lemma A.5 with (4.6) and (4.7) we have

$$\begin{align*} \lim_{m\to\infty}\mathcal{A} (k_m,\mu_j(k_m)) & =\lim_{m\to\infty}\left|{\rm sgn} (\nu_+)\mathcal{M}(\nu_+,k_m)-{\rm sgn} (\nu_-)\mathcal{M}(\nu_-,k_m)\right|\\ & =\left| ({-}1)\left(\frac{\pi}{2}-\tan^{{-}1}\sqrt{\frac{\nu^*_+}{1-\nu^*_+}}\right) \right.\nonumber\\ & \quad \left. -({-}1)\left(\frac{\pi}{2}-\tan^{{-}1}\sqrt{\frac{\nu^*_-}{1-\nu^*_-}}\right) \right|\\ & =\left|\tan^{{-}1}\sqrt{\frac{1+2\sqrt{r^*_j}}{1-2\sqrt{r^*_j}}} -\tan^{{-}1}\sqrt{\frac{1-2\sqrt{r^*_j}}{1+2\sqrt{r^*_j}}}\right|\\ & =2\tan^{{-}1}\sqrt{\frac{1+2\sqrt{r^*_j}}{1-2\sqrt{r^*_j}}}-\frac{\pi}{2}. \end{align*}$$

Solving the equation

\[ 2\tan^{{-}1}\sqrt{\frac{1+2\sqrt{r^*_j}}{1-2\sqrt{r^*_j}}}-\frac{\pi}{2}=\frac{j\pi}{n} \]

with respect to $r^*_j$, we have $r^*_j=\frac {1}{4}\sin ^2(\frac {j\pi }{n})$, which implies

\[ \lim_{k\to 1}\frac{r_j(k)}{1-k^2}=\frac{1}{4}\sin^2\left(\frac{j\pi}{n}\right). \]

By (4.2) we have

(4.8)$$\begin{align} \lambda_j & ={-}3+(\lambda_-{+}3)+r_j(\lambda_-{+}3)\nonumber\\ & ={-}3+3(1-k^2)+\frac{3}{4}\left(\sin^2\left(\frac{j\pi}{n}\right)-5\right)(1-k^2)^2+o((1-k^2)^2) \ \ \textrm{as}\ \ k\to 1. \end{align}$$

Proof of theorem 1.3(ii) We consider the case $(-3+3k^2)/(2-k^2)<\lambda _j<0$, which corresponds to the case $n/2< j< n$. We define $r_j(k)$ by

\[ r_j:=\frac{\lambda_j}{1-k^2}. \]

It is obvious that $\lambda _j\to 0$ as $k\to 1$. Since

\[ \frac{1}{1-k^2}\frac{-3+3k^2}{2-k^2}<\frac{\lambda_j}{1-k^2}<0, \]

there exists $r_j^*\in [-3,0]$ such that

\[ \lim_{m\to\infty}r_j(k_m)=r^*_j,\ \textrm i.e., \lambda_j=(r^*_j+o(1))(1-k_m^2)\ \ \textrm{as}\ \ m\to\infty \]

with a suitable monotonically increasing sequence $\{k_m\}_{m=1}^{\infty }$ satisfying $k_m\to 1$ as $m\to \infty$.

Let $\sigma _{\pm }(k_m)$ be defined by (2.32) with $k_m$ and $\lambda _j(k_m)$. Let $\nu _{\pm }(k_m)$ be defined by (2.28) with $k_m$ and $\mu _j(k_m)$. We calculate the limit $\lim _{m\to \infty }{(1+\nu _+)}/{(1-k_m^2)}$. Since

\[ D_m:=\frac{3k_m^4-(2-k_m^2)^2(\lambda_j^2+2\lambda_j)}{3(2-k_m^2)^2}\to 1\ \ \textrm{as}\ \ m\to\infty, \]

we have

(4.9)$$\begin{align} \frac{1}{(1-k_m^2)\sigma_-} & =\frac{1}{1-k_m^2}\frac{2-k_m^2}{2}\left(\frac{\lambda_j}{3}+1-\sqrt{D_m}\right)\nonumber\\ & =\frac{2-k_m^2}{6}\frac{\lambda_j}{1-k_m^2} \nonumber\\ & \quad +\frac{2-k_m^2}{2(1+\sqrt{D_m})}\left(\frac{4}{(2-k_m^2)^2}+\frac{2}{3}\frac{\lambda_j}{1-k_m^2}+\frac{1}{3}\frac{\lambda_j}{1-k_m^2}\lambda_j\right)\nonumber\\ & \to 1+\frac{r^*_j}{3}\ \ \textrm{as}\ \ m\to\infty. \end{align}$$

Moreover,

(4.10)$$\begin{align} \frac{1}{\sigma_-}=\left(1+\frac{r^*_j}{3}+o(1)\right)(1-k_m^2)\to 0 \ \ \textrm{as}\ \ m\to\infty. \end{align}$$

Hence, by (4.10), (4.9) and (2.34) we have

(4.11)$$\begin{align} \nu^*_+:=\lim_{m\to\infty}\frac{1+\nu_+}{1-k_m^2} =\lim_{m\to\infty}\frac{\frac{1}{(1-k_m^2)\sigma_-}-1}{\frac{1}{\sigma_-}-1} ={-}\frac{r^*_j}{3}. \end{align}$$

Since

\[ \frac{1}{\sigma_+}=\frac{2-k_m^2}{2}\left(\frac{\lambda_j}{3}+1+\sqrt{D_m}\right)\to 1\ \ \textrm{as}\ \ m\to\infty \]

and $\sigma _+<1$, we have

\[ \nu_-{=}\frac{k_m^2\sigma_+}{1-\sigma_+}\to\infty \ \ \textrm{as}\ \ m\to\infty. \]

Since

\[ 0\le\frac{1}{\nu_-} =\frac{2-k_m^2}{2k_m^2} \frac{\frac{2\lambda_j}{3}\left(\frac{2\lambda_j}{3}+\frac{2-2k_m^2}{2-k_m^2}\right)} {\frac{\lambda_j}{3}-\frac{k_m^2}{2-k_m^2}-\sqrt{D_m}} \le C_0(1-k_m^2), \]

by lemma A.3 we have

$$\begin{align*} & 0\le\frac{K(k_m)}{\sqrt{\nu_-{+}1}}\le\frac{K(k_m)}{\sqrt{\nu_-}}\le \sqrt{C_0(1-k_m^2)}\left(\log\frac{1}{\sqrt{1-k_m^2}}+2\log 2+o(1)\right)\\ & \quad \to 0\ \ \textrm{as}\ \ m\to\infty. \end{align*}$$

Let $\mathcal {A}$ and $\mathcal {M}$ be given in lemma 2.3. By lemma A.5 with (4.11) and lemma A.6 we have

$$\begin{align*} \lim_{m\to\infty}\mathcal{A}(k_m,\mu_j(k_m)) & =\lim_{m\to\infty}\left|{\rm sgn}(\nu_+)\mathcal{M}(\nu_+,k_m)-{\rm sgn}(\nu_-)\mathcal{M}(\nu_-,k_m)\right|\\ & =\lim_{m\to\infty}\left|({-}1)\mathcal{M}(\nu_+,k_m)\right.\\ & \left.\quad -\,\sqrt{\frac{k_m^2+\nu_-}{\nu_-}}\left(J(\nu_-,k_m)+\frac{K(k_m)}{\sqrt{\nu_-{+}1}}\right)\right|\\ & =\left|({-}1)\left(\frac{\pi}{2}-\tan^{{-}1}\sqrt{\frac{\nu^*_+}{1-\nu^*_+}}\right)-\frac{\pi}{2}\right|\\ & =\pi-\tan^{{-}1}\sqrt{\frac{-r^*_j}{3+r^*_j}}. \end{align*}$$

Solving the equation

\[ \pi-\tan^{{-}1}\sqrt{\frac{-r^*_j}{3+r^*_j}}=\frac{j\pi}{n} \]

with respect to $r^*_j$, we have $r^*_j=-3\sin ^2(\frac {j\pi }{n})$, which implies

\[ \lim_{k\to 1}r_j(k)={-}3\sin^2\left(\frac{j\pi}{n}\right). \]

Since $\lambda _j/(1-k^2)=r^*_j+o(1)$, we have

\[ \lambda_j={-}3\sin^2\left(\frac{j\pi}{n}\right)(1-k^2)+o(1-k^2) \ \ \textrm{as}\ \ k\to 1. \]

Proof of theorem 1.3(iii) Let $\lambda _{\pm }$ be defined by (2.4). We consider the case $\lambda _j>\lambda _+$, which corresponds to the case $j>n$. Let $\mathcal {A}$ be given in lemma 2.3 and let $\mu _{\pm }$ be defined by (2.22). By the same argument as in the proof of lemma 3.1(vi) we have

\[ \frac{j\pi}{n}=\mathcal{A}(k,\mu)\ge c\sqrt{\mu_j-\mu_+}K(k). \]

This implies that

\[ \lim_{k\to 1}(\mu_j(k)-\mu_+(k))=0, \]

and in particular $\lim _{k\to 1 }\mu _j(k)=1$,

\[ \lim_{k\to 1}a(k,\mu_j(k))=0,\ \ \lim_{k\to 1}b(k,\mu_j(k))=\frac{1}{3}, \]

where $a(k,\mu _j)$ and $b(k,\mu _j)$ are defined by (2.30). Using

\[ a(k,\mu_j)+(b(k,\mu_j)-1)^2=\frac{1}{9k^4}\mu_j(\mu_j+3), \]

by lemma A.8 we have

$$\begin{align*} \sqrt{\mu_j-\mu_+}K(k) & =\left(\frac{9k^4}{\sqrt{\mu_j(\mu_j+3)(\mu_j+3-3k^2)(\mu_j-\mu_-)}}\frac{j\pi}{n}\right.\\ & \quad\left. -\frac{\sqrt{\mu_j-\mu_+}\tilde{J}(a(k,\mu_j),b(k,\mu_j),k)}{\sqrt{a(k,\mu_j)}}\right)\left\{ a(k,\mu_j)+(1-b(k,\mu_j))^2\right\}\\ & =\sqrt{\frac{\mu_j(\mu_j+3)}{(\mu_j+3-3k^2)(\mu_j-\mu_-)}}\frac{j\pi}{n} \\ & \quad -\frac{2\mu_j(\mu_j+3)\tilde{J}(a(k,\mu_j),b(k,\mu_j),k)}{3\sqrt{3}k^2\sqrt{\mu_j-\mu_-}}\\ & \to \frac{j\pi}{n}-\pi\ \ \textrm{as}\ \ k\to 1. \end{align*}$$

Here we used

\[ \lim_{k\to 1}\sqrt{\frac{\mu_j(\mu_j+3)}{(\mu_j+3-3k^2)(\mu_j-\mu_-)}} =\sqrt{\frac{1(1+3)}{(1+3-3)(1-({-}3))}}=1, \]

\[ \lim_{k\to 1}\frac{2\mu_j(\mu_j+3)\tilde{J}(a(k,\mu_j),b(k,\mu_j),k)}{3\sqrt{3}k^2\sqrt{\mu_j-\mu_-}} =\frac{2\cdot 1(1+3)\frac{\pi}{2\sqrt{\frac{1}{3}}\left(1-\frac{1}{3}\right)}}{3\sqrt{3}\cdot1^2\cdot\sqrt{1-({-}3)}}=\pi. \]

By lemmas 2.1 and (1.3) we see that $1-k^2=16\,{\rm e}^{-2/n\varepsilon }(1+o(1))$ and $K(k)=1/(1+o(1))n\varepsilon$. Hence,

\[ 1-k^2=o\left(\frac{1}{K(k)^2}\right). \]

Using this relation, we have

$$\begin{align*} \lambda_j& =\lambda_+{+}\frac{1}{2-k^2}\frac{(j-n)^2\pi^2}{n^2}\frac{1}{K(k)^2}+o\left(\frac{1}{K(k)^{2}}\right)\\ & =1-3(1-k^2)+o(1-k^2)+\frac{(j-n)^2\pi^2}{n^2}\frac{1}{K(k)^2}(1+o(1))+o\left(\frac{1}{K(k)^{2}}\right)\\ & =1+\frac{(j-n)^2\pi^2}{n^2}\frac{1}{K(k)^2}+o\left(\frac{1}{K(k)^{2}}\right). \end{align*}$$

Proof of corollary 1.7 (i) By (4.8) and lemma 2.1 we have

$$\begin{align*} \lambda_j & ={-}3+(1+r_j)(\lambda_-{+}3)\\ & ={-}3+\left\{ 1+\frac{1}{4}\left(\sin^2\frac{j\pi}{n}\right)(1-k^2)+o(1-k^2)\right\}(\lambda_-{+}3)\\ & ={-}3+\left\{ 1+4\left(\sin^2\frac{j\pi}{n}\right){\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}) \right\}(\lambda^{{\pm}}_{0,\varepsilon}+3). \end{align*}$$

(ii) By theorem 1.3(ii) and lemma 2.1 we have

\[ \lambda_j={-}3\left(\sin^2\frac{j\pi}{n}\right)(1-k^2)+o(1-k^2) ={-}48\left(\sin^2\frac{j\pi}{n}\right){\rm e}^{-\frac{2}{n\varepsilon}}+o({\rm e}^{-\frac{2}{n\varepsilon}}). \]

(iii) By theorem 1.3(iii), corollary 1.6(iii) and (1.3) we have

$$\begin{align*} \lambda_j& =\lambda_+{+}\frac{(j-n)^2\pi^2}{(2-k^2)n^2K(k)^2}+o\left(\frac{1}{K(k)^2}\right) = \lambda_+{+}(j-n)^2\pi^2\varepsilon^2+o(\varepsilon^2)\\ & \quad =1+(j-n)^2\pi^2\varepsilon^2+o(\varepsilon^2). \end{align*}$$