## Spectral asymptotics of Pauli operators and orthogonal polynomials in complex domains

### 3 April 2005

Abstract
We consider the spectrum of a two-dimensional Pauli operator with a compactly supported electric potential and a variable magnetic field with a positive mean value. The rate of accumulation of eigenvalues to zero is described in terms of the logarithmic capacity of the support of the electric potential. A connection between these eigenvalues and orthogonal polynomials in complex domains is established.
Keywords: Pauli operator, magnetic field, spectral asymptotics, logarithmic capacity, orthogonal polynomials

1 Introduction

1. The unperturbed Pauli operator. Let $B=B\left(x\right)$  , $x=\left({x}_{1},{x}_{2}\right)\in {\mathbb{R}}^{2}$  , be a real valued function which has the physical meaning of the strength of a magnetic field in ${\mathbb{R}}^{2}$  .
A two-dimensional non-relativistic spin- $1/2$  particle in the external magnetic field $B$  can be described by the Pauli operator $h=\left(\begin{array}{cc}{h}^{+}& 0\\ 0& {h}^{-}\end{array}\right)\text{in}{L}^{2}\left({\mathbb{R}}^{2}\right)\oplus {L}^{2}\left({\mathbb{R}}^{2}\right).$  The standard approach to the definition of the operators ${h}^{±}$  in ${L}^{2}\left({\mathbb{R}}^{2}\right)$  involves introducing the magnetic vector potential $A\left(x\right)=\left({A}_{1}\left(x\right),{A}_{2}\left(x\right)\right)$  such that $B={\partial }_{{x}_{1}}{A}_{2}-{\partial }_{{x}_{2}}{A}_{1}$  and setting
 $\begin{array}{c}{h}^{±}=\left(-i\nabla -A{\right)}^{2}\mp B.\end{array}$ (1.1)
Instead, we adopt the approach advocated in [6, which consists of defining ${h}^{±}$  in terms of a solution $\Psi =\Psi \left(x\right)$  to the differential equation $\Delta \Psi =B$  . Assume that $B$  is such that a solution $\Psi$  can be chosen subject to the condition
 $\begin{array}{c}\Psi \left(x\right)=\frac{{B}_{0}}{4}|x{|}^{2}+{\Psi }_{1}\left(x\right),{\overline{\Psi }}_{1}={\Psi }_{1}\in {L}^{\infty }\left({\mathbb{R}}^{2}\right),{B}_{0}>0.\end{array}$ (1.2)
Important examples of magnetic fields $B$  of this class are periodic fields with mean value ${B}_{0}$  and constant magnetic fields $B\left(x\right)={B}_{0}$  .
Next, denote, as usual, $\partial =\frac{1}{2}\left({\partial }_{{x}_{1}}-i{\partial }_{{x}_{2}}\right)$  and $\overline{\partial }=\frac{1}{2}\left({\partial }_{{x}_{1}}+i{\partial }_{{x}_{2}}\right)$  . Consider the quadratic forms
 $\begin{array}{c}{\mathfrak{h}}^{+}\left[u\right]=4{\int }_{{\mathbb{R}}^{2}}|\overline{\partial }\left({e}^{\Psi \left(x\right)}u\left(x\right)\right){|}^{2}{e}^{-2\Psi \left(x\right)}dx,{\mathfrak{h}}^{-}\left[u\right]=4{\int }_{{\mathbb{R}}^{2}}|\partial \left({e}^{-\Psi \left(x\right)}u\left(x\right)\right){|}^{2}{e}^{2\Psi \left(x\right)}dx,\end{array}$ (1.3)
which are closed on the domains $Dom\left({\mathfrak{h}}^{±}\right)=\left\{u\in {L}^{2}\left({\mathbb{R}}^{2}\right)|{\mathfrak{h}}^{±}\left[u\right]<\infty \right\}$  . Let us define ${h}^{±}$  as the self-adjoint operators in ${L}^{2}\left({\mathbb{R}}^{2}\right)$  , corresponding to the quadratic forms ${\mathfrak{h}}^{±}$  . For a wide class of magnetic fields this definition is equivalent to the standard definition  1.1 with $A=\left(-{\partial }_{{x}_{2}}\Psi ,{\partial }_{{x}_{1}}\Psi \right)$  ; see [6for a detailed analysis of this issue. In fact, the magnetic field $B$  or the magnetic vector potential $A$  do not enter directly either the definition of Pauli operator or any of our considerations; instead, the “potential function” $\Psi$  becomes the main functional parameter. Note that the condition  1.2 is very close to the `admissibility' condition used in [16.
We will denote by ${\mathfrak{h}}_{0}^{±}$  and ${h}_{0}^{±}$  the above defined forms and operators corresponding to the case of the constant magnetic field $B\left(x\right)={B}_{0}>0$  .
2. Zero modes and the spectral gap. It is well known that Pauli operator $h$  has infinite dimensional kernel. More precisely (cf [1) we have $Ker{h}^{-}=\left\{0\right\}$  and
 $\begin{array}{c}Ker{h}^{+}=\left\{u\in {L}^{2}\left({\mathbb{R}}^{2}\right)|u\left(x\right)=f\left(x\right){e}^{-\Psi \left(x\right)},\overline{\partial }f=0\right\},dimKer{h}^{+}=\infty .\end{array}$ (1.4)
Next, the following well known supersymmetric argument (which has appeared in many forms in the literature; see e.g. [9or [16) establishes the existence of a spectral gap $\left(0,m\right)$  , $m>0$  of the operator $h$  . Let ${a}_{0}$  and ${a}_{0}^{*}$  be the annihilation and creation operators in ${L}^{2}\left({\mathbb{R}}^{2}\right)$  , corresponding to the constant component ${B}_{0}>0$  of the magnetic field $B$  :
 $\begin{array}{c}{a}_{0}=-2i{e}^{-{B}_{0}|x{|}^{2}/4}\overline{\partial }{e}^{{B}_{0}|x{|}^{2}/4},{a}_{0}^{*}=-2i{e}^{{B}_{0}|x{|}^{2}/4}\partial {e}^{-{B}_{0}|x{|}^{2}/4}.\end{array}$ (1.5)
Then one can define
 $\begin{array}{c}a={e}^{-{\Psi }_{1}}{a}_{0}{e}^{{\Psi }_{1}}\text{on}Dom\left(a\right)=\left\{{e}^{-{\Psi }_{1}}u|u\in Dom\left({a}_{0}\right)\right\}\text{and}{a}^{*}={e}^{{\Psi }_{1}}{a}_{0}^{*}{e}^{-{\Psi }_{1}}.\end{array}$ (1.6)
In terms of these operators, we have
 $\begin{array}{c}{h}^{+}={a}^{*}a\text{and}{h}^{-}=a{a}^{*}\end{array}$ (1.7)
and therefore $\sigma \left({h}^{+}\\left\{0\right\}\right)=\sigma \left({h}^{-}\\left\{0\right\}\right)$  . Finally, comparing the form ${\mathfrak{h}}^{-}$  with ${\mathfrak{h}}_{0}^{-}$  , one obtains (see e.g. [3or [16,Proposition1.2) ${\mathfrak{h}}^{-}\left[u\right]⩾{e}^{2essinf{\Psi }_{1}}{\mathfrak{h}}_{0}^{-}\left[{e}^{-{\Psi }_{1}}u\right]⩾2{B}_{0}{e}^{2essinf{\Psi }_{1}}\parallel {e}^{-{\Psi }_{1}}u{\parallel }^{2}⩾2{B}_{0}{e}^{-2osc{\Psi }_{1}}\parallel u{\parallel }^{2},u\in Dom\left({\mathfrak{h}}^{-}\right),$  where $osc{\Psi }_{1}=esssup{\Psi }_{1}-essinf{\Psi }_{1}$  . It follows that $\left(0,m\right)$  , $m=2{B}_{0}{e}^{-2osc{\Psi }_{1}}>0$  , is a gap in the spectrum of ${h}^{+}$  and of $h$  . See [5for a different point of view on the issue of existence of the spectral gap and [18for further progress on this topic.
3. Perturbations of the Pauli operator and spectral asymptotics. Let $v\in {L}^{p}\left({\mathbb{R}}^{2}\right)$  , $p>1$  , be a non-negative compactly supported function, which has the physical meaning of the electric potential. The Pauli operator which describes a particle in the external magnetic field $B$  and electric field with the potential $±v$  , is
 $\begin{array}{c}h±vI=\left(\begin{array}{cc}{h}^{+}±v& 0\\ 0& {h}^{-}±v\end{array}\right)\text{in}{L}^{2}\left({\mathbb{R}}^{2}\right)\oplus {L}^{2}\left({\mathbb{R}}^{2}\right).\end{array}$ (1.8)
The main object of interest in this paper is the spectrum of ${h}^{+}±v$  . In order to define ${h}^{+}±v$  as a quadratic form sum, let us establish that $v$  is ${h}^{+}$  -form compact. By  1.6 ,  1.7 and boundedness of ${\Psi }_{1}$  , we see that $v$  is ${h}^{+}$  -form compact if and only if $v{e}^{-2{\Psi }_{1}}$  is ${h}_{0}^{+}$  -form compact. As $v{e}^{-2{\Psi }_{1}}\in {L}^{p}\left({\mathbb{R}}^{2}\right)$  , $p>1$  , we obtain that $v{e}^{-2{\Psi }_{1}}$  is ${h}_{0}^{+}$  -form compact (see [2).
By the above established relative compactness, the essential spectra of ${h}^{+}+v$  , ${h}^{+}-v$  , and ${h}^{+}$  coincide; moreover, due to the assumption $v⩾0$  , the eigenvalues of ${h}^{+}+v$  can accumulate to $0$  only from above, and the eigenvalues of ${h}^{+}-v$  can do so only from below.
Let ${\lambda }_{1}^{-}⩽{\lambda }_{2}^{-}⩽\cdots$  be the negative eigenvalues of ${h}^{+}-v$  , and ${\lambda }_{1}^{+}⩾{\lambda }_{2}^{+}⩾\cdots$  be the eigenvalues of ${h}^{+}+v$  in the spectral gap $\left(0,m\right)$  ; here and in the rest of the paper, we assume eigenvalues to be enumerated with multiplicities taken into account. The main aim of this paper is to describe the rate of convergence ${\lambda }_{n}^{±}\to 0$  as $n\to \infty$  . Roughly speaking, we prove the following asymptotics (precise statements are given in Section  2 ):
 $\begin{array}{c}log\left(±n!{\lambda }_{n}^{±}\right)=nlog\left({B}_{0}/2\right)+2nlogCap\left(suppv\right)+o\left(n\right),n\to \infty ,\end{array}$ (1.9)
where $Cap$  is the logarithmic capacity of a set. The notion of logarithmic capacity is introduced in the framework of potential theory; see e.g. [7, 11. Recall that the logarithmic capacity of compact sets in ${\mathbb{R}}^{2}$  has the following properties:
(i) if ${\Omega }_{1}\subset {\Omega }_{2}$  then $Cap{\Omega }_{1}⩽Cap{\Omega }_{2}$  ; (ii) $Cap\Omega$  coincides with the logarithmic capacity of the outer boundary of $\Omega$  ( $=$  the boundary of the unbounded component of ${\mathbb{R}}^{2}\\Omega$  ); (iii) the logarithmic capacity of a disc of radius $r$  is $r$  ; (iv) if ${\Omega }_{2}=\left\{\alpha x|x\in {\Omega }_{1}\right\}$  , $\alpha >0$  , then $Cap{\Omega }_{2}=\alpha Cap{\Omega }_{1}$  .
We establish  1.9 by means of the following simple chain of equivalent reformulations of the problem. Firstly, a perturbation theory argument reduces the problem to the spectral asymptotics of an auxiliary compact self-adjoint operator ${P}_{0}v{P}_{0}$  , where ${P}_{0}$  is the spectral projection of ${h}^{+}$  , corresponding to the eigenvalue $0$  . Next, we observe that the eigenvalues of ${P}_{0}v{P}_{0}$  coincide with the singular numbers of a certain embedding operator (see  4.1 ).
Using the approach of [13, we relate the singular numbers of this embedding operator to some sequence of orthogonal polynomials in the complex domain (see below). Finally, application of the results of [22concerning the asymptotics of these orthogonal polynomials leads to  1.9 .
Using the same technique, we are also able to treat two similar problems. First, we consider the Pauli operators ${h}_{0}^{+}±v$  in the case of a constant magnetic field and describe the rate of accumulation of the eigenvalues to the higher Landau levels. Secondly, we consider the three-dimensional Pauli Hamiltonian with a constant magnetic field and a compactly supported electric potential and describe the rate of convergence of eigenvalues to 0. These results are presented in Section 2.
The rate of convergence of eigenvalues to zero for Pauli operators in dimensions two and three was investigated before in the case of constant magnetic field $B\left(x\right)={B}_{0}>0$  for various classes of potentials $v$  with power or exponential decay at infinity; see [21, 19, 23, 14, 15, 10, 17.
We refer the reader to the discussion in [17. The case of a constant magnetic field and compactly supported potentials $v$  was considered in [17and [12. The case of a two-dimensional operator with variable magnetic field and potentials $v$  with power or exponential decay and also with compactly supported potentials was treated in [16. The results of [17, 12, 16for the case of compactly supported potentials read as
 $\begin{array}{c}log\left(±{\lambda }_{n}^{±}\right)=-nlogn+O\left(n\right),n\to \infty .\end{array}$ (1.10)
As far as we are aware, a connection between the spectral asymptotics of magnetic operators and logarithmic capacity or orthogonal polynomials has not been made before.
Some physical intuition concerning these problem with constant magnetic field can be gained from [8.
4. Orthogonal and Chebyshev polynomials. We identify ${\mathbb{R}}^{2}$  and $\mathbb{C}$  in a standard way: $z={x}_{1}+i{x}_{2}$  for $\left({x}_{1},{x}_{2}\right)\in {\mathbb{R}}^{2}$  , denote by $dm\left(z\right)$  the Lebesgue measure in $\mathbb{C}$  and consider $v$  as a function of $z$  . It appears that the sequence of polynomials in $z$  , orthogonal with respect to the measure $v\left(z\right)dm\left(z\right)$  , is related to the asymptotics of ${\lambda }_{n}^{±}$  . Here we present necessary facts from the theory of Chebyshev and orthogonal polynomials in complex domains; see e.g. [7and [22for the details. For any $n=0,1,2,...$  , let ${\mathcal{P}}_{n}$  be the set of all monic polynomials in $z$  of degree $n$  :
 $\begin{array}{c}{\mathcal{P}}_{n}=\left\{{z}^{n}+{a}_{n-1}{z}^{n-1}+\cdot \cdot \cdot +{a}_{1}z+{a}_{0}|{a}_{0},...,{a}_{n-1}\in \mathbb{C}\right\}.\end{array}$ (1.11)
Let $\Omega \subset \mathbb{C}$  be a compact set. For a fixed $n$  , consider the problem of minimization of the norm $\parallel t{\parallel }_{C\left(\Omega \right)}\equiv {sup}_{z\in \Omega }|t\left(z\right)|$  on the set $t\in {\mathcal{P}}_{n}$  . It is clear that the minimum is positive and attained at some polynomial ${t}_{n}\in {\mathcal{P}}_{n}$  . The polynomial ${t}_{n}$  is called the $n$  'th Chebyshev polynomial for the set $\Omega$  . (One can prove that such a polynomial is unique, but we will not need this fact). It is well known that all zeros of ${t}_{n}$  lie in the closed convex hull of $\Omega$  . The $n$  'th root asymptotics of ${t}_{n}$  is given by
 $\begin{array}{c}{lim}_{n\to \infty }\parallel {t}_{n}{\parallel }^{1/n}=Cap\Omega .\end{array}$ (1.12)
Next, let $v\in {L}^{1}\left(\mathbb{C},dm\right)$  be a non-negative compactly supported function. Denote
 $\begin{array}{c}{M}_{n}\left(v\right)={inf}_{p\in {\mathcal{P}}_{n}}{\int }_{\mathbb{C}}|p\left(z\right){|}^{2}v\left(z\right)dm\left(z\right)={\int }_{\mathbb{C}}|{p}_{n}\left(z\right){|}^{2}v\left(z\right)dm\left(z\right),\end{array}$ (1.13)
where the sequence $\left\{{p}_{n}{\right\}}_{n=0}^{\infty }$  , ${p}_{n}\in {\mathcal{P}}_{n}$  , is obtained by applying the Gram–Schmidt orthogonalisation process in ${L}^{2}\left(\mathbb{C},v\left(z\right)dm\left(z\right)\right)$  to the sequence $1,z,{z}^{2},...$  . All zeros of ${p}_{n}$  lie in the closed convex hull of $suppv$  .
Regarding the $n$  'th root asymptotics of ${p}_{n}$  , the following facts are known (see [22).
Denote
 $\begin{array}{c}{\rho }_{+}\left(v\right)={limsup}_{n\to \infty }{M}_{n}\left(v{\right)}^{1/n},{\rho }_{-}\left(v\right)={liminf}_{n\to \infty }{M}_{n}\left(v{\right)}^{1/n}.\end{array}$ (1.14)
In general, it can happen that ${\rho }_{-}\left(v\right)<{\rho }_{+}\left(v\right)$  (see the proof of Theorem 1.1.9 in [22). One has the estimates
 $\begin{array}{c}\begin{array}{}\text{(1.15)}& {\rho }_{+}\left(v\right)⩽\left(Capsuppv{\right)}^{2},{\rho }_{-}\left(v\right)⩾\left(Cap{\Omega }_{-}\left(v\right){\right)}^{2},{\Omega }_{-}\left(v\right)=\left\{z\in \mathbb{C}|{limsup}_{r\to +0}\frac{log{\int }_{|z-\zeta |⩽r}v\left(\zeta \right)dm\left(\zeta \right)}{logr}<\infty \right\}.\end{array}\end{array}$ (1.16)
The inequality  1.15 is a part of Corollary 1.1.7 of [22. The inequality  1.16 , although not stated explicitly in [22, follows directly from the proof of Theorem 4.2.1 therein.
Remark 1.1. Let $\Omega \subset \mathbb{C}$  be a compact set with a Lipschitz boundary, and let $v\in {L}^{1}\left(\mathbb{C},dm\right)$  be such that $v\left(z\right)⩾c>0$  for all $z\in \Omega$  and $v\left(z\right)=0$  for all $z\in \mathbb{C}\\Omega$  . Then we easily find that ${\Omega }_{-}\left(v\right)=\Omega =suppv$  and therefore ${\rho }_{+}\left(v\right)={\rho }_{-}\left(v\right)=\left(Cap\Omega {\right)}^{2}$  .

2 Main results

1. Two-dimensional Pauli operators with variable magnetic field. Let ${h}^{+}$  , as in the Introduction, be the Pauli operator defined via  1.3 with $\Psi$  subject to  1.2 . Let $v$  and ${\lambda }_{n}^{±}$  be as in the Introduction.
Theorem 2.1. Let $0⩽v\in {L}^{p}\left({\mathbb{R}}^{2}\right)$  , $p>1$  , be a compactly supported potential and let ${M}_{n}\left(v\right)$  be as defined in  1.13 . Then there exists $k\in \mathbb{N}$  such that
 $\begin{array}{cc}\left({B}_{0}/2\right){M}_{n+k}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right)⩽\left(n!{\lambda }_{n}^{+}{\right)}^{1/n}⩽\left({B}_{0}/2\right){M}_{n-1}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right),& \end{array}$ (2.1)
 $\begin{array}{cc}\left({B}_{0}/2\right){M}_{n-1}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right)⩽\left(-n!{\lambda }_{n}^{-}{\right)}^{1/n}⩽\left({B}_{0}/2\right){M}_{n-k}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right),& \end{array}$ (2.2)
 $\begin{array}{}\end{array}$
as $n\to \infty$  . In particular, ${limsup}_{n\to \infty }\left(±n!{\lambda }_{n}^{±}{\right)}^{1/n}={B}_{0}{\rho }_{+}\left(v\right)/2,{liminf}_{n\to \infty }\left(±n!{\lambda }_{n}^{±}{\right)}^{1/n}={B}_{0}{\rho }_{-}\left(v\right)/2,$  where ${\rho }_{±}\left(v\right)$  are defined by  1.14 . If $v$  is of the class described in Remark  1.1 , then the asymptotics  1.9 holds true.
Remark 2.2. Let $\mu$  be a compactly supported finite measure in ${\mathbb{R}}^{2}$  such that the quadratic form ${\int }_{\mathbb{C}}|u\left(x\right){|}^{2}d\mu \left(x\right)$  is compact with respect to the quadratic form ${\mathfrak{h}}^{+}$  .
Then one can define the self-adjoint operators corresponding to the quadratic forms ${\mathfrak{h}}^{+}\left[u\right]±{\int }_{{\mathbb{R}}^{2}}|u\left(x\right){|}^{2}d\mu \left(x\right).$  All our considerations remain valid for such operators. For example, the case of a measure $\mu$  , supported by a curve, can be interesting.
2. Two-dimensional Pauli operators with constant magnetic field. Let $B\left(x\right)={B}_{0}>0$  ; consider the corresponding operator ${h}_{0}^{+}$  . As it is well known, the spectrum of ${h}_{0}^{+}$  consists of the eigenvalues $\left\{2q{B}_{0}{\right\}}_{q=0}^{\infty }$  of infinite multiplicities; these eigenvalues are known as Landau levels. Consider the problem of accumulation of eigenvalues of ${h}_{0}^{+}±v$  to a fixed higher Landau level $2q{B}_{0}$  , $q⩾1$  . Let ${\lambda }_{q,1}^{-}⩽{\lambda }_{q,2}^{-}⩽\cdots$  be the eigenvalues of ${h}_{0}^{+}-v$  in the interval $\left(2\left(q-1\right){B}_{0},2q{B}_{0}\right)$  , and let ${\lambda }_{q,1}^{+}⩾{\lambda }_{q,2}^{+}⩾\cdots$  be the eigenvalues of ${h}_{0}^{+}+v$  in $\left(2q{B}_{0},2\left(q+1\right){B}_{0}\right)$  .
Theorem 2.3. Let $\Omega \subset {\mathbb{R}}^{2}$  be a compact set with Lipschitz boundary and let $v\in {L}^{p}\left({\mathbb{R}}^{2}\right)$  , $p>1$  , be such that $v\left(x\right)⩾c>0$  for $x\in \Omega$  and $v\left(x\right)=0$  for $x\in {\mathbb{R}}^{2}\\Omega$  .
Then for the corresponding eigenvalues ${\lambda }_{q,n}^{±}$  we have: ${lim}_{n\to \infty }\left(±n!\left({\lambda }_{q,n}^{±}-2q{B}_{0}\right){\right)}^{1/n}=\frac{{B}_{0}}{2}\left(Cap\Omega {\right)}^{2}.$
The rate of convergence of ${\lambda }_{q,n}^{±}\to 2q{B}_{0}$  as $n\to \infty$  , $q⩾1$  , was studied before in [17, 12, where the asymptotics $log\left(±\left({\lambda }_{q,n}^{±}-2q{B}_{0}\right)\right)=-nlogn+O\left(n\right),n\to \infty$  was obtained. Note that if the potential $v$  depends only on $|x|$  , then the result of Theorem  2.3 can be obtained by a direct calculation using separation of variables, see e.g. [17,Proposition3.2.
3. Three-dimensional Pauli operator with a constant magnetic field. Let $H=\left(-i\nabla -\mathbf{A}\left(\mathbf{x}\right){\right)}^{2}-{B}_{0}\text{in}{L}^{2}\left({\mathbb{R}}^{3},d\mathbf{x}\right),\mathbf{x}=\left({x}_{1},{x}_{2},{x}_{3}\right),$  where $\mathbf{A}\left(\mathbf{x}\right)=\left(-\frac{1}{2}{B}_{0}{x}_{2},\frac{1}{2}{B}_{0}{x}_{1},0\right)$  . It is well known that the spectrum of $H$  is absolutely continuous and coincides with the interval $\left[0,\infty \right)$  . The background information concerning the spectral theory of $H$  and its perturbations can be found in [2.
Let $V\in {L}^{3/2}\left({\mathbb{R}}^{3}\right)$  be a non-negative compactly supported potential. The operator of multiplication by $V$  in ${L}^{2}\left({\mathbb{R}}^{3}\right)$  is $H$  -form compact (cf. [2). Thus, one can define the self-adjoint opertor $H-V$  via the corresponding quadratic form; the essential spectrum of $H-V$  is also $\left[0,\infty \right)$  . Let ${\Lambda }_{1}⩽{\Lambda }_{2}⩽\cdots$  be the negative eigenvalues of $H-V$  ; we have ${\Lambda }_{n}\to 0$  as $n\to \infty$  . Below we describe the asymptotic behaviour of ${\Lambda }_{n}$  as $n\to \infty$  in terms of the auxiliary weight function
 $\begin{array}{c}w\left({x}_{1},{x}_{2}\right)={\int }_{-\infty }^{\infty }V\left({x}_{1},{x}_{2},{x}_{3}\right)d{x}_{3},\text{a.e.}\left({x}_{1},{x}_{2}\right)\in {\mathbb{R}}^{2}.\end{array}$ (2.3)
As above, we consider $w$  as a function of $z={x}_{1}+i{x}_{2}$  .
Theorem 2.4. Let $0⩽V\in {L}^{3/2}\left({\mathbb{R}}^{3}\right)$  be a compactly supported potential and $w$  be defined by  2.3 . Then there exists $k\in \mathbb{N}$  such that
 $\begin{array}{c}\left({B}_{0}/2{\right)}^{2}{M}_{n+k}\left(w{\right)}^{2/n}\left(1+o\left(1\right)\right)⩽\left(-\left(n!{\right)}^{2}{\Lambda }_{n}{\right)}^{1/n}⩽\left({B}_{0}/2{\right)}^{2}{M}_{n-k}\left(w{\right)}^{2/n}\left(1+o\left(1\right)\right),\end{array}$ (2.4)
as $n\to \infty$  . In particular, ${limsup}_{n\to \infty }\left(-\left(n!{\right)}^{2}{\Lambda }_{n}{\right)}^{1/n}=\left({B}_{0}{\rho }_{+}\left(w\right)/2{\right)}^{2},{liminf}_{n\to \infty }\left(-\left(n!{\right)}^{2}{\Lambda }_{n}{\right)}^{1/n}=\left({B}_{0}{\rho }_{-}\left(w\right)/2{\right)}^{2},$  where ${\rho }_{±}\left(w\right)$  are defined by  1.14 .
The rate of accumulation ${\Lambda }_{n}\to 0$  for potentials $V$  with power or exponential decay was considered before in [21, 19, 23, 14, 15, 10, 17. For compactly supported potentials, this problem was considered in [17, 12, where the asymptotics $log{\Lambda }_{n}=-2nlogn+O\left(n\right),n\to \infty$  was obtained.
Remark 2.5. Theorem  2.4 remains valid under the following assumptions on $V$  : (i) $V⩾0$  , $V$  is $H$  -form compact; (ii) ${\int }_{{\mathbb{R}}^{3}}V\left(\mathbf{x}\right)\left(1+|{x}_{3}{|}^{2}\right)d\mathbf{x}<\infty$  ; (iii) the function $w$  , defined by  2.3 , is compactly supported.

3 Proof of Theorems  2.1 ,  2.3 and  2.4

Proof of Theorem  2.1 : Let ${\mathcal{ℋ}}_{0}\subset {L}^{2}\left({\mathbb{R}}^{2}\right)$  be the kernel of ${h}^{+}$  , and let ${P}_{0}$  be the corresponding eigenprojection, $Ran{P}_{0}={\mathcal{ℋ}}_{0}$  . Consider the compact self-adjoint operator ${P}_{0}v{P}_{0}$  . The key ingredient in the proof is the following
Lemma 3.1. Let $v\in {L}^{1}\left({\mathbb{R}}^{2}\right)$  be a non-negative compactly supported function and let ${s}_{1}⩾{s}_{2}⩾\cdots >0$  be the eigenvalues of ${P}_{0}v{P}_{0}$  . Then
 $\begin{array}{c}\left(n!{s}_{n+1}{\right)}^{1/n}=\left({B}_{0}/2\right){M}_{n}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right),n\to \infty ,\end{array}$ (3.1)
where ${M}_{n}\left(v\right)$  are defined by  1.13 .
The proof is given in Section  4 . Now it remains to employ a perturbation theory argument (see [16,Proposition3.1or [17,Proposition4.1) based on the Birman-Schwinger principle and on Weyl inequalities for eigenvalues of a sum of compact operators. This argument shows that there exists $k\in \mathbb{N}$  such that for all sufficiently large $n\in \mathbb{N}$  one has ${s}_{n}⩽-{\lambda }_{n}^{-}⩽2{s}_{n-k}.\frac{1}{2}{s}_{n+k}⩽{\lambda }_{n}^{+}⩽{s}_{n}.$  Combining these inequalities with Lemma  3.1 , we obtain the required result.
Proof of Theorem  2.3 : For any $q⩾0$  , denote ${\mathcal{ℋ}}_{q}=Ker\left({h}^{+}-2q{B}_{0}\right)$  and let ${P}_{q}$  be the eigenprojection of ${h}_{0}^{+}$  corresponding to the eigenvalue $2q{B}_{0}$  . Consider the compact self-adjoint operator ${P}_{q}v{P}_{q}$  , and let ${s}_{1}^{\left(q\right)}⩾{s}_{2}^{\left(q\right)}⩾\cdots$  be the eigenvalues of this operator.
As in the proof of Theorem  2.1 , using a perturbation theory argument based on the Birman-Schwinger principle and Weyl inequalities (see [17,Proposition 4.1), one shows that there exists $k\in \mathbb{N}$  such that for all sufficiently large $n\in \mathbb{N}$  ,
 $\begin{array}{c}\frac{1}{2}{s}_{n+k}^{\left(q\right)}⩽±\left({\lambda }_{q,n}^{±}-2qB\right)⩽2{s}_{n-k}^{\left(q\right)}.\end{array}$ (3.2)
Now the proof of Theorem  2.3 reduces to
Lemma 3.2. Let $\Omega \subset {\mathbb{R}}^{2}$  be a compact set with Lipschitz boundary and let $v\in {L}^{1}\left({\mathbb{R}}^{2}\right)$  be such that $v\left(x\right)⩾c>0$  for $x\in \Omega$  and $v\left(x\right)=0$  for $x\in {\mathbb{R}}^{2}\\Omega$  . Fix $q\in \mathbb{N}$  and let ${s}_{1}^{\left(q\right)}⩾{s}_{2}^{\left(q\right)}⩾\cdots$  be the eigenvalues of ${P}_{q}v{P}_{q}$  . Then one has
 $\begin{array}{c}{lim}_{n\to \infty }\left(n!{s}_{n}^{\left(q\right)}{\right)}^{1/n}=\left({B}_{0}/2\right)\left(Cap\Omega {\right)}^{2}.\end{array}$ (3.3)
The proof of Lemma  3.2 is given in Section  5 . From Lemma  3.2 and the estimate  3.2 , we immediately obtain the required result.
Proof of Theorem  2.4 : The proof repeats almost word for word the construction of [19.
According to the Birman-Schwinger principle, for $E>0$  we have:
 $\begin{array}{c}♯\left\{n|{\Lambda }_{n}<-E\right\}={n}_{+}\left(1;\sqrt{V}\left({H}_{0}+E{\right)}^{-1}\sqrt{V}\right).\end{array}$ (3.4)
The operator $\sqrt{V}\left({H}_{0}+E{\right)}^{-1}\sqrt{V}$  can be represented as $\sqrt{V}\left({H}_{0}+E{\right)}^{-1}\sqrt{V}=\frac{1}{2\sqrt{E}}{K}_{1}+{K}_{2}+{K}_{3}.$  Here ${K}_{1}$  , ${K}_{2}$  are the operators in ${L}^{2}\left({\mathbb{R}}^{3}\right)$  with the integral kernels
 $\begin{array}{cc}{K}_{1}\left(\mathbf{x},\mathbf{y}\right)& =\sqrt{V\left(\mathbf{x}\right)}{P}_{0}\left({x}_{\perp },{y}_{\perp }\right)\sqrt{V\left(\mathbf{y}\right)},\end{array}$
 $\begin{array}{cc}{K}_{2}\left(\mathbf{x},\mathbf{y}\right)& =\sqrt{V\left(\mathbf{x}\right)}{P}_{0}\left({x}_{\perp },{y}_{\perp }\right)\frac{{e}^{-\sqrt{E}|{x}_{3}-{y}_{3}|}-1}{2\sqrt{E}}\sqrt{V\left(\mathbf{y}\right)},\end{array}$
 $\begin{array}{}\end{array}$
where the notation ${x}_{\perp }=\left({x}_{1},{x}_{2}\right)$  , ${y}_{\perp }=\left({y}_{1},{y}_{2}\right)$  is used, and ${P}_{0}\left({x}_{\perp },{y}_{\perp }\right)$  is the integral kernel of the operator ${P}_{0}$  in ${L}^{2}\left({\mathbb{R}}^{2}\right)$  . Finally, ${K}_{3}$  is the operator ${K}_{3}=\sqrt{V}{Q}_{0}\left({H}_{0}+E{\right)}^{-1}\sqrt{V},$  where ${Q}_{0}=\left(I-{P}_{0}\right)\otimes I$  in the decomposition ${L}^{2}\left({\mathbb{R}}^{3},d{x}_{1}d{x}_{2}d{x}_{3}\right)={L}^{2}\left({\mathbb{R}}^{2},d{x}_{1}d{x}_{2}\right)\otimes {L}^{2}\left(\mathbb{R},d{x}_{3}\right)$  .
The operators ${K}_{2}$  and ${K}_{3}$  have limits (in the operator norm) as $E\to +0$  ; these limits are compact self-adjoint operators. Thus, by the Weyl's inequalities for eigenvalues (see e.g. [4), we have for $E\to +0$  :
 $\begin{array}{c}\begin{array}{}\text{(3.5)}& {n}_{+}\left(1;\sqrt{V}\left({H}_{0}+E{\right)}^{-1}\sqrt{V}\right)⩽{n}_{+}\left(\frac{1}{2};\frac{1}{2\sqrt{E}}{K}_{1}\right)+{n}_{+}\left(\frac{1}{2};{K}_{2}+{K}_{3}\right)⩽{n}_{+}\left(\sqrt{E};{K}_{1}\right)+O\left(1\right),{n}_{+}\left(1;\sqrt{V}\left({H}_{0}+E{\right)}^{-1}\sqrt{V}\right)⩾{n}_{+}\left(\frac{3}{2};\frac{1}{2\sqrt{E}}{K}_{1}\right)-{n}_{+}\left(\frac{1}{2};-{K}_{2}-{K}_{3}\right)⩾{n}_{+}\left(3\sqrt{E};{K}_{1}\right)-O\left(1\right).\end{array}\end{array}$ (3.6)
Finally, again as in [19, let us prove that the non-zero eigenvalues of ${K}_{1}$  coincide with those of ${P}_{0}w{P}_{0}$  , where $w$  is defined by  2.3 . It suffices to prove this statement for continuous $V$  with compact support; the general case $V\in {L}^{3/2}$  then follows by approximation argument. Let ${N}_{1}:{L}^{2}\left({\mathbb{R}}^{3},d{x}_{1}d{x}_{2}d{x}_{3}\right)\to {L}^{2}\left({\mathbb{R}}^{2},d{x}_{1}d{x}_{2}\right)$  and ${N}_{2}:{L}^{2}\left({\mathbb{R}}^{2},d{x}_{1}d{x}_{2}\right)\to {L}^{2}\left({\mathbb{R}}^{3},d{x}_{1}d{x}_{2}d{x}_{3}\right)$  be the following operators:
$\left({N}_{1}u\right)\left({x}_{1},{x}_{2}\right)={\int }_{-\infty }^{\infty }{V}^{1/2}\left({x}_{1},{x}_{2},{x}_{3}\right)u\left({x}_{1},{x}_{2},{x}_{3}\right)d{x}_{3},\left({N}_{2}u\right)\left({x}_{1},{x}_{2},{x}_{3}\right)={V}^{1/2}\left({x}_{1},{x}_{2},{x}_{3}\right)u\left({x}_{1},{x}_{2}\right).$
Then ${K}_{1}={N}_{2}{P}_{0}{N}_{1}=\left({N}_{2}{P}_{0}\right)\left({P}_{0}{N}_{1}\right)$  and ${P}_{0}w{P}_{0}=\left({P}_{0}{N}_{1}\right)\left({N}_{2}{P}_{0}\right)$  . It follows that the non-zero eigenvalues of ${K}_{1}$  coincide with $\left\{{s}_{n}\right\}$  , the non-zero eigenvalues of ${P}_{0}w{P}_{0}$  , and so ${n}_{+}\left(\sqrt{E};{K}_{1}\right)=♯\left\{n|\left({s}_{n}{\right)}^{2}>E\right\}$  . From here and  3.4 ,  3.5 ,  3.6 it follows that for some $k\in \mathbb{N}$  and all sufficiently large $n\in \mathbb{N}$  , one has
 $\begin{array}{c}\frac{1}{9}\left({s}_{n+k}{\right)}^{2}⩽{\Lambda }_{n}⩽\left({s}_{n-k}{\right)}^{2}.\end{array}$ (3.7)
Combining this with Lemma  3.1 , we get the statement of Theorem  2.4 .

4 Proof of Lemma  3.1

First let us consider the case of a constant magnetic field $B\left(x\right)={B}_{0}>0$  . Let ${F}^{2}$  be the Hilbert space of all entire functions $f$  such that $\parallel f{\parallel }_{{F}^{2}}^{2}={\int }_{\mathbb{C}}|f\left(z\right){|}^{2}{e}^{-{B}_{0}|z{|}^{2}/2}dm\left(z\right)<\infty .$  In the case ${B}_{0}=2$  , the space ${F}^{2}$  is usually called Fock space or Segal-Bargmann space.
By  1.4 , we have an isometry between ${\mathcal{ℋ}}_{0}=Ker{h}^{+}\subset {L}^{2}\left(\mathbb{C},dm\right)$  and ${F}^{2}$  , given by $u\left(z\right)={e}^{-{B}_{0}|z{|}^{2}/4}f\left(z\right)$  , $u\in {\mathcal{ℋ}}_{0}$  , $f\in {F}^{2}$  . Thus, the quadratic form of the operator ${P}_{0}v{P}_{0}{|}_{{\mathcal{ℋ}}_{0}}$  is unitarily equivalent to the quadratic form ${\int }_{\mathbb{C}}|f\left(z\right){|}^{2}v\left(z\right){e}^{-{B}_{0}|z{|}^{2}/2}dm\left(z\right),f\in {F}^{2}.$  It follows that the non-zero eigenvalues ${s}_{n}$  of ${P}_{0}v{P}_{0}$  coincide with the singular values ${\mu }_{n}$  of the embedding operator
 $\begin{array}{c}{F}^{2}\subset {L}^{2}\left(\mathbb{C},v\left(z\right){e}^{-{B}_{0}|z{|}^{2}/2}dm\left(z\right)\right).\end{array}$ (4.1)
The case of a variable magnetic field can be also reduced to the embedding  4.1 . Indeed, using the boundedness of ${\Psi }_{1}$  , one obtains (see [16,Proposition3.2):
${\mu }_{n}{e}^{-2osc{\Psi }_{1}}⩽{s}_{n}⩽{\mu }_{n}{e}^{2osc{\Psi }_{1}},n\in \mathbb{N}.$  Thus, it remains to prove the asymptotic formula
 $\begin{array}{c}\left(n!{\mu }_{n+1}{\right)}^{1/n}=\left({B}_{0}/2\right){M}_{n}\left(v{\right)}^{1/n}\left(1+o\left(1\right)\right),n\to \infty \end{array}$ (4.2)
for the singular values ${\mu }_{n}$  of the embedding  4.1 . We shall assume ${B}_{0}=2$  ; the general case can be reduced to this one by a linear change of coordinates.
Asymptotics of the $n$  -widths of the embedding ${F}^{2}\subset C\left(\Omega \right)$  , where $\Omega$  is a compact set in $\mathbb{C}$  , was studied in [13. Below we repeat the arguments of [13(with trivial modifications) to obtain the required asymptotics.
By the minimax principle, we have the following variational characterisation of ${\mu }_{n}$  :
 $\begin{array}{cc}{\mu }_{n+1}& ={inf}_{{L}_{n}^{+}\subset {F}^{2}}{sup}_{f\in {L}_{n}^{+}\\left\{0\right\}}\frac{{\int }_{\mathbb{C}}|f\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)}{\parallel f{\parallel }_{{F}^{2}}^{2}},codim{L}_{n}^{+}=n,\end{array}$ (4.3)
 $\begin{array}{cc}{\mu }_{n+1}& ={sup}_{{L}_{n}^{-}\subset {F}^{2}}{inf}_{f\in {L}_{n}^{-}\\left\{0\right\}}\frac{{\int }_{\mathbb{C}}|f\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)}{\parallel f{\parallel }_{{F}^{2}}^{2}},dim{L}_{n}^{-}=n+1.\end{array}$ (4.4)
 $\begin{array}{}\end{array}$
1. Upper bound on ${\mu }_{n+1}$  . For the subspaces ${L}_{n}^{+}$  from  4.3 , we will take ${L}_{n}^{+}=\left\{f\in {F}^{2}|f\left(z\right)={p}_{n}\left(z\right)g\left(z\right),\text{}g\text{is entire function}\right\},$  where ${p}_{n}$  is the sequence of monic polynomials orthogonal with respect to the measure $v\left(z\right)dm\left(z\right)$  .
In order to estimate the ratio in  4.3 from above, let us prove the following auxiliary statement. Denote ${R}_{0}={max}_{z\in suppv}|z|$  . We claim that for any $\varepsilon \in \left(0,\frac{1}{3}\right)$  , there exists $N\in \mathbb{N}$  such that for all $n⩾N$  and any $f={p}_{n}g\in {L}_{n}^{+}$  , we have
 $\begin{array}{c}{sup}_{|z|⩽{R}_{0}}|g\left(z\right){|}^{2}⩽\left(1-\varepsilon {\right)}^{-2n}\frac{1}{n!}\parallel {p}_{n}g{\parallel }_{{F}^{2}}^{2}.\end{array}$ (4.5)
Indeed, we have $g\left(z\right)=\frac{1}{2\pi i}{\int }_{|\zeta |=r}\frac{f\left(\zeta \right)}{{p}_{n}\left(\zeta \right)\left(\zeta -z\right)}d\zeta ,r>{R}_{0},$  and therefore ${sup}_{|z|⩽{R}_{0}}|g\left(z\right){|}^{2}⩽\frac{r}{2\pi }{sup}_{|z|⩽{R}_{0}}{\int }_{|\zeta |=r}\frac{|f\left(\zeta \right){|}^{2}}{|{p}_{n}\left(\zeta \right){|}^{2}|\zeta -z{|}^{2}}d|\zeta |$  for any $r>{R}_{0}$  . Denote $R={R}_{0}/\varepsilon$  . Since all zeros of ${p}_{n}$  lie in the closed convex hull of $suppv$  , we obtain:
$|{p}_{n}\left(\zeta \right)||\zeta -z|⩾\left(\left(1-\varepsilon \right)r{\right)}^{n+1},|z|⩽{R}_{0},|\zeta |=r⩾R.$  Thus, we get ${sup}_{|z|⩽{R}_{0}}|g\left(z\right){|}^{2}⩽\frac{{r}^{-2n-1}}{2\pi \left(1-\varepsilon {\right)}^{2n+2}}{\int }_{|\zeta |=r}|f\left(\zeta \right){|}^{2}d|\zeta |,r⩾R.$  Integrating the last inequality over $r$  from $R$  to $\infty$  with the weight ${e}^{-{r}^{2}}{r}^{2n+1}$  , and using the fact that ${\int }_{R}^{\infty }{e}^{-{r}^{2}}{r}^{2n+1}dr=\frac{1}{2}n!-{\int }_{0}^{R}{e}^{-{r}^{2}}{r}^{2n+1}dr⩾\frac{1}{2\pi }\left(1-\varepsilon {\right)}^{-2}n!$  for all sufficiently large $n$  , we obtain  4.5 .
From  4.5 we obtain for any $f={p}_{n}g\in {L}_{n}^{+}$  :
${\int }_{\mathbb{C}}|f\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)⩽{M}_{n}\left(v\right)\parallel g{\parallel }_{C\left(suppv\right)}^{2}⩽{M}_{n}\left(v\right)\left(1-\varepsilon {\right)}^{-2n}\frac{1}{n!}\parallel f{\parallel }_{{F}^{2}}^{2}.$  Together with  4.3 , the last estimate yields
 $\begin{array}{c}\left(n!{\mu }_{n+1}{\right)}^{1/n}⩽\left(1-\varepsilon {\right)}^{-2}{M}_{n}\left(v{\right)}^{1/n}\end{array}$ (4.6)
for all sufficiently large $n$  .
2. Lower bound for ${\mu }_{n+1}$  . Let us use formula  4.4 and take ${L}_{n}^{-}$  to be the set of all polynomials in $z$  of degree $⩽n$  . As in the proof of the upper bound, we denote ${R}_{0}={max}_{z\in suppv}|z|$  , fix $\varepsilon >0$  and set $R={R}_{0}/\varepsilon$  . We shall use the following equivalent norm in ${F}^{2}$  :
 $\begin{array}{c}|||f||{|}_{{F}^{2}}^{2}={\int }_{|z|⩾R}|f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}dm\left(z\right),|||f||{|}_{{F}^{2}}⩽\parallel f{\parallel }_{{F}^{2}}⩽C\left(R\right)|||f||{|}_{{F}^{2}}.\end{array}$ (4.7)
Let ${q}_{n}\in {L}_{n}^{-}\\left\{0\right\}$  be the polynomial which minimizes the ratio
 $\begin{array}{c}\frac{{\int }_{\mathbb{C}}|{q}_{n}\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)}{|||{q}_{n}||{|}_{{F}^{2}}^{2}}\end{array}$ (4.8)
among all polynomials in ${L}_{n}^{-}\\left\{0\right\}$  . The following standard argument shows that all zeros of ${q}_{n}$  are confined to the disk $\left\{z||z|⩽{R}_{0}\right\}$  . Suppose that one of the zeros ${z}_{k}$  is outside the disk; then replace ${q}_{n}\left(z\right)$  by ${q}_{n}\left(z\right)|{z}_{k}|\left(z-{R}_{0}^{2}/{\overline{z}}_{k}\right)/\left({R}_{0}\left(z-{z}_{k}\right)\right)$  . One has $\frac{|{z}_{k}||z-{R}_{0}^{2}/{\overline{z}}_{k}|}{{R}_{0}|z-{z}_{k}|}⩽1\text{for}|z|⩽{R}_{0}\text{and}\frac{|{z}_{k}||z-{R}_{0}^{2}/{\overline{z}}_{k}|}{{R}_{0}|z-{z}_{k}|}⩾1\text{for}|z|⩾{R}_{0}\text{,}$  so this change decreases the ratio  4.8 — contradiction. Next, without the loss of generality, we may assume that ${q}_{n}$  is monic. Denote $m=deg{q}_{n}⩽n$  ; we get the estimate $|||{q}_{n}||{|}_{{F}^{2}}^{2}={\int }_{|z|⩾R}|{q}_{n}\left(z\right){|}^{2}{e}^{-|z{|}^{2}}dm\left(z\right)⩽{\int }_{|z|⩾R}|z{|}^{2m}\left(1+\varepsilon {\right)}^{2m}{e}^{-|z{|}^{2}}dm\left(z\right)⩽\left(1+\varepsilon {\right)}^{2m}\pi m!.$  On the other hand, for the numerator of  4.8 , we have ${\int }_{\mathbb{C}}|{q}_{n}\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)⩾{e}^{-{R}_{0}^{2}}{\int }_{\mathbb{C}}|{q}_{n}\left(z\right){|}^{2}v\left(z\right)dm\left(z\right)⩾{e}^{-{R}_{0}^{2}}{M}_{m}\left(v\right).$  Combining the above estimates, we obtain:
 $\begin{array}{c}{\mu }_{n+1}⩾{inf}_{f\in {L}_{n}^{-}\\left\{0\right\}}\frac{{\int }_{\mathbb{C}}|f\left(z\right){|}^{2}v\left(z\right){e}^{-|z{|}^{2}}dm\left(z\right)}{C\left(R\right)|||f||{|}_{{F}^{2}}^{2}}⩾{min}_{0⩽m⩽n}\frac{{M}_{m}\left(v\right)}{{C}_{1}\left(R\right)\left(1+\varepsilon {\right)}^{2m}m!}.\end{array}$ (4.9)
As $z{p}_{m}\left(z\right)\in {\mathcal{P}}_{m+1}$  , from the definition  1.13 of ${M}_{m}\left(v\right)$  we get a trivial estimate ${M}_{m+1}\left(v\right)⩽{R}_{0}^{2}{M}_{m}\left(v\right)$  . This estimate shows that for a sufficiently large $n$  , the minimum in  4.9 is attained at $m=n$  . Therefore, $\left(n!{\mu }_{n+1}{\right)}^{1/n}⩾{\left(\frac{1}{{C}_{1}\left(R\right)}\right)}^{1/n}\frac{{M}_{n}\left(v{\right)}^{1/n}}{\left(1+\varepsilon {\right)}^{2}}⩾\left(1+\varepsilon {\right)}^{-3}{M}_{n}\left(v{\right)}^{1/n}$  for all sufficiently large $n$  . The latter estimate together with  4.6 completes the proof of the Lemma.

5 Proof of Lemma  3.2

First recall some well known facts concerning the spectral decomposition of the operator ${h}_{0}^{+}$  . The operator ${h}_{0}^{+}$  can be represented in terms of the annihilation and creation operators  1.5 as ${h}_{0}^{+}={a}_{0}^{*}{a}_{0}$  . The operators ${a}_{0}$  , ${a}_{0}^{*}$  obey the commutation relation $\left[{a}_{0},{a}_{0}^{*}\right]=2{B}_{0}$  , wherefrom we get the identity ${a}_{0}^{q}\left({a}_{0}^{*}{\right)}^{q}u=\left(2{B}_{0}{\right)}^{q}q!u$  for all $u\in {\mathcal{ℋ}}_{0}$  and $q\in \mathbb{N}$  . It follows that
 $\begin{array}{c}\left(2{B}_{0}{\right)}^{-q/2}\left(q!{\right)}^{-1/2}\left({a}_{0}^{*}{\right)}^{q}:{\mathcal{ℋ}}_{0}\to {\mathcal{ℋ}}_{q}\text{is an isometry onto}{\mathcal{ℋ}}_{q}\text{.}\end{array}$ (5.1)
Recalling the explicit isomorphism between ${\mathcal{ℋ}}_{0}$  and the space ${F}^{2}$  (see the previous section), we see that the change $u=\left(2{B}_{0}{\right)}^{-q/2}\left(q!{\right)}^{-1/2}\left({a}_{0}^{*}{\right)}^{q}\left({e}^{-{B}_{0}|z{|}^{2}/2}f\left(z\right)\right)$  gives a unitary equivalence between the operator ${P}_{q}v{P}_{q}$  and the operator in ${F}^{2}$  defined by the quadratic form
 $\begin{array}{c}\left(2{B}_{0}{\right)}^{-q}\left(q!{\right)}^{-1}{\int }_{\mathbb{C}}|\left({a}_{0}^{*}{\right)}^{q}\left({e}^{-{B}_{0}|z{|}^{2}/2}f\left(z\right)\right){|}^{2}v\left(z\right)dm\left(z\right).\end{array}$ (5.2)
For simplicity we will consider the case $q=1$  ; the general case can be treated in a similar manner. Also, we will take ${B}_{0}=2$  ; the general case can be reduced to this one by a linear change of variables. With these simplifications, the form  5.2 becomes
 $\begin{array}{c}{\int }_{\mathbb{C}}|{f}^{\prime }\left(z\right)-\overline{z}f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}v\left(z\right)dm\left(z\right),f\in {F}^{2}.\end{array}$ (5.3)
Let us prove the asymptotics  3.3 for the eigenvalues $\left\{{s}_{n}^{\left(1\right)}{\right\}}_{n=1}^{\infty }$  corresponding to the form  5.3 .
1. Upper bound for ${s}_{n}^{\left(1\right)}$  . Let ${\Omega }_{\delta }=\left\{z|dist\left(z,\Omega \right)⩽\delta \right\}$  . By the Cauchy integral formula, we have ${sup}_{\Omega }|{f}^{\prime }|⩽\frac{1}{\delta }{sup}_{{\Omega }_{\delta }}|f|,f\in {F}^{2}.$  Thus, we have the following bound for the form  5.3 :
${\int }_{\mathbb{C}}|{f}^{\prime }\left(z\right)-\overline{z}f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}v\left(z\right)dm\left(z\right)⩽C\parallel f{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2},$  where $C$  depends on $\Omega$  , $\delta$  , $v$  . Let us define ${L}_{n}^{+}=\left\{f\in {F}^{2}|f\left(z\right)={t}_{n}\left(z\right)g\left(z\right),\text{}g\text{entire}\right\},$  where ${t}_{n}$  is the $n$  'th Chebyshev polynomial for the set ${\Omega }_{\delta }$  . Note that the proof of  4.5 uses only the fact that all zeros of ${p}_{n}$  lie in the closed convex hull of $suppv$  . Therefore, the same estimate remains valid with the change ${p}_{n}↦{t}_{n}$  . Thus, we get
${\int }_{\mathbb{C}}|{f}^{\prime }\left(z\right)-\overline{z}f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}v\left(z\right)dm\left(z\right)⩽C\parallel f{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2}⩽C\parallel {t}_{n}{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2}\parallel g{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2}⩽C\left(1-\varepsilon {\right)}^{-2n}\frac{1}{n!}\parallel {t}_{n}{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2}\parallel f{\parallel }_{{F}^{2}}^{2},f\in {L}_{n}^{+},$
for all sufficiently large $n$  . This yields ${limsup}_{n\to \infty }\left(n!{s}_{n}^{\left(1\right)}{\right)}^{1/n}⩽\left(1-\varepsilon {\right)}^{-2}{lim}_{n\to \infty }\parallel {t}_{n}{\parallel }_{C\left({\Omega }_{\delta }\right)}^{2/n}=\left(1-\varepsilon {\right)}^{-2}\left(Cap{\Omega }_{\delta }{\right)}^{2}.$  It remains to note that $\varepsilon$  and $\delta$  can be chosen arbitrary small and that ${lim}_{\delta \to +0}Cap{\Omega }_{\delta }=Cap\Omega$  for any compact set $\Omega$  (see e.g. [7).
2. Lower bound for ${s}_{n}^{\left(1\right)}$  . Due to the compactness of the embedding of the Sobolev space ${W}_{2}^{1}\left(\Omega \right)\subset {L}^{2}\left(\Omega \right)$  , for any $\gamma >0$  there exists a subspace $N\subset {W}_{2}^{1}\left(\Omega \right)$  of a finite codimension such that
 $\begin{array}{c}\parallel u{\parallel }_{{L}^{2}\left(\Omega \right)}⩽\gamma \parallel \nabla u{\parallel }_{{L}^{2}\left(\Omega \right)},\forall u\in N.\end{array}$ (5.4)
Let us choose $\gamma =1/\left(4{R}_{0}\right)$  , ${R}_{0}={max}_{\Omega }|z|$  , consider the corresponding subspace $N$  and denote $codimN=l<\infty$  . Next, let ${L}_{n}^{-}$  , as above, be the set of all polynomials in $z$  of degree $⩽n$  . Consider the subspace ${\stackrel{~}{L}}_{n}^{-}={L}_{n}^{-}\cap N$  ; clearly, $dim{\stackrel{~}{L}}_{n}^{-}⩾n+1-l$  . By  5.4 , for any $f\in {\stackrel{~}{L}}_{n}^{-}$  we have $\parallel f{\parallel }_{{L}^{2}\left(\Omega \right)}⩽\frac{1}{2{R}_{0}}\parallel {f}^{\prime }{\parallel }_{{L}^{2}\left(\Omega \right)}$  . It follows that for any $f\in {\stackrel{~}{L}}_{n}^{-}$  :
$\parallel {f}^{\prime }{\parallel }_{{L}^{2}\left(\Omega \right)}⩽\parallel {f}^{\prime }-\overline{z}f{\parallel }_{{L}^{2}\left(\Omega \right)}+\parallel \overline{z}f{\parallel }_{{L}^{2}\left(\Omega \right)}⩽\parallel {f}^{\prime }-\overline{z}f{\parallel }_{{L}^{2}\left(\Omega \right)}+\frac{1}{2}\parallel {f}^{\prime }{\parallel }_{{L}^{2}\left(\Omega \right)},$  and so $\parallel {f}^{\prime }-\overline{z}f{\parallel }_{{L}^{2}\left(\Omega \right)}⩾\frac{1}{2}\parallel {f}^{\prime }{\parallel }_{{L}^{2}\left(\Omega \right)}⩾{R}_{0}\parallel f{\parallel }_{{L}^{2}\left(\Omega \right)}.$  Thus, for the quadratic form  5.3 we have ${\int }_{\mathbb{C}}|{f}^{\prime }\left(z\right)-\overline{z}f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}v\left(z\right)dm\left(z\right)⩾C{\int }_{\Omega }|f\left(z\right){|}^{2}dm\left(z\right),f\in {\stackrel{~}{L}}_{n}^{-}.$  According to the second inequality in  4.9 with $v={\chi }_{\Omega }$  (we denote by ${\chi }_{\Omega }$  the characteristic function of $\Omega$  ) we have therefore $\frac{{\int }_{\mathbb{C}}|{f}^{\prime }\left(z\right)-\overline{z}f\left(z\right){|}^{2}{e}^{-|z{|}^{2}}v\left(z\right)dm\left(z\right)}{\parallel f{\parallel }_{{F}^{2}}^{2}}⩾C\frac{{\int }_{\Omega }|f\left(z\right){|}^{2}dm\left(z\right)}{\parallel f{\parallel }_{{F}^{2}}^{2}}⩾{C}^{\prime }\frac{{M}_{n}\left({\chi }_{\Omega }\right)}{\left(1+\varepsilon {\right)}^{2n}n!},f\in {\stackrel{~}{L}}_{n}^{-}.$  It follows that ${s}_{n+1-l}^{\left(1\right)}⩾{C}^{\prime }\frac{{M}_{n}\left({\chi }_{\Omega }\right)}{\left(1+\varepsilon {\right)}^{2n}n!},$  for all sufficiently large $n$  . As stated in Remark  1.1 , ${lim}_{n\to \infty }{M}_{n}\left({\chi }_{\Omega }{\right)}^{1/n}=\left(Cap\Omega {\right)}^{2}$  . Thus, ${liminf}_{n\to \infty }\left(n!{s}_{n}^{\left(1\right)}{\right)}^{1/n}⩾\frac{\left(Cap\Omega {\right)}^{2}}{\left(1+\varepsilon {\right)}^{2}}$  for any $\varepsilon >0$  .
Acknowledgements We are indebted to M. Sh. Birman, Yu. Netrusov, G. Raikov, G. Rozenblum, A. Sobolev and H. Stahl for useful discussions. The work was supported by the Royal Society grant 2004/R1-FS. The authors are grateful to the Mathematisches Forschungsinstitut Oberwolfach for hospitality and financial support. The first named author is also grateful to Loughborough University for hospitality.
References

1. Y. Aharonov, A. Casher, Ground state of a spin $\frac{1}{2}$  charged particle in a two-dimensional magnetic field. Phys. Rev. A (3) 19 (1979), no. 6, 2461–2462.
2. J. Avron, I. Herbst and B. Simon, Schrödinger operators with magnetic fields. I. General interactions, Duke Math. J. 45 (1978), no. 4, 847–883.
3. A. Besch, Eigenvalues in spectral gaps of the two-dimensional Pauli operator. J. Math. Phys. 41, 7918–7931.
4. M. Sh. Birman and M. Z. Solomyak, Spectral theory of self-adjoint operators in Hilbert space, Dordrecht, D.Reidel P.C., 1987
5. B. A. Dubrovin, S. P. Novikov, Fundamental states in a periodic field. Magnetic Bloch functions and vector bundles. (Russian) Dokl. Akad. Nauk SSSR 253 (1980), no. 6, 1293–1297.
6. L. Erdös, V. Vougalter, Pauli operator and Aharonov-Casher theorem for measure valued magnetic fields, Comm. Math. Phys. 225 (2002), no. 2, 399–421.
7. E. Hille, Analytic function theory. Vol. II, Ginn and Co., Boston, Mass., 1962
8. K. Hornberger and U. Smilansky, Magnetic edge states, Physics Reports 367 (2002), 249–385.
9. A. Iwatsuka The essential spectrum of two-dimensional Schrödinger operators with perturbed constant magnetic fields, J. Math. Kyoto Univ. 23 no. 3 (1983), 475–480.
10. V. Ivrii, Microlocal analysis and precise spectral asymptotics, Springer, Berlin, 1998
11. N. S. Landkof, Foundations of modern potential theory, Springer, New York, 1972
12. M. Melgaard and G. Rozenblum, Eigenvalue asymptotics for weakly perturbed Dirac and Schrödinger operators with constant magnetic fields of full rank, Comm. Partial Differential Equations 28 (2003), no. 3-4, 697–736.
13. O. G. Parfënov, The widths of some classes of entire functions, Mat. Sb. 190 (1999), no. 4, 87–94; translation in Sb. Math. 190 (1999), no. 3-4, 561–568.
14. G. D. Raikov, Eigenvalue asymptotics for the Schrödinger operator with homogeneous magnetic potential and decreasing electric potential. I. Behaviour near the essential spectrum tips, Comm. Partial Differential Equations 15 (1990), no. 3, 407–434; Errata: Comm. Partial Differential Equations 18 (1993), no. 11, 1977–1979.
15. G. D. Raikov, Border-line eigenvalue asymptotics for the Schrödinger operator with electromagnetic potential, Integral Equations Operator Theory 14 (1991), no. 6, 875–888.
16. G. D. Raikov, Spectral asymptotics for the perturbed 2D Pauli operator with oscillating magnetic fields. I. Non-zero mean value of the magnetic field, Markov Processes Relat. Fields 9, 775-794 (2003).
17. G. D. Raikov and S. Warzel, Quasi-classical versus non-classical spectral asymptotics for magnetic Schrödinger operators with decreasing electric potentials, Rev. Math. Phys. 14 (2002), no. 10, 1051–1072.
18. G. Rozenblum and N. Shirokov, Infiniteness of zero modes for the Pauli operator with singular magnetic field, preprint 2005, http://lanl.arxiv.org/abs/math-ph/0501059
19. A. V. Sobolev, Asymptotic behavior of energy levels of a quantum particle in a homogeneous magnetic field perturbed by an attenuating electric field. I. (Russian) Probl. Mat. Anal., 9, 67–84, Leningrad. Univ., Leningrad, 1984. English translation in: J. Sov. Math. 35 (1986), 2201–2212.
20. A. V. Sobolev, On the Lieb-Thirring estimates for the Pauli operator, Duke Math. J. 82 (1996), no. 3, 607–635.
21. S. N. Solnyshkin, Asymptotic behavior of the energy of bound states of the Schrödinger operator in the presence of electric and homogeneous magnetic fields. (Russian) Probl. Mat. Fiz., 10, 266–278, Leningrad. Univ., Leningrad, 1982.
22. H. Stahl and V. Totik, General orthogonal polynomials, Cambridge Univ. Press, Cambridge, 1992.
23. H. Tamura, Asymptotic distribution of eigenvalues for Schrödinger operators with homogeneous magnetic fields, Osaka J. Math. 25 (1988), no. 3, 633–647.