^{-2};q)_{\infty}
      [(1+z_1)(1+z_1^{-1})]^{j-1}
      [(1+z_2)(1+z_2^{-1})]^{L-j} }
     {(az_1,a z_1^{-1},q z_1 z_2, q z_1^{-1} z_2^{-1},
       q z_1 z_2^{-1},q z_1^{-1} z_2^{-1}, bz_2,b z_2^{-1};q)_{\infty}
      (z_2 + z_2^{-1} -z_1-z_1^{-1})}.\end{aligned}$$ When $0<a,b<1$, the contours of $z_1$ and $z_2$ are both unit circles. For other valures of the parameters, the contours are deformed as poles in the integrands move in and out of the unit circles.

The evaluation of the integrals $I_1$ and $I_2$ in the thermodynamic limit for each phase is relegated to Appendices. Basically, the integral $I_1$ gives the density at bulk region whereas the integral $I_2$ gives the density near the right boundary. The results are summerized in the following. We thus obtain the phase diagram shown in Fig. 5. Notice that this phase diagram was corrrectly predicted in [@me99]. Therein the correlation length for each phase was assumed to be given by the logarithm of the ratio of the largest and the second largest eigenvalue of the matrix $C$.

-   Phase $C$ ($\tilde{\alpha} > 1/2$ and $\tilde{\beta}>1/2$; $0<a,b<1$)

    In this phase, the average density at bulk is $1/2$. The average density decays near the right bounary as $$\label{density-max}
      \langle n_j \rangle_L
      =
      \frac12 
      -
      \frac{1}{2\sqrt{\pi} l^{\frac12}}.$$ Here and in the following, we set $l=L-j+1$. The density decays algebraically and hence the correlation length is infinite. This average density profile is exactly the same as that for the totally asymmetric case. The density decay at the left boundary can
^{-2};q)_{\infty } 
       [ (1+z_1)(1+z_1^{-1})]^{j-1 } 
       [ (1+z_2)(1+z_2^{-1})]^{L - j } } 
      { (az_1,a z_1^{-1},q z_1 z_2, q z_1^{-1 } z_2^{-1 }, 
        q z_1 z_2^{-1},q z_1^{-1 } z_2^{-1 }, bz_2,b z_2^{-1};q)_{\infty } 
       (z_2 + z_2^{-1 } -z_1 - z_1^{-1})}.\end{aligned}$$ When $ 0 < a, b<1 $, the contours of $ z_1 $ and $ z_2 $ are both unit circles. For early valures of the argument, the contours are deformed as pole in the integrands go in and out of the unit circles. 

 The evaluation of the integral $ I_1 $ and $ I_2 $ in the thermodynamic limit for each phase is relegated to Appendices. Basically, the built-in $ I_1 $ give the density at bulk area whereas the integral $ I_2 $ gives the density near the right limit. The results are summerized in the pursuit. We thus receive the phase diagram read in Fig. 5. Notice that this phase diagram was corrrectly predicted in [ @me99 ]. Therein the correlation distance for each phase was assumed to be given by the logarithm of the ratio of the largest and the second large eigenvalue of the matrix $ C$. 

 -    Phase $ C$ ($ \tilde{\alpha } > 1/2 $ and $ \tilde{\beta}>1/2 $; $ 0 < a, b<1 $) 

     In this phase, the average density at bulk is $ 1/2$. The average density decay near the right bounary as $ $ \label{density - max } 
       \langle n_j \rangle_L 
       = 
       \frac12 
       - 
       \frac{1}{2\sqrt{\pi } l^{\frac12}}.$$ Here and in the following, we set $ l = L - j+1$. The concentration disintegrate algebraically and hence the correlation length is infinite. This average density profile is exactly the same as that for the wholly asymmetrical case. The density decay at the left limit can
^{-2};q)_{\invty}
      [(1+z_1)(1+z_1^{-1})]^{j-1}
      [(1+z_2)(1+z_2^{-1})]^{L-j} }
     {(az_1,a z_1^{-1},q z_1 z_2, q z_1^{-1} z_2^{-1},
       q z_1 z_2^{-1},q z_1^{-1} x_2^{-1}, bz_2,b z_2^{-1};q)_{\infty}
      (z_2 + z_2^{-1} -z_1-z_1^{-1})}.\env{alitned}$$ Qhen $0<a,b<1$, the contours ow $z_1$ and $z_2$ are boty unmt circles. For ovger valmxes or the 'arameters, the gontours ara deformed as [oued in the integrands move in and out of the ujit circles.

The evaktatikn of the integrals $I_1$ and $I_2$ in the thermovynamic limit fpr each phase is relegated to Wppendices. Basicalpy, the intetral $U_1$ gives the aensity at bulk region whereas the integral $I_2$ gives tfe deusity near jkw rlcht boundarb. The gesults are smkmerizad in tne following. Wc thuv ovtain the phase diagrem shown in Fig. 5. Notyce that dhns phase diagram was xoerrecjly psedizred in [@mx99]. Tgerein thx correlatikn length fir each phase was axstned to be givsn by ehq logarithm of the ratio of the largest ans the second largest eitenvalue of the matrid $C$.

-   Phafe $C$ ($\tilde{\alpha} > 1/2$ and $\tilde{\beta}>1/2$; $0<a,b<1$)

    In this phave, thx xvexqge ddbslty at bulk is $1/2$. The average density decays newd uhe right bounary as $$\label{demsltu-iax}
      \langue n_j \xznfle_L
      =
      \frwc12 
      -
      \frac{1}{2\sqwt{\pi} l^{\frac12}}.$$ Here and in the folliwing, we set $o=L-j+1$. The density deeays algebrancally and nence the correlation lzngth js infinite. This avedxge density profkle iv exactly the same as that for the votalky asymoetroc casq. The denslty dcway at the left bokndari can
^{-2};q)_{\infty} [(1+z_1)(1+z_1^{-1})]^{j-1} [(1+z_2)(1+z_2^{-1})]^{L-j} } {(az_1,a z_1^{-1},q z_1 z_1^{-1} q z_1 z_1^{-1} z_2^{-1}, bz_2,b When the contours of and $z_2$ are unit circles. For other valures of parameters, the contours are deformed as poles in the integrands move in and of the unit circles. The evaluation of the integrals $I_1$ and $I_2$ in thermodynamic for phase relegated to Appendices. Basically, the integral $I_1$ gives the density at bulk region whereas the integral gives the density near the right boundary. The are summerized in the We thus obtain the phase shown Fig. 5. that phase was corrrectly predicted [@me99]. Therein the correlation length for each phase was assumed to be given by the logarithm of ratio of and the largest of matrix $C$. - ($\tilde{\alpha} > 1/2$ and $\tilde{\beta}>1/2$; $0<a,b<1$) the average density at bulk is $1/2$. The density decays the right bounary as $$\label{density-max} \langle \rangle_L = \frac12 - \frac{1}{2\sqrt{\pi} l^{\frac12}}.$$ Here and the following, we set $l=L-j+1$. The density decays algebraically and hence the correlation length is average density profile is the same as for totally case. density decay the left boundary can
^{-2};q)_{\infty}
      [(1+z_1)(1+z_1^{-1})]^{j-1}
      [(1+z_2)(1+z_2^{-1})]^{L-j} }
     {(az_1,a z_1^{-1},q z_1 z_2, q z_1^{-1} z_2^{-1},
       q z_1 Z_2^{-1},q z_1^{-1} z_2^{-1}, bz_2,b z_2^{-1};q)_{\iNfty}
      (z_2 + Z_2^{-1} -z_1-z_1^{-1})}.\End{AlIgneD}$$ WheN $0<a,b<1$, the contours OF $z_1$ anD $z_2$ are both unit circles. FoR otheR vALureS Of The paRameterS, ThE COntOuRs Are DeFOrMed as PolEs in the Integrands MovE iN and out of the UNiT circles.

ThE evAluation of thE inTegralS $I_1$ And $i_2$ In the TheRmodyNamic lIMit for Each phase Is RElegatED to AppeNDIcEs. BaSically, the integraL $i_1$ gIVes the density aT bulk rEgIOn WHEreAs tHe integral $i_2$ gIves tHE densitY NeAR THe rIGht boundary. ThE results are SUmmErized In The FOllowiNg. We tHuS ObtAin the phase DiagRam shown iN Fig. 5. NoTIce that THis phasE diagrAm wAs cOrrrECtLy PreDiCTed IN [@mE99]. ThEReiN the corrElAtIon leNgth FOR EAch pHasE was AssumEd to be given by The LogaRIthM of thE ratiO of tHe LargeSt and tHe secOnD largest eigenvaLue oF the matriX $C$.

-   PHaSe $C$ ($\TiLde{\alPHa} > 1/2$ and $\tIldE{\beTa}>1/2$; $0<a,b<1$)

    In tHis phasE, The AvERAGe Density at bulk is $1/2$. The AvERAgE density Decays NEaR tHE right boUnAry As $$\laBEL{densIty-mAX}
      \lAngle n_j \rAngle_L
      =
      \FRaC12 
      -
      \fRac{1}{2\sqrt{\Pi} L^{\frac12}}.$$ HErE anD in The foLLowiNg, we seT $l=L-j+1$. The dEnsitY Decays algebraiCAlly and hence tHE cORReLAtioN leNgth is infinIte. THIs avEragE DeNsiTY profIle is ExACtLY the same as that for thE tOtally AsymmEtric case. The dEnsity decaY AT The left bOundARy CAn
^{-2};q)_{\infty}
      [( 1+z_1)(1+z _1^{- 1}) ]^{ j- 1}
    [(1+z_2)(1+z_ 2 ^{-1 })]^{L-j} }
     {(az_ 1,a z _1 ^ {-1} , qz_1 z _2, q z _ 1^ { - 1}z_ 2^ {-1 },      q z _1 z_2^ {-1},q z_1 ^{- 1}  z_2^{-1}, b z _2 ,b z_2^{-1 };q )_{\infty}
    (z_2  +  z_ 2 ^{-1}  -z _1-z_ 1^{-1} ) }.\end {aligned} $$ When $ 0 <a,b<1$ , th e co ntours of $z_1$ a n d$ z_2$ are bothunit c ir c le s .  Fo r o ther valur es  of t h e param e te r s ,  th e  contours are  deformed a s  po les in  t hei ntegra nds m ov e  in  and out of  the  unit cir cles.The eva l uationof the  in teg rals $I _1 $ a nd $I_ 2 $int hethermody na mi c lim it f o r e achpha se i s rel egated to App end ices .  Ba sical ly, t he i nt egral  $I_1$  give sthe density atbulk  region w her ea s t he  inte g ral $I _2$  gi ves the  densit y  ne ar t h eright boundary. Th er e su lts aresummer i ze di n the fo ll owi ng.W e  thus  obt a in  the pha se dia g ra mshown i nFig. 5 .Not ice  that this  phase  diagram  wasc orrrectly pred i cted in [@me9 9 ]. T he r einthe  correlatio n le n gthfore ac h p h ase w as as su m ed to be given by thelo garith m ofthe ratio ofthe larges t a nd the s econ d  l a rgest eigenval ue of  the matri x  $C$.

-    Ph ase $C$($\tilde{ \ a lpha} >1/2 $ a nd$\t i l de {\beta}>1/2$; $ 0<a, b< 1$)

   I n thispha se,  th e a ve rage dens ity at b ul kis  $ 1/2 $. Th e  average  d ens it y d ecays near t he ri ghtbo un a ryas $$\l a be l { dens it y- max} 
    \la ngle n_j  \rangl e_L
     =    \ fr ac12 
     -
     \f rac{1}{2\s qr t{\ pi} l^ { \ frac12}} .$$ Here and in the fol l owing,weset $ l=L- j+1$. The  de nsitydec a ys alg ebraic allyan d h e n ce th e co rre la tion lengt h isinfin it e. T his ave rage density profi l e i s exactly the  sa me a s th atf or the  t o tal l y  asymmetric cas e. The den si t ydecay at t h e l ef t bound ary can 
^{-2};q)_{\infty}
 _ _  [(1+z_1)(1+z_1^{-1})]^{j-1}
 _ __ [(1+z_2)(1+z_2^{-1})]^{L-j}_}
_  _ {(az_1,a z_1^{-1},q_z_1 z_2, q z_1^{-1}_z_2^{-1},
  __   q z_1 z_2^{-1},q z_1^{-1} z_2^{-1}, bz_2,b z_2^{-1};q)_{\infty}
    __(z_2 +_z_2^{-1}_-z_1-z_1^{-1})}.\end{aligned}$$_When $0<a,b<1$, the contours of_$z_1$ and $z_2$ are both_unit circles._For other valures of the parameters, the contours_are_deformed as poles_in the integrands move in and out of the_unit circles.

The evaluation of the integrals_$I_1$ and $I_2$_in_the_thermodynamic limit for each_phase is relegated to Appendices. Basically,_the integral $I_1$ gives the density_at bulk region whereas the integral $I_2$_gives the density near the right_boundary. The results are summerized_in the_following. We thus obtain the_phase diagram shown_in Fig._5. Notice that_this phase diagram was corrrectly predicted_in [@me99]. Therein_the correlation length for each phase_was_assumed to be_given_by_the logarithm_of the ratio_of_the largest_and_the second largest eigenvalue of the_matrix_$C$.

-   Phase $C$ ($\tilde{\alpha} >_1/2$ and $\tilde{\beta}>1/2$; $0<a,b<1$)

__  In this_phase, the average density at_bulk is $1/2$. The average density_decays near_the right_bounary as $$\label{density-max}
      \langle n_j \rangle_L
_     =
_   _ \frac12_
_  __ -
_     \frac{1}{2\sqrt{\pi} l^{\frac12}}.$$_Here and_in the following, we set $l=L-j+1$._The density decays algebraically_and_hence the correlation length is infinite._This average density profile is exactly_the same as that for_the_totally_asymmetric case. The density decay_at the left boundary can