two - dimensional ( 2d ) superconducting arrays , weakly coupled by josephson junctions , have attracted much attention  @xcite . in the classical array , where the charging energy is neglected , the relevant variable is the phase of the superconducting order parameter at each grain , and the logarithmic interaction between vortices is well known to lead to the berezinskii - kosterlitz - thouless ( bkt ) transition  @xcite .
recent development of the fabrication techniques , on the other hand , allows to make regular arrays composed of very small superconducting grains . in such arrays the
charging energy can not be neglected any longer  @xcite , and studies considering both the self - capacitance @xmath0 and the junction capacitance @xmath1 are required . in the presence of ohmic dissipation ,
the charges , which are conjugate to the phase variables , take continuous values and the corresponding phase diagrams have been obtained through the use of a variational method  @xcite .
this has revealed the existence of a low - temperature reentrant transition , which appears consistent with the recent quantum monte carlo study  @xcite . without dissipation , in contrast
, the charges change in discrete quanta @xmath2 with @xmath3 being the electron charge . in this case
the phase diagrams have been investigated via the coarse - graining approximation , which is mean - field - like in nature  @xcite . in general , the mean - field approximation does not give reliable results in two dimensions , where fluctuations are too strong to be neglected .
it is thus desirable to study the 2d superconducting arrays , with general capacitance matrices and discrete charge states , beyond the mean - field approximation .
the quantum phase transitions in the superconducting arrays have drawn much interest also in relation to the bose - hubbard model , which describes strongly interacting bosons under the competition between the kinetic - energy and the potential - energy effects . here the kinetic energy and the potential energy correspond to the josephson energy and the charging energy in the superconducting array , respectively . in the presence of the on - site potential ,
the system at zero temperature exhibits a transition between the mott insulating phase and the superfluid phase , displaying lobe - like structures in the resulting phase diagram  @xcite .
in addition , nearest - neighbor interactions and next nearest - neighbor interactions produce various phases such as the checker - board - type solid , the striped solid , and the supersolid .
in particular , the supersolid phase , which is characterized by superfluidity in the presence of underlying commensurability , has been studied extensively  @xcite .
this paper investigates the phase transitions in 2d superconducting arrays with emphasis on the effects of charge frustration and on the competition between the self- and junction capacitances . at zero temperature
, we use the perturbation expansion and the simulated annealing method to obtain various insulating lobes . here
rational charge frustration introduces commensurability effects to the array , leading to the solid phases with various charge densities .
the detailed phase boundaries depend crucially on the ratio of the junction to self - capacitances , @xmath4 : for example , the insulating lobes become narrower and their central positions approach the charge frustration values , as @xmath4 is increased .
further , as the system size is increased , more and more insulating lobes are expected to be observed in the phase diagram , suggesting the interesting possibility of the lobe at every rational frustration in the thermodynamic limit . to check the robustness of the lobes against finite temperatures , we apply the monte carlo ( mc ) method to the system without the josephson term and find the persistence of the insulating lobes with low - order rational values of the charge density .
we also use a variational method to obtain an effective classical hamiltonian describing the system at finite temperatures .
the corresponding phase diagrams are found approximately in a self - consistent manner , which again display lobes with the half - filled charge density . in the absence of charge frustration ,
the effective classical hamiltonian has only real terms , which allows us to use the mc method : we obtain the critical strengths of the josephson coupling and compare them with ones obtained from the self - consistent approximation . in the self - charging limit
, the possibility of a reentrant transition is also pointed out .
this paper is organized as follows : section  [ sec_model ] introduces the quantum mechanical hamiltonian for the superconducting arrays with general capacitance matrices and discrete charge states .
such a superconducting array is described by quantum mechanical conjugate variables , charges and phases , and charge frustration is induced by applying an external voltage between the array and the substrate . in sec .
[ sec_zero ] , we study the zero - temperature phase transitions in the system via the perturbation expansion together with the simulated annealing method . the phase transitions at vanishing charge density is investigated in sec .
[ subsec_0 ] , where the results in the nearest - neighbor charging limit are also compared with those of the previous works . in sec .
[ subsec_pq ] , phase transitions at finite charge densities are studied , and the nature as well as the existence of the supersolid phase is discussed .
section  [ sec_finite ] is devoted to the investigation of the finite - temperature phase transitions by means of a variational method and mc simulations .
the system without the josephson coupling is described by only charge variables , and we present the corresponding mc results in sec .
[ subsec_mc ] , where the charge densities are obtained as functions of the charge frustration at various temperatures .
the resulting step structures at zero temperature are compared with those obtained in sec .
[ subsec_pq ] . whereas the phase diagrams in the presence of the josephson coupling are obtained in a self - consistent manner in sec .
[ subsec_var ] , the comparison with mc results reveals that the self - consistent approximation is accurate only at sufficiently high temperatures .
finally , a summary of the paper is given in sec .
[ sec_conclusion ] .
we consider an @xmath5 ( @xmath6 ) square array of josephson junctions .
the extra charge @xmath7 on the superconducting grain at site @xmath8 can be written as    ' '' ''    ' '' ''    @xmath9 where @xmath0 and @xmath1 are the self- and the junction capacitances , respectively , and @xmath10 is the electric potential of the grain @xmath8 with respect to the substrate . if we apply an external voltage @xmath11 between the substrate and the ground , eq .
( [ eq_q ] ) is changed to @xmath12+@xmath13 , with the `` gauge charge '' @xmath14 .
the hamiltonian describing such a superconducting array consists of two parts : @xmath15    ' '' ''    2 where @xmath16 represents the phase of the superconducting order parameter at site @xmath8 .
note that the charge @xmath7 and the phase @xmath16 are quantum mechanical conjugate variables satisfying the commutation relation @xmath17 = 2ei$ ] ( @xmath18 ) . in the charging energy term @xmath19 ,
the charges at sites @xmath8 and @xmath20 interact via the inverse of the capacitance matrix @xmath21 whereas the summation in the josephson energy term @xmath22 is to be done over all the nearest neighboring pairs .
when the lattice constant of the array approaches zero , we can write eq .
( [ eq_q ] ) in the continuum form @xmath23 and obtain @xmath24 due to a single positive charge at the origin  @xcite : @xmath25 with the modified bessel function @xmath26 , where the screening length @xmath27 ( @xmath28 ) measures the interaction range between charges . the above hamiltonian in eq .
( [ eqh1 ] ) is closely related to the bose - hubbard hamiltonian and to the spin-1/2 @xmath29 hamiltonian  @xcite .
the charging energy term in eq .
( [ eqh1 ] ) drives the system to arrange its discrete charges such that the average charge per site is as close to @xmath13 as possible , and the josephson energy corresponds to the hopping energy in the bose - hubbard model  @xcite .
consequently , if the charging energy is sufficiently larger than the josephson energy , hopping is suppressed and all the charges are arranged to form a periodic lattice , leading to the insulating phase at zero temperature . in particular , when @xmath30 , the energy eigenstate can be chosen also as a simultaneous eigenstate of the charge operator since @xmath31=0 $ ] .
the uncertainty relation between the charge and the phase then asserts that phase coherence and consequently , superconductivity can not be attained in this limit . in the opposite limit of the strong josephson coupling ,
all the charges are in the extended state owing to the dominant hopping term , and the superconducting phase is favored . for convenience ,
we write @xmath32 , @xmath33 , @xmath34 , and @xmath35 , and obtain eq .
( [ eqh1 ] ) in the form    ' '' ''    ' '' ''    @xmath36    ' '' ''    2 where the inverse of the dimensionless capacitance matrix is given by @xmath37
in general the state which minimizes the energy determines the zero - temperature phase of the system .
suppose that we have obtained the energy levels of the system which correspond to the ground state , the first excited state , etc .
if we introduce a small change in the control parameters , the energy levels shift by a small amount , and it is possible that the ground state and the first excited state cross each other as we pass through a set of values of the control parameters ; this `` level crossing '' leads to a zero - temperature phase transition in the parameter space  @xcite .
it should be noted here that consideration of only the two low - lying states , the ground state and the first excited one , may not be sufficient .
figure  [ fig_sch ] shows schematic dependence of the energy levels @xmath38 on the control parameter @xmath39 in two possible cases . in fig .
[ fig_sch ]  ( a ) , the ground state and the first excited state cross at the critical value @xmath40 of the control parameter , which corresponds to the true transition point . in fig .
[ fig_sch ]  ( b ) , on the other hand , the true transition point is given by @xmath41 , where the ground state and the second excited one cross each other . in this case consideration of only the two low - lying states would lead to @xmath40 , which is not the true transition point . to find the true transition point
, we thus need to consider all the energy levels , which is not possible in practice . despite this
, considering the level crossing of the ground state and the first excited state is very useful and convenient to obtain the upper limit of the true transition point .    in terms of charges
the system can be either in the incommensurate conducting phase or in the commensurate insulating phase ; in the latter charges are expected to form a superlattice which covers periodically the whole system . if the commensurate state has energy lower / higher than that of the incommensurate state for given values of parameters , the system should be in the insulating / conducting phase .
( here @xmath42 , @xmath43 , and @xmath4 constitute the control parameters . ) accordingly , for given total number of charges , it is needed to find the charge configurations of the commensurate state and of the low - lying incommensurate states .
the symmetry of the hamiltonian under the transformations @xmath44 allows us to consider @xmath43 only in the range @xmath45 . in this restricted range @xmath46 ,
the commensurate states of the system may be classified in terms of the average number of charges , @xmath47 @xmath48 . for @xmath49 with relatively prime integers @xmath50 and @xmath51 , the charges are expected to form a @xmath52 superlattice , on the unit cell of which @xmath53 charges , each with the unit charge , are located .
we first investigate the phase transitions in the system with @xmath54 .
the energy of the @xmath54 commensurate state is computed via the second - order perturbation expansion with the josephson term as a perturbation in eq .
( [ eqh2 ] ) , and is written as @xmath55 without the perturbation , the hamiltonian of the system contains only charge operators and yields the eigenfunction    ' '' ''    ' '' ''    @xmath56    2 [ ] which is simply the plane - wave state of the usual free particle system , with the interpretation of @xmath16 and @xmath57 as the position and the momentum , respectively . for the commensurate state with @xmath58 , all the charge eigenvalues
@xmath57 s are zero and it is easy to compute the zeroth - order contribution @xmath59 and the first - order contribution @xmath60 : @xmath59 = @xmath61@xmath62 = @xmath63 and @xmath60 = @xmath64 = @xmath65 @xmath66 the second - order contribution is given by @xmath67 where @xmath68 is the zeroth - order energy for the charge configuration @xmath69 , and    ' '' ''    ' '' ''    @xmath70 \nonumber \\ & = & -\frac{e_j}{2}\sum_{\langle i , j\rangle}\left[\prod_{k\ne i , j }    \delta(q'_k,0)\right ] \left[\delta(q'_i,1 ) \delta(q'_j,-1 ) +      \delta(q'_j,1 ) \delta(q'_i,-1 ) \right ] .\end{aligned}\ ] ]    ' '' ''    2 [ ] inserting eq .
( [ eq_v ] ) into eq .
( [ eq_ec2 ] ) , we obtain the energy of the @xmath54 commensurate state to the 2nd order in @xmath71 : @xmath72    we next consider the incommensurate state with @xmath73 . among the various excitations present in the system with the constraint @xmath54
, the lowest excitation is expected to be point - like and we consider two types of excitations : one is the creation of a single charge ( sc ) at one site  @xcite , and the other is the creation of a charge dipole ( cd ) at one bond [ @xmath74 , @xmath75 with @xmath76 being the nearest - neighboring pair ] .
both excitations do not change the value of @xmath77 in the thermodynamic limit ( @xmath78 ) . in the self - charging limit ( @xmath79 )
, the interaction matrix @xmath21 takes a diagonal form and the first excited state should be of the sc type , since the cd type excitation requires the creation energy of two charges . on the other hand , in the nearest - neighbor ( nn ) charging limit ( @xmath80 ) , only the excitations which does not change the total number of charges are allowed and the first excited state is expected to be of the cd type .
accordingly , as the value of @xmath4 is increased , the first excited state should change from the sc type to the cd type at a critical value of @xmath4 .
in particular , the system in the nn charging limit displays the charge - anticharge unbinding transition  @xcite , the bkt nature of which is manifested by the square - root cusp in the resistance  @xcite .
this charge bkt transition is suppressed by the presence of the self - capacitance , making the screening length finite , and it seems likely that the change of the first excited state according to the value of @xmath4 in turn leads to the change in the nature of the transition .    between the above two types of excitations ,
the lowest one is determined from the comparison of the zeroth - order energies of both types : if @xmath81 given by    ' '' ''    ' '' ''
@xmath82 is negative / positive , the creation of a single charge / charge dipole is the lowest excitation .
figure  [ fig_1st ] shows the boundary for @xmath83 in the limit of @xmath84 .
it is found that for @xmath85 the lowest excitation is of the sc type regardless of the value of @xmath43 , and that the region of the cd type excitation ( the region below the curve in fig .
[ fig_1st ] ) is quite small , allowed only for @xmath86 .
therefore we below focus on the sc type excitation for @xmath87 and @xmath88 , and then consider the cd type excitation only in the nn charging limit ( @xmath89 ) with @xmath90 .
we write the energy of the first excited state of the sc type as @xmath91 and the eigenstate as @xmath92 which is expected to be a superconducting state since the extra charge at site @xmath93 can hop to any site without energy cost . inserting eq .
( [ eqql ] ) to @xmath19 in eq .
( [ eqh2 ] ) , we obtain the zeroth - order contribution : @xmath94 = @xmath95 , where we have used the translational symmetry @xmath96 and @xmath97 . since the above excited state in eq .
( [ eqql ] ) has an @xmath98-fold degeneracy ( @xmath93 can take values @xmath99 ) , we use the degenerate perturbation theory to calculate the first - order contribution @xmath100 .
we thus need to diagonalize the matrix @xmath22 whose components are given by @xmath101 for that purpose , it is convenient to use the fourier transform : @xmath102 which gives straightforwardly the eigenvalues @xmath103 . the first - order contribution
is then given by the lowest eigenvalue : @xmath104 .
for this lowest - lying excited state with @xmath105 , it is also straightforward to obtain the second - order contribution , which leads to the energy of the incommensurate state @xmath106 here @xmath107 denotes the summation over @xmath8 and its four nearest - neighbors @xmath20 with the restriction @xmath108 .
the zero - temperature transition occurs when the energy of the commensurate state given by eq .
( [ eq_ec0 ] ) and that of the incommensurate state given by eq .
( [ eqeic0 ] ) becomes equal .
this reveals the zero - temperature phase transition in the system with @xmath54 , as @xmath43 is varied .
namely , the system becomes superconducting as @xmath43 is increased beyond the critical value @xmath109 , which is given by @xmath110 the resulting phase boundary between the insulating phase and the superconducting phase is shown in fig .
[ q0_bound ] , obtained for a @xmath111 array .
the junction capacitance here tends to enhance superconductivity , which is in accord with the result in ref .  .
it is also obvious that as the system approaches the nn charging limit , the insulating region shrinks to zero ; _ this apparently cures the well - known artifact of the coarse - graining approximation _
, where the system remains in the insulating phase for arbitrarily large values of @xmath71 in the nn charging limit  @xcite .    to check the validity of the second - order perturbation expansion ,
we have computed the third - order contribution in the self - charging limit , and obtain the critical value @xmath109 : @xmath112    2 [ ] the three curves in fig .
[ 3rd_bound ] represent the phase boundaries obtained from eq .
( [ eq_3rd ] ) up to the first- , the second- , and the third - order contributions , respectively .
it is shown that the third - order contribution does not change significantly the boundary obtained from the second - order calculation , inducing 8.3% of the difference in the critical value of @xmath42 at @xmath90 .
we next consider the nn charging limit in the absence of charge frustration ( @xmath90 ) . in this case , eq .  ( [ eq_ec0 ] ) gives the energy of the commensurate state in the form @xmath113 where @xmath114 and @xmath115 in the thermodynamic limit , we replace the discrete sum by the integral and obtain    ' '' ''    ' '' ''    @xmath116    ' '' ''    2 [ ] leading to @xmath117 .
the first - excited incommensurate state in the nn charging limit is of the cd type ( see fig .  [ fig_1st ] ) with the zeroth - order energy @xmath118=@xmath119-@xmath120 = @xmath121 and the null first - order energy @xmath122 .
it is complicated but straightforward to calculate @xmath123 from the second - order degenerate perturbation theory , which leads to a @xmath124 matrix .
the fourier transform then yields @xmath98 @xmath125 hermitian matrices , each of which can be diagonalized .
the lowest eigenvalue then gives the energy of the cd type incommensurate state , up to the second order : @xmath126 where the number 62.000 has been obtained numerically for a sufficiently large array ( @xmath127 ) ; at this size finite - size effects have been confirmed to be negligible .
the critical value of @xmath128 , below which the array is in the insulating phase , is obtained from the condition @xmath129 : @xmath130 note that this value is somewhat larger than the critical value @xmath131 obtained in refs .  .
one possible explanation of the above discrepancy is the failure of the second - order perturbation expansion .
namely , consideration of higher - order terms may eliminate the discrepancy .
the other possibility is that the cd type excitation does not give the true transition point : see fig .
[ fig_sch ]  ( b ) for a schematic picture , where @xmath40 may correspond to the obtained transition point @xmath132 and @xmath41 the true transition point . in experiment ,
the even larger critical value @xmath133 has been obtained  , via the curve - fit to the square - root cusp form of the resistance  @xcite .
recent quantum monte carlo simulations , on the other hand , have yielded the result @xmath134 .      in the nn charging limit
, there exists an interesting duality between the charge and the vortex  @xcite , which concludes that in the absence of the josephson term , the charges form a superlattice determined by the value of charge frustration @xmath43 .
this is a counterpart of the vortex superlattice in the classical array without the charging energy , formed in the presence of magnetic frustration  @xcite . in the system considered here , on the other hand , the non - vanishing self - capacitance gives the interaction range between charges @xmath136 in units of the lattice constant , thus making it finite  @xcite .
accordingly , we expect that the formation of the @xmath52 charge superlattice with @xmath137 is improbable in the system with the rational charge frustration @xmath138 .
this implies that @xmath77 is not necessarily equal to @xmath43 in the array with non - vanishing @xmath0 , which reflects that a short - ranged potential in general makes the commensurate phase extend over a finite range of the chemical potential .
( note that the charge frustration in this work plays the role of the chemical potential . )
we now investigate the phase transitions for @xmath139 @xmath140 , using the canonical ensemble with the total number of charges fixed .
to find the zero - temperature phase diagrams we need to know the configurations of the commensurate states and of the lowest - lying incommensurate states , which is a highly nontrivial problem in the 2d system considered here , especially for large @xmath51 .
for the 1d josephson junction array with @xmath141 , the commensurate state , the first excited state , and the resulting zero - temperature phase boundaries have been obtained via the perturbation expansion  @xcite .
we write the energy of the commensurate state with @xmath142 in the form @xmath143 where the relations @xmath144 and @xmath97 give @xmath145 . in general , when @xmath146 , the commensurate ground states are degenerate , which makes it necessary to consider the matrix element :    ' '' ''    ' '' ''    @xmath147 & &   \left [ \delta\left(q_i^{(l)},q_i^{(m ) } - 1\right )   \delta\left(q_j^{(l)},q_j^{(m ) } + 1\right ) \right .
\nonumber \\   & & + \left .
\delta\left(q_i^{(l)},q_i^{(m ) } + 1\right )   \delta\left(q_j^{(l)},q_j^{(m ) } - 1\right ) \right ] \end{aligned}\ ] ] with @xmath148 being the @xmath93th degenerate commensurate state .
note that the configuration of the superlattice unit cell in a degenerate commensurate state is in general different from that in a different ( degenerate ) state .
therefore , for @xmath149 , @xmath150 differs from @xmath151 at the total of @xmath152 sites for finite @xmath51 , leading to @xmath153 .
further , it is obvious that @xmath154 .
it is thus concluded that all the components of the matrix @xmath22 vanish , which results in @xmath155 .
the second - order contribution is given by @xmath156 where @xmath157 and @xmath158 is the zeroth - order energy for given charge configuration @xmath69 and @xmath159 , respectively , and the summation @xmath160 is over the states @xmath159 making the denominator nonzero .
the superscript @xmath93 on @xmath69 has been dropped since all the degenerate commensurate states give the same contribution to @xmath161 .
namely , the degeneracy is not removed in a finite order since @xmath150 and @xmath151 can not be connected by a finite number of charge hoppings . with eq .
( [ eq_h ] ) , eq .
( [ eq_i ] ) reads @xmath162 where @xmath163 @xmath164 @xmath165 is obtained from @xmath69 by a charge hopping between sites @xmath8 and @xmath20 .
it is easy to compute the denominator in eq .
( [ eq_j ] ) : @xmath166 \nonumber \\ & = & 8e_0\left [ \sum_l\left ( { \widetilde c}_{jl}^{-1 } -{\widetilde c}_{il}^{-1 } \right)q_l - \left ( { \widetilde c}_{00}^{-1 } - { \widetilde c}_{10}^{-1}\right ) \right ] , \end{aligned}\ ] ] which leads to the expression for the energy of the commensurate state with @xmath167 : @xmath168 here @xmath57 is the charge at site @xmath8 in the commensurate state with @xmath169 .
since the ( commensurate ) ground state configurations are not known for general @xmath170 , @xmath57 in this work is obtained via the simulated annealing method .    to compute the energies of the incommensurate states with @xmath171
, we classify excitations according to the number @xmath172 of extra charges ( with the sign taken into account ) , and let @xmath173 denote the energy of this excited state which has the total number of charges @xmath174 . to the second order in @xmath71 , it is straightforward to obtain the energy of the @xmath172-charge excited state : @xmath175    ' '' ''    2 where @xmath176 is the charge at site @xmath8 in the @xmath172-charge excited state , and the first - order contribution @xmath177 has been observed to vanish in numerical calculations , regardless of the values of @xmath170 @xmath178 and @xmath172 considered in this work . further , in eq .
( [ eq_eic ] ) , we have assumed that degeneracy is not removed in the second order although there are some exceptions , e.g. , for @xmath170=1/2 and @xmath179 ( see below ) . to find the phase boundary , we need to know the charge configuration which corresponds to the first excited state for given values of @xmath4 and @xmath43 . for this purpose
, we use the following numerical procedure :    1 .    distribute @xmath180 charges randomly .
2 .   use the simulated annealing method to find the charge configuration @xmath176 , which minimizes the zeroth - order energy @xmath181 .
3 .    change the value of @xmath172 , and repeat the steps 1 - 3 .
4 .    among the obtained charge configurations with their zeroth - order energies ,
find one corresponding to the lowest - lying incommensurate state .    in the numerical simulations , we have used the conventional simulated annealing method adopting the metropolis algorithm for the arrays of sizes @xmath182 and @xmath183 .
starting from high temperatures ( @xmath184 ) , we decrease the temperature slowly with the decrements @xmath185 ( at high temperatures ) and 0.002 ( at low temperatures ) , and perform 10000 mc steps per site at each temperature step .
we have found that the charge configuration corresponding to the lowest energy depends crucially on the value of @xmath4 , and thus used three different values : @xmath4 = 0.1 , 1.0 , and 5.0 . through 20 independent runs ,
we have obtained the charge configurations corresponding to the lowest energy , which have also been checked via the entropic annealing algorithm  @xcite . although we have allowed double occupancy ( @xmath186 ) as well , in the lowest - lying configurations @xmath57
is always found to have the value either 0 or 1 .
once the charge configuration @xmath176 which minimizes the zeroth - order energy in eq .
( [ eq_eic ] ) is found , the phase boundary is determined from the condition @xmath187 .
this leads to the critical value @xmath109 separating the superconducting phase from the insulating phase :    ' '' ''    ' '' ''    @xmath188 \nonumber \\     & + & \frac{(e_j / e_0)^2}{256q_e}\left [        { \sum_{(i , j)}}^\prime         \frac{1}{\sum_l\left ( { \widetilde c}_{jl}^{-1 } -{\widetilde c}_{il}^{-1 }         \right)q^{(q_e)}_l - \left ( { \widetilde c}_{00}^{-1 } -          { \widetilde c}_{10}^{-1}\right ) } \right
.   \nonumber \\      & & \left .
-{\sum_{(i , j)}}^\prime         \frac{1}{\sum_l\left ( { \widetilde c}_{jl}^{-1 } -{\widetilde c}_{il}^{-1 }         \right)q_l - \left ( { \widetilde c}_{00}^{-1 } -          { \widetilde c}_{10}^{-1}\right ) } \right ] , \end{aligned}\ ] ]    ' '' ''    2 where @xmath189 again describes the charge configuration of the commensurate state .
when @xmath190 , the system , displaying superconductivity for @xmath191 , becomes insulating as @xmath43 is increased beyond @xmath109 .
this behavior according to whether @xmath43 is larger or smaller than @xmath109 is reversed when @xmath192 .
table  [ tab_qe ] shows the values of @xmath172 in the lowest excitations , breaking the commensurability of the states with @xmath77= 1/2 , 1/3 , 1/4 , and 1/5 , in the arrays of sizes @xmath193 10 and 12 .
except for @xmath194 , there are two values of @xmath172 in each case , one for the positive - charge excitation ( @xmath195 ) and the other for the negative - charge excitation ( @xmath196 ) .
since the value of the charge frustration has been restricted in the range @xmath45 , we need to consider only the negative - charge excitations for @xmath197 .
it displays that as @xmath1 is increased , the absolute value of @xmath172 decreases , manifesting that the excitations become more point - like . as an example
, we present the charge configurations for @xmath198 in fig .
[ q2_conf ] , where ( a ) , ( b ) , and ( c ) correspond to the commensurate state , the excited state with @xmath179 , and the excited state with @xmath199 , respectively , in a @xmath200 array .
it is of interest to note that the configuration ( c ) , which has 60 charges , is precisely the same as the ground state configuration of vortices in the presence of the magnetic frustration 5/12  @xcite .
when @xmath179 is the lowest excitation , there exists @xmath201-fold degeneracy since the vacant site , represented by a small square in ( b ) , can be located on any of the total @xmath201 sites . in this case , the degenerate states are coupled in @xmath202 , and we should use the second - order degenerate perturbation theory to calculate the energy . it is again straightforward to diagonalize the second - order matrix by means of the fourier transform , which yields    ' '' ''    ' '' ''    @xmath203 , \end{aligned}\ ] ]    ' '' ''    2 in place of eq .
( [ eq_eic ] ) . here
@xmath204 denotes the position of the vacant site .
figure  [ fig_zerot ] shows the resulting phase diagrams for @xmath4= ( a ) 0.1 , ( b ) 1.0 , and ( c ) 5.0 .
the various lobe - like structures , in which the array is in insulating phase , are obtained .
since the screening length @xmath136 , beyond which the magnitude of the interaction between charges decreases exponentially , is much smaller than the system size , the finite - size effects on the phase boundaries are presumably negligible .
we have checked the phase boundary for @xmath197 , and indeed found that it hardly changes with the system size considered here ( see below ) .
accordingly , the phase boundaries in fig .
[ fig_zerot ] are expected to be qualitatively correct even in the thermodynamic limit .
as we increase the array size , more insulating - lobes should be observed , since the array with @xmath30 should be insulating , regardless of the value of @xmath43 . at this stage , however , it is still unclear whether the insulating lobe exists at every rational value of @xmath43 . in the nn charging limit , the charging - energy term takes the form @xmath205 with @xmath206 given by eq .
( [ eq_cbar ] ) . here
the divergence of @xmath207 leads to @xmath208 , which is simply the counterpart of the vortex - neutrality condition in the classical array .
indeed fig .
[ fig_zerot ] shows that the position of the insulating region approaches @xmath208 , as @xmath4 is increased .    to study the effects of the junction capacitance versus the self - capacitance in detail
, we now concentrate on the @xmath197 insulating phase , whose commensurate charge configuration displays the @xmath209 superlattice structure .
[ see fig .
[ q2_conf ]  ( a ) . ]
similarly to the vortex superlattice in the fully - frustrated classical array  @xcite , the commensurate state has the two - fold degeneracy .
if we remove a single positive charge from the system , commensurability is destroyed near the resulting point defect .
the charge configuration of the corresponding lowest - lying incommensurate state has been shown in fig .
[ q2_conf ]  ( b ) , which describes the first excited state for @xmath4 = 1.0 or 5.0 . when @xmath4 = 0.1 , in contrast , the first excited state has the configuration shown in fig .  [ q2_conf ]  ( c ) , with the number of extra charge @xmath199 .
this observation that the nature of the lowest excited state changes with @xmath4 has an important implication with regard to the supersolid phase , where a commensurate charge - density wave is believed to coexist with superfluidity  @xcite : the charge configuration in fig .
[ q2_conf ]  ( b ) regarded as a snapshot of the charge configuration in the supersolid phase demonstrates that the underlying global commensurability still holds despite of the local point defect . furthermore , since the point defect can hop to any site without energy cost , it is delocalized in the energy eigenstate , leading to superfluidity .
accordingly , we conclude that the existence of the supersolid phase depends on the value of @xmath4 . in particular , when @xmath4 is sufficiently small
, the supersolid phase does not exist .
it is of interest to compare this argument pointing out the importance of the vacant defect in supersolidity with that of ref .  , emphasizing the role of the double occupancy . in the numerical investigations to find the lowest - energy charge configurations , we have found that all the charges have the value either 0 or 1 at zero temperature , and that those configurations with @xmath210 or 2 at some sites have much higher energies .
the @xmath197 phase boundaries in fig .
[ fig_zerot ]  ( b ) and ( c ) may be considered to divide the insulating commensurate phase and the supersolid phase . as we go further into the supersolid phase starting from the phase boundary , the more point defects are generated ; the commensurability should be completely destroyed when we cross another phase boundary separating the supersolid phase from the superfluid phase , although within the formalism adopted here , we are unable to find this new boundary . to obtain the phase diagram for @xmath197 , we have considered @xmath5 arrays up to @xmath211 , and found that the finite - size effects are negligible for @xmath212 , which corresponds to the interaction range @xmath213 , and that the qualitative features of the phase diagram do not change even when @xmath214 .
we have further found that the difference between the phase boundary obtained from the excitation @xmath179 and that from @xmath199 is insignificant in the @xmath200 array with @xmath215 .
therefore , for simplicity , we consider only the excitation @xmath179 and compute the phase boundary from eqs .
( [ eq_o ] ) and ( [ eq_deg ] ) .
the resulting phase diagram obtained for the array of size @xmath211 is shown in fig .
[ q2_bound ] .
the insulating region is shown to expand with @xmath4 when @xmath216 and then to shrink as @xmath4 is increased further .
this is to be compared with the case @xmath54 ( shown in fig .
[ q0_bound ] ) , where superconductivity is enhanced monotonically with @xmath4 , even for small @xmath4 .
in this section , we investigate the phase transitions of the superconducting arrays at finite temperatures , where thermal fluctuations as well as quantum fluctuations suppress the ordering of phases . while at high temperatures the system should be in the normal ( disordered ) state due to large thermal fluctuations , strong charging effects at low temperatures produce large quantum fluctuations , tending to destroy superconductivity of the system .
we thus expect rich behaviors resulting from the interplay of thermal and quantum fluctuations , which may be controlled by the temperature and charging energy .
we first study the system in the absence of the josephson coupling via mc simulations , and then investigate the phase transitions in the presence of the josephson coupling via a variational method together with mc simulations .      in the absence of the josephson coupling ( @xmath217 ) ,
the hamiltonian is diagonal in the charge representation and written in the form @xmath218 where @xmath57 is the eigenvalue of the charge operator at site @xmath8 .
we perform mc simulations on the above hamiltonian to investigate how the insulating phase corresponding to the regions inside the lobes in fig .
[ fig_zerot ] changes with the temperature . using the conventional mc method with the metropolis algorithm
, we start from a sufficiently high temperature ( @xmath219 ) , and decrease the temperature slowly with the decrements @xmath185 ( at high temperatures ) and 0.002 ( at low temperatures ) , during the simulations , we allow creation and annihilation of a charge at a site , as well as hopping to other sites . at each temperature step , we compute the thermal average @xmath220 using 10000 mc steps per site , after 1000 mc steps of equilibration .
figure  [ fig_mc ] displays the results for arrays of sizes @xmath221 ( left ) and 12 ( right ) , with @xmath222 ( a ) 0.1 , ( b ) 1.0 , and ( c ) 5.0 . at zero temperature , the curves representing @xmath77 versus @xmath43 exhibit steps at rational values giving the ( commensurate ) insulating phase . such steps structures , which correspond to the lobe structures in fig
.  [ fig_zerot ] , are common in systems displaying commensurate - incommensurate transitions , and signal the commensurate locking  @xcite .
it is shown that , as expected , the width of a step decreases with @xmath4 , and the curve approaches the straight line @xmath223 in the limit @xmath224 , manifesting the suppression of the insulating phase . comparing the results for the @xmath221 array with those for the @xmath225 array
, we also expect more rational steps in larger arrays , which is apparently suggestive of the devil s staircase structure  @xcite in the thermodynamic limit . for conclusive results , however ,
more detailed study is needed .
table  [ tab_mc ] presents the widths of the rational steps of @xmath226 1/3 , and 1/2 in the @xmath225 array at zero temperature .
one can see here that the values in table  [ tab_mc ] are in good agreement with the @xmath217 results in fig .
[ fig_zerot ] , which supports the validity of the approach in sec .
[ sec_zero ] .
figure  [ fig_mc ] also shows that steps tend to disappear at finite temperatures , reflecting that thermal fluctuations suppress the charge ordering due to quantum coherence as well as the phase ordering .
accordingly , the insulating commensurate phase in general turns into the normal ( disordered ) phase as the temperature is increased from zero .
note , however , that some steps survive weak thermal fluctuations . namely , whereas steps at higher - order rationals easily disappear , low - order steps such as @xmath54 , 2 , and 1/3 persist at low but nonzero temperatures . these features presumably can also be observed in the systems with nonzero josephson couplings .
thus as we increase the temperature from zero , higher - order insulating lobes in fig .
[ fig_zerot ] disappear , but some low - order lobes are expected to remain at low temperatures .      in the presence of the josephson coupling , both the charge and the phase variables should be treated quantum - mechanically , which prohibits exact analytical treatment . here , we use the giachetti - tognetti - feynman - kleinert ( gtfk ) variational method  @xcite to evaluate the path integral and to obtain the effective classical hamiltonian , which is convenient for studying the phase transitions at finite temperatures .
the gtfk method , which remains reliable in known cases even at zero temperature , has been successfully applied to the superconducting arrays with continuous charge states in the presence of ohmic dissipation  @xcite .
we begin with the hamiltonian given by eq .
( [ eqh2 ] ) and write the partition function of the system in terms of the path integral  @xcite : @xmath227 \right\}\ ] ] with the euclidean action @xmath228 where @xmath229 is the inverse temperature , and the planck constant has been set equal to unity ( @xmath230 ) .
the euclidean lagrangian @xmath231 can be obtained from the hamiltonian via the wick rotation @xmath232 and the legendre transformation .
we use the representation @xmath233 obeying the commutation relation @xmath234 = i$ ] , and get    ' '' ''    ' '' ''    @xmath235 with @xmath236 , where the charge frustration enters through a purely imaginary term .    in the presence of ohmic dissipation ,
the charge at each site takes continuous values , and only the paths satisfying @xmath237 contribute to the partition function , with the dissipative action term added  @xcite . in this case , the charge frustration does not appear in the euclidean action since @xmath238 , and the resulting phase diagram does not depend on @xmath43 .
if there is no ohmic dissipation present , in contrast , only discrete charge states are allowed , and all the paths satisfying @xmath239 with integer @xmath240 contribute to the path integral in eq .
( [ eqz ] )  @xcite .
therefore the path integral should include the summation over the `` winding numbers '' @xmath241 , yielding interesting behaviors associated with the charge frustration  @xcite .
accordingly , the allowed charge states , i.e. , continuous or discrete , are crucial in the resulting phase diagrams .    from the boundary condition
@xmath242 , we decompose @xmath16 into the periodic variable @xmath243 satisfying @xmath244 and @xmath240 , according to @xmath245 , and obtain the euclidean action in the form @xmath246 where @xmath247 \right\ } .\ ] ] following the gtfk method , it is straightforward to obtain the effective classical hamiltonian @xmath248 where @xmath249 is determined by @xmath250^{-1 } \right\ } \nonumber\\ & = & \label{eq : g0 } \frac{1}{8 g_0 \beta e_j } ( x_1 \coth x_1 + x_2 \coth x_2 - 2)\end{aligned}\ ] ] with @xmath251 and @xmath252 .
here we have used identity @xmath253 when eq .
( [ eq : g0 ] ) has more than one zeros , we should choose the largest one according to the extremal principle .
the detailed procedure is entirely similar to that in ref .  , and will not be repeated here . as a result ,
the quantum fluctuations renormalize the josephson coupling @xmath71 to @xmath254 . furthermore
, in contrast to the continuous charge case  @xcite , the charge frustration appears explicitly in the hamiltonian , and is expected to play an important role in the phase transition of the system .
unfortunately , the kronecker delta @xmath255 between winding numbers of nearest neighboring pairs in eq .
( [ eqhcl ] ) makes further analysis very difficult . for simplicity ,
we thus replace it by its self - consistent average @xmath256 and obtain the approximate hamiltonian @xmath257 where @xmath258 is the ensemble average with respect to the approximate hamiltonian @xmath259 .
the poisson summation formula then allows us to write @xmath260 }      { \sum_{\{m_i\ } }      \exp \left [ - 4\beta e_0      \sum_{i , j } ( m_i - q_x ) { \widetilde c}_{ij}^{-1 } ( m_j - q_x ) \right]}\ ] ]    ' '' ''    2 with @xmath261 , which manifests the absence of dependence on @xmath71 . once @xmath256 is determined from eq .
( [ eqde ] ) , the partition function of the system obtains the form @xmath262 where the charge - part @xmath263\ ] ] has been obtained from the poisson summation formula and the phase - part is given by @xmath264.\ ] ] unless the ratio @xmath4 is infinite , the interaction @xmath265 between charges is short - ranged and does not have any singularity .
it is then reasonable to assume that the criticality of the system is governed by the phase - part , which implies that the system exhibits a vortex - unbinding bkt transition at the critical temperature given by  @xcite @xmath266    we first consider the self - charging limit ( @xmath79 ) . in this limit , @xmath256
can easily be evaluated since the terms in eq .
( [ eqde ] ) factorize .
the resulting phase diagrams determined from eq .
( [ eqcritical ] ) at temperatures @xmath267 , and @xmath268 ( in units of @xmath269 ) are shown in fig .
[ fig_selfbd ] .
the region in the right - hand side of each line corresponds to the superconducting phase whereas the left region represents the normal ( disordered ) phase . in this phase diagram
, we observe that the superconducting region expands as @xmath43 is increased at all temperatures .
( recall that the symmetry of the system allows to restrict within the range [ 0,1/2 ] . )
therefore it is concluded that the charge frustration in general enhances superconductivity .
the qualitative features are largely similar to those of the results from the coarse - graining approach  @xcite at temperatures @xmath270 . unlike the latter , however , fig .
[ fig_selfbd ] also shows reentrant behavior at low temperatures , for small @xmath43 . at zero temperature ,
unfortunately , we obtain the unphysical result that the system is in the superconducting state for @xmath271 , even in the limit of the vanishing josephson coupling .
this is presumably an artifact of the self - consistent approximation ( sca ) replacing of the kronecker delta by its average , with correlations neglected .    when @xmath272 , on the other hand , the terms in eq .
( [ eqde ] ) do not factorize , and the average @xmath256 can not be evaluated analytically , which makes it inevitable to resort to a numerical method . at extremely low temperatures , @xmath273 in eq .
( [ eqde ] ) can have the value either 0 or 1 , which allows us to calculate @xmath256 by direct summation of @xmath274 terms for small @xmath275 . at high temperatures , other values of @xmath273
are allowed ; this makes the direct summation impractical and we use mc simulations .
figure  [ fig_1.0bd ] shows the resulting phase diagram obtained for a @xmath125 array with @xmath276 at temperature @xmath277 , 0.2 , 0.3 , and 0.4 . at @xmath278
, @xmath256 has been calculated by means of both the direct summation with @xmath279 and @xmath280 and mc simulations , which give results in good agreement with each other . at higher temperatures ,
@xmath256 has been obtained via mc simulations , yielding the enhancement of superconductivity by charge frustration in a similar manner to that in the self - charging limit . at low temperatures ( @xmath277 ) , on the other hand , there exists a small lobe near @xmath281 , which indicates the expansion of the insulating region due to the commensurability effects on the ordering of charges .    to check the validity of the sca adopted here , replacement of @xmath282 to its self - consistent value @xmath256 ,
we have performed mc simulations with @xmath90 , at which the hamiltonian in eq .
( [ eqhcl ] ) has only real terms :    ' '' ''    ' '' ''    @xmath283 we have considered a @xmath284 array with various values of @xmath42 , and measured the helicity modulus @xmath285 with 100000 mc steps : @xmath286 ,      \end{aligned}\ ] ]    ' '' ''    2 where @xmath287 with @xmath288 denoting the @xmath289 coordinate of the @xmath8th grain in units of the lattice constant .
the universal jump condition of the helicity modulus  @xcite @xmath290 then determines the critical coupling @xmath291 , the values of which are displayed in table  [ tab_j ] , for @xmath292 and 1 at various temperatures .
table  [ tab_j ] also displays the results from the sca replacing the kronecker delta by its average , which has yielded figs .
[ fig_selfbd ] and [ fig_1.0bd ] : the two results indeed exhibit behaviors qualitatively in accord with each other . in particular , for @xmath293 , both results give the value of @xmath291 larger at @xmath277 than at @xmath294 , strongly suggesting the reentrant transition .    from a quantitative viewpoint , on the other hand , the agreement between the two results is not satisfactory , especially at low temperatures , which seems to indicate that the sca fails at low temperatures . to examine this ,
we have also measured in the mc simulations the average value of @xmath255 for @xmath79 and display the results in fig .
[ fig_delta ] . at high temperatures ,
the mc data indeed agree well with the results from eq .
( [ eqde ] ) which is represented by the solid line in fig .
[ fig_delta ] . however , at low temperatures the difference is apparent unless @xmath42 is sufficiently small .
in particular , for @xmath295 , the mc simulations give the value approaching unity as the temperature is decreased to zero , whereas the value from eq .
( [ eqde ] ) approaches zero .
this is the case that the second term in eq .
( [ eq_hcl ] ) is dominant over the first term , which leads to the ordering of both @xmath296 and @xmath297 .
this reflects strong correlations between the two , which have not been taken into account in the sca .
it is thus concluded that the _ quantitative _ behaviors of the phase diagrams in fig .
[ fig_selfbd ] are accurate only at sufficiently high temperatures as already shown in table  [ tab_j ] .
we have investigated the quantum phase transitions in two - dimensional superconducting arrays with various values of external charge frustration at both zero and finite temperatures .
we have considered general capacitance matrices , allowing both the self- and the junction capacitances . in the absence of dissipation ,
the charge variable in the array changes in discrete quanta of @xmath2 , and the phase variable becomes compact in the interval @xmath298 .
we have studied in detail zero - temperature phase transitions at various values of the charge density and also considered the system at finite temperatures with and without the josephson energy .
the states of the system at zero temperature have been determined from the investigation of the crossing of the ground - state level and the first - excited - state level . at vanishing charge density ( @xmath54 ) ,
we have considered the commensurate insulating state , where @xmath299 at all sites , and two types of the incommensurate superconducting state , i.e. , the single - charge ( sc ) type and the charge - dipole ( cd ) type .
unless the self - capacitance vanishes , the sc type incommensurate state has been found to have lower energy than the cd type one except at very small values of the charge frustration . from the comparison of the energies of the commensurate insulating state and the sc type superconducting state , computed via the second - order perturbation expansions
, the phase boundary has been obtained .
the validity of the second - order expansion has been checked for the array in the self - charging limit , which reveals that the inclusion of the third - order term does not change significantly the phase diagram . at finite charge densities ,
identification of the commensurate state and of the first - excited incommensurate state is a highly nontrivial problem .
we have thus adopted the simulated annealing method and applied the second - order perturbation theory to obtain a number of lobe - like insulating phases such as @xmath197 , 1/3 , 1/4 , and 1/5 in @xmath182 and @xmath183 arrays . as the system size
is increased , more insulating lobes are expected to be observed .
the effects of the junction capacitance have been investigated in detail in the insulating phases of @xmath54 and 1/2 : as @xmath4 is increased , the insulating region shrinks monotonically for @xmath54 , while the insulating region for @xmath197 expands with @xmath1 and then shrinks , with its maximum area at @xmath300 .
the mc simulations have displayed that the lowest excited state changes as the junction capacitance is increased , which in turn indicates that the existence of the supersolid phase with @xmath197 depends on the interaction range of charges : the supersolid phase exists for sufficiently large values of the ratio of the junction capacitance to the self - capacitance . in the supersolid phase , it has also been shown that the point defect in the charge configuration plays an important role .    at finite temperatures , thermal fluctuations as well as quantum fluctuations
play an important role in the phase transitions of the system . in the limit of vanishing josephson coupling strength ,
we have performed the mc simulations to obtain @xmath77 as functions of the charge frustration at various temperatures , and pointed out the possibility of the devil s staircase structure in the thermodynamic limit .
as the temperature is increased , higher - order insulating states disappear but some lower - order ones still remain at sufficiently low temperatures . in the presence of the josephson coupling , we have obtained the effective classical hamiltonian via a variational method in the path integral formalism .
a self - consistent approximation has been applied to the resulting hamiltonian and phase boundaries have been computed for @xmath293 and 1 .
a reentrant transition at @xmath301 has been observed for the array in the self - charging limit ( @xmath302 ) , while an insulating lobe - like structure has been found near @xmath303 in the case @xmath304 .
the validity of the formalism used here , the self - consistent approximation ( sca ) on the effective classical hamiltonian , has been checked using mc simulations in the absence of charge frustration . although the reentrant behavior of the array in the self - charging limit has also been observed in mc simulations , the sca appears to be accurate only at sufficiently high temperatures .
the quantum phase transitions of the superconducting arrays have been studied theoretically by many authors via mean - field approximations such as the coarse - graining approximation  @xcite and via the semiclassical methods including the variational methods  @xcite and the wkb approximation  @xcite .
although the mean - field approximations are in general expected to give better results at zero temperature , where the dimension of the system is effectively increased by one , than at finite temperatures , it is well known that the coarse - graining approximation gives unphysical results even at zero temperature , for the arrays without self - capacitances . on the other hand ,
the semiclassical methods should be applied with great care at zero temperature , where quantum fluctuations are strong . in comparison with these methods , the formalism adopted here to study the zero - temperature phase transitions , perturbation expansion combined with simulated annealing , appears to be more reliable and
can be applied to other systems systematically .
we thank g.s .
jeon for useful discussions .
this work was supported in part by the basic science research institute program , ministry of education , and in part by the korea science and engineering foundation through the src program .
99 for a list of references , see , e.g. , _ proceedings of the 2nd ctp workshop on statistical physics : kt transition and superconducting arrays _ , edited by d. kim , j.s .
chung , and m.y .
choi ( min - eum sa , seoul , 1993 ) ; physica b * 222 * 253 - 406 ( 1996 ) .
berezinskii , zh .
. fiz . * 59 * , 907 ( 1970 ) [ sov .
jetp * 32 * , 493 ( 1971 ) ] ; j.m .
kosterlitz and d.j .
thouless , j. phys .
c * 6 * , 1181 ( 1973 ) ; j.m .
kosterlitz , j. phys .
c * 7 * , 1046 ( 1974 ) ; j.v .
jos , l.p .
kadanoff , s. kirkpatrick , and d.r .
nelson , phys .
b * 16 * , 1217 ( 1977 ) .
geerligs , m. peters , l.e.m .
de groot , a. verbruggen , and j.e .
mooij , phys .
lett . * 63 * , 326 ( 1989 ) ; h.s.j .
van der zant , l.j .
geerligs , and j. e. mooij , europhys .
* 19 * , 541 ( 1992 ) .
kim and m.y .
choi , phys . rev .
b * 52 * , 3624 ( 1995 ) . c. rojas and v. jos , phys . rev .
b * 54 * , 12361 ( 1996 ) . c. bruder , r. fazio , a. kampf , a. van otterlo , and g. schn , phys . scr .
t * 42 * , 159 ( 1992 ) ; a. van otterlo , k .- h .
wagenblast , r. fazio , and g. schn , phys .
b * 48 * , 3316 ( 1993 ) . m.p.a .
fisher , p.b .
weichman , g. grinstein , and d.s .
fisher , phys .
b * 40 * , 546 ( 1989 ) ; g.g .
batrouni , r.t .
scalettar , and g.t .
zimanyi , phys .
lett . * 65 * , 1765 ( 1990 ) .
freericks and h. monien , europhys .
26 * , 545 ( 1994 ) .
meisel , physica b * 178 * , 121 ( 1992 ) ; g.g .
batrouni , r.t .
scalettar , g.t .
zimanyi , and a.p .
kampf , phys .
lett . * 74 * , 2527 ( 1995 ) ; a. van otterlo and k .- h .
wagenblast , _ ibid . _
* 72 * , 3598 ( 1994 ) ; e. roddick and d. stroud , phys .
b * 51 * , 8672 ( 1995 ) .
a. van otterlo , k .- h .
wagenblast , r. baltin , c. bruder , r. fazio , and g. schn , phys .
b * 52 * , 16176 ( 1995 ) .
keldysh , pisma zh .
* 29 * , 716 ( 1979 ) [ jetp lett .
* 29 * , 658 ( 1979 ) ] ; r. fazio and g. schn , phys .
b * 43 * , 5307 ( 1991 ) ; j.e .
mooij , b.j .
van wees , l.j .
geerligs , m. peters , r. fazio , and g. schn , phys .
rev . lett . * 65 * , 645 ( 1990 ) .
cha , m.p.a .
fisher , s.m .
girvin , m. wallin , and a.p .
young , phys .
b * 44 * , 6883 ( 1991 ) .
n. goldenfeld , _ lectures on phase transitions and the renormalization group _
( addison - wesley , 1992 ) , p. 44 . h.s.j .
van der zant , l.j .
geerligs , and j.e .
mooij , europhys . lett . * 19 * , 541 ( 1992 ) .
a. kanda and s. kobayashi , j. phys .
. jpn . * 64 * , 3172 ( 1995 ) ; a. kanda , s. katsumoto , and s. kobayashi , _ ibid . _
* 63 * , 4306 ( 1994 ) ; r. yagi , t. yamaguchi , h. kazawa , and s. kobayashi , czech .
* 46 * , suppl .
s2 , 693 ( 1996 ) .
lobb , d.w .
abraham , and m. tinkham , phys .
b * 27 * , 150 ( 1983 ) .
straley and g.m .
barnett , phys .
b * 48 * , 3309 ( 1993 ) .
a.a . odintsov , phys .
b * 54 * , 1228 ( 1996 ) .
j. lee , phys .
lett . * 71 * , 211 ( 1993 ) ; j. lee and m.y .
choi , phys .
e * 50 * , r651 ( 1994 ) .
choi and s. doniach , phys .
b * 31 * , 4516 ( 1985 ) .
p. bak , rep .
. phys . * 45 * , 587 ( 1982 ) .
r. giachetti and v. tognetti , phys .
lett . * 55 * , 912 ( 1985 ) ; r.p . feynman and h. kleinert ,
a * 34 * , 5080 ( 1986 ) .
s. kim and m.y .
choi , phys .
b * 41 * , 111 ( 1990 ) ; * 42 * , 80 ( 1990 ) .
r.p . feynman and a.r .
hibbs , _ quantum mechanics and path integrals _ ( mcgraw - hill , new york , 1965 ) . g. schn and a.d .
zaikin , phys .
rep . * 198 * , 237 ( 1990 ) .
see , e.g. , s. oh , s. lee , and k .- c .
lee , j. korean phys .
* 27 * , 365 ( 1994 ) and references therein .
nelson and j.m .
kosterlitz , phys .
lett . * 39 * , 1201 ( 1977 ) .
.values of @xmath172 in the lowest excited states of the arrays with size ( a ) @xmath221 and ( b ) @xmath225 .
it is shown that @xmath305 tends to decrease as @xmath1 is increased . [
cols=">,^,^,^",options="header " , ]