statistical mechanics is a bridge between mechanics and thermodynamics .
indeed , it may be assumed as the microscopic mechanical basis for the macroscopic thermodynamical properties of a physical system .
the key concept which links these two arenas is the notion of the phase space , that is the space of all distinct possible microstates .
the dynamical evolution is nothing but a smooth transition from one microstate to another and statistical mechanics is about how to count the number of these microstates . in a more technical language
, the phase space of a dynamical system is a symplectic manifold which naturally is equipped with a symplectic structure @xcite .
the symplectic structure covariantly determines the poisson brackets and liouville volume and the kinematics is then defined in a coordinate - independent manner on the phase space . the dynamics will be also specified upon taking the hamiltonian function as the generator of time evolution in this setup .
therefore , the hamiltonian determines how the system evolves through the microstates and the liouville volume specifies how one can count and measure them .
consider , for instance , a mechanical system consisting of @xmath0 particles which are subjected to some forces .
each particle moves in euclidean three - dimensional space @xmath1 and the corresponding configuration space is then @xmath2 .
therefore , the phase space will be the cotangent bundle of this space , that is , @xmath3 , which naturally admits a symplectic structure . in terms of the positions and momenta of the particles , the associated symplectic structure takes the standard canonical form .
such a canonical symplectic structure , which is also dynamically invariant , may be used to determine the fundamental canonical poisson brackets and the kinematics for the system under consideration .
subsequently , based on the canonical symplectic structure one leads to the canonical measure which can be viewed as a mathematical expression for the principle of equal probabilities in standard statistical mechanics @xcite .
although it might be thought that the symplectic formulation is nothing but a covariant reformulation of the standard statistical mechanics , its advantages may become more clear when one is dealing with the hamiltonian systems with nontrivial structures ( for instance when the standard @xmath4 topology of the phase space is replaced with a more complicated topology or the standard canonical poisson algebra is replaced with a deformed noncanonical one ) .
such systems , are investigated in the context of doubly special relativity theories , as a flat limit of quantum gravity , where the four - momentum space has a curved ( de sitter or anti - de sitter ) geometry with nontrivial topology @xcite .
the space - time structure then turns out to be noncommutative and the associated poisson algebra between the positions and momenta become noncanonical @xcite .
the space - time with noncommutative coordinates was first proposed by snyder in ref .
@xcite in order to regularize quantum field theory in the ultraviolet regime @xcite . in more recent times , it was also used by the stability theory of the lie algebras to be a minimal candidate for the relativistic quantum mechanics @xcite and by string theory at a fundamental level @xcite .
the duality between curved momentum space and noncommutative structure was first pointed out by majid in ref .
the direct consequence of space - time with noncommutative coordinates is the modification to the dispersion relation and the associated density of states @xcite . in symplectic formulation ,
the density of states is properly determined by the liouville measure which will be constructed from the symplectic structure @xcite . in noncommutative space - time , however , the corresponding symplectic structure seems to be noncanonical which induces a nonuniform probability distribution over the set of microstates through a deformed liouville volume element .
this dynamically invariant liouville volume enters in the definition of the gibbs entropy and partition function from which all the thermodynamical properties of a system can be extracted .
therefore , the topology of the phase space will play an important roll since all integrals are evaluated over the entire phase space .
for instance , let us consider a system which has an ultraviolet cutoff due to the existence of a maximal momentum .
it is easy to see that the topology of the spatial sector of the momentum space may be compact .
moreover , systems with a totally compact phase space , which can be considered as the classical limit of quantum systems with finite dimensional hilbert space ( e.g , angular momentum for systems with fixed total angular momentum ) were recently studied in the context of loop quantum gravity @xcite .
for such a compact phase space , the liouville volume is finite and consequently a finite number of microstates naturally arises @xcite .
this shows the correspondence between the topology of the phase space and the number of microstates which can be realized when the symplectic geometry is implemented . the black hole physics is a good example in which the system obeys the thermodynamical laws and its statistical origin is not thoroughly formulated yet @xcite in the framework of hamiltonian formalism ( see also ref .
@xcite ) .    in this
regard , formulating the statistical mechanics in its most fundamental form is important , at least , for two reasons : ( i ) finding a more precise and fundamental interpretation for the basis of the statistical mechanics , and ( ii ) studying the statistical mechanics for the hamiltonian systems with nontrivial structure , such as those we have mentioned above .
motivated by the stated issues , in this paper we are going to formulate statistical mechanics for a general hamiltonian system in a covariant ( coordinate - independent ) manner .
the paper is organized as follows . in section slowromancap2@
, we deal with the kinematics and dynamics of a many - particle system with general hamiltonian structure in the context of symplectic geometry . the statistical mechanics for such systems
is formulated in section slowromancap3@. in section slowromancap4@ , by means of the constructed setup we study the statistical mechanics in the snyder noncommutative space - time with curved energy - momentum space as an explicit example of a physical system with nontrivial ( noncommutative ) hamiltonian structure . as a case study ,
the thermodynamical properties of a harmonic oscillator system in snyder space is also presented at the end of this section .
section slowromancap5@ is devoted to the summary and conclusions .
in this section , in order to study the statistical mechanics in the context of symplectic geometry , we consider the kinematics and dynamics of an @xmath0-particle system .
let us consider a system of @xmath0 particles similar to conventional systems in statistical mechanics .
the corresponding @xmath5-dimensional kinematical phase space may be obtained by coupling @xmath0 single - particle phase spaces @xmath6 with @xmath7 as @xmath8 the phase space @xmath9 naturally admits a symplectic structure @xmath10 which is a nondegenerate closed 2-form @xcite .
the classical state of this @xmath0-particle system is determined by a point in @xmath9 ( phase point ) .
the evolution of phase points is then determined by the hamiltonian vector field @xmath11 , whose integral curves are trajectories of the phase points ( phase trajectories ) . from the fact that @xmath10 is nondegenerate , one can assign a vector field to a function @xmath12 as @xmath13 and the poisson bracket between two observables ( real - valued functions on the phase space ) can be defined as @xmath14 these functions with the above structure constitute a lie algebra under poisson brackets on @xmath9 .
we also have the so - called liouville volume on @xmath9 with its standard definition as @xmath15 in terms of positions and momenta @xmath16 of the particles , the liouville volume takes the local form @xmath17 where @xmath18 and @xmath19 are the @xmath20th component of the positions and momenta of the @xmath21th particle , respectively , and @xmath22 is the matrix representation of the symplectic structure @xmath10 with components @xmath23 .
now , one may be able to decompose the symplectic structure as @xmath24 where @xmath25 is the symplectic structure on @xmath26 .
although this assumption has no mathematical basis , from a physical point of view it is acceptable since the particles are assumed to be kinematically separated . in other words ,
the lie algebras of the functions defined by @xmath27 over @xmath28 are separately closed .
therefore , the liouville volume over @xmath9 takes the form @xmath29 where the wedge products take place @xmath30 times .
the volume element then works out to be @xmath31 where @xmath32 is the volume element of the corresponding @xmath33 , and we have also used the fact that @xmath34 , for @xmath35 . clearly , @xmath36 , that is , a function on @xmath9 , factorizes into products of @xmath37 .
this can not be realized in the most general situation , where the symplectic structure can not be written as ( [ decomposition ] ) .    as a special case ,
consider a system consisting of @xmath0 particles subject to some forces moving in euclidean three - dimensional space with standard @xmath38 topology .
therefore , the phase space will be @xmath39 which is indeed the cotangent bundle of the configuration space , that is , @xmath40 .
geometrically , the cotangent bundle is a symplectic manifold endowed with a canonical symplectic structure @xmath41 in which @xmath42 and @xmath43 are interpreted as the positions and momenta of particles , respectively .
the summation is over both @xmath44 and @xmath45 .
this is the canonical representation of symplectic structure and coordinates @xmath46 are known as the canonical coordinates . substituting the canonical structure ( [ canonical structure ] ) into the definition ( [ pbd1 ] )
, the poisson bracket will be @xmath47 which leads to the standard canonical poisson algebra @xmath48 also , with the help of ( [ canonical structure ] ) we get from ( [ vol - d ] ) the liouville volume as @xmath49 which is nothing but the standard volume element on @xmath50 .
thus , the kinematics of the standard statistical mechanics can be recovered as a special case of this setup .
taking into account the fact that @xmath51 , the measure associated to the standard volume element ( [ canonical measure ] ) assigns a uniform probability distribution to the set of microstates .
this is justifiable in the light of laplace s principle of indifference , which states that in the absence of any further information , all outcomes are equally likely @xcite .
this is the fundamental basis of the statistical mechanics .    in the case where the space - time has a noncommutative structure ,
apart from the details of the different models , there is always a deformation parameter which can be linked to a minimal length associated with the system under consideration .
the direct consequence of this setup is the modification of the dispersion relation and the density of states ( see ref .
@xcite for the special case of the stable noncommutative algebra for the relativistic statistical mechanics ) .
the density of states is determined by the symplectic structure and the associated liouville volume .
since the symplectic structure takes a noncanonical form in terms of the physical positions and momenta @xmath52 in noncommutative spaces @xcite , we have @xmath53 .
this result immediately shows that laplace s indifference principle is no longer applicable for a noncommutative space - time .
this is because of the fact that there is extra information in these setups ( minimal length or maximal momentum ) and the liouville measure then assigns a nonuniform probability distribution over the set of microstates by means of the relation ( [ vol - d2 ] ) . this simple consideration of noncommutative space - time as an example of a kinematically deformed system shows the advantage of the coordinate - independent symplectic formulation of the statistical mechanics ( see refs .
@xcite for explicit examples ) .
we will explicitly consider such an example in section iv .
it also should be noted that , according to the darboux theorem , there is always a local chart on @xmath9 in which any symplectic structure takes the canonical form with @xmath54 . in the case of kinematically deformed systems such as the noncommutative space - time , however , the new canonical coordinates can not be interpreted as the positions and momenta of the particles .
moreover , in the case of standard statistical mechanics , one can also find a local chart in which the standard canonical structure @xmath55 takes a noncanonical form with @xmath56 .
but again , these new noncanonical coordinates can not be interpreted as positions and momenta of particles .
in all of the above cases , as we will see in the next section , the partition function and consequently resultant thermodynamical properties are independent of the chart in which the system is considered .
as we have mentioned above , the time evolution of the system is determined by the hamiltonian vector field @xmath11 .
it can be obtained by solving the dynamical equation @xmath57 where @xmath58 is the hamiltonian function of the system that is a real - valued function over @xmath9 and @xmath59 is an interior product of @xmath10 and @xmath11 .
additionally , @xmath11 can be expanded as @xmath60 where @xmath61 is the projection of @xmath62 on @xmath26 .
note that , in general , the domain of components of @xmath61 is the whole phase space @xmath9 and consequently the lie brackets between @xmath63 will not vanish : @xmath64\neq0\,.\ ] ] however , as a special case , one can consider a system consisting of noninteracting particles .
therefore , the hamiltonian function can be written as the sum of individual hamiltonians , that is , @xmath65 so , from decomposition ( [ decomposition ] ) for the symplectic structure , the dynamical equation ( [ can - dy - eq ] ) becomes @xmath66 using the fact that @xmath67 , for @xmath68 , one is led to a set of @xmath0 independent differential equations @xmath69 now , it is clear that the domain of the components of @xmath70 will be @xmath26 and the lie brackets between @xmath71 and @xmath72 vanish @xmath73=0\,.\ ] ] this can be seen as a criterion for the _ kinematical and dynamical separability _ of the particles . in the appendix , we will see that this benchmark shows itself as a criterion for the statistical independence of the macroscopic subsystems which in turn is a definition of equilibrium .
the other important property of the hamiltonian systems is the existence of an invariant measure that enables one to construct the equilibrium statistical mechanics through the well - known liouville theorem as @xmath74 where @xmath75 is the total time derivative and @xmath76 denotes the lie derivative with respect to @xmath62 .
the conservation of @xmath77 can be traced back to the fact that @xmath78 and @xmath79 , where we have used @xmath80 and the closure of the symplectic structure @xmath81 .
it is obvious that @xmath77 is not explicitly time dependent .
note also that the liouville theorem arises from the particular dynamics of a hamiltonian system .
consider the standard statistical mechanics as a special case , where the kinematics is given by canonical symplectic structure ( [ canonical structure ] ) .
the dynamical equation ( [ can - dy - eq ] ) then leads to the familiar form of the hamilton s equations in terms of the physical positions and momenta as @xmath82 note that when the symplectic structure has a noncanonical form , the form of hamilton s equations would deviate from the relation ( [ h equation ] ) .
the triple @xmath83 constitutes a general hamiltonian system and we explore the statistical physics of this general system in the next section by means of an invariant volume element @xmath77 which is covariantly defined in the relation ( [ vol - d ] ) .
statistical mechanics links macroscopic properties of a system to microscopic laws that are quantum mechanical at the fundamental level .
the basic concept is the von neumann entropy that is defined as @xmath84 where @xmath85 is the density operator @xcite . using the principle of maximum entropy and considering specific relevant statistical constraints
, we can find the density operator which contains all the statistical information about the system at equilibrium .
for example , for a system in contact with a thermal bath at temperature @xmath86 and therefore with a fixed mean energy , the resultant normalized density operator works out to be @xmath87\,,\end{aligned}\ ] ] where @xmath88 is the hamiltonian operator of the system and @xmath89 is called the canonical partition function which is defined as @xmath90=\sum_{i}\exp\big[-\varepsilon_i / t \big]\,,\end{aligned}\ ] ] where @xmath91 are the energy eigenvalues of @xmath88 and the summation is taken over all the accessible microstates for the system . at the fundamental level ,
the microscopic laws are quantum mechanical and the microstates are determined by quantum mechanics . on the other hand , in the classical limit , the microstate of a system represents itself as a point on the corresponding phase space and due to the infinite resolution of phase space points , the summation over microstates is really not well defined .
consequently , one can not construct a full classical statistical mechanics in essence
. however , quantum mechanics provides the semiclassical approximation in the spirit of the uncertainty principle .
indeed , a phase space with finite resolution and well - defined number of classical microstates can be achieved through the approximation @xcite @xmath92 where @xmath93 is due to the fact that the particles may be indistinguishable and @xmath77 is the dynamically invariant liouville volume element of the @xmath5-dimensional symplectic manifold @xmath9 which is defined in the relation ( [ vol - d ] ) .
it is important to note that the semiclassical approximation ( [ approx ] ) coincides with the full quantum consideration at high temperature regime . in this limit
, the von neumann entropy ( [ von entropy ] ) is replaced with the gibbs entropy @xmath94 where @xmath95 is now interpreted as a probability density over the space of microstates @xmath9 .
the maximum entropy principle for a system in heat bath at temperature @xmath86 leads to the gibbs state @xmath96\,,\end{aligned}\ ] ] where @xmath58 is the hamiltonian function of the system and @xmath89 is the semiclassical canonical partition function @xmath97\omega^{3n}\,.\end{aligned}\ ] ] the result of the integral over the phase space for the total partition function ( [ partition f ] ) is independent of a chart in which the system is considered when the dynamical equation ( [ can - dy - eq ] ) is satisfied in a chart - independent manner .
while the local form of the symplectic structure determines the probability distribution over the set of microstates , from the statistical point of view , we expect that the resultant thermodynamical properties extracted from a partition function are independent of a chart in which the physical system is considered ( see , for instance , ref .
@xcite in which the partition function in two different charts is calculated ) .    by considering the decompositions ( [ decomposition ] ) and ( [ decomposition h ] ) for the symplectic structure @xmath10 and hamiltonian function @xmath58 , respectively , the probability density will also be decomposed as @xmath98 where @xmath99 is the gibbs state defined on @xmath6 .
hence the integral over @xmath9 can be factorized into the products of the integrals over the phase spaces of the particles as @xmath100\,\omega^3_{\alpha}\bigg)\,.\end{aligned}\ ] ] now , if the particles are influenced by the same kinematics and dynamics , the partition function can be written in the well - known form @xmath101 where we have defined the single - particle state partition function @xmath102\,\omega_1 ^ 3\,.\end{aligned}\ ] ] it is usually claimed that the relation ( [ partition fghi ] ) is applicable when one ignores the quantum correlations between particles .
this means that decompositions such as ( [ decomposition ] ) and ( [ decomposition h ] ) may be considered as conditions for the absence of quantum correlations .
moreover , the general formalism which is formulated in this section , can also support a hamiltonian system with very nontrivial structure , _
e.g. _ a kinematically classical entangled system in which the decomposition ( [ decomposition ] ) is no longer applicable .
existence of a minimal length is suggested by any quantum gravity candidate such as string theory and loop quantum gravity @xcite .
it is then widely believed that a nongravitational theory which supports the existence of a minimal length would arise at the flat limit of quantum gravity .
such a theory , would be reduced to the standard relativistic quantum mechanics at low energy regime where the effects of a minimal length are negligible . evidently , taking a minimal length into account
naturally leads to the space - time with noncommutative coordinates . the first attempt in this direction
was done by snyder in 1947 who formulated a lorentz - invariant discrete space - time with noncommutative coordinates @xcite .
the space - time with noncommutative structure is also suggested by the theory of stability of the lie algebras @xcite . on the other hand , existence of a universal minimal length
can not be supported by the special relativity since any length scale in one inertial frame may be different in another observer s frame through the well - known lorentz - fitzgerald contraction .
thus , the doubly special relativity theories are formulated in order to take into account an observer - independent minimum length as well as the velocity of light @xcite .
the curved four - momentum space then naturally appears to be a suitable framework to formulate the doubly special relativity theories @xcite .
furthermore , the different doubly special relativity theories can be obtained from the different basis of the @xmath103-poincar algebra on the associated @xmath103-minkowski noncommutative space - time @xcite . in this respect
, it is also shown that the snyder noncommutative algebra which is proposed in ref . @xcite and the stable algebra of relativistic quantum mechanics that is obtained in ref .
@xcite can be obtained from a ten - dimensional phase space with curved geometry for the four - momentum space through the symplectic reduction process @xcite .
thus , it seems that the space - time with curved four - momentum space and noncommutative coordinates is a fundamental framework for taking into account a minimal length scale in the flat limit of quantum gravity @xcite .
it would also be claimed that the lorentz invariance is an approximate symmetry and will be broken at high energy regime .
therefore , the lorentz - violating noncommutative algebra was proposed by camelia in ref .
@xcite which supports the existence of a minimal observer - independent length scale .
magueijo and smolin showed that the lorentz invariance can also be preserved by a nonlinear action of the lorentz group on the momentum space @xcite . in order to preserve the lorentz invariance ,
the following commutation relations for the generators will be held @xmath104 where @xmath105 is the minkowski metric and @xmath106 .
the commutation relations between @xmath107 then classify the different lorentz - invariant algebras such as the snyder noncommutative algebra @xcite , the stable algebra @xcite and also the standard lorentz algebra .
while the relativistic algebras are written in the eight - dimensional relativistic phase space , to formulate the statistical mechanics one needs to work in a six - dimensional nonrelativistic subalgebra of the deformed relativistic algebras in which the time parameter is fixed and the corresponding poisson brackets take a noncanonical form . by means of the resulting deformed hamiltonian system , one can study the statistical mechanics with the help of the constructed setup in the pervious sections . in recent years , many works have gone in this direction in order to study the effects of an ultraviolet cutoff ( minimal length and maximal momentum ) on the thermodynamical properties of the physical systems .
for instance , thermodynamics of some physical systems in noncommutative spaces is studied in ref .
for the special case of the snyder noncommutative space , see ref .
@xcite , in which the thermodynamics of the early universe is explored .
it is shown that the liouville theorem for the deformed phase space with general noncanonical poisson algebra is satisfied . however , as we have shown in section slowromancap3@ , the liouville theorem can be justified in a very simple manner in the covariant formalism , which also shows the advantage of the setup .
furthermore , thermodynamical ideal gases are studied in the context of doubly special relativity theories @xcite .
also , motivated by the string theory , the generalized uncertainty principle is suggested which supports the existence of a minimal length as a nonzero uncertainty in position measurement @xcite . for the statistical mechanics in this setup ,
inspired by loop quantum gravity , the polymer quantum mechanics is formulated which supports the existence of a minimal length known as the polymer length scale @xcite .
thermodynamical properties of the ideal gases and harmonic oscillator are also studied in refs .
@xcite . to see the effects of minimal length on the thermodynamics of the black holes in noncommutative space , the generalized uncertainty principle framework , and the polymer quantization scheme ,
@xcite . apart from the details of the above - mentioned effective approaches to quantum gravity phenomenology , all of them try to consider the effects of a minimal length in a relevant manner , and implementing the covariant formalism can clarify the consequences and applications in its most fundamental form in the language of symplectic geometry .
as we have mentioned above , the advantages of the covariant formulation may be revealed when a hamiltonian system with nontrivial structure is considered .
therefore , in this subsection we study the statistical mechanics in the snyder noncommutative space as a well - known example of a kinematically deformed hamiltonian system .
the snyder relativistic algebra preserves the lorentz invariance and therefore the commutation relations ( [ algebra ] ) are satisfied by the corresponding generators .
the commutation relations between @xmath107 which define the snyder algebra are given by @xcite @xmath108 where @xmath103 is the deformation parameter with the dimension of the inverse of length which plays the role of the universal quantum gravity scale in this setup .
the four - momentum space of the eight - dimensional relativistic phase space of a test particle which moves on the space - time with noncommutative algebra ( [ snyder - algebra ] ) , is a de sitter space with topology @xmath109 @xcite . identifying the energy space with @xmath110 , the topology for the space of spatial momenta will be @xmath111 . in order to study the statistical mechanics in this setup
, one should consider the nonrelativistic subalgebra of ( [ snyder - algebra ] ) and replace the standard canonical poisson algebra with them .
the nonrelativistic snyder algebra for a particle moves in three - dimensional euclidean space @xmath38 and is then given by @xcite @xmath112 it is straightforward to show that the above commutation relations can be realized from the symplectic structure @xcite @xmath113 through the covariant definition ( [ pbd1 ] ) for the poisson brackets .
the associated phase space of the particle is then simply the cotangent bundle of the configuration space which now has the nontrivial @xmath114 topology .
indeed , the space of the spatial momenta has compact @xmath115 topology rather than the standard @xmath38 one ( see refs .
@xcite for more details ) .
since the quantum gravity cutoff ( which is defined by the deformation parameter @xmath103 ) is universal , it will be the same for all of the particles .
thus , for the kinematically deformed system consisting of @xmath0 particles which obeys the snyder noncommutative algebra ( [ snyder - algebra - nr ] ) , the total symplectic structure can be obtained by the relation ( [ decomposition ] ) .
the phase space @xmath6 for all of the particles is the same and has the nontrivial @xmath116 topology .
the associated total phase space @xmath117 can be easily obtained through the coupling ( [ gamma split n ] ) , which has the following nontrivial topology : @xmath118 substituting the symplectic structure ( [ snyder - symplectic ] ) into the definition ( [ vol - d ] ) , the corresponding liouville volume takes the form @xmath119 where @xmath18 and @xmath19 are the @xmath20th component of the positions and momenta of the @xmath21th particle respectively . from this relation
it is clear that the probability distribution is nonuniform in snyder noncommutative space and the standard uniform one can be recovered only at low energy ( low temperature ) regime @xmath120 .
indeed , the existence of a minimal length , as an extra information , changes the probability distribution at the high energy regime .
thus , in all of the kinematically deformed noncommutative phase spaces , the local form of the symplectic structure ( in terms of the positions and momenta of the particles ) is noncanonical and therefore the probability distribution will be nonuniform such that the microstates with higher momenta are less probable .
this is the advantage of the symplectic covariant formulation of the statistical mechanics which shows that the probability distribution will be nonuniform on any phase space with noncanonical poisson algebra like the snyder algebra ( [ snyder - algebra - nr ] ) .
however , in light of the darboux theorem , one can always find a local chart through the darboux transformation @xmath121 , in which the symplectic structure ( [ snyder - symplectic ] ) takes the canonical form @xmath122 for all of the particles .
thus , the associated liouville volume becomes @xmath123 , which induces a uniform probability distribution over the set of microstates .
but , it is important to note that @xmath124 are different from the standard canonical ones @xmath125 which are defined in relations ( [ canonical structure ] ) and ( [ canonical measure ] ) . indeed , while @xmath125 are the positions and momenta of the particles by definition , @xmath126 can not be interpreted as the positions and momenta of the particles .
also , the noncanonical snyder liouville volume ( [ snyder - vol ] ) reduces to the standard canonical one ( [ canonical measure ] ) in the limit of @xmath127 and the noncanonical variables @xmath128 then coincide with standard canonical variables @xmath125 in this limit .
however , the canonical variables @xmath126 are obtained through a darboux transformation and the functional form of the hamiltonian function will be modified ( compared with the standard functional form ) in this chart as @xmath129 since the dynamical equation ( [ can - dy - eq ] ) will be satisfied in a chart - independent manner .    substituting the liouville volume ( [ snyder - vol ] ) in the relation ( [ partition f ] ) , the total partition function will be @xmath130}{\prod_{\alpha=1}^{n}\big(1+(p^\alpha/\kappa)^2 \big)}d^{3n}qd^{3n}p\,.\end{aligned}\ ] ] since the snyder liouville volume ( [ snyder - vol ] ) is decomposed similar to the relation ( [ decomposition 4 ] ) , for the noninteracting systems in which the total hamiltonian function also decomposes as the relation ( [ decomposition h ] ) , one can rewrite the total partition function ( [ snyder - tpf ] ) in terms of the single - particle partition function . taking the decomposition ( [ decomposition h ] ) into account and substituting the volume ( [ snyder - vol ] ) into relation ( [ snyder - tpf ] )
, the total partition function can be rewritten in the form of the relation ( [ partition fgh ] ) as @xmath131}{\big[1+(p^\alpha /\kappa)^2\big]}\,d^3q_{\alpha}d^3p^\alpha\bigg)\nonumber
\\=\frac{1}{n ! } { \mathcal z_1}^n,\hspace{6cm}\end{aligned}\ ] ] where we have defined the single - particle partition function as @xmath132}{\big(1+(p^1/\kappa)^2\big)}\,,\ ] ] and we have also assumed that the hamiltonian functions are the same for all of the particles .
then , all the thermodynamical properties of a system can be deduced from the total partition function ( [ snyder - tpf ] ) or ( [ snyder - d - tpf ] ) in which , as we have mentioned in the pervious section , the gibbs factor @xmath133 is considered for the nonlocalized systems such as the ideal gas and it should then be removed for the localized systems such as the harmonic oscillator which is the subject of the next subsection .
now let us consider a system of @xmath0 independent and similar three - dimensional harmonic oscillators .
since the oscillators are assumed to be independent , the relations ( [ decomposition ] ) and ( [ decomposition h ] ) will be satisfied and one then can implement the relation ( [ snyder - d - tpf ] ) in order to obtain the partition function of the system .
the hamiltonian of the three - dimensional simple harmonic oscillator is given by @xmath134 , where @xmath135 is the mass of the oscillator , @xmath136 is the frequency , @xmath137 and @xmath138 . substituting this hamiltonian into the relation ( [ snyder - pf1 ] ) , gives the following single - particle partition function @xmath139\exp\bigg[\frac{\kappa^2}{2mt}\bigg]}\right)\,.\end{aligned}\ ] ] at the high temperature regime , the second term in the right - hand side of the above relation is negligible . under this condition
the dominant contribution is due to the first term and the partition function is approximated by @xmath140 .
this shows that in the domain of the high temperature , the number of degrees of freedom will be reduced since the partition function is proportional to @xmath141 rather than the standard one which is proportional to @xmath142 .
we will show this feature in a more precise manner in this subsection .
the total partition function can be obtained from the relation ( [ snyder - d - tpf ] ) as @xmath143 where the gibbs factor @xmath133 is dropped since the oscillators are localized . substituting the single - particle partition function ( [ ho - pf ] ) into this relation , the total partition function for @xmath0 three - dimensional harmonic oscillators can be obtained from which one is able to study the thermodynamical behavior of the system in snyder space . however , before doing this task , it is useful to consider the limit @xmath144 which is related to the first order corrections that would arise from the minimal length effects .
this limit also allows us to compare our result with the full quantum partition function . in this limit
, we may use the relation @xmath145 , to arrive at the following asymptotic expression for the partition function @xmath146 in which we have defined the dimensionless variable @xmath147 . in order to obtain the quantum partition function
, one should solve the hamiltonian eigenvalue problem in snyder space which is not an easy task at all since the schrdinger equation takes a complicated form in this setup .
however , in the context of the generalized uncertainty principle the problem of a @xmath148-dimensional harmonic is exactly solved @xcite and the results are applicable also for the snyder algebra ( [ snyder - algebra - nr ] ) as a particular case .
the energy eigenvalues for a three - dimensional harmonic oscillator in snyder space are given by @xmath149,\hspace{0.5cm}\end{aligned}\ ] ] where @xmath150 , and @xmath151 and also we have considered relevant parameter identification to rewrite the result of the ref .
@xcite in our notation and units .
the associated quantum partition function then will be obtained from the definition ( [ qpf ] ) as @xmath152\,.\ ] ] the above summation can not be evaluated analytically . however , in ref .
@xcite , it is shown that up to the first order of approximation , in the limit of @xmath153 , it leads to @xmath154 $ ] , which is nothing other than the first term in the summation ( [ ho - pf - asym ] ) .
this coincidence shows that the quantum partition function ( [ ho - qpf ] ) correctly leads to the semiclassical partition function ( [ ho - pf ] ) up to the first order of approximation at the high temperature regime .
it also reveals the advantages of semiclassical approximation ( [ approx ] ) in the sense that although there is not an analytical solution for the case of full quantum partition function ( [ ho - qpf ] ) , the relation ( [ ho - pf ] ) provides an analytical expression for the high temperatures .
therefore , in dealing with the statistical considerations of the physical systems in the presence of a minimal length , since one is usually interested in the high temperature regime ( such as early universe thermodynamics ) , the usage of the semiclassical approximation seems to be quite reasonable .    substituting the single - particle partition function ( [ ho - pf ] ) into the relation ( [ ho - tpf ] ) gives the total partition function .
the internal energy then can be obtained from the standard definition @xmath155 as @xmath156\,{e}^{\theta^2}\right)^ { -1}-\theta^2\right].\hspace{.5cm}\end{aligned}\ ] ] in figure [ fig:1 ] we have plotted the internal energy versus temperature in comparison with its nondeformed counterpart .
the specific heat can also be deduced from the definition @xmath157 which leads to a somehow cumbersome expression .
+ 3em internal energy versus temperature is plotted .
the solid and dashed lines represent the internal energy for the snyder - deformed and nondeformed cases , respectively .
the deviation from the standard case arises at the high temperature regime when the minimal length effects dominate .
the figure is plotted for @xmath158.,title="fig:",width=288 ]    an interesting feature here is the reduction of the number of degrees of freedom in snyder space . according to the well - known equipartition theorem of energy , for the hamiltonians of the form @xmath159 , each of the 6 degrees of freedom makes a contribution of @xmath160 towards the internal energy .
thus one may consider @xmath161 as the number of degrees of freedom for each three - dimensional harmonic oscillator . in figure [ fig:2 ] the number of degrees of freedom for each three - dimensional harmonic oscillator is plotted . as
this figure shows the number of degrees of freedom will be reduced from 6 to 4 in snyder space .
this may be compared with the statistical mechanics of the ideal gases considered in ref .
@xcite in which the same result is obtained .
however , we would like to emphasize that the reduction of the number of degrees of freedom for the case of ideal gas can be interpreted as an effective dimensional reduction of the space at the high temperature regime which is a common feature of quantum gravity proposals @xcite and also phenomenological approaches to the minimal length conjecture @xcite .
+ 3em the number of degrees of freedom versus temperature .
as it is clear from the figure , in snyder space the number of degrees of freedom for a three - dimensional harmonic oscillator is reduced from 6 to 4 at the high temperature regime.,title="fig:",width=288 ]    considering the limit @xmath162 is useful at least for two reasons .
first , we see that the modified thermodynamical relations recover the standard results at the sufficiently low temperature regime when the minimal length effect are negligible ( correspondence principle ) .
second , one can estimate the magnitude of the order of the quantum gravity corrections to the thermodynamical quantities which may potentially lead to observable effects . in this limit ,
the series in the relation ( [ ho - pf - asym ] ) converges and one can use it to obtain all the minimal length ( quantum gravity ) corrections to the internal energy and specific heat . substituting the relation ( [ ho - pf - asym ] ) into ( [ ho - tpf ] ) and then using the definition @xmath163 ,
it is easy to show that the internal energy in the limit of @xmath164 will be @xmath165 where @xmath166 is the internal energy for the standard nondeformed @xmath0 independent three - dimensional harmonic oscillators . in this limit
, the specific heat can be also obtained by substituting the above relation into the definition @xmath167 which gives @xmath168 where @xmath169 is the specific heat of the corresponding nondeformed case .    as is clear from the relations ( [ ho - ie - asym ] ) and ( [ ho - sh - asym ] ) , the first minimal length ( quantum gravity ) corrections to the internal energy and specific heat of the harmonic oscillator are of the order of @xmath170 . to estimate the magnitude of the order of these corrections ,
consider the vibrational oscillations of a carbon monoxide molecule which may be modeled with the oscillators of mass @xmath171 and frequency @xmath172 . with this numerical value for the mass and also considering @xmath173 and @xmath174
, we get @xmath175 . therefore , the first minimal length corrections to the internal energy and specific heat of the harmonic oscillator in snyder space are of the order of @xmath176 which is too small to be experimentally detected for the accessible energy scales .
appearance of the hamiltonian systems with nontrivial structure in the context of phenomenological quantum gravity candidates , such as the doubly special relativity theories with deformed noncommutative phase spaces , naively suggests the revision of the statistical mechanics formalism . in the first step ,
we have formulated the statistical mechanics of a general hamiltonian system in a covariant ( chart - independent ) manner by means of the symplectic geometry .
the results show that the two properties of the phase space , as a symplectic manifold , play distinguished roles in statistical consideration of a system : the topology of the phase space as a global property and the local form of the symplectic structure that determines the geometry of the phase space . for a topologically trivial phase space with the canonical symplectic structure ,
the standard statistical mechanics emerges as it is expected .
the subtleties arise when the phase space has nontrivial topology or geometry such as the phase spaces with curved momentum space and noncommutative structure .
the topology of the phase space turns out to be related to the total number of microstates , for instance , the number of accessible microstates will be finite for a system with compact phase space .
the symplectic structure also affects the probability distribution of the microstates .
the canonical form of the symplectic structure in terms of the positions and momenta of the particles plays an important role in the standard statistical mechanics .
it induces a uniform probability distribution over the microstates by defining a uniform measure on the corresponding phase space .
however , the noncommutative phase space always provides a noncanonical symplectic structure in terms of the positions and momenta of the particles and the probability distribution then will be nonuniform .
this is a general result for the phase spaces with noncanonical ( noncommutative ) structure which can be immediately realized from the covariant formulation of the phase space .
we implemented the covariant formalism in order to study the statistical mechanics in snyder noncommutative space as an explicit example of a phase space with nontrivial topology and geometry .
we obtained the associated partition function from which all the thermodynamical properties of the system can be obtained .
as a particular example , we obtained the partition function for the three - dimensional harmonic oscillator and we have shown that our result is in good agreement with that which arises from the full quantum consideration at the high temperature regime . while there is no analytical expression for the full quantum partition function of the harmonic oscillator in snyder space , the semiclassical approximation provides an analytical partition function which is applicable at the high temperature regime . using the obtained partition function ,
we have studied the thermodynamical properties of the system of harmonic oscillators in snyder space and our analysis shows that the number of microstates will be drastically reduced in this setup due to the existence of a minimal length .
this result justifies our general claim for the deformed spaces which take into account a minimal length scale .
apart from the details of the models which deal with a universal minimal length such as the noncommutative spaces , doubly special relativity theories , the generalized uncertainty principle , and polymer quantization , the covariant formalism reveals the main role of the minimal length in statistical mechanics : the microstates with higher energy or momenta are less probable when there is a minimal length scale for the system . also , we calculated the quantum gravity corrections to the internal energy and specific heat of a system of harmonic oscillators .
we estimated the order of magnitude of the first minimal length corrections which are of the order of @xmath177 for both of the internal energy and specific heat .
the equilibrium as an essential criterion for the statistical system is also considered in the appendix by means of the dynamically invariant liouville volume which enters in the definition of the gibbs entropy and partition function of the statistical system .
relation ( [ decomposition 2 ] ) can be seen as a particular case in which the system is divided into statically independent subsystems and subsystems were taken to be the particles themselves .
it is clear that the statistical independence between particles will be spoiled for a system in which the particles interact with each other .
however , when the system is in equilibrium , one can still divide it into @xmath178 macroscopic subsystems ( @xmath179 ) which are statistically independent . in this respect ,
the total phase space @xmath9 and the associated symplectic structure @xmath10 will be @xmath180 the corresponding liouville volume is then given by @xmath181 where @xmath182 is the number of particles of the subsystems and is assumed to be sufficiently large @xmath183 to guarantee that the subsystems are macroscopic .
the statistical independence of the subsystems is then defined as @xmath184 this criterion is usually employed as an intrinsic characterization of an equilibrium probability distribution @xcite .
indeed , the system is in equilibrium if components of the system are in relative equilibrium with respect to each other .
based on this internal definition of the equilibrium , the gibbs state can be understood as an equilibrium state as follows .
consider the total hamiltonian function of the system as @xmath185 where @xmath186 are the hamiltonian functions of the @xmath178 macroscopic subsystems and we also disregard the interaction between subsystems based on the fact that the relative number of particles which take part in the interactions is negligible compared to @xmath182 .
so , using the relations ( [ can rho class ] ) , ( [ gamma split m ] ) , and ( [ decomposition 8 ] ) , the gibbs state turns out to be an equilibrium state based on ( [ decomposition 5 ] ) .
additionally , it is straightforward to show that @xmath187=0\,,\ ] ] which is the criterion for the kinematical and dynamical separability of the subsystems .
the relation ( [ equilibrium ] ) then can be viewed as the equilibrium criterion when the system is in a gibbs state .
note that the decomposition ( [ decomposition 8 ] ) holds only over not - too - long intervals of time since the effects of interaction of the subsystems will eventually be dominated even if such interactions are too weak such as the universal gravitational effects @xcite .
another more common definition for the equilibrium statistical state of the system is @xmath188 in which @xmath189 denotes the density that corresponds to the equilibrium state .
the fact that the gibbs state is an equilibrium state then will be justified through the liouville equation @xmath190 in which the liouville theorem ( [ liouvillee ] ) is used .
the liouville equation should be satisfied by any statistical state as well as gibbs one and it is , indeed , the necessary requirement for the consistent probabilistic interpretation of @xmath95 . using the fact that @xmath191 , the liouville equation ( [ liouville e ] ) can be rewritten in the more well - known form @xmath192 from the relation ( [ can rho class ] ) , it is clear that @xmath193 , and the gibbs state then satisfies the relation @xmath194 , which shows that it is indeed an equilibrium state through the definition ( [ equili ] ) .
if we use the liouville equation ( [ liouville e ] ) and consider the fact that the liouville volume @xmath77 does not change explicitly with time , it turns out that @xmath195 which shows that the fine - grained entropy ( [ gibbs entropy ] ) does not change with time as expected for closed systems . the above relation and the liouville equation ( [ liouville e ] ) guarantee that the number of microstates does not change through the dynamical evolution of the system , which makes the statistical formulation a consistent setup to give the thermodynamical properties in any ensemble .
j. kowalski - glikman , phys .
b * 547 * ( 2002 ) 291 + j. kowalski - glikman and s. nowak , class .
quantum grav . * 20 * ( 2003 ) 4799 + f. girelli , e. r. livine and d. oriti , nucl . phys .
b * 708 * ( 2005 ) 411 + j. kowalski - glikman , int .
phys . a * 28 * ( 2013 ) 1330014 + s. meljanac , a. pachol , a. samsarov and k. s. gupta , phys .
d * 87 * ( 2013 ) 125009 .
j. kowalski - glikman , phys .
lett . a * 299 * ( 2002 ) 454 .
g. amelino - camelia , l. freidel , j. kowalski - glikman and l. smolin , phys .
d * 84 * ( 2011 ) 084010 + g. amelino - camelia , l. freidel , j. kowalski - glikman and l. smolin , gen .
. grav . * 43 * ( 2011 ) 2547 .
j. kowalski - glikman , phys .
lett . a * 299 * ( 2002 )
454 + x. zhang , l. shao and b. -q .
ma , astropart .
* 34 * ( 2011 )
840 + s. das , s. ghosh and d. roychowdhury , phys . rev .
d * 80 * ( 2009 ) 125036 + s. das and d. roychowdhury , phys .
d * 81 * ( 2010 ) 085039 + n. chandra and s. chatterjee , phys .
d * 85 * ( 2012 ) 045012 + g. amelino - camelia , n. loret , g. mandanici and f. mercati , int . j. mod .
d * 21 * ( 2012 ) 1250052 .
g. t hooft , nucl .
b * 335 * ( 1990 ) 138 + e. witten , phys .
d * 44 * ( 1991 ) 314 + l. susskind and j. uglum , phys .
d * 50 * ( 1994 ) 2700 + a. strominger and c. vafa , phys .
b * 379 * ( 1996 ) 99 + g. t. horowitz and a. strominger , phys .
* 77 * ( 1996 ) 2368 + j. m. maldacena , [ arxiv : hep - th/9607235 ] + j. m. maldacena , a. strominger and e. witten , jhep * 9712 * ( 1997 ) 002 .    c. rovelli , phys .
* 77 * ( 1996 ) 3288 + a. ashtekar , j. baez , a. corichi and k. krasnov , phys .
* 80 * ( 1998 ) 904 + m. domagala and j. lewandowski , class .
quantum grav .
* 21 * ( 2004 ) 5233 + k. a. meissner , class .
quantum grav .
* 21 * ( 2004 ) 5245 .
d. j. gross and p. f. mende , nucl .
b * 303 * ( 1988 ) 407 + d. amati , m. ciafaloni and g. veneziano , phys .
b * 216 * ( 1989 ) 41 + m. kato , phys .
b * 245 * ( 1990 ) 43 + l. garay , int .
a * 10 * ( 1995 ) 145 .    c. rovelli and l. smolin , nucl .
b * 331 * ( 1990 ) 80 + c. rovelli and l. smolin , nucl .
b * 442 * ( 1995 ) 593 + a. ashtekar and j. lewandowski , class .
quantum grav .
* 14 * ( 1997 ) a55 + a. ashtekar and j. lewandowski , class .
quantum grav .
* 21 * ( 2004 ) r53 .
o. f. dayi and a. jellal , phys .
lett . a * 287 * ( 2001 ) 349 + f. g. scholtz and j. govaerts , j. phys .
a * 41 * ( 2008 ) 505003 + w. -h .
huang and k. -w .
huang phys .
b * 670 * ( 2009 ) 416 + a. shariati , m. khorrami and a. h. fatollahi , j. phys .
a * 43 * ( 2010 ) 285001 .
m. maggiore , phys .
b * 319 * ( 1993 ) 83 + m. maggiore , phys .
d * 49 * ( 1994 ) 5182 + a. kempf , g. mangano and r. b. mann , phys .
d * 52 * ( 1995 ) 1108 + a. kempf and g. mangano , phys .
d * 55 * ( 1997 ) 7909 + k. nozari and a. etemadi , phys .
d * 85 * ( 2012 ) 104029 .
s. k. rama , phys .
b * 519 * ( 2001 ) 103 + l. n. chang , d. minic , n. okamura and t. takeuchi , phys .
d * 65 * ( 2002 ) 125028 + m. lubo , phys .
d * 68 * ( 2003 ) 125004 + k. nozari and b. fazlpour , gen .
reltiv . grav .
* 38 * ( 2006 ) 1661 + p.
wang , h. yang and x. zhang , jhep * 08 * ( 2010 ) 043 + d. mania and m. maziashvili , phys . lett .
b * 705 * ( 2011 ) 521 + p. wang , h. yang and x. zhang , phys .
b * 718 * ( 2012 ) 265 + b. vakili and m. a. gorji , j. stat .
( 2012 ) p10013 .
a. ashtekar , s. fairhurst and j. willis , class .
quantum grav . * 20 * ( 2003 ) 1031 + k. fredenhagen and f. reszewski , class .
quantum grav .
* 23 * ( 2006 ) 6577 + a. corichi , t. vukainac and j. a. zapata , phys .
d * 76 * ( 2007 ) 044016 .
g. m. hossain , v. husain and s. s. seahra , class .
quantum grav .
* 27 * ( 2010 ) 165013 + v. husain , s. s. seahra and e. j. webster , phys .
d * 88 * ( 2013 ) 024014 + g. chacn - acosta , e. manrique , l. dagdug and h. a. morales - tcotl , _ sigma _ * 7 * ( 2011 ) 110 + e. castellanos and g. chacn - acosta , phys .
b * 722 * ( 2013 ) 119 + g. chacn - acosta and h. hernandez - hernandez , int .
d * 24 * ( 2015 ) 1550033 .
p. nicolini , a. smailagic and e. spallucci , phys .
b * 632 * ( 2006 ) 547 + s. ansoldi , p. nicolini , a. smailagic and e. spallucci , phys . lett .
b * 645 * ( 2007 ) 261 + y. s. myung , y .- w .
kim and y. -j .
park , jhep * 0702 * ( 2007 ) 012 + k. nozari and s. h. mehdipour , class .
quantum grav . * 25 * ( 2008 ) 175015 + r. banerjee , b. r. majhi and s. samanta , phys .
d * 77 * ( 2008 ) 124035 + r. banerjee , b. r. majhi and s. k. modak , class . quantum grav .
* 26 * ( 2009 ) 085010 + c. bastos , o. bertolami , n. c. dias and j. n. prata , phys .
d * 80 * ( 2009 ) 124038 .
f. scardigli , phys .
b * 452 * ( 1999 ) 39 + a. j. m. medved and e. c. vagenas , phys .
d * 70 * ( 2004 ) 124021 + g. amelino - camelia , m. arzano , y. ling and g. mandanici , class .
quantum grav .
* 23 * ( 2006 ) 2585 + y. s. myung , y. -w .
kim and y. -j .
park , phys .
b * 645 * ( 2007 ) 393 + k. nozari and a. s. sefiedgar , phys .
b * 635 * ( 2006 ) 156 + a. bina , s. jalalzadeh and a. moslehi , phys . rev .
d * 81 * ( 2010 ) 023528 + l. xiang and x. q. wen , jhep * 0910 * ( 2009 ) 046 .


we will consider the ginzburg
landau model of superconductivity .
if a @xmath1-dimensional superconducting sample with ginzburg  landau parameter @xmath0 is submitted to a uniform magnetic field of strength @xmath2 , then ( by a theorem of giorgi and phillips  @xcite ) there exists a field strength @xmath3 such that if @xmath4 , then the sample will be in its normal state , i.e. superconductivity is lost altogether .
it is at first sight natural to expect this phenomenon to mark a monotone transition , i.e. to expect that the material is in its superconducting ( possibly mixed ) state for all @xmath5 .
indeed , such a monotonicity result has been proved recently in a number of geometric situations and in both @xmath1 and @xmath6 dimensional settings  @xcite in the case where the ginzburg
landau parameter @xmath0 is large ( it also follows from asymptotic expansions obtained in other works such as @xcite ) .
however , nature does not support this monotonicity in general .
the famous little  parks effect  @xcite shows that for narrow cylinders ( or @xmath1d annuli ) one has an oscillatory behavior instead of monotonicity .
in this paper we will establish such ` oscillatory ' effects rigorously in different geometric settings .
the lack of monotonicity comes from the topology / geometry of the annulus .
it is natural to ask whether one can get such an oscillatory effect for ( non - vanishing ) magnetic fields defined on domains without topology . from the previous investigations
@xcite we know this to be impossible for a uniform magnetic field , but how about more general fields ?
the analysis of constant magnetic fields tells us that this question is linked to a purely spectral problem , namely whether the first eigenvalue of the schrdinger operator @xmath7 is monotone increasing in the parameter ( strength of the magnetic field ) @xmath8 for sufficiently large values of @xmath8 .
this property has been called ` strong diamagnetism ' and has been proved for large classes of magnetic fields  it is even ` generically ' satisfied @xcite
. however , we produce counterexamples in the general case .
the ginzburg  landau theory of superconductivity is based on the energy functional @xmath9 here @xmath10 is a material parameter ( the ginzburg
landau parameter ) , @xmath11 is a parameter measuring the intensity of the external magnetic field .
the domain @xmath12 is the part of space occupied by the superconducting material . for @xmath13
there are two natural choices .
one can take @xmath14 .
that will not be our choice here because for reasons of simplicity we want to avoid an unnecessary technical complication connected with unbounded domains in @xmath15 ( for details on how to handle this issue see @xcite ) .
one can also  and that will be our convention here
take @xmath13 to be the smallest simply connected domain containing @xmath16 , i.e. the union of @xmath16 and all the ` holes ' in @xmath16 .
the function @xmath17 is the profile of the external magnetic field .    in the setting of bounded @xmath18 the functional @xmath19
is naturally defined on @xmath20 .
the functional is immediately seen to be gauge invariant , @xmath21 . the vector field @xmath22 models the induced magnetic vector potential .
the function @xmath23 measures the superconducting properties of the material , with @xmath24 being a measure of the local density of cooper pairs .
we say that a minimizer @xmath25 of the ginzburg
landau functional is trivial if @xmath26 and @xmath27 . in each of the situations we will encounter
, the notation @xmath28 will be reserved for a fixed choice of vector potential with @xmath29 . for trivial minimizers
we clearly have @xmath30 . for a nontrivial minimizer
the functional must be negative , since one gets from the euler - lagrange equations of a minimizer that @xmath31 if @xmath25 is a minimizer .
we define the set @xmath32 following @xcite one typically defines the third critical field to be given by @xmath33 , which is finite by @xcite .
however , unless @xmath34 is an interval , this definition is not the only natural one to take  see @xcite for a discussion .
we will see below that this is not always the case .
let @xmath35 denote the annulus with inner radius @xmath36 and outer radius @xmath37 , let @xmath38 .
in this case we will write @xmath39 i.e. the disc of radius @xmath37 .
[ thm : main3 ] there exists an annulus @xmath16 and a @xmath40 such the set @xmath41 is not an interval .
the mechanism behind this result is a convergence of the magnetic quadratic form on the annulus to the corresponding form on the circle .
this convergence was already noticed in the works @xcite , where also ` annuli ' of non - uniform width were considered .
it is likely that one could deduce theorem  [ thm : main3 ] from these works , however , we prefer to give a simple independent proof which also emphasizes the connection to the bohm  aharonov - effect .
[ rem : largerkappa ] by shrinking the inner radius @xmath36 of the annulus , we can get @xmath42 as large as we want , since the eigenvalues of the limiting problem will then cross at a level @xmath43 .
in particular it is possible to have @xmath44 , which means that theorem  [ thm : main3 ] also applies to superconductors of type ii .
one may criticize the result of theorem  [ thm : main3 ] on two accounts .
one could desire not to have the topology fixed a priori , but rather have it generated by localization properties of the minimizer .
also most previous mathematical analysis has considered the limit of large values of @xmath0 .
one can show that for sufficiently large values of @xmath0 the set @xmath34 of a superconducting sample in the shape of an annulus will behave as the one of the disc with the same outer radius , and it is known that for the disc and with constant magnetic field  for sufficiently large values of @xmath0@xmath34 is indeed an interval @xcite .
our next theorem remedies these defects .
[ thm : hc3-oscillation ] let @xmath16 be the unit disc in @xmath45 .
there exists an everywhere positive magnetic field @xmath46 such that for all @xmath40 there exists @xmath47 satisfying that @xmath34 is not an interval .
in fact , the magnetic field can be chosen as @xmath48 , where @xmath49 is some sufficiently small constant .
theorem  [ thm : hc3-oscillation ] follows directly from theorem  [ thm : main ] ( or theorem  [ thm : intasymptot ] ) below using ( * ? ? ?
actually , it easily follows from theorem  [ thm : intasymptot ] below , that for all integers @xmath50 we can choose @xmath51 so small that @xmath52 will consist of at least @xmath53 intervals for all @xmath0 sufficiently large .      for easy reference
we collect the notation and assumptions concerning the magnetic fields that we will treat .
we will work on an open set @xmath16 being one the following three cases @xmath54 .
[ ass : magnetic ] suppose that @xmath55 , is a non - negative , radial magnetic field , possessing five continuous derivatives in an open neighborhood @xmath56 of the unit circle @xmath57 .
define @xmath58 and assume that @xmath59 and write @xmath60 when @xmath61 , we assume that @xmath62 where @xmath63 is the spectral constant recalled in appendix  [ sec : dg ] .
when @xmath64 , we impose the stronger assumption that @xmath65 has a unique , non - degenerate minimum at @xmath66 and that @xmath67    the assumptions assure that ground state eigenfunctions will be localized near @xmath66 . for @xmath64
, we have @xmath68 by assumption , but that is not necessarily true in the cases with boundary .
we define @xmath69 i.e. @xmath70 denotes the magnetic flux through the unit disc .    for a magnetic field satisfying assumption  [ ass : magnetic ] and @xmath71
, we study the lowest eigenvalue @xmath72 of the self - adjoint magnetic schrdinger operator @xmath73 in @xmath74 . here
@xmath28 is a magnetic vector potential associated with the magnetic field @xmath75 .
we refer the reader to section  [ sec : prel ] for a more complete definition of this operator and the eigenvalue .
we will study this eigenvalue problem in three cases , namely for @xmath16 the unit disc , the complement of the unit disc and the whole plane @xmath45 .
if @xmath16 has a non - empty boundary we impose a magnetic neumann boundary condition .
the next theorem states that if @xmath16 is the unit disc or its complement , then special choices of magnetic fields satisfying assumption  [ ass : magnetic ] will give that the function @xmath76 is _ not _ monotonically increasing for large @xmath8 . before stating the theorems , we remind the reader that @xmath77 are universal ( spectral ) constants coming from the de gennes model operator  this is recalled in appendix  [ sec : dg ] .    [
thm : main ] let @xmath16 be the unit disc or its complement .
suppose that @xmath75 satisfies assumption  [ ass : magnetic ] .
assume that @xmath49 and @xmath78 then for all @xmath79 there exist @xmath80 and @xmath81 , with @xmath82 , such that @xmath83 on the other hand , if @xmath84 then there exists @xmath79 such that @xmath85 is monotone increasing on @xmath86 .    in particular , holds for the magnetic field @xmath87 for all @xmath49 sufficiently small  the flux in this case is @xmath88 .
therefore , this magnetic field will not display monotonicity for large field strength .
theorem  [ thm : main ] is a consequence of the following precise asymptotic formulas for the ground state eigenvalue given as theorem  [ thm : extasymptot ] and theorem  [ thm : intasymptot ] .
[ thm : extasymptot ] suppose that @xmath16 is the complement of the unit disc , that @xmath75 satisfies assumption  [ ass : magnetic ] and that @xmath49
. then there are constants @xmath89 and @xmath90 such that if @xmath91 then , as @xmath92 , @xmath93    by a careful reading of the proof , one will realize that the constant @xmath89 is independent of @xmath51 but that @xmath90 depends on @xmath51 . however ,
for our purposes this extra information is irrelevant .
a similar expansion holds in the interior of the unit disc .
[ thm : intasymptot ] suppose that @xmath16 is the unit disc , that @xmath75 satisfies assumption  [ ass : magnetic ] and that @xmath49
. then there exist constants @xmath94 and @xmath95 such that if @xmath96 then , as @xmath92 , @xmath97    notice that for the disc or its complement , the constant magnetic field @xmath98 satisfies assumption  [ ass : magnetic ] , so theorems  [ thm : extasymptot ] and  [ thm : intasymptot ] imply this special case .
this agrees with the calculations in  @xcite ( see also  @xcite ) . in the case of constant field
is not satisfied , and one does get monotonicity of the ground state energy for large magnetic field ( this is discussed in detail in  @xcite ) .
we continue with @xmath99 . here , we are only able to destroy monotonicity in the case @xmath100 .
[ thm : wholeplane ] let @xmath99 .
then , for all @xmath101 and all magnetic fields satisfying assumption  [ ass : magnetic ] there exists a @xmath79 such that @xmath72 is monotonically increasing for @xmath102 . however , if @xmath103 , then @xmath85 is not monotone increasing on any unbounded half - interval .    as for the disc and the exterior of the disc , the proof of this result goes via asymptotic expansions .
[ thm : r2deltapos ] suppose that @xmath99 , and that @xmath75 satisfies assumption  [ ass : magnetic ] with @xmath49 . then , as @xmath92 , @xmath104    [ thm : r2deltazero ] let @xmath105 and @xmath106 be the spectral constants from   and   respectively .
suppose that @xmath99 , and that and that @xmath75 satisfies assumption  [ ass : magnetic ] with @xmath100 .
there exist constants @xmath107 and @xmath108 such that if @xmath109 then , as @xmath92 , @xmath110    in all of the results above the ground state has angular momentum @xmath111 ( to leading order in @xmath8 ) .
we recall that @xmath112 is the total flux through the unit disc  the bounded domain enclosed by the curve where we have localization .
the possibility to obtain non - monotonicity comes from the condition that @xmath113 must be an integer , which leads to frustration .
this is similar to examples in  @xcite .
theorem  [ thm : wholeplane ] raises the question whether one can break strong diamagnetism with a strictly positive magnetic field on the whole plane .      in the next section
we define the operators involved and perform the fourier decomposition reducing the study to a family of ordinary differential operators .    in section  [ sec : annulus ]
we prove a non - monotonicity result for an annulus and use that to prove theorem  [ thm : main3 ] . in section  [ sec : section4 ]
we work in the exterior of the unit disc and prove theorem  [ thm : extasymptot ] .
we indicate in section  [ sec : disc ] how the proof of theorem  [ thm : extasymptot ] can be modified to give the proof of theorem  [ thm : intasymptot ] . in section  [ sec : nonmondiscanddiscext ] we see how theorem  [ thm : extasymptot ] and theorem  [ thm : intasymptot ] imply theorem  [ thm : main ] .    in section  [ sec : planeg0 ]
we prove theorem  [ thm : r2deltapos ] and in section  [ sec : plane0 ] we prove theorem  [ thm : r2deltazero ] .
these two results are used to prove theorem  [ thm : wholeplane ] .
we consider the self - adjoint magnetic neumann schrdinger operator @xmath114 with domain @xmath115    here @xmath116 is the interior unit normal to @xmath117 , @xmath118 and @xmath119 is the strength of the magnetic field .    in general , for a self - adjoint operator @xmath120 that is semi - bounded from below we will write @xmath121 for the lowest point of the spectrum of @xmath120 .    in the case of the disc or
if @xmath122 as @xmath123 the operator has compact resolvent ( see  @xcite ) . if @xmath16 is unbounded and if @xmath124 , then the essential spectrum will be bounded below by @xmath125 ( see  @xcite for the case of @xmath45 and @xcite for the case of the exterior of the disc ) . in any case , as it will follow by the results below , @xmath72 will be an eigenvalue .
we will work in domains @xmath16 that are rotationally symmetric .
for that reason , we will often work in polar coordinates @xmath126 here @xmath127 will be an interval .    moreover
, we will work with magnetic fields that depends only on @xmath128 .
for a radial magnetic field @xmath129 we will work with the gauge @xmath130 where , so @xmath131 has an immediate interpretation in terms of the flux through the disc of radius @xmath132 . ]
@xmath133 in calculations , we will often meet the expression @xmath134 .
this we can write as @xmath135 where @xmath136 thus , under assumption  [ ass : magnetic ] , as @xmath137 ,    @xmath138    with @xmath139 .
the expression for the operator @xmath140 in polar coordinates becomes @xmath141 we decompose the hilbert space as ( here @xmath142 denotes any of the intervals @xmath143 , @xmath144 , @xmath145 or @xmath146 ) @xmath147 that is , for a function @xmath148 , we write @xmath149 where @xmath150 .
next , we write the operator @xmath140 corresponding to this decomposition as @xmath151 where @xmath152 is the self - adjoint operator acting in @xmath153 , given by @xmath154 with neumann boundary conditions at the endpoints of @xmath142 .
the quadratic form corresponding to @xmath152 is given by @xmath155=\int_i \bigl[|\psi'(r)|^2 + \bigl(\frac{m}{r}-ba(r)\bigr)^2|\psi(r)|^2   \bigr]r\ , dr.\ ] ] it holds that @xmath156
in this section we will let @xmath157 we aim to prove theorem  [ thm : main3 ] .
we first notice the non - monotonicity of the function @xmath76 .
[ thm : mainlin ] let @xmath158 and @xmath159 . then the operator @xmath140 in the annulus @xmath16 satisfies @xmath160 in particular , the function @xmath76 is monotonically decreasing around @xmath161 .
one might suspect that some properties of @xmath140 are carried over to some model problem on the circle , as @xmath162 .
let @xmath163 be the self - adjoint operator @xmath164 in @xmath165 with periodic boundary conditions .
its spectrum is easily seen to consist of eigenvalues @xmath166 . in particular
@xmath167    our next theorem states that @xmath72 will tend to @xmath168 as @xmath162 .
[ thm : main2 ] let @xmath71
. then @xmath169    as a direct consequence of theorem  [ thm : main2 ] it is possible to find an annulus such that the function @xmath76 is monotonically increasing and decreasing alternatively as many times as desired .
another direct consequence of theorem  [ thm : main2 ] is that , although the diamagnetic inequality tells us that @xmath170 for all @xmath71 we can actually get @xmath72 to be arbitrary close to zero if @xmath171 , @xmath172 by choosing @xmath37 close enough to @xmath36 .
theorem  [ thm : main2 ] can easily be extended to thin cylinders in three dimensions , since the third variable then separates .      in this section
we will prove the spectral results theorem  [ thm : mainlin ] and theorem  [ thm : main2 ] .
we will work in polar coordinates .
we recall that here @xmath158 .
let @xmath173 denote the potential in the quadratic form @xmath174 in  .
we start by showing that if @xmath175 and @xmath176 then @xmath177    the function @xmath178 is positive for @xmath179 . indeed ,
@xmath180 . if @xmath181 then @xmath182 is decreasing , and @xmath183 . if @xmath184 then @xmath185 which is clearly positive for all @xmath179 . the inequality   follows by a comparison of quadratic forms .
next , we show that if @xmath186 then @xmath187    by perturbation theory it holds that @xmath188 where @xmath189 denotes the eigenfunction corresponding to @xmath190 .
moreover the factor @xmath191 is negative for all @xmath192 if @xmath193 .
inserting this into   gives     it is now easy to finish the proof of theorem  [ thm : mainlin ] .
let @xmath159 .
inequality   and analytic perturbation theory imply that @xmath194 for @xmath8 in a neighborhood of @xmath195 .
since , by  , it holds that the derivative of @xmath196 is negative at @xmath197 the same is true for the derivative of @xmath72 . by continuity of the derivative this holds in a neighborhood of @xmath197 .
in particular we conclude that the function @xmath198 is strictly decreasing for these values of @xmath8 .    since @xmath199 theorem  [ thm : main2 ] is a direct consequence of the fact that , for @xmath200 and @xmath119 , @xmath201    to get an upper bound we use a trial state .
in fact , we use the simplest possible one . let @xmath202
. then @xmath203 .
a simple calculation shows that @xmath204   & = \lim_{\ro\searrow\ri } \biggl ( \frac{2m^2}{\ro+\ri}\frac{\log\ro-\log\ri}{\ro-\ri } -bm + \frac{b^2}{8}\bigl(\ri^2+\ro^2\bigr ) \biggr)\\ & = \bigl(\frac{m}{\ri}-\frac{b\ri}{2}\bigr)^2 . \end{aligned}\ ] ] hence @xmath205 .
the lower bound is obtained by using the potential @xmath206 .
let @xmath189 be a normalized eigenfunction corresponding to @xmath190 .
then @xmath207 \geq   \int_{\ri}^{\ro } \bigl(\frac{br}{2}-\frac{m}{r}\bigr)^2|u|^2 r\,dr \geq \min_{\ri\leq r\leq \ro}\bigl(\frac{br}{2}-\frac{m}{r}\bigr)^2.\ ] ] since @xmath208 as @xmath162 we conclude that @xmath209 this completes the proof of  , and thus finishes the proof of theorem  [ thm : main2 ] .      in this section
we prove theorem  [ thm : main3 ] .
we recall the reader that @xmath210 below denotes the disc with radius @xmath37 , centered at the origin .
we need the following lemma , and refer to  @xcite for its proof .
[ lem : curl ] let @xmath36 be fixed and let @xmath211 .
there exists a constant @xmath212 ( independent of @xmath37 ) such that for all @xmath213 we have @xmath214 combining this with the sobolev embedding we get the existence of a constant @xmath215 ( independent of @xmath216 $ ] ) such that for all @xmath213 @xmath217    given @xmath218 , the cauchy inequality implies that @xmath219 and so @xmath220 here we used   and @xmath221 to get the last inequality .
if we choose @xmath222 , then we see that if @xmath223 , then for all @xmath37 sufficiently close to @xmath36 and all @xmath25 , @xmath224    on the other hand , if @xmath225 , then we have ( with @xmath226 and @xmath189 the normalized eigenfunction corresponding to @xmath227 ) @xmath228 for sufficiently small values of @xmath229 .
therefore , by the explicit spectrum of @xmath163 we can choose @xmath40 and @xmath82 such that @xmath230 define @xmath231 . by the convergence of the spectrum given in theorem  [ thm : main2 ] and   we find the existence of @xmath232 such that @xmath233 has a non - trivial minimizer for all @xmath234 and @xmath235 .
on the other hand , it follows from that the minimizer of @xmath236 is trivial for all @xmath237 sufficiently close to @xmath36 .
we conclude the existence of @xmath237 such that there exist non - trivial minimizers when @xmath238 and @xmath239 but not when @xmath240 . since @xmath241
it is clear that @xmath41 is not an interval .
in this section we consider the case @xmath242 and assume that the magnetic field satisfies assumption  [ ass : magnetic ] with @xmath49 .
our aim is to prove theorem  [ thm : extasymptot ] .      before continuing we give an agmon estimate for the lowest eigenfunction .
[ prop : agmon_outside ] assume that @xmath75 satisfies assumption  [ ass : magnetic ] with @xmath49 .
let @xmath243 . then there exist positive constants @xmath244 , @xmath245 and @xmath246 such that if @xmath102 , and if @xmath23 is an eigenfunction of @xmath140 corresponding to an eigenvalue @xmath247 .
then @xmath248    theorem  8.2.4 in  @xcite gives the same estimate with the restriction that the domain should be bounded .
however , since we give a similar agmon estimate in section  [ sec : planeg0 ] with proof we omit the proof here .
we recall that the quadratic form after decomposition is given by ( with @xmath249 from ) @xmath250=\int_{1}^{+\infty}\bigl(|u'(r)|^2 + \bigl(\frac{m}{r}-b a(r)\bigr)^2|u(r)|^2\bigr)r\,dr.\ ] ] notice that at @xmath66 the potential takes the value @xmath251 this suggests that we will find the lowest energy for @xmath252 . that this is the case is the content of the following lemma .
[ lem : restrictiononm ] let @xmath243 .
suppose @xmath253 is an eigenfunction of @xmath140 with eigenvalue @xmath247 .
then @xmath254    we neglect the kinetic energy in the expression for @xmath174 .
recall the calculation  . for @xmath255 , we get @xmath256 so , estimating the quadratic form with the potential , combining   and  , and using proposition  [ prop : agmon_outside ] , we get @xmath257 & \geq   \int_1 ^ 2 \frac{1}{r^2}\bigl(m - b r a(r)\bigr)^2 |u_m(r)|^2 r\ , dr\\   & \geq \int_1 ^ 2 \frac{1}{r^2}\bigl[\frac{1}{2}(m-\phi b)^2-(cb)^2(r-1)^2\bigr]|u_m(r)|^2 r\ , dr\\   & \geq \frac{1}{8}(m-\phi b)^2\bigl(1+\mathcal{o}(b^{-\infty})\bigr ) -\widetilde{c}b,\end{aligned}\ ] ] from which the lemma follows .    [
lem : uniqueeigenvalue ] let @xmath243 .
there exists @xmath79 such that if @xmath258 and @xmath259 , then @xmath152 admits at most one eigenvalue below @xmath260 .
fix @xmath261 with @xmath262 . by the lower bound , we see that there exist @xmath263 such that if @xmath264 , then @xmath265 .
so we will restrict attention to @xmath113 s such that @xmath266 , with @xmath267 .
suppose , to get a contradiction , that @xmath268 are eigenfunctions of @xmath174 corresponding to eigenvalues below @xmath260 .
we write @xmath269 with @xmath270 as @xmath271 .
so @xmath272 using the agmon estimates , this yields the following bound on normalized functions @xmath273 in @xmath274 .
@xmath275 & \geq\int_1^{\infty } \bigl ( |v'(r)|^2 +   \frac{1}{r^2}\bigl(\delta m - b \delta ( r-1 )
\bigr)^2 |v(r)|^2\bigr ) r\,dr + { \mathcal o}(b^{1/2})\\ & = \widetilde{\q}_m[v ] + { \mathcal o}(b^{1/2 } ) , \end{aligned}\ ] ] with @xmath276 = \int_1^{\infty }    |v'(r)|^2 + \bigl(\delta m - b \delta ( r-1 ) \bigr)^2 |v(r)|^2\,dr.\ ] ] by translation and scaling @xmath277 is unitarily equivalent to ( the quadratic form of ) a de gennes operator ( see appendix  [ sec : dg ] ) and therefore has spectrum given by @xmath278 only the first of these @xmath279counted with multiplicity  is below @xmath195 ( for some values of @xmath280 ) , so we reach a contradiction if we have a subspace of dimension @xmath1 on which the quadratic form is small .    [
lem : furtherrestrictiononm ] let @xmath281 .
suppose @xmath152 admits an eigenvalue below @xmath282 .
then there exists a constant @xmath283 such that @xmath284    by lemma  [ lem : restrictiononm ] , @xmath285 .
assuming that @xmath189 is the eigenfunction corresponding to the unique ( by lemma  [ lem : uniqueeigenvalue ] ) eigenvalue @xmath286 below @xmath282 we can use the estimate in   to find that @xmath250\geq \widetilde{\q}_m[u ] + \mathcal{o}(b^{1/2}),\ ] ] with @xmath277 as in  . implementing the change of variable @xmath287 , we get ( here we write @xmath288 ) @xmath289=\delta b\int_0^{+\infty }     \biggl(\rho-\xi_0+\xi_0-\frac{m-\phi b}{(\delta b)^{1/2}}\biggr)^2     de  gennes operator ( see appendix  [ sec : dg ] ) . by noticing that the first eigenvalue @xmath290 has a quadratic minimum @xmath291 at @xmath292 ( and using the bound on @xmath293 )
we find that there exists a positive constant @xmath294 such that @xmath289\geq \biggl[\theta_0\delta b   + c_0\delta b\biggl(\xi_0-\frac{m-\phi b}{(\delta b)^{1/2}}\biggr)^2\,\biggr ] \| v \|^2.\ ] ] the second term above is bounded by some constant times @xmath295 according to the assumption .
this in turn gives the existence of a positive constant @xmath244 such that   holds .    in the remainder of this section
we will always restrict our attention to @xmath113 s satisfying the conclusion of lemma  [ lem : furtherrestrictiononm ] .
the strategy of the rest of the proof is as follows .
we will construct an explicit trial state for the operator @xmath296 ( here we suppress the dependence on @xmath113 and @xmath8 for the simplicity of notation ) .
this trial state will be constructed as the first terms of a formal expansion . by taking only finitely many terms ( for our purposes @xmath6 terms suffice ) and performing a localization
one gets a well - defined trial state . in terms of the objects calculated below our explicit trial state
will be as follows .
let @xmath297 let furthermore , @xmath298 , with @xmath299 , and define ( with suitable @xmath300 , say @xmath301 ) @xmath302 then @xmath303 and @xmath304 by self - adjointness of @xmath305 we get that @xmath306 . since we by lemma  [ lem : uniqueeigenvalue ] know that @xmath305 has at most one eigenvalue near @xmath307
, we can conclude that @xmath286 gives the first terms of the asymptotic expansion of that lowest eigenvalue of  @xmath305 .
we proceed to the termwise construction of the trial state .
since ( by proposition  [ prop : agmon_outside ] ) we have localization around @xmath66 , we implement unitarily the change of variables @xmath308 here , the @xmath51 is included for convenience .
then    @xmath309    here the estimate on the remainder should be understood in the following sense : we will only act with our operator on the function @xmath310 from which is localized near @xmath66 on the scale @xmath311 .
so we may consider @xmath312 as a quantity of order @xmath195 .    by lemma  [ lem : furtherrestrictiononm ] the constant term @xmath313 vanishes to leading order .
for reasons of expositions we will write @xmath314 and not insert the choice @xmath315 until later . recall that @xmath316 is bounded
. integrating by parts , we find ( with @xmath288 ) @xmath317 we expand our operator @xmath305 as @xmath318 and obtain @xmath319 we make the ansatz @xmath320 equating order by order in the relation @xmath321 gives :    [ [ order - b0 ] ] * order @xmath322 : * + + + + + + + + + + + + + + + + + + + + + + + + + +    to leading order we find @xmath323 which is the eigenvalue problem for the de gennes operator discussed in appendix  [ sec : dg ] .
the optimal eigenvalue @xmath307 is attained for @xmath324 and @xmath325 .
[ [ order - b-12 ] ] * order @xmath311 :
* + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    here we get @xmath326 by taking scalar product ( with measure @xmath327 ) , we find @xmath328 via the formulas   we find @xmath329    in particular @xmath330 is independent of @xmath331 . moreover , since we can choose @xmath332 , we can let @xmath333 be the regularized resolvent @xmath334 of @xmath335 .
this regularized resolvent is defined as the inverse of the operator @xmath336 restricted to the space @xmath337 .
so we have , @xmath338.\end{aligned}\ ] ]    [ [ order - b-1 ] ] * order @xmath339 : * + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we get @xmath340 taking scalar product with @xmath341 again gives @xmath342 we will not calculate this expression in all detail .
we are only interested in the dependence on @xmath331 .
an inspection gives that it will be a polynomial of degree two .
we will calculate the coefficient in front of @xmath343 to see that it is positive so that @xmath344 has a unique minimum with respect to @xmath331 .    the term @xmath345 is easily calculated since @xmath346 contains one @xmath343 only .    for the term @xmath347
we find one @xmath331 in @xmath348 and therefore also one in @xmath333 .
the coefficient in front of @xmath331 in that term becomes @xmath349\bigr\rangle\\ = -4\bigl\langle ( \rho-\xi_0)\phi_{\xi_0 } , \bigl(\dg(\xi_0)-\theta_0\bigr)^{-1}_{\text{reg } } \bigl[(\rho-\xi_0)\phi_{\xi_0}\bigr]\bigr\rangle.\end{gathered}\ ] ] so , the coefficient in front of @xmath343 in @xmath344 will be ( see  ) @xmath350\bigr\rangle   = \xi_0\phi_{\xi_0}(0)^2>0.\ ] ] this means that we can write @xmath351 where @xmath89 and @xmath90 depend only on @xmath352 , @xmath51 , @xmath292 and @xmath353 ( but not on @xmath70 ) .
we summarize these findings in a lemma .
[ lem : expansion ] suppose @xmath354 with @xmath355 .
then @xmath356    we have to control the asymptotic expansion in @xmath331 subject to the bound @xmath357 .
define @xmath358 with @xmath359 being the constants from above and @xmath344 being the quadratic function of @xmath331 from .
we also define an approximate eigenfunction by @xmath360 with @xmath361 , @xmath333 given by and @xmath362 being given by solving , i.e. @xmath363.\end{aligned}\ ] ] notice from the explicit form of the operators that @xmath333 depends linearly on @xmath331 and @xmath362 depends quadratically , so @xmath273 is normalized to leading order .
also , by the mapping properties of @xmath334 each @xmath364 is a smooth , rapidly decreasing function ( see lemma  3.2.9 in  @xcite ) .
we can now estimate as follows @xmath365 ) v \|.\end{aligned}\ ] ] by the decay properties of @xmath273 , the last term is bounded by @xmath366 .
our choice of @xmath273 gives that the first term is equal to @xmath367 + b^{-2 } ( \h_2 - \lambda_2 ) v_2 \| , \end{aligned}\ ] ] which is easily seen to be bounded by @xmath368 .
using lemma  [ lem : furtherrestrictiononm ] , theorem  [ thm : extasymptot ] follows from lemma  [ lem : expansion ] by the following argument .
notice that the positive quadratic term in @xmath369 dominates the error term @xmath370 unless @xmath331 is bounded in which case the dependence on @xmath331 in the error term disappears .
this finishes the proof of theorem  [ thm : extasymptot ] .
in this section we will indicate a similar calculation of the ground state eigenvalue in the case of the unit disc , thereby proving theorem  [ thm : intasymptot ] , i.e. we work on @xmath371 and for a magnetic field satisfying assumption  [ ass : magnetic ] .
we mainly give the results of the calculations referring to the exterior case for details .
we will have exponential localization estimate like the one of proposition  [ prop : agmon_outside ] ( with domain of integration being @xmath372 , of course ) .
therefore , also the rough ` localization ' of the relevant angular momenta
lemma  [ lem : restrictiononm]will hold in this case as well .
so we can proceed to make a change of variable to the region near ( on the scale @xmath373 as before ) the boundary .
the leading order terms in the expansion of the operator become very similar to the case of the exterior of the disc :    @xmath374    the same calculations ( using the same ansatz ) as in the previous section show that ( with @xmath375 ) @xmath376 for some constants @xmath94 and @xmath95 , depending only on the spectral parameters and @xmath51 .
thus , theorem  [ thm : intasymptot ] follows from calculations / arguments completely analogous to the ones in section  [ sec : section4 ] and we omit the details .
using the results of theorem  [ thm : extasymptot ] and  [ thm : intasymptot ] it is now easy to prove theorem  [ thm : main ] .
we only consider the case of the disc , the complement of the disc being similar ( using theorem  [ thm : extasymptot ] instead of theorem  [ thm : intasymptot ] ) .
assume first that @xmath377 denote by @xmath182 the function @xmath378 notice that @xmath379 is increasing for all large vaues of @xmath8 .
choose a sequence @xmath380 such that @xmath381 , i.e. is a half - integer .
let @xmath382 .
choose @xmath383 . then , for all sufficiently large @xmath53 , @xmath384 .
so @xmath385 and @xmath386 so we get from the eigenvalue asymptotics that @xmath387 \nonumber \\ & \quad+ \xi_0\phi_{\xi_0}(0)^2[(1/2 - \phi \varepsilon)^2 - 1/4 ] + o(1 ) \nonumber \\ & = \theta_0 \delta \varepsilon - \xi_0\phi_{\xi_0}(0)^2[\phi \varepsilon - \phi^2 \varepsilon^2 ] + o(1),\end{aligned}\ ] ] which is negative for small @xmath388 ( and for all sufficiently large @xmath53 ) since @xmath389 by assumption .
suppose now that @xmath390 we restrict attention to the interval near infinity on which @xmath391 is increasing .
here we can calculate the right - hand derivative @xmath392 so we see that for any @xmath393 there exists @xmath79 such that for all @xmath394 and all @xmath395 , @xmath396 we aim to prove monotonicity of @xmath72 , so it suffices to prove a positive lower bound on its right hand derivative @xmath397 , which exists by perturbation theory .
perturbation theory yields , for any @xmath398 , @xmath399 here we completed the square and used the variational characterization of the eigenvalue in order to get the inequality .    since @xmath400
, for some constant @xmath401 independent of @xmath8 , we can estimate , using the eigenvalue asymptotics and @xmath402 since @xmath403 were arbitrary , we get that @xmath404 in particular , @xmath72 is monotone increasing for large value of @xmath8 if   is satisfied .
in this section we will consider the case @xmath99 and a magnetic field @xmath75 satisfying assumption  [ ass : magnetic ] with @xmath49 . we aim to prove theorem  [ thm : wholeplane ] for @xmath49 .
this , however , follows directly once the asymptotic expansion in theorem  [ thm : r2deltapos ] is obtained , since then it follows that ( see  ( * ? ? ?
* section  2.3 ) ) @xmath405 the proof of theorem  [ thm : r2deltapos ] follows the same idea as the proof of theorem  [ thm : extasymptot ] .
we use a localization of the ground state to restrict the situation to certain values of the angular momentum .
then we show that if we find a trial state with low enough energy , it must be related to the ground state energy .
finally we expand our operator formally and construct a trial state that has the correct energy .
we start with a localization estimate valid for @xmath407 . for @xmath100
it gives the right length scale of the localization .
[ prop : firstagmon ] suppose @xmath75 satisfies assumption  [ ass : magnetic ] with @xmath407 .
let @xmath23 be an eigenfunction of @xmath140 corresponding to an eigenvalue @xmath408 for some @xmath409 .
then there exist positive constants @xmath244 and @xmath246 such that @xmath410 and @xmath411 if @xmath102 .    by the localization estimates of proposition  [ prop : firstagmon ] , the quadratic forms @xmath174 are well approximated by harmonic oscillators , whose ground state eigenvalues are simple .
this implies simplicity of the low - lying eigenvalues of @xmath152 .
[ lem : uniqueeigenvalue2 ] let @xmath101 .
let @xmath412 .
there exists @xmath79 such that if @xmath258 and @xmath259 , then @xmath152 admits at most one eigenvalue below @xmath413 .
the proof of lemma  [ lem : uniqueeigenvalue2 ] is similar to that of lemma  [ lem : uniqueeigenvalue ] and will be omitted .
let @xmath414 be a smooth cut - off function of the real variable @xmath415 satisfying @xmath416 and such that @xmath417 for all @xmath415 , and @xmath418 .
next , let @xmath419 and @xmath420 be positive ( to determined below ) real numbers and define in @xmath45 the functions @xmath421 and @xmath422 via @xmath423 and @xmath424
. then there exists a constant @xmath107 such that @xmath425 next , for @xmath426 , let @xmath427 .
then , pointwise in @xmath428 , it holds that @xmath429 as @xmath430 .
moreover , @xmath431 is differentiable almost everywhere and if @xmath432 its gradient satisfies @xmath433 moreover , @xmath431 is bounded for all @xmath426 , so the function @xmath434 belongs to the form - domain of @xmath140 .    with the ims formula
, we find that @xmath435+\q[\chi_2\tpsi ] & \leq \bigl(2c_1m^2b^{2\alpha}+\lambda+16b^{2\sigma}\bigr)\|\tpsi\|^2\\ & \leq \bigl(2c_1m^2b^{2\alpha } + \delta b+\omega b^{1/2}+16b^{2\sigma}\bigr)\|\tpsi\|^2 . \end{aligned}\ ] ] using that the smallest dirichlet eigenvalue is greater than the smallest value of the magnetic field ( again , see  @xcite ) , we find that @xmath436\geq \delta b\|\chi_1\tpsi\|^2\ ] ] and @xmath437\geq   \bigl(\delta b+\frac{kb^{1 - 2\alpha}}{4m^2}\bigr)\|\chi_2\tpsi\|^2.\ ] ] inserting this into   we find that @xmath438 which can be written @xmath439 choosing @xmath440 we find that all @xmath8s factor out , and hence @xmath441 with @xmath419 so small that the left parenthesis above becomes positive , we find that there exists a constant @xmath108 such that @xmath442 on the support of @xmath421 it holds that @xmath443 , and hence @xmath444 inserting this in   above yields @xmath445 using monotone convergence we find that @xmath446 on the other hand , since @xmath443 on the support of @xmath421 it is clear that @xmath447 combining these two last inequalities we find   with @xmath448 .    to prove we essentially only have to reinsert the @xmath449-estimate in the previous calculations . by monotone convergence and the ims - formula ,
we have @xmath450   - \int |\nabla \phi_\ell|^2 |\psi|^2 \,dx.\end{aligned}\ ] ] the last term is negative , and we can estimate the first term using again the ims - formula and as @xmath451 \leq \q[\chi_1\psi ] + \q[\chi_2 \psi ]   \leq ( \delta b + c_2 b^{1/2 } ) \| \psi \|^2\end{aligned}\ ] ] ( with @xmath452 and using @xmath453 ) . now   follows from  .    with the help of proposition  [ prop :
firstagmon ] , we now get a first control of the involved angular momenta .
[ lem : locangmomfirst ] let @xmath407 .
suppose @xmath454 is an eigenfunction of @xmath140 with eigenvalue below @xmath413 .
then @xmath455    the proof of lemma  [ lem : locangmomfirst ] is similar to the one of lemma  [ lem : restrictiononm]taking into account the weaker localization given by proposition  [ prop : firstagmon]and will be omitted .      by lemma  [ lem : uniqueeigenvalue2 ]
there is at most one eigenvalue of @xmath152 for sufficiently low energy .
so it suffices to construct a trial state .
the trial function ( and all its derivatives ) will be localized on the length scale @xmath311 near @xmath66 ( see for the explicit choice of trial state ) . also the function has support away from @xmath457 .
the calculation is slightly different in different regimes of angular momenta @xmath113 . in this subsection
, we consider angular momenta satisfying that latexmath:[\ ] ] in particular this expression is positive .
for @xmath584 , we define the montgomery type operator @xmath585 as a self - adjoint operator on @xmath586 .
let us denote by @xmath587 the eigenvalues of @xmath536 , with corresponding eigenfunctions @xmath588 , @xmath589 and so on .
[ thm : non - degenerate ] the function @xmath590 has a unique minimum at @xmath551 .
furthermore , the minimum is non - degenerate , i.e. @xmath591    we introduce the following notation , @xmath592    by the estimates in  @xcite we know that @xmath593 a numerical value of @xmath106 , calculated by v. bonnaillie - nol , is @xmath594 .
[ prop : secdiffzeropos ] it holds that @xmath595    perturbation theory .
the eigenvalues of @xmath152 can be solved explicitly in terms of confluent hypergeometric functions , and plotted by the computer .
below we include a figure with the lowest eigenvalue of the limit operator @xmath163 and the lowest eigenvalue of the annulus of inner radius @xmath158 and outer radius @xmath596 .
( dotted ) and the lowest eigenvalue @xmath168 ( solid ) for @xmath597 and @xmath158 .
right : the lowest eigenvalue @xmath72 plotted for @xmath598 .
the dotted lines are the lowest eigenvalue of @xmath152 for @xmath599 , @xmath158 and @xmath596 . ]
( dotted ) and the lowest eigenvalue @xmath168 ( solid ) for @xmath597 and @xmath158 .
right : the lowest eigenvalue @xmath72 plotted for @xmath598 .
the dotted lines are the lowest eigenvalue of @xmath152 for @xmath599 , @xmath158 and @xmath596 . ]
sf was partially supported by the lundbeck foundation , the danish natural science research council and the european research council under the european community s seventh framework program ( fp7/20072013)/erc grant agreement 202859 .