the model of quantum lattice bose gas has a long history and has been suggested initially for conventional superconductors @xcite , quantum crystals such as @xmath4he where superfluidity coexists with a crystalline order.@xcite afterwards , the bose - hubbard ( bh ) model has been studied as a model of the superconductor - insulator transition in materials with local bosons , bipolarons , or preformed cooper pairs.@xcite two - dimensional bh models have been addressed as relevant to describe the superconducting films and josephson junction arrays .
the most recent interest to the system of hard - core bosons comes from the delightful results on bose - einstein ( be ) condensed atomic systems produced by trapping bosonic neutral atoms in an optical lattice .
@xcite the progress in boson physics generates an especial interest around nonlinear topological excitations ( skyrmions , vortices ) which play an increasingly significant role in physics .
one of the fundamental hot debated problems in bosonic physics concerns the evolution of the charge ordered ( co ) ground state of 2d hard - core bh model ( hc - bh ) with a doping away from half - filling .
numerous model studies steadily confirmed the emergence of `` supersolid '' phases with simultaneous diagonal and off - diagonal long range order while penrose and onsager @xcite were the first showing as early as 1956 that supersolid phases do not occur .
the most recent quantum monte - carlo ( qmc ) simulations @xcite found two significant features of the 2d bose - hubbard model with a screened coulomb repulsion : the absence of supersolid phase at half - filling , and a growing tendency to phase separation ( co+bs ) upon doping away from half - filling .
moreover , batrouni and scalettar @xcite studied quantum phase transitions in the ground state of the 2d hard - core boson hubbard hamiltonian and have shown numerically that , contrary to the generally held belief , the most commonly discussed `` checkerboard '' supersolid is thermodynamically unstable and phase separates into solid and superfluid phases .
the physics of the co+bs phase separation in bose - hubbard model is associated with a rapid increase of the energy of a homogeneous co state with doping away from half - filling due to a large `` pseudo - spin - flip '' energy cost .
hence , it appears to be energetically more favorable to `` extract '' extra bosons ( holes ) from the co state and arrange them into finite clusters with a relatively small number of particles .
such a droplet scenario is believed to minimize the long - range coulomb repulsion .
however , immediately there arise several questions .
whether a simple mean - field approximation ( mfa ) and classical continuum model can predict such a behavior ?
what is the detailed structure of the co+bs phase separated state ?
what are the main factors governing the essential low - energy and long - wavelength physics ? is it possible to make use of simple toy models ?
are there the analogies with fermion hubbard model ?
the behavior of the latter under doping away from half - filling is extensively studied in last decade in frames of stripe scenario which implies that doping may proceed through the formation of stripes , or charged domain walls being specific topological solitons in nel antiferromagnets . in particular , emery and
kivelson @xcite argued that the phase separation reflects a universal tendency of the correlated antiferromagnet to expel the doped holes .
white and scalapino @xcite showed that the pure 2d t - j model , in the small - j / t regime and with dopings near @xmath5 , has a striped ground state .
in general , the topological phase separation is believed to be a generic property of 2d fermion hubbard model .    in the paper
we present a topological scenario of co+bs phase separation in 2d hc - bh model with inter - site repulsion .
the extra bosons / holes doped to a checkerboard co phase of 2d boson system are believed to be confined in the ring - shaped domain wall of the skyrmion - like topological defect which looks like a bubble domain in an easy - axis antiferromagnet .
this allows us to propose the mechanism of 2d _ topological _ co+bs _ phase separation _ when the doping of the bare checkerboard co phase results in the formation of a multi - center topological defect , which simplified pseudo - spin pattern looks like a system of bubble co domains with bose superfluid confined in charged ring - shaped domain walls .
the rest of the paper is organized as follows . in sec.ii
we give a short overview of the conventional hc - bh model in frames of a pseudo - spin formalism and mfa description .
in sec.iii we show that the doping in 2d hc - bh system can be accompanied by the formation of a topological defect like a single or multi - center skyrmion .
we present the scenario of the essential low - energy physics for the topologically doped hc - bh system . in sec.iv
we address a simple model of nanoscopic bubble - like domain in a checkerboard co phase for a discrete square lattice that allows us to demonstrate a subtle microscopic structure of such a center with a multi - component order parameter .
the hamiltonian of hard - core bose gas on a lattice can be written in a standard form as follows : @xmath6 where @xmath7 is the projection operator which removes double occupancy of any site , @xmath8 , @xmath9 the chemical potential determined from the condition of fixed full number of bosons @xmath10 or concentration @xmath11 $ ] .
the @xmath12 denotes an effective transfer integral , @xmath13 is an intersite interaction between the bosons . here
@xmath14 are the pauli creation ( annihilation ) operators which are bose - like commuting for different sites @xmath15=0,$ ] if @xmath16 @xmath17=1 - 2n_i$ ] , @xmath18 ; @xmath19 is a full number of sites .
it is worth noting that near half - filling ( @xmath20 ) one might introduce the renormalization @xmath21 , or neutralizing background , that immediately provides the particle - hole symmetry .
the model of hard - core bose - gas with an intersite repulsion is equivalent to a system of spins @xmath22 exposed to an external magnetic field in the @xmath23-direction . for the system with neutralizing background we arrive at an effective pseudo - spin hamiltonian @xmath24 where @xmath25 , @xmath26 , @xmath27 , @xmath28 .
in a linear approximation the hamiltonian for the coupling with an electromagnetic field reads as follows @xmath29_{z}-\frac{1}{2}(\phi _
{ j}-\phi _ { i})^{2}({\hat { \bf s}_{i}}\cdot { \hat { \bf s}_{j}}));\ ] ] @xmath30 where @xmath31 is the vector - potential , and integration runs over line binding the @xmath32 and @xmath33 sites ; @xmath34_{z}=\frac{i}{2}(\hat b_{i}^{\dagger}\hat b_{j}-\hat b_{i}\hat b_{j}^{\dagger}),\ ] ] however , the pseudo - spins are assumed to lie in @xmath35 plane .
then the boson current density operator one may represent to be a sum of the paramagnetic @xmath36_{z},\end{aligned}\ ] ] and diamagnetic @xmath37 terms , respectively .      below we make use of a conventional two - sublattice mfa approach .
for the description of the pseudospin ordering to be more physically clear one may introduce two vectors , the ferromagnetic and antiferromagnetic ones : @xmath38 where @xmath39 .
for the plane where these vectors lie one may introduce two - parametric angular description : @xmath40 . the hard - core boson system in a two - sublattice approximation
is described by two diagonal order parameters @xmath41 and two complex off - diagonal order parameters @xmath42 and @xmath43 .
the complex superfluid order parameter @xmath44 is determined by the in - plane components of ferromagnetic vector : @xmath45 , @xmath46 being the length of the in - plane component of ferromagnetic vector .
so , for a local condensate density we get @xmath47 .
it is of interest to note that in fact all the conventional uniform @xmath48 states with nonzero @xmath49 imply simultaneous long - range order both for modulus @xmath50 and phase @xmath51 .
the in - plane components of antiferromagnetic vector @xmath52 describe a staggered off - diagonal order .
it is worth noting that by default one usually considers the negative sign of the transfer integral @xmath12 , that implies the ferromagnetic in - plane pseudospin ordering .
the model exhibits many fascinating quantum phases and phase transitions .
early investigations predict at @xmath48 charge order ( co ) , bose superfluid ( bs ) and mixed ( co+bs ) supersolid uniform phases with an ising - type melting transition ( co - no ) and kosterlitz - thouless - type ( bs - no ) phase transitions to a non - ordered normal fluid ( no).@xcite the detailed mean - field and spin - wave analysis of the uniform phases of 2d hc - bh model is given by pich and frey.@xcite mfa yields for the conventional uniform supersolid phase @xcite @xmath53 with a constant chemical potential @xmath54 .
it should be noted that the supersolid phase appears to be unstable with regard to small perturbations in the hamiltonian .
the mean - field energy per site of the uniform supersolid phase is written as follows : @xmath55 where @xmath56 .
the cost of doping both for co and co+bs phase appears to be rather high as compared with an easy - plane bs phase at half - filling where the chemical potential turns into zero.@xcite
magnetic analogy allows us to make unambiguous predictions as regards the doping of bh system away from half - filling .
indeed , the boson / hole doping of checkerboard co phase corresponds to the magnetization of an antiferromagnet in @xmath23-direction . in the uniform easy - axis @xmath57-phase of anisotropic antiferromagnet
the local spin - flip energy cost is rather big . in other words ,
the energy cost for boson / hole doping into checkerboard co phase appears to be big due to a large contribution of boson - boson repulsion .
however , the magnetization of the anisotropic antiferromagnet in an easy axis direction may proceed as a first order phase transition with a `` topological phase separation '' due to the existence of antiphase domains.@xcite the antiphase domain walls provide the natural nucleation centers for a spin - flop phase having enhanced magnetic susceptibility as compared with small if any longitudinal susceptibility thus providing the advantage of the field energy .
namely domain walls would specify the inhomogeneous magnetization pattern for such an anisotropic easy - axis antiferromagnet in relatively weak external magnetic field . as concerns the domain type in quasi - two - dimensional antiferromagnet
one should emphasize the specific role played by the cylindrical , or bubble domains which have finite energy and size .
these topological solitons have the vortex - like in - plane spin structure and often named `` skyrmions '' .
the classical , or belavin - polyakov ( bp ) skyrmion @xcite describes the solutions of a non - linear @xmath58-model with a classical isotropic 2d heisenberg hamiltonian and represents one of the generic toy model spin textures .
it is of primary importance to note that skyrmion - like bubble domains in easy - axis layered antiferromagnets were actually observed in the experiments performed by waldner,@xcite that were supported later by different authors ( see , e.g. refs . ) .
although some questions were not completely clarified and remain open until now , @xcite the classical and quantum @xcite skyrmion - like topological defects are believed to be a genuine element of essential physics both of ferro- and antiferromagnetic 2d easy - axis systems .
the magnetic analogy seems to be a little bit naive , however , it catches the essential physics of doping the hc - bh system .
as regards the checkerboard co phase of such a system , namely a finite energy skyrmion - like bubble domain @xcite seems to be the most preferable candidate for the domain with antiphase domain wall providing the natural reservoir for extra bosons .
the classical description of nonlinear excitations in hc - bh model includes the skyrmionic solution given @xmath59 .
@xcite the skyrmion spin texture consists of a vortex - like arrangement of the in - plane components of ferromagnetic @xmath60 vector with the @xmath57-component reversed in the centre of the skyrmion and gradually increasing to match the homogeneous background at infinity .
the simplest spin distribution within classical skyrmion is given as follows @xmath61 @xmath62 where @xmath63 is a global phase ( @xmath64 order parameter ) , @xmath65 skyrmion radius .
the skyrmion spin texture describes the coexistence and competition of the staggered charge order parameter @xmath57 and bs order parameter @xmath46 ( see fig.[fig1 ] ) that reflects a complex spatial interplay of potential and kinetic energies .    ;
b ) the radial distribution of the modulus of the superfluid order parameter @xmath66 ; c ) the radial distribution of the charge density @xmath67 for the charged skyrmion ; d ) the orientation of the superfluid order parameter @xmath66 given two values of the phase order parameter @xmath68 : @xmath69 , respectively . rings in all the pictures correspond to the skyrmion ( bubble domain ) wall , the dashed circle being its center.,width=321 ]    the skyrmion looks like a bubble domain in an easy - axis magnet . it should be noted that the domain wall in such a bubble domain somehow created in the checkerboard co phase of 2d hc - bh system represents an effective ring - shaped reservoir both for bose - superfluid and extra boson / hole . indeed , the soliton energy depends quadratically on the number @xmath70 of bosons bound in domain wall,@xcite similarly to that of homogeneous bs phase.@xcite in other words , one might say about a zero value of the effective boson / hole chemical potential for the co bubble domain configuration , provided it were a ground state .
the numerical calculations performed in frames of a classical continuous model @xcite show that the doped bosons appear to be trapped inside the bubble domain wall .
the spatial distribution of the doped boson / hole density ( @xmath71 ) in the charged skyrmion is shown schematically in fig .
skyrmions are characterized by the magnitude and sign of its topological charge , by its size ( radius ) , and by the global orientation of the spin .
the scale invariance of classical bp skyrmionic solution reflects in that its energy does not depend on radius , and global phase .
an interesting example of topological inhomogeneity is provided by a multi - center bp skyrmion @xcite which energy does not depend on the position of the centers .
the latter are believed to be addressed as an additional degree of freedom , or positional order parameter .    in the continuous model
the classical bp skyrmion is a topological excitation and can not dissipate .
however , the classical static skyrmion is unstable with regard to an external field or anisotropic interactions both of easy - plane and easy - axis type .
small easy - axis anisotropy or external field are sufficient to shrink skyrmion to a nanoscopic size when magnetic length @xmath72 : @xmath73 is of the order of several lattice parameters , and the continuous approximation fails to correctly describe excitations .
nonetheless , abanov and pokrovsky @xcite have shown that the easy - axis anisotropy together with fourth - order exchange term can stabilize skyrmion with radius @xmath74      a skyrmionic scenario in hc - bh model allows us to make several important predictions .
away from half - filling one may anticipate the nucleation of a topological defect in the unconventional form of the multi - center skyrmion - like object with ring - shaped bose superfluid regions positioned in an antiphase domain wall separating the co core and co outside of the single skyrmion .
the specific spatial separation of bs and co order parameters that avoid each other reflects the competition of kinetic and potential energy .
such a _ topological _ ( co+bs ) _ phase separation _ is believed to provide a minimization of the total energy as compared with its uniform supersolid counterpart .
thus , the parent checkerboard co phase may gradually lose its stability under boson / hole doping , while a novel topological self - organized texture is believed to become stable .
the most probable possibility is that every domain wall accumulates single boson , or boson hole .
then , the number of centers in a multi - center skyrmion nucleated with doping has to be equal to the number of bosons / holes .
in such a case , we anticipate the near - linear dependence of the total bs volume fraction on the doping . generally speaking
, one may assume scenario when the nucleation of a multi - center skyrmion occurs spontaneously with no doping .
in such a case we should anticipate the existence of neutral multi - center skyrmion - like object with equal number of positively and negatively charged single skyrmions .
however , in practice , namely the boson / hole doping is likely to be a physically clear driving force for a nucleation of a single , or multi - center skyrmion - like self - organized collective mode in the form of multi - center charged topological defect which may be ( not strictly correctly ) referred to as multi - skyrmion system akin in a quantum hall ferromagnetic state of a two - dimensional electron gas.@xcite in such a case , we may characterize an individual skyrmion by its position ( i.e. , the center of skyrmionic texture ) , its size ( i.e. , the radius of domain wall ) , and the orientation of the in - plane components of pseudo - spin ( u(1 ) degree of freedom ) .
an isolated skyrmion is described by the inhomogeneous distribution of the co parameter , or staggered boson density @xmath57 , order parameter @xmath75 characterizing the deviation from the half - filling , and @xmath46 that corresponds only to the modulus of the superfluid order parameter .
it seems likely that for a light doping any doped particle ( boson / holes ) results in a nucleation of a new single - skyrmion state , hence its density changes gradually with particle doping .
therefore , as long as the separation between skyrmionic centers is sufficiently large so that the inter - skyrmion interaction is negligible , the energy of the system per particle remains almost constant .
this means that the chemical potential of a boson or hole remains unchanged with doping and hence apparently remains fixed .
the multi - skyrmionic system in contrast with uniform ones can form the structures with inhomogeneous long - range ordering of the modulus of the superfluid order parameter accompanied by the non - ordered global phases of single skyrmions .
such a situation resembles in part that of granular superconductivity .    in the long - wavelength limit
the off - diagonal ordering can be described by an effective hamiltonian in terms of u(1 ) ( phase ) degree of freedom associated with each skyrmion .
such a hamiltonian contains a repulsive , long - range coulomb part and a short - range contribution related to the phase degree of freedom .
the latter term can be written out in the standard for the @xmath76 model form of a so - called josephson coupling @xmath77 where @xmath78 are global phases for skyrmions centered at points @xmath79 , respectively , @xmath80 josephson coupling parameter . namely the josephson coupling gives rise to the long - range ordering of the phase of the superfluid order parameter in a multi - center skyrmion .
such a hamiltonian represents a starting point for the analysis of disordered superconductors , granular superconductivity , insulator - superconductor transition with @xmath81 array of superconducting islands with phases @xmath82 . calculating the phase - dependent part of skyrmion - skyrmion interaction timm _
et al._@xcite arrived at @xmath83 sign of @xmath80 that favors antiparallel alignment of the u(1 ) pseudospins .
in other words , two skyrmions are believed to form a peculiar josephson @xmath84 micro - junction .
there are a number of interesting implications that follow directly from this result:@xcite the spontaneous breaking of time - reversal symmetry with non - zero supercurrents and magnetic fluxes in the ground state , long - time tails in the dynamics of the system , unconventional aharonov - bohm period @xmath85 , negative magnetoresistence .    to account for coulomb interaction and allow for quantum corrections we should introduce into effective hamiltonian the charging energy @xcite @xmath86 where @xmath87 is a boson number operator for bosons bound in @xmath32-th skyrmion ; it is quantum - mechanically conjugated to @xmath51 : @xmath88 , @xmath89 stands for the capacitance matrix , @xmath90 for bosonic charge .
the classical ground state energy of skyrmion lattice for all reasonable two - dimensional lattice structures was minimized by timm _
_ taking the u(1 ) order into account .
@xcite besides the expected triangular lattice with frustrated antiferromagnetic u(1 ) order and square lattice with nel u(1 ) order the authors have also obtained ground - state energies for all 2d bravais lattices .
such a system appears to reveal a tremendously rich quantum - critical structure . in the absence of disorder ,
the @xmath48 phase diagram of the multi - skyrmion system implies either triangular , or square crystalline arrangements ( skyrmion , or bubble crystal ) with possible melting transition to a skyrmion ( bubble ) liquid .
the melting of skyrmion lattice is successfully described applying the berezinsky - kosterlitz - thouless ( bkt ) theory to dislocations and disclinations of the lattice and proceeds in two steps .
the first implies the transition to a liquid - crystal phase with short - range translational order , the second does the transition to isotropic liquid .
disorder pins the skyrmion lattice and also causes the crystalline order to have a finite correlation length .
for such a system provided the skyrmion positions fixed at all temperatures , the long - wave - length physics would be described by an @xmath76 model with expectable bkt transition and gapless @xmath76 spin - wave mode .
the classical phase diagram of skyrmion lattice is quite rich .
depending on the relation between coulomb and josephson coupling , and density one may arrive at different lattice structures with continuous or first order phase transitions . as regards the superfluid properties the skyrmionic system reveals unconventional behavior with two critical temperatures @xmath91 and @xmath92 , @xmath93 being the temperature of the ordering of the modulus and @xmath94 that of the phase of order parameter @xmath95 .
the low temperature physics in skyrmion crystals is governed by an interplay of two bkt transitions , for the u(1 ) phase and positional degrees of freedom , respectively .
@xcite dislocations in most skyrmion lattice types lead to a mismatch in the u(1 ) degree of freedom , which makes the dislocations bind fractional vortices and lead to a coupling of translational and phase excitations .
both bkt temperatures either coincide ( square lattice ) or the melting one is higher ( triangular lattice).@xcite quantum fluctuations can substantially affect these results .
quantum melting can destroy u(1 ) order at sufficiently low densities where the josephson coupling becomes exponentially small .
similar situation is expected to take place in the vicinity of structural transitions in skyrmion crystal . with increasing the skyrmion density
the quantum effects result in a significant lowering of the melting temperature as compared with classical square - root dependence .
the resulting melting temperature can reveal an oscilating behavior as a function of particle density with zeros at the critical ( magic ) densities associated with structural phase transitions .    in terms of our model
, the positional order corresponds to an incommensurate charge density wave , while the u(1 ) order does to a superconductivity . in other words , we arrive at a subtle interplay between two orders .
the superconducting state evolves from a charge order with @xmath96 , where @xmath97 is the temperature of a melting transition which could be termed as a temperature of the opening of the insulating gap ( pseudo - gap ! ? ) .
the normal modes of a dilute skyrmion system ( multi - center skyrmion ) include the pseudo - spin waves propagating in - between the skyrmions , the positional fluctuations , or phonon modes of the skyrmions which are gapless in a pure system , but gapped when the lattice is pinned , and , finally , fluctuations in the skyrmionic in - plane orientation and size .
the latter two types of fluctuation are intimitely connected , since the @xmath23-component of spin and orientation are conjugate coordinates because of commutation relations of quantum angular momentum operators .
so , rotating a skyrmion changes its size .
the orientational fluctuations of the multi - skyrmion system are governed by the gapless @xmath76 model.@xcite the relevant model description is most familiar as an effective theory of the josephson junction array .
an important feature of the model is that it displays a quantum - critical point .
the low - energy collective excitations of skyrmion liquid includes an usual longitudinal acoustic phonon branch .
the liquid crystal phases differ from the isotropic liquid in that they have massive topological excitations , _
i.e. _ , the disclinations .
one should note that the liquids do not support transverse modes , these could survive in a liquid state only as overdamped modes . so that it is reasonable to assume that solidification of the skyrmion lattice would be accompanied by a stabilization of transverse modes with its sharpening below melting transition .
in other words an instability of transverse phonon modes signals the onset of melting .
a generic property of the positionally ordered skyrmion configuration is the sliding mode which is usually pinned by the disorder .
the depinning of sliding mode(s ) can be detected in a low - frequency and low - temperature optical response .
our speculations as regards the topological phase separation and long wave - length physics of 2d hc - bh model are summarized in a qualitative @xmath98-t ( @xmath99 ) phase diagram in fig.[fig3 ] .
first of all , the phase diagram implies a scenario in which the topological phase separation state evolves with a minimal doping , though it is worth noting that the ideal checkerboard co phase can not be a nominal @xmath100 limit of any topologically phase separated phase .    despite the qualitative character of phase diagram we took into account some quantitative results of quantum monte - carlo calculations for the 2d hc - bh model with nearest neighbor coulomb repulsion.@xcite first
, it concerns the doping range of additive , or near - linear concentration behavior of bose - condensate density @xmath101 and co structure factor @xmath102 .
the superfluid density increases approximately linearly with doping except for the most overdoped point @xmath103 , where it turns down again .
diagonal long - range charge order is characterized by the equal time structure factor at the ordering vector @xmath104 , @xmath105 in the presence of long - range order @xmath102 will diverge with the system size @xmath106 for a given ordering momentum @xmath107 , which characterizes the ordered phase . for the checkerboard order @xmath108 . in the concentration range @xmath109 , where both quantities vary linearly with @xmath110,@xcite we may approximate the topological defect to be a system of @xmath111 interacting single - charged skyrmions .
it should be noted that the both co and bs order parameters coexist in a rather narrow doping concentration interval : @xmath112 . beyond this `` overdoping '' region
we deal with an inhomogeneous boson liquid which pseudo - spin picture implies a strongly frustrated singlet - triplet system that resembles the spin glass , and maybe termed as a dynamical `` singlet - triplet '' pseudo - spin glass .
such a system can be characterized by a dynamical short - range diagonal and off - diagonal ordering with a wide distribution of respective correlation lengths and relaxation rates .
interestingly , that in frames of a continuous model this phase is still described to be a strongly correlated multi - center topological defect .
however , such a model fails obviously to describe the real system where the inter - center spacing is comparable with the lattice parameter .
the temperature evolution of hc - bh system with large inter - site boson repulsion implies the highest critical temperature @xmath113 separating the high - temperature non - ordered no phase ( boson liquid ) and a low - temperature quasi - co phase with a disordered system of skyrmions .
the next critical temperature @xmath114 points to the first order phase transition with a formation of inhomogeneous bose condensate in a single skyrmion with the vortex - like texture of the quasi - local order parameter @xmath115 . in other words , it is a temperature of the topological co+bs phase separation ( tps ) . in frames of our scenario
such tps state emerges with a minimal doping , and @xmath116 is likely to be nearly constant in a linear doping regime ( @xmath117 ) .
the low - temperature part of the phase diagram , which is schematically shown in fig . [ fig3 ] , describes the gradual skyrmionic solidification and may include the isotropic skyrmion(bubble ) liquid phase , the skyrmion liquid - crystal phase , the incommensurate skyrmion crystal phase , and u(1 ) ordering , or superfluid phase , restricted by the temperatures of the proper bkt transitions . for a small deviation from half - filling ( `` under - doping '' ) the temperatures of bubble crystallization / melting are seemingly to obey the square root concentration law @xmath118 with a strongly developed quantum melting effect when approaching an `` overdoped '' regime or concentrations limiting the linear regime .
the superfluid phase in fig.[fig3 ] is arbitrarily chosen to lie inside the skyrmion solid phase .
one should emphasize the specific role played by quantum fluctuations : these lead to the melting of the bubble crystal at high densities and orientational disordering @xcite at low densities .
both effects are of primary importance in the overdoped and heavily underdoped regions of the phase diagram , respectively
. moreover , quantum melting effect may strongly affect the phase diagram near the `` magic '' doping level , where the skyrmion lattice undergoes the structural phase transition . for illustration , in the right hand side panel in fig .
[ fig3 ] we present the possible phase diagram for a hc - bh system with a quantum melting effect near the `` magic '' doping level @xmath119 . in all the cases , the critical doping level for the superfluid formation is determined by the magnitude of the josephson coupling constant .
it is worth noting that the bubble crystallization is accompanied by different ( pseudo)gap effects .
our interpretation of the phase transition at @xmath120 differs from that in ref .. this temperature is governed by the magnitude of transfer integral , and believed to be a temperature of the emergence of the nonzero magnitude of the modulus of the superconducting order parameter rather than a critical temperature for an insulator to superconductor transition as it is stated in ref .. this conclusion seems to be a result of finite size effects and boundary conditions in qmc calculations despite the authors make use of the most efficient qmc strategy .
such problems seem to be typical for a finite - size simulation of many phase transitions .
in addition , we should note that mc calculations need substantial increase in lattice size to reproduce quantitatively the long - wavelength physics because the size of skyrmion and skyrmion lattice parameter appear to be new characteristic lengths .
the close inspection of the phase diagram in fig.[fig3 ] , where @xmath121 , does not provide the optimistic expectations as regards the high - temperature superconductivity in 2d hc - bh systems , even if the magnitude of the local boson transfer integral were as large as @xmath122 k. nevertheless , the attractively large temperatures @xmath120 of the topological phase separation with the emergence of nonzero local condensate density engender different reasonable speculations as regards its practical realization .
the making use of the mean - field approximation together with simplified classical continuous models can hardly provide the quantitative description of a quantum lattice boson system . both quantum effects and
the discreteness of skyrmion texture can result in substantial deviations from the predictions of classical model .
the continuous model is relevant for discrete lattices only if we deal with long - wave length inhomogeneities when their size is much bigger than the lattice spacing . in the discrete lattice
the very notion of topological excitation seems to be inconsistent . at the same time , the discreteness of the lattice itself does not prohibit from considering the nanoscale ( pseudo)spin textures whose topology and spin arrangement is that of a skyrmion.@xcite naturally , the typical continuous models fail to describe properly such short - wave length nanoscopic inhomogeneities .
hereafter , we discuss a simple model which seems to catch the main effects of discreteness and quantization .      what is the lattice counterpart of the small skyrmion - like bubble domain
?        fig .
[ fig4 ] presents the schematic view of such a smallest skyrmion - like bubble domain in a checkerboard co phase for 2d square lattice with an effective size of 4 lattice spacings .
the domain wall is believed to include as a minimum eight sites forming a ring - shaped system of four dimers each composed of two sites .
there are two types of such domains which differ by a rotation by @xmath123 .
the formation of such a center seems not to require a big energy .
indeed , the simple estimate of a change in potential energy : @xmath124 points to a near cancelation of @xmath0 and @xmath1 contributions . the remarkable feature of the domain is in a rather small magnitude both of scalar potential and electric field inside the 8-site domain wall ( see fig .
[ fig4 ] ) .
the flip - energy for the dimer dipole moment is estimated to be @xmath125 , compared with that for bare checkerboard co phase ( @xmath126 ) that implies rather subtle competition between @xmath0 and @xmath1 couplings .
however , the difference @xmath127 is believed to be always positive and large . in other words , the flip - energy for
the dimer dipole may be relatively small .
respective dipole fluctuations would result in an effective screening of electrostatic repulsion energy thus providing a stabilization of a bubble - like defect .
an important additional mechanism of the domain stabilization in the hc - bh systems with local bosons composed from electron pairs may arise from the electron and lattice polarization effects.@xcite as is well known , the respective energies are comparable in value with the intersite coulomb interaction . probably , namely both these effects might strongly contribute to the domain stabilization energy . on the other hand ,
the bubble geometry implies the formation of the electrostatic potential well inside the domain wall both for positive and negative charges .
it means that the doping into a domain wall stabilizes the domain configuration .
such a doping maybe energy costless , while the energy cost of the pseudo - spin - flip in a checkerboard co phase is rather high : @xmath128 , if the neutralizing background is taken into account .      with taking into account the kinetic energy ( quantum tunneling ) only inside the 8-site domain wall we shall consider it as a quantum system in an external electrostatic potential field assuming a rigid checkerboard co ordering overall beyond the domain wall .
further , taking into account the reasonable relation : @xmath129 , we shall consider the domain wall to be a system of four dimers , or pairs of nearest neighboring sites forming a quantum cluster .
transfer integral .
the arrows mark the dipole - allowed transitions ( see text for details).,width=321 ]    we anticipate the model though being simplified would be very instructive analyzing the role played both by quantum and discreteness effects .
let first address a simple two - site , or a dimer system .
the effective pseudo - spin - hamiltonian for such a cluster , or a single dimer hamiltonian can be written as follows @xmath130\,+d { \hat s}_{z}^2\ ,
-\,h_s\,{\hat s}_{z}-\,h_t\,{\hat t}_{z } , \label{h}\ ] ] where @xmath131 is the total pseudo - spin momentum , @xmath132 is an operator that changes pseudo - spin multiplicity , @xmath133 , @xmath134 , @xmath135 .
the first two terms in this effective pseudo - spin hamiltonian describe the intra - dimer interactions , while the last two ones describe its coupling with the off - domain - wall surroundings .
it is worth noting that the condition @xmath134 means that the effective magnetic field produced by this surroundings turns into zero that provides the particle - hole symmetry of the dimer physics , in particular , for the domain wall doping .
actually , we arrive at an effective singlet - triplet model.@xcite both @xmath136 and @xmath137 operators have a rather simple physical sense : the former corresponds to the total `` quantum '' charge of the dimer , while the latter does to the total `` quantum '' dipole moment . strictly speaking , the diagonal order parameters @xmath138 and @xmath139 describe the charge and dipole density , while the off - diagonal order parameters @xmath140 and @xmath141 describe the corresponding phase ordering .
it is of primary importance to note that these order parameters are not independent because the respective operators obey the kinematic constraint : @xmath142 stemmed from a simple spin algebra .
it should be noted that there are two operators : @xmath143 and @xmath144 $ ] that change the pseudo - spin multiplicity with their matrices being symmetric and antisymmetric , respectively : @xmath145 @xmath146 the hamiltonian ( [ h ] ) points to the competition of s- and t - orders in the ring - shaped domain wall .
fig.[fig5 ] shows a step - by - step formation of the energy spectrum of such a two - site cluster , or dimer .
the arrows mark the dipole - allowed transitions that could be revealed in optical spectra .
the tunnel states @xmath147 and @xmath148 describing purely `` half - filled '' dimer states are mixed due to the electric field with the mixing level governed by the ratio @xmath149 .
if @xmath150 , we arrive at classical nel - like `` up - and - down '' dimer state .
hereafter , to describe our singlet - triplet quantum pseudo - spin system we start with trial functions @xmath151 where the spin functions @xmath152 in cartesian basis are used : @xmath153 and @xmath154 .
the conventional spin operator is represented on this basis by a simple matrix : @xmath155 and for the order parameters one easily obtains : @xmath156 , \langle\{\hat{s}_i\hat{s}_j\}\rangle=2(\delta_{ij}-a_ia_j - b_ib_j ) \label{med}\ ] ] given the normalization constraint @xmath157 .
thus , for the case of spin-1 system the order parameters are determined by two classical vectors ( two real components of one complex vector @xmath158 from ( [ wf1 ] ) ) .
the two vectors are coupled , so the minimal number of dynamic variables describing the @xmath159 spin system appears to be equal to four .        for the dimer as the singlet - triplet center
we have additional unconventional t , j - order parameters . for the respective averages we can easily obtain @xmath160
given the real value of @xmath161 parameter we arrive at a simple form:@xmath162 , and @xmath163 .
now we may proceed with the inter - dimer coupling in the domain wall . introducing the symmetrized superpositions of s and t vectors , we can write the pseudo - spin hamiltonian of the inter - dimer interaction as follows @xmath164 @xmath165 + \frac{1}{2}(v_{nnn}-t_{nnn})[s_{z}^{2}(a_{1g})- s_{z}^{2}(b_{1 g } ) + \ ] ] @xmath166,\ ] ] where @xmath167 @xmath168 @xmath169 @xmath170 @xmath171 @xmath172 @xmath173 @xmath174 ( see fig.[fig4 ] ) .
the chemical potential @xmath175 for the system of four dimers in a domain wall is determined now by the condition : @xmath176 for undoped and singly - doped ( boson / hole ) domain wall , respectively .
the respective mean values we can address to be order parameters that describe the subtle structure of a domain wall : as regards the diagonal order @xmath177 specifies the full charge ; @xmath178 and @xmath179 do the electric dipole moments , @xmath180 the component of the quadrupole momentum , and one might introduce the higher order multipole moments . as regards the off - diagonal order we should , in general , proceed with the three types ( @xmath181 , or @xmath182 ) of the s - order , and three types ( @xmath183 , or @xmath184 ) of the t - order , which can be defined as follows : @xmath185 noteworthy to mention the kinematic constraint that couples different order parameters . in the absense of an external magnetic field the energy
does not depend on @xmath186 vector , it is restricted only to lie in @xmath35-plane . in other words , we deal with the uncertainty of j - order parameter .
the mean - field domain - wall ground state corresponds to the coexistent @xmath187 modes given the negative sign of the transfer integral @xmath188 or @xmath189 modes given the positive sign of the transfer integral @xmath188 .
the relation between s- and t - mode weight is specified by the relationship between the @xmath0 transfer integral @xmath190 and electric field @xmath191 .
if we assume @xmath192 , then given @xmath193 the dimers have purely @xmath159 ground state , and we arrive at s - type off - diagonal order . in contrast , given @xmath194 we arrive at t - type off - diagonal order . in general , neglecting the st - mixing term leads to two independent phase order parameters . generally speaking , the interference dipole - dipole st - term would result in st - mixing accompanied by the constraint on the phase order parameters with an appearance of a noncollinearity effect .
thus we arrive at the conclusion that the symmetry of the order parameter distribution in the domain wall would be specified only by the sign of the transfer integral .
in addition , we see that the problem of the order parameter associated with our bubble domain is much more complicated than in conventional bcs - like approach @xcite due to its multicomponent nature .
moreover , we deal with a system with different symmetry of low - lying excited states and competing order parameters that implies their possible ambiguous manifestation in either properties .
the low - symmetry crystalline electric fields or crystal distortions would result in a mixing of the order parameters with different symmetry .
the bubble domain ( fig.[fig4 ] ) yields a simple and instructive toy model to describe such effects .
it is worth noting that the flux quantization effects for the bubble domain , in particular , the localizing effect of the magnetic field on the moving bosons , are expected to be observed only for rather large fields . indeed
, nanoscopic atomic systems such as a square plaquette with a size of around 10 @xmath195 require a field of @xmath196 tesla for a half flux quantum per plaquette .    above we have concerned a simple one - center topological defect .
however , in general one has to make use of more complicated topological excitation like a multi - center bp skyrmion.@xcite the question arises , whether such an entity may be described as a system of weakly coupled individual one - center defects ? in practice , it seems to be a rather reasonable approach only for slight deviation from half - filling , when the mean separation between doped particles is much larger than the effective domain size .
the bubble domain may be addressed to be well isolated only if its environment does not contain another domain(s ) which could be involved to a `` dangerous '' nearest - neighbor inter - dimer coupling .
each domain in fig.[fig4 ] has @xmath197 of such neighbors .
hence the concentration of well isolated domains can be written as follows : @xmath198 where @xmath23 is a number of `` dangerous '' neighbors , @xmath199 the boson concentration .
the @xmath200 maximum is reached at @xmath201 . in our case @xmath202 , or @xmath203 . with increasing doping the deviation of @xmath200 from the linear law rises . on the other hand , knowing the effective domain area @xmath204 we can roughly estimate the limiting concentration of the single - domain model description to be @xmath205 , @xmath206 .
the inter - domain coupling includes both long - range ss- and short - range tt - terms .
we should emphasize a strong anisotropy of this coupling .
for nearest neighbors along [ 1,0 ] and [ 0,1 ] directions namely the short - range ( dipole - dipole ) tt - coupling would result in an effective suppression of the main s ordering ( see fig.[fig4 ] ) , in contrast with [ 11 ] direction .
thus we may conclude that the above mentioned simple continuous model for josephson junction arrays should be strongly modified to describe the multi - domain configurations in 2d quantum hc - bh lattices , both as regards the subtle internal domain structure , the competition of different order parameter , and the anisotropy of josephson coupling .
in conclusion , the boson / hole doping of the hard - core boson system away from half - filling is assumed to be a driving force for a nucleation of a multi - center skyrmion - like self - organized collective mode that resembles a system of co bubble domains with a bose superfluid and extra bosons both confined in domain walls .
such a _ topological _ co+bs _ phase separation _ , rather than an uniform mixed co+bs supersolid phase , is believed to describe the evolution of hc - bh model away from half - filling .
starting from the classical model we predict the properties of the respective quantum system . in frames of our scenario
we may anticipate for the hc - bh model the emergence of an inhomogeneous bs condensate for superhigh temperatures @xmath207 , and 3d superconductivity for rather high temperatures @xmath208 .
the system is believed to reveal many properties typical for granular superconductors , cdw materials , wigner crystals , and multi - skyrmion system akin in a quantum hall ferromagnetic state of a 2d electron gas .
topological inhomogeneity is believed to be a generic property of 2d hard - core boson systems away from half - filling .
such a behavior represents a boson counterpart of the so - called _ topological doping _ being a general feature of mott - insulator or 2d fermion hubbard model.@xcite    despite all shortcomings , mfa and continuous approximation are expected to provide a physically clear semiquantitative picture of rather complex transformations taking place in bare co system with doping , and can be instructive as a starting point to analyze possible scenarios .
first of all , the mfa analysis allow us to consider the antiphase domain wall in co phase to be a very efficient ring - shaped potential well for the localization of a single extra boson ( hole ) thus forming a novel type of a topological defect with a single - charged domain wall .
such a defect can be addressed as a charged skyrmion - like quasiparticle which energy can be approximated by its classical value for co bubble domain .
it is of great importance to note that domain wall simultaneously represents a ring - shaped reservoir for bose superfluid .
unfortunately , we have no experience to deal with multi - center skyrmions as regards its structure , energetics , and stability .
it should be noted that such a texture with strongly polarizable centers is believed to provide an effective screening of long - range boson - boson repulsion thus resulting in an additional self - stabilization .
nucleation of topological phase is likely to proceed in the way typical for the first order phase transitions .
the role played by quantum effects and lattice discreteness has been illustrated in frames of the simplest nanoscopic counterpart of the bubble domain in a checkerboard co phase of 2d hc - bh square lattice .
it is shown that the relative magnitude and symmetry of multi - component order parameter are mainly determined by the sign of the @xmath0 and @xmath1 transfer integrals .
the topologically inhomogeneous phase of the hc - bh system away from the half - filling can exhibit the signatures both of @xmath2 , and @xmath3 symmetry of the off - diagonal order .
the model allows us to study the subtle microscopic details of the order parameter distribution including its symmetry in a real rather than momentum space , though the problem of the structure and stability of nanoscale domain configurations remains to be solved .
the present paper establishes only the framework for analyzing the subtleties of the phase separation in a lattice hc - bh model away from half - filling .
much work remains to be done both in a macroscopic and microscopic approaches .
i acknowledge stimulating discussions with c. timm , s .- l .
drechsler , t. mishonov , and the support by smwk grant , intas grant no .
01 - 0654 , crdf grant no .
rec-005 , rme grant no .
e 02 - 3.4 - 392 and no .
ur.01.01.062 , rfbr grant no .
04 - 02 - 96077 .
i would like to thank leibniz - institut fr festkrper- und werkstoffforschung dresden where part of this work was made for hospitality .
special thanks are to referees for their remarks improving the manuscript .
it is worth noting that in frames of bcs - like scenario the symmetry of the order parameter is strictly defined in a momentum space albeit the discussion of different experimental data has usually been performed with a real - space distributed order parameter .


ever since landau@xcite suggested that the electron can be trapped by the deformable lattice , strongly coupled electron - phonon systems have been the subject of intensive examination . besides the investigations of those systems in which the lattice is coupled to the whole electronic band , there has been significant interest in the physics of a single polaron , in which the electron and the associated lattice deformation form a quasi - particle , spatially and spectrally decoupled from the rest of the system .
lattice degrees of freedom make even a single polaron problem a many - body one .
the analytical and numerical examinations of most electron - phonon models are thus difficult . for this reason
, even the simple holstein  model@xcite ( suggested in 1959 ) is still being investigated in recent works .
various methods have been proposed in order to calculate the polaron ground state of the holstein  model .
almost exact results ( except for the adiabatic limit ) have been obtained with the quantum monte carlo calculations,@xcite the global - local method,@xcite the density matrix renormalization group method,@xcite and some exact diagonalization methods.@xcite    the main goal of this paper is to determine the elements important for the qualitative description of the polarons in the whole range of electron - phonon coupling .
some of them , although already known , are found to be better understood when supplemented with additional details .
the ground state of the holstein  system changes from the delocalized polaron state , in which the electron is nearly free , to the small and self - trapped polaron state , as the electron - phonon coupling @xmath0 increases .
these two opposite limits are usually identified as the weak and strong coupling regime , respectively . in the weak coupling regime
the influence of the small lattice deformation on polaron dynamics is very small , which makes the energies of polaron and electron hopping to neighboring sites similar , @xmath1 .
the exact ground state is an eigenstate of the system momentum , regardless of the coupling.@xcite therefore the polaron ground state is delocalized for all parameters .
however , in the strong coupling regime @xmath2 becomes negligible , leading to _ self - trapped _
polaron states . according to ref .
, the dynamics of the small self - trapped polaron can be separated into two time scales . on the short time scale ,
the lattice deformation is centered at some lattice site , and the electron can virtually hop among the neighboring lattice sites . only after a certain number of such events ( of the order @xmath3 ) , the polaron as a whole tunnels to a new central site .
the localized polaron states thus may lead to a very accurate estimation of the polaron ground state energy in the strong coupling regime .
accordingly , for the self - trapped polarons , an exact diagonalization method of calculating the localized polaron states , rather than the translationally invariant ones , can be used . by comparison to the results of other methods
, it is shown that this approach introduces only minor errors in the ground state energy of the self - trapped polaron .
moreover , the localized polaron functions permit , unlike the translationally invariant ones , a separate analysis of the electron and phonon properties of the polaron . the local electron density , the mean lattice deformation , and the on - site zero point motions of the self - trapped polarons can be calculated in this way .    in the crossover ( intermediate ) regime , which is between the weak and the strong coupling regime
, no known perturbation calculation converges , which complicates the discussion of the polaron nature .
although it has been proved that the change of the polaron ground state with @xmath0 is smooth,@xcite the physics of the rapid crossover between two opposite limits of the electron - phonon coupling , in the small interval of @xmath0 s , is not completely clear . in ref .
it has been claimed that the phonon excitation associated with the first excited state of the system is uncorrelated to the electron in the weak coupling regime , while it is confined to the electron in the strong coupling regime .
the transition occurs in the crossover regime in which , in addition , the energy difference between the ground and excited states becomes small . in the present paper ,
the excited polaron states are treated by the new method which uses variational approach in order to define and solve the generalized eigenvalue problem .
even if this method does not converge systematically , it does provide new results concerning the nature of the polaron ground and first excited state .
the holstein  hamiltonian  reads    @xmath4    it describes the tight - binding electrons in the nearest - neighbor approximation , coupled to one branch of dispersionless optical phonons .
@xmath5 is the creation operator for the electron of spin @xmath6 at lattice site @xmath7 , and @xmath8 is the creation operator for the phonon .
@xmath9 is the transfer ( hopping ) integral of the electron .
@xmath10 and @xmath0 are the phonon and the electron - phonon coupling energies , respectively .
the holstein  hamiltonian  depends only on two ratios of relevant energy parameters : @xmath11 and @xmath12 , i.e. , the results will use @xmath10 as the energy unit .
it is often convenient to express the lattice vibrations in terms of the nuclei space and momentum coordinates ,    @xmath13    where @xmath14 and @xmath15 are space and momentum uncertainties of the harmonic oscillator ground state .
@xmath16 is the mass of a nucleus and @xmath17 denotes the spring constant , @xmath18 . by noting that the electron - lattice displacement coupling constant @xmath19 , in @xmath20 , is independent of @xmath16 , an alternative set of holstein  hamiltonian  parameters can be introduced ,    @xmath21    which is convenient for the discussion of the adiabatic regime @xmath22 . here , @xmath23 is the electron effective mass and @xmath24 is the lattice constant .
it is worth noting that @xmath9 and @xmath25 ( the binding energy of the polaron for @xmath26 ) are independent of @xmath16 and thus they are the only parameters relevant in the adiabatic limit .    by using standard conventions for @xmath27 and @xmath28 ,    @xmath29    eq .
( [ h0 ] ) can be rewritten in momentum space and separated in two mutually commuting parts , @xmath30 ,    @xmath31    where @xmath32 , @xmath33 , and @xmath34 .
the part which involves only the @xmath35 phonon mode , @xmath36 , can easily be transformed into the diagonal form , by using the unitary operator @xmath37 ,    @xmath38    where @xmath39 .
any eigenstate of @xmath40 , @xmath41 , has the same mean total lattice deformation , @xmath42 ,    @xmath43    which is independent of @xmath9 .
since the @xmath35 part of the hamiltoniancommutes with @xmath44 , it can be concluded that eq .
( [ xtot ] ) is also valid for all eigenstates of the total holstein  hamiltonian .
the phonon part of these eigenstates can be represented in the form of a direct product of two groups of states , the first one includes @xmath45 phonon modes , while the second includes only the @xmath35 phonon mode .
this is useful because one can always check approximate computations by calculating @xmath46 , or include this property in the computation itself .
this specific property of the @xmath35 mode is not restricted to the particular dimension of the system , nor to the number of electrons .
moreover , it can also be found in some other models in which the electron - phonon coupling consists of the lattice deformation linearly coupled to the local electron density .
a hint in this direction was reported for the first time in ref .  .
the total momentum of the system , @xmath47 , is the sum of the electron and phonon momenta ,    @xmath48    and it commutes with the hamiltonian . in the present treatment ,
only the low energy polaron states ( the low lying states of the system for which the electron and lattice part of the wave function are spatially bound ) , are explicitly calculated . in this case , the total momentum @xmath49 of the system is also the polaron momentum .
the eigenstate computations reported here are based on the variational approach .
still , from the physical and mathematical point of view , there are significant differences among them .
each method is therefore described separately in the present section .
it is important to notice that two qualitatively different kinds of polaron functions are employed , the _ localized _ ones and the _ translationally invariant _ ones .
the primary objective of the methods with localized functions are self - trapped polaron states , i.e. , the strong coupling regime . on the other hand , translationally
invariant states are devised in the first place with the weak coupling and crossover regime in mind .
the methods with a small number of variational parameters are meant to help in understanding the basic properties of the polarons in different regimes of the holstein  model .
again , the methods with a very large number of variational parameters are necessary to obtain accurate results for the polaron states .
let us start with a simple localized polaron wave function formed as a product of the electron and the lattice part , centered at the lattice site @xmath51 ,    @xmath52    here , @xmath53 is the normalized electron function at site @xmath54 , @xmath55 , while @xmath56 denotes the coherent state operator acting on site @xmath57 , with a complex amplitude , @xmath58 ,    @xmath59    it is easy to see that operator @xmath56 shifts the space and momentum coordinates of lattice vibration at a site @xmath57 by @xmath60 and @xmath61 , respectively ,    @xmath62 } e^{[-i2\re(\xi_m)x_0\hat p_{j+m}/\hbar]}\\ & \times&e^{[i\re(\xi_m)\im(\xi_m)]}\;.\end{aligned}\ ] ]    the variational energy of the state ( [ iafunction ] ) , @xmath63 , is independent of @xmath51 , and is given by    @xmath64    the minimization of the energy with respect to @xmath65 establishes a simple relationship between the lattice mean deformation and electron density @xmath66 ,    @xmath67    so that only the equation for @xmath53 has to be solved .
the well known approximate solution to this problem is the large holsteinpolaron@xcite valid in the long wave limit ,    @xmath68    the numerical scheme suggested in ref .  , and denoted here by @xmath50 ( @xmath50 for localized ) , is not restricted to long waves , and is used here in order to obtain the exact minimum of eq .
( [ iaenergy ] ) .
the energy , henceforth referred to as @xmath69 , depends only on two relevant hamiltonianparameters , @xmath9 and @xmath70 , and therefore @xmath69 and @xmath71 are both independent of @xmath16 . for this reason eq.([locking ] ) is sometimes referred to as the adiabatic locking of electron and lattice coordinates .
it is noteworthy that the lattice part of the function @xmath72 is a simple product of coherent states with real amplitudes .
consequently , the state of the lattice at some site is defined by the ground state of the displaced harmonic oscillator .
the mean lattice deformation @xmath71 corresponds to the equilibrium position of that oscillator , while the lattice zero point motion is approximated in eq .
( [ iafunction ] ) by that of the free lattice .
next , we shall study a translationally invariant solution composed of a linear superposition of the localized states ( [ iafunction ] ) ,    @xmath75    @xmath76 describes the polaron state with the momentum  @xmath49 .
a similar type of function was first proposed by toyozawa.@xcite in the present work , @xmath77 is used so that the mean total deformation of function ( [ tivfunction ] ) satisfies eq .
( [ xtot ] ) ,    @xmath78    the expectation value of the polaron energy , @xmath79 , may be written in terms of @xmath80 ,    @xmath81    @xmath82 and @xmath83 are given in the appendix . a simple method for calculating the minimum of the energy @xmath79
has not yet been proposed , but accurate results have been obtained in ref .   by using the toyozawa method , which includes a very large number of variational parameters .
some additional approximations may be found in refs . and .
the approximation used here simplifies the general expression in eq.([tivfunction ] ) by introducing the exponential form for functions @xmath53 and @xmath84,@xcite    @xmath85    equation  ( [ ve ] ) defines a polaron function @xmath86 which is completely determined by two parameters , @xmath87 and @xmath88 .    in what follows ,
two different approaches are presented .
the first , denoted by the index @xmath73 ( @xmath73 for translational ) , treats @xmath87 and @xmath88 as the variational parameters for which the energy minimum @xmath89 has to be found , and its corresponding polaron function is @xmath90 .
the @xmath73 method gives good results in the weak and the strong coupling regime .
namely , in both of these limits , the function @xmath90 becomes similar to the polaron function obtained by the appropriate perturbative calculations .
in the second approach , presented here for the first time , the variational method is used to define a generalized eigenvalue problem as follows .
the polaron wave function , denoted by the index @xmath74 ( @xmath74 for combination of translational functions ) , is rewritten as a linear combination of functions @xmath91 ,    @xmath92    it is understood here that the functions @xmath93 form a set of @xmath94 generally nonorthogonal functions , defined by @xmath94 different pairs of parameters @xmath95 .
again , the coefficients @xmath96 should be determined from the requirement that the expectation value of the energy ,    @xmath97    is minimal ,    @xmath98    or ,    @xmath99    the solution of this generalized eigenvalue problem is a set of @xmath94 orthogonal polaron functions , @xmath100 , with corresponding energies @xmath101 .
the ground state energy is @xmath102 .
one may always include the function @xmath103 in the sum ( [ lvfunction ] ) in order to ensure that the energy @xmath102 is the same or better than @xmath89 , the energy computed by the @xmath73 method .
moreover , by paying some further attention to the starting set of functions @xmath91 in eq.([lvfunction ] ) at the outset , one is able to investigate the first excited state @xmath104 of the system , when this state corresponds to an excited polaron .
the best results for the @xmath74 method are obtained when the number @xmath94 of @xmath91 functions in eq .
( [ tivfunction ] ) changes with the hamiltonian  parameters .
the special case where the @xmath74 method is used with constant @xmath105 is denoted by the index @xmath106 .
finally , this paper presents the results of two numerical exact diagonalization methods.@xcite in order to compute the low energy polaron states , one approximates the infinite dimensional hamiltonian  matrix with a finite one , and proceeds with the exact diagonalization of this matrix .
the lowest eigenvalue and eigenvector in such a reduced hilbert  space correspond to the polaron energy and wave function , respectively .
for a large sparse matrix , the energy and the wave function can be calculated very accurately , by using an appropriate numerical scheme , in the present case the lanczos algorithm .
the two exact diagonalization methods used here differ in the choice of the basis of the hilbert  space . in the first method,@xcite denoted by @xmath107 ( @xmath109 stands for exact diagonalization and @xmath73 for translational ) ,
the general orthonormal state is given by    @xmath110    which describes an eigenstate of the system with momentum @xmath49 .
@xmath111 is the number of phonons at the @xmath112th lattice site away from the electron .
for example , at the site @xmath51 , and at the nearest neighbor lattice sites left and right from it , there are @xmath113 , @xmath114 , and @xmath115 phonons , respectively . the hamiltonian  ( [ h0 ] ) does not mix states ( [ tbstates ] ) with different momenta .
therefore , the polaron function obtained by the @xmath107 method has the same @xmath49 momentum as the basis states .
the current implementation of the @xmath107 method is highly accurate , and from a practical point of view may be treated as exact.@xcite for this reason , the @xmath107 results can be used to determine the numerical errors present in the other methods .
the minimal number of states of the reduced basis necessary to obtain accurate results depends on the hamiltonian  parameters .
this number for the @xmath107 method increases very rapidly for @xmath116 , @xmath9 , which prevents its use for both large @xmath0 and @xmath9 .    in the second exact diagonalization method , denoted by @xmath108 ( @xmath109 stands for exact diagonalization and @xmath50 for localized ) ,
the general orthonormal state of the chosen basis is more complicated than for the @xmath107 method ,    @xmath117\ ; \mid n_0 , n_{-1 } , n_1, ...
,n_m, ... \rangle_j\;.\label{lbstates}\end{aligned}\ ] ]    here the @xmath118 and @xmath112 indices are given with respect to the center of polaron , which is placed at the site @xmath51 .
thus @xmath119 creates an electron at the @xmath118th site from the polaron center at @xmath51 , @xmath111 is the number of extra phonons at the @xmath112th site away of the polaron center , when the lattice is already distorted by the coherent state operators @xmath56 , i.e. , @xmath120 .
the @xmath108 method calculates localized polaron wave functions .
namely , eq .  ( [ lbstates ] ) describes a localized state , with the polaron center at site @xmath51 kept constant .
the @xmath108 method is therefore accurate only in the strong coupling regime , in which the effects of polaron delocalization are negligible ( self - trapped polarons ) .    for a given set of hamiltonian  parameters , @xmath84 in eq .
( [ lbstates ] ) are determined by the use of the @xmath50 method , i.e. , by minimizing the energy ( [ iaenergy ] ) , @xmath121 .
if only those states ( [ lbstates ] ) with all @xmath122 are used in calculations , the @xmath108 method gives the same polaron wave function as the @xmath50 method .
the additional states ( [ lbstates ] ) , with @xmath111 phonon excitations , are necessary to obtain the actual equilibrium positions of the lattice in the exact localized polaron state and the zero point motion of the renormalized lattice vibrations .
it is worth noting that the electron and phonon parts of the @xmath108 polaron wave function can not be completely separated , as in the case of the @xmath50 function .    in the case of @xmath108 method ,
the mean lattice deformation of the localized polaron is approximately taken care of by the product of the coherent states operators , @xmath123 , which keeps the necessary number of states ( [ lbstates ] ) , in the @xmath108 method , relatively small . in order to reduce the basis of the hilbert  space , the maximal allowed distance of the electron and phonons from the polaron center
has been limited here by the choice @xmath124 .
the distance @xmath125 has been determined from the condition that @xmath126 if @xmath127 .
the maximal total number of phonons has been kept limited , retaining the states with @xmath128 . as the sum @xmath129 does not include the phonons associated with the coherent state operators @xmath56
, the small value of @xmath129 is not a restriction on the overall amplitude of the lattice deformation .
the accuracy of the results obtained by the @xmath108 method , supplemented by the two abovementioned criteria , depends of course on the values of parameters .    for the purpose of clarity
, it seems appropriate at this point to review briefly the notation @xmath50 , @xmath73 , @xmath74 , @xmath107 , and @xmath108 of all five presented methods .
all methods based on the localized polaron function are denoted by the letter @xmath50 ( @xmath50 and @xmath108 methods ) , whereas all methods based on the translational function are denoted by the letter @xmath73 ( @xmath73 , @xmath74 , and @xmath107 methods ) . the letter @xmath109 denotes an exact diagonalization method ( @xmath107 and @xmath108 methods ) , while a single letter notation ( @xmath50 and @xmath73 ) suggests the simplest form of the method .
as the variational methods of the preceding section themselves , the results of the corresponding calculations may be best understood in terms of the weak and strong coupling limits and the crossover between them .
it has been shown previously,@xcite on the basis of the global - local method for @xmath130 , that the empirical relation    @xmath131    describes well the values of parameters for which the variation of the effective polaron mass with @xmath0 is the fastest .
it will be also argued here that @xmath132 of eq .
( [ romero ] ) describes accurately the crossover from the weak to the strong coupling limit with respect to the nature of the @xmath133 ground state .
the latter changes continuously from the light delocalized state in the weak coupling limit to the heavy self - trapped state in the strong coupling limit , with the anticrossing of the two states at @xmath134 .
the physical content of eq .
( [ romero ] ) is best understood by considering the limits of small and large @xmath9 with respect to @xmath10 , when , respectively , @xmath135 and @xmath136 .
both these conditions were qualitatively explained in ref .  .
in the present section we first discuss the strong and the weak coupling limit , and then devote most of our attention to the crossover between the two .
the nature of the self - trapped polaron may be discussed by analyzing the properties of the two polaron wave functions given by eqs .
( [ iafunction ] ) and ( [ tivfunction ] ) . when the localized polaron functions @xmath80 at different lattice sites are not orthogonal ,    @xmath137    the local properties and delocalization effects of @xmath76 in eq .
( [ tivfunction ] ) are a complex mixture . however , in both the weak and strong coupling regimes , the translational polaron function can be approximately written in terms of orthogonal localized polaron functions . in the weak coupling regime ,
the orthogonality follows from the electron part of the wave function ,    @xmath138    while in the strong coupling regime , it follows from the lattice part ,    @xmath139 } \approx\delta_{0,\delta}\;.\label{phorth}\ ] ]    in eq .
( [ phorth ] ) , @xmath140 is the debye - waller factor .
the condition ( [ phorth ] ) corresponds to the regime of self - trapped polarons .
the negligible contribution of the lattice part to the overlap of any two localized polaron functions at different lattice sites results then in a negligible polaron hopping energy @xmath2 .    one may notice that in the limit @xmath141 the translational form of the polaron function ( [ tivfunction ] ) has no consequences on polaron energy , and the minimal values of the variational energies ( [ iaenergy ] ) and ( [ tivenergy ] ) coincide , i.e. , @xmath82 , @xmath142 in eq .
( [ tivenergy ] ) .
the hopping of the self - trapped polaron occurs only at the time scale which is much larger than the scale relevant for the local interplay between the electron and the lattice deformation .
therefore an accurate description of the local polaron properties may be obtained even if the hopping of the polaron is completely omitted by using only localized polaron functions .
the translational invariance of the polaron may be , however , always restored in the same way as the function ( [ iafunction ] ) is used to obtain ( [ tivfunction ] ) .
the @xmath108 method provides very accurate results for the ground state energy of the self - trapped polaron . except that it neglects the polaron hopping ,
the local polaron function is calculated exactly .
the error of the @xmath108 method may be estimated by using the @xmath107 method , that is , by subtracting the exact energy of zero momentum polaron state @xmath143 from @xmath144 . for the electron - phonon coupling which is greater than the critical electron - phonon coupling in eq.([romero ] ) just by @xmath10 , i.e. , @xmath145 , the maximal error of the @xmath108 method is @xmath146 .
this is shown in fig .
[ fig1 ] , in which two qualitative estimates of @xmath2 in the strong coupling regime are plotted as well .
the first estimate gives @xmath2 as one fourth of the polaron bandwidth computed exactly by the @xmath107 method ,    @xmath147\;,\label{bandw}\ ] ]    while the second estimate , based on the @xmath50 method , multiplies the electron hopping energy with the small debye - waller factor ,    @xmath148}\;.\label{overlap}\ ] ]    the estimation of the debye - waller factor in eq .
( [ overlap ] ) is based on the evaluation of the lattice part of the overlap of the two neighboring @xmath149 functions ( the @xmath50 method gives good results for the mean lattice deformation of the self - trapped polarons ) .
all three curves in fig .
[ fig1 ] are similar functions of @xmath0 , which is not surprising since all plotted quantities are related to the polaron hopping energy @xmath2 in the strong coupling regime .
moreover , it may be seen from fig .
[ fig1 ] that the error of the @xmath108 method is almost equal to one half of the polaron bandwidth ,    @xmath150    one of the advantages of the @xmath108 method is that it permits separate calculations of the electron and the lattice properties , in spite of the fact that the electron and the phonon part of the @xmath108 wave function are not separable .
for instance , for the polaron centered at the origin , the associated mean lattice deformation @xmath151 is given simply by the expectation value of @xmath152 of eq .
( [ zpm ] ) . besides @xmath151 , in the present paper the mean real space uncertainty of the on - site lattice vibration @xmath153 and the product @xmath154 , where @xmath155 is the mean momentum uncertainty of the on - site lattice vibration , are calculated .
figure [ fig2 ] shows the data for two sets of parameter , corresponding to a small self - trapped polaron and a self - trapped polaron extended over few lattice sites , respectively . for the second set of parameter ( large @xmath0 and @xmath9 )
, the @xmath107 method is not accurate enough , and the question arises of whether the @xmath108 polaron state is really a good approximation of the system ground state .
it is difficult to prove that a polaron state , with non - negligible hopping energy @xmath2 and an energy close or lower than @xmath108 energy , does not exist .
this question is currently under investigation .
differences between the results of the @xmath50 and the @xmath108 methods are also analyzed .
it is found that they are more pronounced for small @xmath0 .
we may see from the results in fig .
[ fig2 ] , that the mean lattice deformation differs between these two methods . in the case of the @xmath108 method ,
the mean lattice deformation is more extended and the width of the polaron is slightly larger .
additionally , the electron density of the @xmath108 method remains approximately proportional to the mean lattice deformation ( as in eq .
( [ locking ] ) , valid for the @xmath50 method ) , since    @xmath156/x_n^{\;el}<1\%\;.\ ] ]    as has already been mentioned , the lattice part of the @xmath50 function describes a set of displaced , but unrenormalized , harmonic oscillators , so @xmath157 , @xmath158 , where @xmath159 and @xmath160 are defined in eq .
( [ zpm ] ) .
one may notice from fig .
[ fig2 ] that real space uncertainties of the lattice vibrations on the sites occupied by the polaron are larger in the @xmath108 case than in the @xmath50 case , @xmath161 , but uncertainties of on - site momentum lattice vibrations are smaller , @xmath162
.    it may be concluded , from the first plot in fig .
[ fig2 ] that the electron affects mostly the central site of the small polaron .
the product of uncertainties for this site stays close to the free lattice value , @xmath163 .
therefore the phonon mode at the central site of the small polaron may be treated , in a good approximation , as harmonic . in particular , the renormalized frequency of this mode , @xmath164 , can be roughly estimated from the relation    @xmath165    since @xmath166 , it follows that @xmath167 .
therefore in the strong coupling regime the first excited state should correspond to the excitation of the renormalized phonon mode , rather than to the excitation of the phonon of energy @xmath10 , which is uncorrelated to the polaron .
one may also notice that the energy of the mean lattice deformation is larger for the @xmath108 method than the @xmath50 method ( for the @xmath50 method this energy is minimal ) .
this is compensated for , however , with the lower energy associated with the zero point motion of the @xmath164 phonon mode , which makes the total polaron energy of the @xmath108 method lower .    for the more extended self - trapped polarons ,
the renormalized normal phonon modes are expected to be spread over a number of lattice sites .
consequently , a number of different phonon modes contribute to the lattice displacement at the lattice sites occupied by the polaron .
thus the analysis of the on - site vibrations can not give direct information on the renormalized lattice modes .
nevertheless , from the second plot in fig .
[ fig2 ] one may notice that the product of uncertainties @xmath168 shows a minor deviation from that of the harmonic vibration .
this suggests that the renormalized lattice modes of the extended self - trapped polaron are harmonic as well in a good approximation .      in the weak coupling regime
, the @xmath73 method gives results close to the @xmath107 results .
since the form of the function @xmath73 is quite simple , it will be used in this section as a basis for further discussion . in the weak coupling regime
, the minimum of the variational energy @xmath89 corresponds to small values of the variational parameters @xmath87 and  @xmath169 .
the standard perturbative ground state of the system with momentum @xmath49 in terms of the @xmath73 method polaron function may be written as follows :    @xmath170 the wave function ( [ wc1function ] ) has two parts .
the main part corresponds to the free electron of momentum @xmath49 , and the smaller part , proportional to @xmath169 , corresponds to the electron dressed by one spatially correlated virtual phonon . at the threshold @xmath171 , the energy of such a polaron state intersects with the energy of the system consisting of the zero momentum polaron and one extra phonon with momentum @xmath171 , eq .
( [ wc2function ] ) .
so , for @xmath172 the ground state is achieved with one real phonon in the system which carries the system momentum and which is spatially uncorrelated with the polaron.@xcite for @xmath173 this state becomes the first excited state of the system .
the difference between the energies of the ground state and the first excited state is the largest for @xmath133 , and is equal to @xmath10 .
the validity of the perturbative treatment requires that the weight of the second term in eqs .
( [ wc1function ] ) and ( [ wc2function ] ) is small , i.e. , the mean number of phonons associated with the lattice deformation has to satisfy    @xmath174    here , @xmath169 has been eliminated by using eq .
( [ xtot ] ) .
there are two ways to satisfy condition ( [ nphpol ] ) .
either the electron - phonon coupling is small , @xmath175 , or the lattice deformation is spread to a large number of lattice sites , @xmath176 . in the latter case , the total mean polaron deformation does not have to be small , @xmath177 , since @xmath0 can be larger than @xmath10 .
the translationally invariant form of the wave function for the @xmath73 method , given by eq .
( [ tivfunction ] ) , provides an energy gain due to the polaron delocalization . at the same time
, the spatial correlation between the electron and the lattice deformation has a finite length .
for instance , in eq.([tivfunction ] ) , this length is of the same order for @xmath76 and @xmath80 .
the perturbative calculation for @xmath88 in eq.([wc1function ] ) gives    @xmath178 ^ 2 - 1}\;.\ ] ]    @xmath88 measures the electron - lattice correlation length .
@xmath88 is independent of @xmath0 , which makes the correlation length finite , even in the limit @xmath179 for which the lattice deformation vanishes , @xmath180 .    on the other hand ,
the electron - lattice deformation correlation length and the polaron delocalization range are of the same order for the localized functions .
this may be easily seen from eq .
( [ iafunction ] ) . in the limit
@xmath179 they both become infinite .
thus in the weak coupling regime one obtains a localized polaron state @xmath181 of very large width , but associated with a tiny lattice deformation .
this is specific for one dimensional systems in which an attractive symmetric potential has always a bound electron state .
therefore the corresponding polaron energy is less than the free electron energy @xmath182 . in higher dimensions ,
an arbitrary attractive symmetric potential does not have a bound electron state for sufficiently small g ( or the electron binding energy is too small to balance the lattice deformation energy ) , and the total polaron energy is larger than @xmath183 , where @xmath184 is the dimension of the system .
this explains why , in the weak coupling regime , the adiabatic ( localized ) polaron functions@xcite fail to have an energy lower than the free electron energy , for the dimension of the system greater than one .
a detailed parallel perturbative analysis of polarons in one , two , and three dimensions is given in ref .  .
in the crossover regime the polaron hopping energy @xmath2 is not negligible , and the translationally invariant polaron functions should be used in order to obtain the full physical picture of the polaron . in order to calculate the numerical errors of different methods accurately ,
the present discussion of the crossover regime is restricted to the values of @xmath9 smaller than @xmath185 , for which results from the @xmath107 method are available .    in order to examine the crossover regime
it is instructive to calculate the energy difference between the ground and first excited state , @xmath186 .
let @xmath187 denote the value of @xmath0 for which this difference is minimal ,    @xmath188    our results indicate that @xmath187 is very close to the value of @xmath132 given by eq .
( [ romero ] ) .
for example , for @xmath189 , @xmath190 , while @xmath191 .
for smaller @xmath9 , @xmath187 and @xmath132 coincide even better .
the analysis of the effective mass,@xcite variational energy of the polaron ground state,@xcite polaron size,@xcite as well as the behavior of the first excitation energy , suggest that a dramatic change in the nature of the polaron ground state and the first excited state are intimately related in the crossover regime .
this can be well understood by considering the properties of the wave function in the @xmath73 method . near @xmath187 , this polaron wave function has two separate energy minima in the @xmath87-@xmath88 parameter space defined by eq .
( [ ve ] ) , which become degenerated for @xmath192 .
let the symbol @xmath193 denotes the lower minimum at @xmath194 , and the symbol @xmath195 the lower minimum at @xmath196 . @xmath197 and @xmath198 are the corresponding polaron wave functions .
even if they are not mutually orthogonal , they are still physically quite different .
the numerical data show that the translational invariance of @xmath197 contributes strongly to the polaron energy . on the other hand ,
the translational invariance of @xmath198 has almost negligible energy contributions , i.e. , @xmath198 describes an almost self - trapped polaron .
it is worth noting that the degenerate nature of the variational energy which has been reported for the toyozawa method@xcite is of the same kind as the one of the @xmath73 method discussed here .
namely , both @xmath73 and toyozawa method are based on the same polaron function ( [ tivfunction ] ) .
@xmath197 and @xmath198 can be combined to form new polaron functions ,    @xmath199    it should be mentioned that a similar combination of two states has been already used in ref .   in order to calculate the polaron ground state .
@xmath200 corresponds to eq.([lvfunction ] ) with @xmath105 , which means that the @xmath106 method is implied . from this treatment the improved ground state @xmath201 and the approximate first excited
polaron state @xmath202 are obtained .
it may be seen from fig .
[ fig3 ] that @xmath106 method describes an anticrossing of @xmath198 and @xmath197 states , which yields two orthogonal states , @xmath201 and @xmath203 . in order to have a better illustration of that anticrossing in fig .
[ fig3 ] , the exact ground state energy @xmath143 is subtracted from all plotted energy curves . for @xmath194
the _ light _ state @xmath197 is lower in energy , and participates in the ground state more than the _ heavy _ state @xmath198
. however , as @xmath0 increases this balance changes continuously in favor of @xmath198 .
the opposite trend is observed for the first excited state , which is heavier than the ground state for smaller @xmath0 and lighter for larger @xmath0 .
the @xmath74 method gives better results for the ground and first excited state when a large number of functions ( large @xmath94 ) in eq .
( [ lvfunction ] ) is used . from fig .
[ fig3 ] , one may estimate that @xmath204 at @xmath192 .
moreover , the energy @xmath205 , unlike @xmath206 , satisfies @xmath207 for all @xmath196 , which is an important improvement over the @xmath105 result .
figure [ fig4 ] shows the results for the ground state energy and the energy of the first excited state with total system momentum @xmath133 , as functions of @xmath0 , obtained by several different methods .
two plots correspond to @xmath208 and @xmath209 , respectively .
the ground state energy obtained by the @xmath107 method is subtracted from all the other results . @xmath210 and @xmath211
are used to mark three different polaron regimes with respect to the strength of the electron - phonon coupling .
for @xmath212 the mean number of phonons of the lattice deformation is less than one .
thus for @xmath213 we recognize the weak coupling regime .
for @xmath214 the @xmath108 polaron energy has a negligible error ( see fig .
[ fig1 ] ) , which means that for @xmath215 polarons are self - trapped , and we recognize the strong coupling regime .
the crossover regime is found in the interval @xmath216 , with @xmath217 .
it may be noted from fig .
[ fig4 ] that in the weak coupling regime , the energy obtained by the @xmath50 method ( @xmath69 ) is close to the free electron energy @xmath182 ( for @xmath218 , the absolute error of the @xmath50 method becomes greater than @xmath10 in some intervals of @xmath0 ) . in the strong coupling regime , for large @xmath0
, @xmath69 approaches the exact polaron energy .
the error of the @xmath73 method , as may be seen from fig .
[ fig4 ] , is the largest in the crossover regime in which the results can be improved by a better choice of @xmath53 and @xmath84 in eq .
( [ ve ] ) . in the strong coupling regime , the translationally invariant form of the function @xmath73 has no effect on the polaron energy , i.e. , both the @xmath73 and @xmath50 methods give very similar results .
the energy of the first excited polaron state , @xmath219 , intersects the energy of the ground state plus one phonon , @xmath220 , for @xmath221 . after the minimum of @xmath186 is reached in the crossover regime at @xmath222 , @xmath219 approaches @xmath220 asymptotically for @xmath223 . as has already been pointed out in sec .
[ f06 ] , the excitation of the renormalized phonon mode explains the nature of the first excited state in the strong coupling limit . in ref .
perturbation theory was used to calculate the frequency of this local phonon mode , @xmath224 , to the lowest order in @xmath225 ,    @xmath226    the obtained perturbative correction to the phonon frequency is adiabatic , i.e. , the square root in eq .
( [ kabanov ] ) is independent of mass @xmath16 .
it is worth noting that for large @xmath0 , @xmath219 calculated here shows approximately the same behavior as the energy of the ground state with one additional phonon of frequency ( [ kabanov ] ) ,    @xmath227    the lattice part of the @xmath74 function is spatially symmetric with respect to the electron .
therefore in the strong coupling regime the first excited state , which has been identified here as a local renormalized phonon mode of the self - trapped polaron , should be basically a symmetric oscillation of the lattice deformation around the central polaron site . in the vicinity of @xmath132 , on the other hand ,
the nature of the first excited state is currently explained by the anticrossing of the self - trapped and the delocalized polaron state .
the question of how the excited self - trapped polaron state in the strong coupling regime and the first excited state near @xmath132 may be linked together is a matter of further considerations .
one may speculate that the anticrossing of the excited self - trapped polaron state and the delocalized state will give the answer .
the present paper discusses the ground and first excited states of the polaron for three different regimes of the electron - phonon coupling parameter @xmath0 .
the results can be briefly summarized as follows . in the strong coupling regime
the polaron hopping energy to neighboring sites is negligible , and the self - trapped polaron states are obtained .
the results of the @xmath108 method suggest that the adiabatic picture of the localized polaron state is valid , in which some of the local lattice vibrations are renormalized by the presence of the electron .
the numerical data show no significant deviation from the adiabatic locking relation ( [ locking ] ) of the electron site density and the mean lattice deformation . in the small polaron case , the predominant effect of the electron is the lowering of the frequency of the vibration at the central polaron site .
the excitation of the renormalized phonon mode corresponds to the first excited state of the small self - trapped polaron .
the nature of the polaron ground state in the crossover regime has been discussed in a number of papers , and its rapid change with @xmath0 has been well established numerically . the difference between the energies of the first excited state and the ground state , as a function of @xmath0 , has a minimum for @xmath192 .
it is shown , by using the @xmath106 method , that the anticrossing of the self - trapped and the delocalized polaron state can link the behavior of the ground and first excited polaron state . according to the @xmath74 method , for @xmath196
the effective mass of the ground state is larger than the effective mass of the first excited state , while for @xmath194 the opposite is true .
in addition , it is found that @xmath187 , which characterizes the first excited state , and @xmath132 , obtained from ground state analysis ( eq.([romero ] ) ) almost coincide for @xmath228 .    upon further reduction of @xmath0 the total mean number of phonons bound by
the polaron becomes smaller than one , and the weak coupling regime is reached .
the nearly free electron is dressed by a cloud of virtual phonons , and its mass is slightly renormalized .
the first excited state of the system , with momentum @xmath229 , can be viewed as the ground state of the polaron plus one additional uncorrelated phonon , rather than as an excited polaron .
finally , it is worth noting that there is a simple sum rule for the mean total lattice deformation , eq .
( [ xtot ] ) , which is valid for any number of electrons and is independent of the system dimension .
this important sum rule may be extended to some other models in which the lattice deformation is linearly coupled to the local electron density .
@xmath230}\;,\\ s_\delta&=&y_\delta\;\sum_n\eta^*_n\eta_{n+\delta}\;,\\ e_\delta&=&-t\;y_\delta\;\sum_n\eta^*_n(\eta_{n+\delta+1}+\eta_{n+\delta-1 } ) \\&+&\hbar\omega\;s_\delta\;\sum_m\xi_m^*\xi_{m+\delta}\\ & -&g\;y_\delta\;\sum_n\eta^*_n\eta_{n+\delta}(\xi_n^*+\xi_{n+\delta})\;.\end{aligned}\ ] ]      l. d. landau , z. phys . * 3 * , 644 ( 1933 ) .
t. holstein , ann .
( n.y . ) * 8 * , 325 ( 1959 ) . h. de raedt and a. lagendijk , phys .
* 49 * , 1522 ( 1982 ) ; phys . rev .
b * 27 * , 6097 ( 1983 ) ; * 30 * , 1671 ( 1984 ) ; p. e. kornilovitch and e. r. pike , phys .
b * 55 * , r8634 ( 1997 ) ; d. w. brown , k. lindenberg , and y. zhao , j. chem . phys . *
107 * , 3179 ( 1997 ) .
c. zhang , e. jeckelmann , and s. r. white , phys .
lett . * 80 * , 2661 ( 1998 ) ; phys .
b * 60 * , 14 092 ( 1999 ) ; e.  jeckelmann and s.  r.  white , _ ibid . _ * 57 * , 6376 ( 1998 ) .
p. gosar , j. phys .
c * 8 * , 3584 ( 1975 ) .
a. s. alexandrov , phys .
b * 61 * , 12315 ( 2000 ) .
b. gerlach and h. lwen , phys .
b * 35 * , 4291 ( 1987 ) ; rev . mod . phys . * 63 * , 63 ( 1991 ) ; h.  lwen , phys . rev .
b * 37 * , 8661 ( 1988 ) .
d. feinberg , s. ciuchi , and f. de pasquale , int .
b * 4 * , 1395 ( 1990 ) .
g. kalosakas , s. aubry , and g. p. tsironis , phys .
b * 58 * , 3094 ( 1998 ) .
y. toyozawa , prog .
* 26 * , 29 ( 1961 ) .
y. zhao , d. w. brown , and k. lindenberg , j. chem . phys . * 107 * , 3159 ( 1997 ) . h. b. shore and l. m. sander , phys .
b * 7 * , 4573 ( 1973 ) . k. h. hock , h. nickisch , and h. thomas , helv .
acta * 56 * , 237 ( 1983 ) .
g. venzl and s. f. fischer , phys .
b * 32 * , 6437 ( 1985 ) .
a. klamt , j. phys .
c * 21 * , 1953 ( 1988 ) .
a. la magna and r. pucci , phys .
b * 53 * , 8449 ( 1996 ) .
j. ranninger and u. thibblin , phys .
b * 45 * , 7730 ( 1992 ) .
f. marsiglio , phys .
a * 180 * , 280 ( 1993 ) ; physica c * 244 * , 21 ( 1995 ) .
a. s. alexandrov , v. v. kabanov , and d. k. ray , phys .
b * 49 * , 9915 ( 1994 ) . g. wellein and h. fehske , phys . rev .
b * 56 * , 4513 ( 1997 ) ; * 58 * , 6208 ( 1998 ) . j. bona , s. a. trugman , and i. batisti ' c , phys .
b * 60 * , 1633 ( 1999 ) .
a. h. romero , d. w. brown , and k. lindenberg , phys .
b * 59 * , 13 728 ( 1999 ) .
m. capone , w. stephan , and m. grilli , phys .
b * 56 * , 4484 ( 1997 ) .
j. m. robin , phys .
b * 58 * , 14 335 ( 1998 ) .
d. emin and t. holstein , phys .
lett . * 36 * , 323 ( 1976 ) . v. v. kabanov and o. yu .
mashtakov , phys .
b * 47 * , 6060 ( 1993 )
. a. h. romero , d. w. brown , and k. lindenberg , phys .
a * 254 * , 287 ( 1999 ) .
a. h. romero , d. w. brown , and k. lindenberg , phys .
b * 60 * , 4618 ( 1999 ) .
a. h. romero , d. w. brown , and k. lindenberg , phys .
a * 266 * , 414 ( 2000 ) . v. cataudella , g. de filippis , and g. iadonisi , phys .
b * 60 * , 15 163 ( 1999 ) ; * 62 * , 1496 ( 2000 ) .