Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: during the past decade it has become abundantly clear that there are several sub - classes of two - dimensional ( 2d ) spin liquids . @xcite perhaps best understood are the topological spin liquids that support gapped excitations which carry fractional quantum numbers . @xcite another possibility are gapless spin liquids with no topological structure . these gapless spin liquids will generically exhibit spin correlations that decay as a power law in space and which can oscillate at particular wavevectors . in such  algebraic " or  critical " spin liquids @xcite these wavevectors are limited to a finite discrete set , and their effective field theories will often exhibit a relativistic structure  as is the case when the spinons in a slave particle construction are massless dirac fermions coupled to a u(1 ) gauge field . however , it is possible that the spin correlations exhibit singularities along surfaces in momentum space . these are analogous to the fermi surfaces in a fermi liquid but in such  quantum spin - metals " a quasiparticle picture will presumably be inapplicable . spin liquids with a spinon fermi surface are examples of particular interest , and have been studied in a number of rvb works @xcite  most recently as a candidate for a spin liquid phase @xcite in the organic compound @xmath1-(et)@xmath2cu@xmath2(cn)@xmath3 .    developing a theoretical framework for itinerant 2d non - fermi liquids is arguably even more challenging than for spin liquids . towards this end recent work@xcite explored the possibility of an uncondensed quantum phase of itinerant bosons in two dimensions which is a conducting fluid but not superfluid . because of the characteristic nodal structure observed by moving one boson around another , ref .   called such phases `` d - wave correlated bose liquids '' ( dbl ) . various physical correlators were found to be singular across surfaces in momentum space in the dbl . for example , the boson occupation number @xmath4 is singular at some `` bose surfaces '' like those illustrated in fig . [ fig : qycuts ] . other examples with critical surfaces have been studied recently in ref .  . the central idea underlying this paper is that 2d phases of quantum spins and itinerant fermions or bosons with singular surfaces in momentum space should have definite signatures when restricted to a quasi one - dimensional ( 1d ) geometry , e.g.  when the system is placed onto an @xmath0-leg ladder . to be concrete , suppose we have a singular 2d surface in momentum space in the ground state of some square lattice quantum hamiltonian . if we put the hamiltonian on an @xmath0-leg ladder , the transverse momentum @xmath5 becomes discrete . with large @xmath0 one would expect that the short - range energetics that stabilizes the 2d phase will still be present on the ladder and will lead to a set of 1d gapless modes where the discrete momenta cut through the singular surface , as illustrated in fig . [ fig : qycuts ] . the long wavelength physics of these 1d modes will be described in terms of some multi - channel luttinger liquid , behavior different than the parent 2d state . but the number of gapless 1d modes and their respective momenta characterizing the oscillatory power law decays would provide a distinctive  fingerprint " of the parent 2d liquid . for example , the number of gapless 1d modes would grow linearly with @xmath0 and the sum of their momenta could satisfy some generalized luttinger volume . moreover , these signatures could be potentially exploited to approach the 2d phases from controlled numerical and analytical studies in quasi-1d . indeed , as we demonstrate explicitly in this paper , the ( putative ) singular surface of the 2d dbl phase is already manifest on the 2-leg ladder !     cuts through bose surfaces in some anisotropic 2d system placed on an @xmath0-leg ladder ; here for illustration @xmath6 . such bose surfaces appear in a dbl construction described in the text ; to avoid clutter , we do not show the fermi surfaces of the @xmath7 and @xmath8 slave particles . ] -2 mm    the itinerant uncondensed dbl of ref .   is constructed by writing a boson in terms of fermionic partons , @xmath9 and considering un - coupled mean field states of the @xmath7 and @xmath8 . the dbl phase on a square lattice is obtained by taking the @xmath7 fermions to hop preferentially along the @xmath10-axis and the @xmath8 fermions to hop preferentially along the @xmath11-axis . one can enforce the square lattice symmetry in the boson state by requiring that the @xmath7 and @xmath8 hopping patterns are related by a 90 degree rotation . to recover the physical hilbert space requires the condition @xmath12 at each site . within a gutzwiller projection approach the corresponding boson wavefunction is taken as product of two fermion slater determinants , @xmath13 in the dbl the @xmath7 and @xmath8 fermi surfaces are different , compressed in the @xmath10 and @xmath11-directions , respectively . this results in a characteristic d - wave - like nodal structure observed when moving one boson around another in @xmath14 . the boson momentum distribution function @xmath4 in the dbl is singular on two surfaces that are constructed from the @xmath7 and @xmath8 fermi surfaces as enveloping @xmath15 surfaces . more details are in ref .  , while fig . [ fig : qycuts ] shows one example of the @xmath15 loci for an open @xmath7 fermi surface and a closed @xmath8 fermi surface ( the fermi surfaces themselves are not shown ; this dbl state would be anisotropic on the square lattice but is closer to the ladder states considered here where there is no 90 degree rotation symmetry ) . as an alternative to the gutzwiller construction one can project into the physical hilbert space within a gauge theory approach , coupling @xmath7 and @xmath8 with opposite gauge charges to an emergent u(1 ) gauge field .    ref .   also proposed a frustrated boson model with ring exchanges that can potentially realize such dbl phases : @xmath16 with @xmath17 . in addition to the usual boson hopping term this hamiltonian contains a four - site ring term acting on each square plaquette . with @xmath18 positive this term is  frustrating " , violating the marshall sign rule of the hopping term  making the system intractable by quantum monte carlo simulations . this boson hamiltonian was constructed@xcite by taking the strong coupling limit of the lattice gauge theory for the @xmath7 and @xmath8 fermions . increasing the disparity in the fermion hopping along the @xmath10 and @xmath11 axes corresponds to an increase in @xmath18 .    in this paper , we make a first step in the proposed program of ladder studies by exploring the ring model eq . ( [ hring ] ) on a 2-leg ladder . of course , the picture of discrete @xmath5 cuts through two - dimensional surfaces is pushed to the extreme here . nevertheless and quite remarkably , in a wide parameter regime of intermediate to large ring couplings , we find an unusual phase which can be understood via the dbl construction eq . ( [ psidbl ] ) . this is a strong - coupling phase with no perturbative picture in terms of the original bosons . the slave - particle approach provides a new starting point , and the resulting gauge theory of the dbl can be solved by the conventional 1d bosonization tools , @xcite providing a consistent picture of the new phase . the paper is organized as follows . section  [ sec : dbl_2leg ] sets the stage by describing dbl states on the 2-leg ladder . the main section  [ sec : jk ] presents results of exact numerical studies of the @xmath19 model . the focus here is on the dbl phase that emerges when the conventional `` @xmath20 liquid '' is destroyed by the ring exchanges . however , the phase diagram is even more rich and contains two more unusual phases in which bosons are paired ; one of these phases can be accessed analytically as an instability of yet a different dbl phase . appendix  [ app : gt ] summarizes the technical gauge theory description of the dbl phases while appendix  [ app : oddends ] presents simple analytical results for the ring model in some limiting cases . section  [ sec : concl ] concludes with an outline of future program of ladder studies . the dbl construction eq . ( [ psidbl ] ) proceeds by taking distinct hopping ground states for the @xmath7 and @xmath8 fermions . the hopping patterns of the two fermion species can be independent and need to only respect the ladder symmetries . there is a bonding and an antibonding band for each species . we take the @xmath7 fermions to hop more strongly along the chains and always assume that the chemical potential crosses both bands as illustrated in fig . [ fig : bands ] . we take the @xmath8 fermions to hop more strongly on the rungs and consider two cases shown in the bottom panels in fig . [ fig : bands ] : the  dbl[2,1 ] phase " when only the bonding band of the @xmath8 fermions crosses the fermi level , and the  dbl[2,2 ] phase " when both @xmath8 bands are partially occupied  the first and second argument in the square brackets refer to the number of involved @xmath7 and @xmath8 bands respectively .     and @xmath8 bands for the dbl construction on the 2-leg ladder . top panel : @xmath7 fermions partially occupy both the bonding and anti - bonding bands . bottom left : in the dbl[2,1 ] phase , @xmath8 fermions occupy only the bonding band . bottom right : in the dbl[2,2 ] phase , @xmath8 fermions partially occupy both bands ; this state is likely unstable towards a boson - paired phase , as discussed in sec . [ subapp : dbl22instab ] . in the wavefunction construction , only the occupied orbitals and not the band dispersions are important . ] -2 mm    let us first discuss properties of the dbl[2,1 ] phase . we denote the fermi momenta of the @xmath7 bonding and antibonding bands as @xmath21 and @xmath22 , corresponding to the transverse momenta @xmath23 and @xmath24 respectively ( which is a convenient alternative to the quantum number for the leg interchange symmetry ) . we denote the fermi momentum of the @xmath8 bonding band as @xmath25 ( there is no @xmath8 antibonding occupation in the dbl[2,1 ] ) . the fermi momenta satisfy @xmath26 , where @xmath27 is the boson density per site and remembering that there are two sites in each rung while the lattice spacing is set to one .    in a mean field treatment , the boson green s function is a product of two fermion green s functions and therefore oscillates at different possible @xmath28 wavevectors and decays as @xmath29 . going beyond the mean field as described in appendix  [ app : gt ] , we expect that the oscillations at wavevectors @xmath30 have slower power law decay while the other two wavevectors @xmath31 have faster power law decay than the mean field . the complete result is , dropping all amplitudes , @xmath32}{x^{2 - \eta^{(0 ) } } } + e^{i\pi y } \frac{\cos[(k^\pd_{f2 } - k_{f1}^{(\pi ) } ) x ] }                    { x^{2 - \eta^{(\pi)}}}\\ + \frac{\cos[(k^\pd_{f2 } + k_{f1}^{(0 ) } ) x]}{x^{2 + \dots } } + e^{i\pi y } \frac{\cos[(k^\pd_{f2 } + k_{f1}^{(\pi ) } ) x]}{x^{2 + \dots}}.\end{aligned}\ ] ] the first line shows the enhanced correlators . these can be guessed using an `` amperean rule '' mnemonic inherited from the fermion - gauge studies in 2d and verified by the bosonization solution in appendix  [ app : gt ] : composites that involve fermion bi - linears with opposite group velocities but that produce parallel gauge currents are enhanced by the gauge fluctuations . thus , low - energy `` bosons '' @xmath33 carrying @xmath34 are expected to be enhanced . the contributions in the second line , on the other hand , can be shown to be always suppressed beyond the mean field . appendix  [ app : gt ] presents a complete low - energy theory of the dbl[2,1 ] which  once the gauge fluctuations are treated  contains two phonon modes and in principle allows calculating all exponents . since the power laws are affected also by short - range interactions , and there are many parameters allowed in the theory , we do not pursue this calculation in general , but mainly rely on the amperean rules that capture the gauge fluctuation effects ; this is reasonable if the gauge interactions dominate . the boson singularities at non - zero wavevectors become manifest upon fourier transformation of @xmath35 to obtain the boson occupation number @xmath4 . for example , at the enhanced wavevectors @xmath36 and @xmath37 , we have @xmath38    the boson density correlations are analyzed similarly . in a naive mean field , @xmath39 , which oscillates at various `` @xmath40 '' vectors and decays as @xmath29 . beyond the mean field , the power laws get modified : @xmath41}{x^{2 - \dots } } + \frac{\cos[2 k_{f1}^{(\pi ) } x]}{x^{2 - \dots } } \\ & & + e^{i\pi y}\frac{\cos[(k_{f1}^{(0 ) } + k_{f1}^{(\pi ) } ) x]}{x^{2 - \dots } } + e^{i\pi y}\frac{\cos[(k_{f1}^{(0 ) } - k_{f1}^{(\pi ) } ) x]}{x^{2 + \dots } } ~ , \nonumber \\ & & d_{d2}({\bf r } ) \sim \frac{1}{x^2 }   + \frac{\cos[2 k_{f2 } x]}{x^{2 - \dots } } ~.\end{aligned}\ ] ] the distinct momenta are : @xmath42 , @xmath43 , @xmath44 , @xmath45 , @xmath46 , and @xmath47 . the first four involve a particle and a hole of the same species but with opposite group velocities and are therefore expected to be enhanced by the gauge fluctuations ; in the above equations , this is indicated schematically with @xmath48 . however , these correlators can also be affected  either enhanced or suppressed  by short - range interactions , and without knowing all parameters in the theory , we can not calculate the exponents reliably . the oscillation at wavevector @xmath46 can be shown to be always suppressed beyond the mean field . finally , the zero momentum power law remains unchanged . more details on the dbl[2,1]-phase can be found in appendix  [ app : gt ] . in particular , we show that the two - boson green s function exhibits some internal d - wave character , which originates from the non - trivial wavefunction signs . it is these signs that prompted the name `` d - wave '' boson liquid ( dbl ) in the 2d continuum setting,@xcite where the wavefunction goes through a sequence of signs @xmath49 upon taking one particle around another . however , such correlations do not necessarily mean that the system is near a d - wave - paired phase . for example , we also examine a potential instability of the dbl[2,1 ] driven by an allowed non - marginal four - fermion term and find that the resulting phase is a boson - paired liquid with an internal s - wave character . let us briefly mention the dbl[2,2 ] case , where in the mean field both @xmath7 and @xmath8 have partial occupations of both the bonding and antibonding bands . one can calculate various correlations in the mean field as before , and the structure is more rich since there is an additional band present . some details are given in appendix  [ subapp : dbl22instab ] . beyond the mean field but ignoring non - marginal four - fermion terms , this phase would have three gapless modes . however , an analysis of allowed interactions suggests that this phase has a strong instability towards a boson - paired phase with an s - wave character . remarkably , as we will see in the next section , a boson - paired phase with similar properties is found in the ring model in the regime of small interchain coupling @xmath50 where we initially hoped to find the dbl[2,2 ] phase . thus , the slave - particle formulation and the gauge theory description solved via subsequent bosonization open a non - perturbative access to this interesting boson - paired phase . the symbols represent dmrg estimates with shaded grey circles indicating the superfluid phase , the filled blue circles the d - wave correlated bose liquid ( dbl ) , and the open red diamonds the s - wave paired state . the schematic line boundaries are obtained from vmc calculations . ] we study the ring model eq . ( [ hring ] ) on a 2-leg ladder with a boson hopping term @xmath51 that allows different amplitudes along and perpendicular to the chains : @xmath52 where @xmath53 represent integer lattice sites . we take @xmath54 with @xmath55 and @xmath56 with @xmath57 the length of the 2-leg system . the ring terms are associated with the square plaquettes of the 2-leg ladder . we use exact diagonalization ( ed ) and density matrix renormalization group ( dmrg)@xcite methods supplemented with trial wavefunction variational monte carlo ( vmc)@xcite to determine the nature of the ground state of the hamiltonian eq . ( [ hring ] ) at boson filling number @xmath58 , where @xmath59 is the total number of bosons in the system . the obtained phase diagrams at densities @xmath60 and @xmath61 are shown in fig . [ fig : phased1 ] and fig . [ fig : phased2 ] respectively . these are typical results for the boson ring model away from half - filling .    using ed and dmrg , we find four different quantum phases upon increasing the ring - exchange coupling @xmath18 and varying the interchain coupling @xmath50 . at small @xmath18 , the `` superfluid '' phase ( sf ) with @xmath20 quasi - long - range order ( qlro ) is stable for both filling numbers . the dbl phase , which is produced by the ring exchanges , dominates the intermediate parameter space at higher filling number @xmath60 , while it occupies a substantially smaller region at @xmath62 . interestingly , at @xmath60 and small @xmath50 , an s - wave paired phase emerges between the sf and dbl as a consequence of the competition between the boson hopping and ring - exchange terms . at @xmath62 and small @xmath50 , a d - wave paired phase is found adjacent to the dbl . finally , for strong ring exchanges phase separation eventually wins near @xmath63 for both fillings , separating into a region with @xmath64 and an empty region @xmath65 . the characteristic features of each phase will be discussed in the following sections .    . the symbols indicate dmrg estimates for the various phases with shaded grey circles for the superfluid phase , filled blue circles for the dbl , and open green diamonds for the d - wave paired state . the schematic line boundaries are obtained from vmc calculations . ]    our complementary vmc study is centered around the dbl phase . the dbl wavefunction eq . ( [ psidbl ] ) is constructed from the slave - particle formulation and thus allows us to compare between the exact ed / dmrg results and the gauge theory description . jastrow sf and boson - paired trial wavefunctions ( described in appendix  [ app : psipaired ] ) are also considered in order to better understand the exact numerical results . we note that this variational approach has been rather successful in one dimension , e.g. , in the original studies of the 1d @xmath66 model . @xcite we borrow some ideas from these studies such as allowing the fermionic determinants to also play a jastrow role via @xmath67 , where @xmath68 is a variational parameter . in our ring model , the phase boundary estimated from vmc almost coincides with the one determined from dmrg as indicated in fig . [ fig : phased1 ] for @xmath60 filling . at smaller filling number , dmrg and vmc also find the same topology of the phase diagram as given in fig . [ fig : phased2 ] . we first describe measurements in the ground state . imposing periodic boundary conditions in the @xmath10-direction in our numerical calculations , we can characterize the ed states by a momentum @xmath69 . in our dmrg calculations we work with real - valued wavefunctions which gives no ambiguity when the system exhibits a unique ground state that caries zero momentum . if the ground state carries nontrivial momentum @xmath69 , its time - reversed partner carries @xmath70 , and the dmrg state is some real - valued combination of these . while the lattice - space measurements may depend slightly on the specific combination ( but vanishingly in the large system limit ) , the momentum space measurements @xmath4 and @xmath71 described below are unique and insensitive to this . in our dmrg calculations , more than 1500 states were kept in each block@xcite to ensure accurate results , and the density matrix truncation error is of the order of @xmath72 . the typical error for the ground - state energy is of the order of @xmath73 for all systems we have studied . the relative error in the correlation functions varies between @xmath74 , depending on the type of the correlations and the spatial distance between operators . we measure boson correlator in the ground state @xmath75 and its fourier transform @xmath76 which is interpreted as mode @xmath77 occupation number . we also measure the boson density correlator @xmath78 where @xmath79 is the average boson density . the density structure factor is @xmath80 for both the boson or boson density correlations , a power law @xmath81 envelope in real space gives rise to a singularity @xmath82 in momentum space .    finally , we measure two - boson ( pairing ) correlator , @xmath83 which creates a pair of bosons at @xmath84 and removes a pair at @xmath85 . this is useful for detecting boson - paired phases where the single boson correlator decays exponentially . in the pairing correlator figures below , we present @xmath86 for a fixed created @xmath87 diagonal pair @xmath88 near the origin , while the removed pair is either a @xmath87 diagonal @xmath89 or a @xmath90 diagonal pair @xmath91 . we have measured the correlator for other pair orientations as well , but find that the diagonal - diagonal correlations are the most distinguishing ones in the two paired phases . when the amplitudes for the @xmath92 diagonals have the same sign , we refer to this as s - wave character ; on the other hand , when the signs are opposite , we call this d - wave . an extended discussion can be found in appendix  [ app : gt ] following eq . ( [ phipair2d ] ) .       and ( b ) the density - density structure factor @xmath71 for the 2-leg system at @xmath60 , @xmath93 , @xmath94 ; the system length is @xmath95 . the results are representative for the superfluid phase . ]     and ( b ) the density - density structure factor @xmath71 for the 2-leg system at @xmath60 , @xmath93 , @xmath94 ; the system length is @xmath95 . the results are representative for the superfluid phase . ]    . we show correlations for pairs on the diagonals @xmath96 and @xmath97 , but roughly similar values are obtained for all pair orientations . ] let us first discuss the phase at very small @xmath18 . the system is not frustrated in the @xmath98 limit and the ground - state wavefunction ( wf ) components in the boson occupation number basis should all be of the same sign ( positive ) to gain the lowest hopping energy . interestingly , the ed calculations for small @xmath99 and @xmath100 systems find that this feature of the wf remains robust up to a finite value of @xmath101 , close to the stability boundary of the sf phase . this is illustrated in the medium panel of fig . [ fig : ed_results ] where the average sign of coefficients in the ground - state wavefunction is plotted . the boson occupation number @xmath102 has a sharp peak at the wavevector @xmath103 . this is shown in fig . [ fig : sf_nq ] ( a ) for a system at @xmath60 , @xmath104 , @xmath94 , and length @xmath95 , calculated using the dmrg method . for this system size , the peak value at @xmath105 is about 100 times larger than the value at wavevector @xmath106 . this is consistent with the expected singular behavior @xmath107 and non - singular @xmath108 in the sf phase with 1d qlro ; the boson correlator decays as @xmath109 in real space and we estimate @xmath110 .    as shown in fig . [ fig : sf_nq ] ( b ) for the same system , the boson density - density structure factor @xmath111 has a @xmath112 dependence around @xmath113 reflecting the global boson number conservation , while @xmath114 is nonzero and has no visible singularities . the pairing correlations in real space show power - law behavior , with nearly equal magnitudes for all pair orientations . this is shown in fig . [ fig : sf_pair ] where we create a @xmath87 diagonal pair near origin and remove either a @xmath87 diagonal pair or a @xmath90 diagonal pair near @xmath115 . furthermore , the exponent @xmath116 of the power - law @xmath117 in the pairing correlator is found to be @xmath118 , which is about four times that for the single boson correlator , @xmath119 . such behavior of the pairing correlations is expected in the sf phase , which has qlro in the single boson correlations . these features characterize the sf phase . it also has a finite @xmath120 , where @xmath121 is the `` superfluid stiffness '' , which we can measure in our numerical calculations by imposing a small twist in the boundary conditions . results from ed of a small system with @xmath99 are given in the top panel of fig . [ fig : ed_results ] and related dmrg measurements for larger systems are in fig . [ fig : delta_e ] ; these are discussed in more detail in sec . [ subsubsec : rhos ] .     at filling @xmath60 and @xmath122 . shown are the superfluid stiffness ( top panel ) , the average sign of the ground - state wavefunction ( medium panel ) and the average double occupancy of rungs ( lower panel ) as a function of the ring exchange @xmath18 . ]     at filling @xmath60 and @xmath122 . shown are the superfluid stiffness ( top panel ) , the average sign of the ground - state wavefunction ( medium panel ) and the average double occupancy of rungs ( lower panel ) as a function of the ring exchange @xmath18 . ]     at filling @xmath60 and @xmath122 . shown are the superfluid stiffness ( top panel ) , the average sign of the ground - state wavefunction ( medium panel ) and the average double occupancy of rungs ( lower panel ) as a function of the ring exchange @xmath18 . ]       and ( b ) the density - density structure factor @xmath71 for the 2-leg system at @xmath60 , @xmath93 , @xmath123 ; the system length is @xmath95 . the results are representative for the dbl phase . ]     and ( b ) the density - density structure factor @xmath71 for the 2-leg system at @xmath60 , @xmath93 , @xmath123 ; the system length is @xmath95 . the results are representative for the dbl phase . ] as we further increase the strength of the ring - exchange @xmath18 , the system undergoes a ( first order ) phase transition where the stiffness @xmath121 drops steeply by a factor of two , see figs . [ fig : ed_results],[fig : delta_e ] and the subsection below . there is also an abrupt change in the character of the wavefunction sign which we can monitor in the ed analysis for small system sizes by simple counting and shown in the middle panel in fig . [ fig : ed_results ] . , there is also an abrupt change inside the dbl phase where the ground - state momentum jumps from @xmath124 to @xmath125 at around @xmath126 . we can understand these quantum numbers from specific orbital occupations in the dbl wavefunction approach , and the change of the total ground - state momentum arises from a change of orbital occupations in the @xmath8 bands . near @xmath127 vmc finds the optimal dbl state to be @xmath128_{\rm abc } \times    { \rm det}_2 [ n_2^{(0 ) } = 6 , n_2^{(\pi ) } = 2]_{\rm abc}$ ] . the label `` @xmath129 '' stands for antiperiodic boundary conditions when specifying the mean field orbitals ; the total bosonic wavefunction satisfies periodic boundary conditions and carries momentum @xmath130 . for the larger @xmath18 , vmc finds the optimal dbl state @xmath128 \times   { \rm det}_2 [ n_2^{(0 ) } = 5 , n_2^{(\pi ) } = 3]$ ] , where now the orbitals are specified with periodic boundary conditions ; this state carries momentum @xmath125 . ] this new phase has unusual phenomenology and can be identified as the dbl[2,1 ] phase : there is no boson peak at @xmath103 , but instead the boson occupation number @xmath102 shows singular points at nonzero wavevectors @xmath131 and @xmath132 , as if there is a `` condensate '' of bonding and antibonding boson modes at these @xmath133 . [ fig : dbl_nq ] presents characteristic plots of the @xmath102 and the density - density structure factor @xmath134 deep in this phase for filling @xmath60 taken at @xmath122 , @xmath123 , and a system length @xmath95 .     for a system with @xmath122 and filling @xmath60 . data for systems with length @xmath135 and @xmath136 are given . the momenta are determined from the peak values of @xmath102 such as shown in the top panel of fig . [ fig : dbl_nq ] . ]    the singular boson wavevectors @xmath137 and @xmath138 in general vary with parameters @xmath18 and @xmath50 . an example is shown in fig . [ fig : dbl_peak ] for the system at @xmath60 with @xmath93 and varying @xmath18 . just into the dbl phase near @xmath139 , we find a small @xmath137 and a large @xmath138 . in the wavefunction construction , this corresponds to rather different @xmath21 and @xmath22 , with significantly more @xmath7 fermions in the bonding than antibonding bands , see  fig . [ fig : bands ] . when we increase the ring coupling , @xmath137 decreases while @xmath138 increases , corresponding to moving the @xmath7 fermions from the bonding to antibonding levels and making the two bands more equally populated . interestingly , we find that the relation @xmath140 is satisfied , consistent with the prediction from gauge theory of a luttinger - type theorem in the correlated boson phase . the identification with the dbl is bolstered by a trial wavefunction study . the optimal such state for the system in fig . [ fig : dbl_nq ] is @xmath141_{\rm abc }   \right)^{0.75 } \\ & & \times \left ( { \rm det}_2 [ n_2^{(0 ) } = 32 , n_2^{(\pi ) } = 0]_{\rm abc } \right)^{0.1 } ~.\end{aligned}\ ] ] the length of the system is @xmath95 and the total boson number is @xmath142 . the @xmath7 determinant is constructed from @xmath143 bonding and @xmath144 antibonding orbitals , while the @xmath8 determinant is constructed from bonding orbitals only . to define the orbitals , we use antiperiodic boundary conditions ( `` abc '' ) along the @xmath10-direction for both @xmath7 and @xmath8 , which gives closed shells for the above occupation numbers , while the physical boson wavefunction respects periodic boundary conditions and carries zero total momentum . the determinant powers are understood as @xmath145 . the @xmath8 determinant can be written explicitly as @xmath146 \propto   \prod_{i < j } \sin\frac{\pi ( x_i - x_j)}{l_x } ~ , \label{detjw}\ ] ] and in particular prevents two bosons from sharing the same rung . in fact , the exact ground states in the ed / dmrg have very small double occupancy of rungs deep in this phase , as illustrated in the lower panel of fig .  [ fig : ed_results ] . the characteristic boson wavevectors in the above dbl state are @xmath147 and @xmath148 . these bosons are created by occupying @xmath7 and @xmath8 orbitals on the opposite sides , which is motivated by the amperean enhancement of such composites as described earlier . the overall match of the trial wavefunction results with the exact dmrg correlations is striking , reproducing the singular features as illustrated in fig . [ fig : dbl_nq ] .    turning to the density structure factor @xmath134 , there are many features but they are also weaker than in the @xmath102 . for @xmath149 , we expect a linear @xmath112-behavior near @xmath113 and signatures at @xmath150 , @xmath151 , which we see , and also at @xmath152 , which we do not see in either dmrg or vmc ( we can understand this in the wavefunction because the power for the @xmath153 is small and suppresses density correlations coming from this piece ) . at @xmath154 , we expect signatures at @xmath155 and at @xmath156 , both of which we see . as discussed in appendix  [ app : gt ] , it is a complicated matter to predict the exponents of the different singularities , and given our limited system size we do not attempt measuring these . it is remarkable that the dbl wavefunction can reproduce the singular features , but we do not know at this point whether it can reproduce all exponents in general .     calculated by dmrg and vmc . we show correlations for pairs on the diagonals @xmath96 and @xmath97 , while values for other pair orientations are much smaller . ]    finally , in fig . [ fig : dbl_pair ] we show the diagonal - diagonal pairing correlator . it is more short - ranged than in the sf phase as well as the two paired phases described below , but is characteristic of the dbl phase . first , @xmath157 has robust sign oscillations at wavevector @xmath158 . second , the @xmath87 and @xmath90 diagonal pairs have opposite signs ; the nontrivial sign structure in the wavefunction gives some `` d - wave''-correlated character to the boson liquid , which is observed in the two - boson correlator .       to the periodic boundary connection along the @xmath10-direction . the energy difference @xmath159 , multiplied by @xmath57 , ( top panel ) for a uniform system with @xmath122 for all bonds ; ( lower panel ) for a specially designed system with @xmath122 in one half and @xmath160 in the other half of the system . ]     to the periodic boundary connection along the @xmath10-direction . the energy difference @xmath159 , multiplied by @xmath57 , ( top panel ) for a uniform system with @xmath122 for all bonds ; ( lower panel ) for a specially designed system with @xmath122 in one half and @xmath160 in the other half of the system . ]    in this subsection we want to further highlight some characteristic difference between the sf and dbl phases . in the sf state , the wavefunction has qlro at @xmath20 . imposing a twist boundary phase @xmath161 in the boson hopping term in the hamiltonian ( [ hring ] ) , one can measure the superfluid stiffness as @xmath162 , which is found to scale as @xmath163 in the sf phase . this stiffness calculated from ed is shown for the small @xmath164 system in the top panel of fig . [ fig : ed_results ] . for larger system lengths , the stiffness can alternatively be estimated by dmrg calculations of the energy difference for a @xmath165 twist , i.e. , @xmath166 .    in the top panel of fig . [ fig : delta_e ] , we show @xmath167 as a function of the ring exchange @xmath168 for several system sizes @xmath169 . for isotropic hopping @xmath122 , we find that @xmath167 is essentially constant in the sf phase , but drops abruptly at the sf to dbl transition . in the dbl phase , @xmath167 shows rather irregular behavior including sign changes when varying either @xmath18 or @xmath57 . for a fixed system length @xmath57 , the more abrupt changes can be matched with the step - like changes of the singular momenta positions shown , e.g. , in fig . [ fig : dbl_peak ] for @xmath170 . for fixed ring exchange @xmath18 , the abrupt changes with varying @xmath57 are indicative of strong incommensurate correlations , where different @xmath57 offer varying degree of matching at the ( twisted ) periodic boundary : @xmath171 means that @xmath172 boundary conditions provide better matching than @xmath165 , while the situation is reversed if @xmath173 . a finite @xmath174 in the large system limit in the above setting implies that there is a propagating `` bonding '' mode ( i.e. , with transverse momentum @xmath175 ) . based on our picture of the dbl phase , there must also be a propagating `` antibonding '' mode with @xmath176 , which we would like to detect directly and contrast with the sf phase where there is no such mode .    to this end , we design a system composed of two halves connected at both ends , where the left half is a @xmath177-hamiltonian while the right half is a @xmath178-hamiltonian . first , we note that a model with @xmath160 can be related to a model with @xmath122 by a transformation @xmath179 . this ensures that the boson density is the same in the two parts . crucially , this transformation interchanges the notions of bonding and antibonding modes . if the designed system is in the regime of the regular sf phase , the gapless bonding bosons in the half with @xmath122 can not penetrate into the other half where the gapless mode is anti - bonding , and vice versa . since the bosons are not able to propagate around the loop , a @xmath165 twist should not produce significant change @xmath180 which will vanish exponentially with @xmath57 . on the other hand , in the dbl phase there are gapless modes of both bonding and antibonding type in either region , so that bosons can propagate around the loop , and the @xmath181-twist should change the energy by an amount proportional to @xmath163 . the  bonding / antibonding propagation stiffness " @xmath167 for the designed two - half system is shown in fig . [ fig : delta_e](b ) as a function of @xmath168 . indeed , we find a vanishingly small @xmath167 in the sf phase and a rather irregular but nevertheless finite @xmath167 in the dbl phase . in fig . [ fig : delta_e](b ) maintains a positive sign throughout the dbl phase , which can be understood roughly as follows : a boson propagates via a bonding mode ( with wavevector @xmath137 ) in one half and via an antibonding mode ( @xmath138 ) in the other half . the total phase accumulated by going around the system is @xmath182 for @xmath60 and the shown @xmath57 multiples of 6 . thus , the zero twist boundary condition provides the best ( crude ) matching for all @xmath18 . when we tried different @xmath57 ( not shown ) , the @xmath180 could be of either sign . ]       and ( b ) the density - density structure factor @xmath71 for a 2-leg system at @xmath60 , @xmath183 , @xmath184 ; the system length is @xmath95 . the results are representative for the s - wave paired phase . ]     and ( b ) the density - density structure factor @xmath71 for a 2-leg system at @xmath60 , @xmath183 , @xmath184 ; the system length is @xmath95 . the results are representative for the s - wave paired phase . ]    for our 2-leg ladder system one might expect to find the dbl[2,2 ] phase if the interchain coupling @xmath50 is much smaller than the coupling @xmath51 along the chains : when the boson hopping between the chains is reduced , one expects that the fermion hopping between the chains is also reduced for both species , which in turn could bring the antibonding band of the @xmath8 fermions to cross the chemical potential as shown in the bottom right panel of fig . [ fig : bands ] . dmrg studies of the ring model at filling @xmath60 and small @xmath185 indeed find a phase that is distinct from both the dbl[2,1 ] and sf as shown in the phase diagram fig . [ fig : phased1 ] around @xmath186 . however , our measurements indicate that this phase is not a dbl[2,2 ] , but a paired phase . as shown in fig . [ fig : sw_nq ] for a system at @xmath60 with @xmath183 , @xmath184 , and length @xmath95 , this phase is characterized by a quite flat region in the boson occupation number @xmath187 near @xmath113 without any singular points , which is different from both the sf and dbl phases . on the other hand , the density - density structure factor @xmath134 has a singular behavior at a wavevector @xmath188 . this singular behavior can be approximated by @xmath189 , with @xmath190 between @xmath191 providing the best fit to the data , which corresponds to power law envelope @xmath192 in real space . [ note that when @xmath193 the number of bosons in each chain is independently conserved and @xmath194 must vanish at @xmath113 , which explains the small value of @xmath195 extracted by dmrg when @xmath196 in fig . [ fig : sw_nq ] . ] the observed phase can be identified as an s - wave boson paired phase with qlro pairing correlation as shown in fig . [ fig : sw_pair ] . the pairing correlation @xmath157 in this phase is strongly enhanced at large distances @xmath115 when compared to the correlation in the dbl phase ; the values are roughly similar for pairs straddling the chains but smaller for pairs inside the chains . the decay of @xmath157 can be fitted by a power law behavior @xmath197 with @xmath118 . we have also checked the single boson correlation in the real space : it is much reduced compared to the nearby sf and dbl phases and decays very quickly  essentially exponentially ; in particular , it is much smaller than the boson - pair correlation beyond few lattice spacings . we then conclude that there is a gap to single - boson excitations , but not to pairs . this identification of the phase is also supported by a vmc study using paired boson wavefunctions described in appendix  [ subapp : psisw ] , which improve the energetics over that of the trial dbl states and reproduce qualitatively the dmrg correlations .    . we show diagonal pairs but roughly similar correlations are obtained also for pairs on the rungs , while correlations for pairs in the chains are significantly smaller . ] interestingly , the s - wave paired phase can be regarded as an instability of the dbl[2,2 ] phase . in appendix [ subapp : dbl22instab ] , we argue that without specially adjusted short - range interactions , the dbl[2,2 ] is unstable once the gauge fluctuations are included , and the most natural outcome is a paired phase with the same phenomenology as described above . in particular , the pairs are s - wave , carry zero momentum , and show qlro with power law @xmath198 , while we also predict power law charge correlations at wavevector @xmath188 with @xmath192 . in appendix [ subapp : dbl22instab ] we also predict that the two exponents are inverses of each other , @xmath199 . the dmrg estimates of the power laws differ somewhat from this , but they are hard to make accurate given the slow decay of the correlations and our limited system sizes .       and ( b ) the density - density structure factor @xmath71 for a 2-leg system at @xmath62 , @xmath196 , @xmath200 ; the system length is @xmath201 . the results are representative for the d - wave paired phase . ]     and ( b ) the density - density structure factor @xmath71 for a 2-leg system at @xmath62 , @xmath196 , @xmath200 ; the system length is @xmath201 . the results are representative for the d - wave paired phase . ]    consider now bosons at small filling @xmath62 , fig . [ fig : phased2 ] . as we further increase @xmath18 , the dbl phase in the intermediate range becomes unstable and a new phase emerges . as shown in fig .  [ fig : dw_nq ] for a system at filling @xmath62 with @xmath202 , @xmath200 , and length @xmath201 , this phase is characterized by @xmath203 , which is satisfied very accurately and is indicative of strongly suppressed boson correlator between the chains . we think that the peaks in @xmath4 are non - singular and represent short - range boson correlations , while @xmath35 decays exponentially on long distances : while we can not ascertain the exponential decay with our system sizes , the single boson correlator is already more than an order of magnitude smaller and decays much faster than the pair correlator described below . this phase is further identified as a d - wave boson paired phase with power - law pairing correlation as shown in fig . [ fig : dw_pair ] . the diagonal pairing correlations @xmath157 show opposite signs depending on the relative orientation of the two pairs : positive when both are @xmath87 diagonals and negative when one is @xmath87 and the other @xmath90 , which is consistent with the d - wave pairing picture . we find power law correlation @xmath204 with exponent @xmath205    the identification of this as d - wave paired phase is supported by the reasoning at low densities described in appendix  [ subapp : boundstate ] . thus , for a pair of bosons on an otherwise empty ladder , there is a transition at @xmath206 where it becomes favorable for bosons to form a molecule with internal d - wave character . if we now have a system of bosons at small density , it is natural that there will be a parameter range where the bosons are paired and the resulting molecules form a luttinger liquid .    to bolster this identification , we also studied such states variationally . as described in appendix  [ subapp : psisw ] , we construct a trial wavefunction as a product of a pfaffian and a determinant , since these are simple to work with in vmc . optimal trial state for the system in fig . [ fig : dw_nq ] is @xmath207 \right)^{0.3 } ~.\end{aligned}\ ] ] the system length is @xmath208 and the boson number is @xmath209 . the pfaffian is for an antisymmetric ( spinless fermion ) pair - function @xmath210 that connects only sites on different chains ; the optimal `` size of the pair '' is @xmath211 . this is a `` @xmath212 '' pairing and when composed with a `` @xmath213 '' character present in the determinant this gives a d - wave character to the boson pairs . strictly speaking , the above state is valid only in the @xmath193 limit , since it prevents fluctuation of the boson number in each chain . still , as we can see from fig . [ fig : dw_nq ] , the variational wavefunction reproduces the dmrg correlations surprisingly well in this phase . the boson momentum distribution function for the trial state is ( by construction ) non - singular , so that the close agreement with the dmrg found for @xmath208 strengthens the identification of this phase as a d - wave paired .    . ]    finally we mention the signatures in the boson density shown in the lower panel in fig . [ fig : dw_nq ] . @xmath111 shows the familiar @xmath112 at small momenta and also a singularity at @xmath158 . the latter gives a wavelength equal to the spacing between the pairs , which is expected in a luttinger liquid of such molecules . the density singularity is not too strong and is in line with the slow decay of the pair correlator . as we further increase the ring exchange towards @xmath214 , the uniform phase becomes unstable to phase separation . for such large @xmath18 , the boson density will separate into half - filled and empty regions , which is energetically favored by the ring exchange term . in the ed of the small @xmath99 system , phase separation is identified to happen at @xmath215 where all the lowest energy states from different momentum sectors become nearly degenerate . in the dmrg calculation , the phase separation is indicated by a nonuniform boson density in the obtained ground state , which breaks the translational symmetry due to small perturbations from the cut - off of the hilbert space in dmrg . the phase separation occurs around @xmath214 as observed from dmrg . finally , the phase separation can also be studied in vmc by measuring trial energies for different boson densities and performing maxwell construction ( this is described in more detail , e.g. , in ref . vii for the 2d ring model ) . the tendency towards phase separation for large @xmath18 is expected since the ring term provides effective attraction between bosons , see appendix  [ subapp : boundstate ] and ref . also , in the large @xmath18 limit , we can solve the ladder hamiltonian exactly and find phase separation  this is shown in appendix  [ subapp : konly ] . nevertheless , for intermediate @xmath18 the various phases discussed above are stabilized by the boson hopping . an obvious virtue of searching for the fingerprint of putative 2d gapless quantum phases in quasi-1d models is the relative numerical and analytic tractability of the latter . an important requirement of this approach is that the 2d quantum phase has singular correlations along surfaces in momentum space . in this case the number of 1d gapless modes present on an @xmath0-leg ladder will necessarily scale linearly with @xmath0 . while the particular fingerprint of the 2d quantum phase can already be present on the 2-leg ladder , as was the case for the dbl phase studied exhaustively in this paper , a compelling approach to 2d will certainly require analyzing models with larger @xmath0 . already when @xmath216 exact diagonalization studies become problematic since finding the ground states of boson or spin models with more than @xmath217 sites becomes virtually impossible . thus a @xmath218-leg ladder of length @xmath219 is already near the limit . dmrg is much more promising , and past studies on spin models reveal that convergence to the ground state is possible for lengths @xmath220 of order of a @xmath221 with @xmath222 , respectively . variational monte carlo for gutzwiller - type wavefunctions constructed as products of determinants or pfaffians as considered in this paper is possible for significantly larger @xmath0 , and equal time correlation functions can be computed with a fair degree of accuracy . but vmc is only as good as the variational wavefunction being employed , and the quality of the wavefunction is hard to assess in the absence of other more accurate methods . the ground - state energy of all physically reasonable hamiltonians is sensitive almost entirely to very short - range correlations , whereas for gapless quantum systems it is usually the longer range correlations which are necessary in order to distinguish between various competing phases . quantum states in 2d with gapless critical surfaces in momentum space are somewhat of an exception since the location of such surfaces is presumably due to short - range correlations in the signs and amplitudes of the wavefunction . but the precise singular structure on such surfaces reflects longer range correlations in the wavefunction . gauge theory@xcite seems to offer one of the few analytic approaches to access putative gapless spin - liquid type quantum phases in 2d that have no quasi - particle description . and for mean field states having a spinon fermi sea , a controlled inclusion of the gauge fluctuations is problematic , even in the simplest case of a compact @xmath223 gauge field . accounting for the instanton ( monopole ) events which reflect the compact nature of the gauge field is exceedingly difficult . even for a non - compact @xmath223 gauge field coupled to a fermi sea , while there does exist an rpa - type approach , @xcite it involves uncontrolled and worrisome approximations . but for the ladder systems , gapless fermions ( spinons ) coupled to a @xmath223 gauge field is eminently tractable , although this fact has been exploited rather infrequently.@xcite in 1d the gauge field has no transverse component and the bosonization method can be employed , leading to a gaussian effective field theory for the low energy excitations . this is another distinctive advantage of studying ladder systems relative to their 2d counterparts . one extension of this paper involves studying the boson ring hamiltonian on ladders with @xmath224 and larger . already for @xmath216 there are a number of new features that are not present on the 2-leg ladder . first , one can consider either periodic or anti - periodic boundary conditions for the bosons in the rung - direction , and it is quite plausible that the phase diagram and the presence / absence of the dbl phase will depend on this choice . second , while our search for a dbl phase with three gapless modes on the 2-leg ladder was unsuccessful due to instabilities of the dbl[2,2 ] variational state , it seems likely that for @xmath225 a dbl phase with more gapless modes than @xmath0 should be accessible . this is interesting since in such a case at least two of the 1d modes would have the same transverse momentum . for a general @xmath0-leg ladder one can construct variational states of the form dbl[@xmath226 , @xmath227 with integer @xmath228 satisfying @xmath229 ( with @xmath230 the number of partially filled bands for the @xmath231 partons , respectively ) . the number of 1d gapless modes in such a phase is @xmath232 , so that dbl[3,2 ] for @xmath216 , say , would have four 1d modes . if @xmath233 the @xmath8 fermion has only one band partially occupied , and assuming the band has @xmath23 , the corresponding variational wavefunction vanishes whenever two or more bosons occupy a rung . in this case of no - double rung occupancy , the determinant for the @xmath8 fermion can be viewed as a ( strictly ) 1d jordan - wigner string multiplying the @xmath7 determinant . whenever @xmath234 this is no longer the case , and the @xmath8 determinant would have a more subtle effect on the signs of the boson wavefunction . it would thus be desirable to access a dbl phase with @xmath234 . as detailed in appendix  [ subapp : dbl22instab ] , the gauge theory suggested an instability of the dbl[2,2 ] state for @xmath235 due to a rather special  nesting " condition present for @xmath236 . for @xmath237 this will not generally be the case so we expect that such dbl phases should be more easily accessible for @xmath238 . it would also be interesting to use the dmrg to access some dynamical information about the dbl , which should be possible at least for @xmath239 for @xmath240 the ladder model recovers the full symmetry of the 2d square lattice , and in future work it would be desirable to study how to approach this limit .      it would be most interesting to search for possible 2d spin liquid phases with singular surfaces in momentum space in models possessing su(2 ) spin symmetry . recently , several authors @xcite have suggested that the putative spin liquid observed in @xmath1-(et)@xmath2cu@xmath2(cn)@xmath3 is of such character , and have proposed that a mean field state with a fermi surface of spinons is an appropriate starting point . variational energetics on the corresponding gutzwiller projected spinon fermi sea have been performed on the triangular lattice heisenberg antiferromagnet with a four - site cyclic ring exchange term , a model argued to be relevant for a mott insulator such as @xmath1-(et)@xmath2cu@xmath2(cn)@xmath3 that has a comparatively small charge gap . it should be possible to study the heisenberg plus ring exchange spin hamiltonian on the triangular strip to search for signatures of the proposed parent 2d spin liquid . in the absence of ring exchange this model is equivalent to the 1d heisenberg model with both first and second neighbor interactions , @xmath241 , and has been explored using dmrg in earlier work.@xcite the phase diagram as a function of @xmath242 appears to have two phases : a bethe - chain phase stable for small @xmath242 with one gapless mode and a fully gapped dimerized state for intermediate and large @xmath242 ; the gap decreases exponentially for large @xmath242 where we have nearly decoupled legs ( @xmath243 ) of the triangular ladder ( with interchain @xmath244 ) .    with a 4-site ring exchange term ( coupling @xmath18 ) the triangular strip model is characterized by two dimensionless parameters , @xmath245 . are there any additional phases in this phase diagram ? the spin liquid phase on the triangular strip that descends from the ( putative ) 2d spinon fermi surface state should have 3 gapless modes ( 2 for spin , 2 for the two legs , minus one from the gauge constraint ) and corresponding incommensurate spin correlations . variational monte carlo together with a gauge theory analysis can provide a detailed characterization of this state , which could then be compared with dmrg , precisely as was done in this paper for the boson ring model . if present , one expects the spin liquid phase to appear for intermediate values of both @xmath242 and @xmath246 . in addition , one might study a half - filled hubbard model on the triangular strip , which should exhibit a metal insulator transition at intermediate coupling , @xmath247 . of particular interest is the nature of the quantum state just on the insulating side of the mott transition . if the mott transition is weakly first order , there will be substantial charge fluctuations and ring exchange interactions in this part of the phase diagram , and possibly a spin liquid phase , in addition to phases that have already been identified ( see refs .   and references therein ) . what are the prospects of using ladders to approach non - fermi liquid phases of 2d itinerant electrons ? one complication is as follows . imagine a weak coupling hubbard model on the square lattice ( at densities well away from special commensurate values ) which has a conventional 2d fermi liquid ground state . on an @xmath0-leg ladder the free fermi points present for @xmath248 would be converted into luttinger liquids . the jump discontinuity in the momentum distribution function at each fermi point would be lost , but singularities would remain , characterized by some luttinger liquid exponents . moreover , the location of the singularities would still satisfy the luttinger sum rule . now imagine some strong coupling 2d electron hamiltonian with a non - fermi liquid phase that has a residual luttinger fermi surface but with @xmath249 , analogous to a 1d luttinger liquid . the corresponding @xmath0-leg ladder descendant would be qualitatively indistinguishable from the phase of the weak coupling hamiltonian . conversely , the presence of a luttinger - satisfying fermi surface on an @xmath0-leg ladder would not enable one to distinguish between the two different 2d phases , one a fermi liquid and the other not .    on the other hand , the @xmath0-leg ladder descendant of a 2d non - fermi liquid ground state with momentum space singularities that violate the luttinger sum rule ( surfaces with the  wrong " volume , or perhaps even arcs ) would have qualitatively distinct signatures . several recent papers@xcite have proposed such 2d non - fermi liquid phases , and it would be extremely interesting to find evidence for their ladder descendants . the wavefunction for one of these phases was constructed@xcite by taking the product of a free fermion determinant and the dbl wavefunction in eq . ( [ psidbl ] ) . this phase , which inherits the d - wave sign structure from the bosonic dbl wavefunction , was called a d - wave metal and would exhibit distinctive signatures if present on a ladder system . a possible hamiltonian that might possess the d - wave metal can be expressed by adding an  itinerant - electron - ring " term to the usual square lattice @xmath66 model : @xmath250 with , @xmath251 where the sites @xmath252 are taken to run clockwise around an elementary square plaquettes and the summation runs over all packets . here we have defined an electron singlet creation operator : @xmath253 the ring term rotates a singlet on the @xmath254 diagonal into a singlet on the @xmath255 diagonal . it would be interesting to study this or other such strong coupling hamiltonians on the @xmath0-leg ladder .      to conclude , in this paper we have initiated a study of @xmath0-leg ladder systems that are descendants of candidate 2d quantum phases with low - energy excitations residing on singular surfaces . motivated by one such proposal for the so - called dbl phases of un - condensed itinerant bosons,@xcite we have studied the 2-leg model with frustrating ring exchanges using exact numerical approaches . we have indeed found the dbl[2,1 ] ladder version arising prominently in the phase diagram . we have also searched for the dbl[2,2 ] version but from the gauge theory analysis concluded that it is likely unstable due to special kinematics on the 2-leg ladder and instead gives rise to a boson - paired phase with extended s - wave character in the pairs . this paired phase is in fact realized in our model , while the dbl theory gives us tools to understand its properties . while our focus has been on the dbl ideas , we have explored the full phase diagram of the specific ring model in fair detail and found also other strong - coupling phases such as the above s - wave paired phase and the d - wave paired phase . the latter is not directly accessible from the dbl theory but is characteristic of the binding tendencies in our ring terms , which eventually cause the bosons to phase - separate for large ring exchanges . we hope to pursue similar ideas in the @xmath225 ladders and also in su(2 ) spin and itinerant electronic models with ring exchanges , which are particularly exciting . we would like to thank l. balents , t.  senthil , a.  vishwanath for useful discussions . this work was supported by doe grant de - fg02 - 06er46305 ( dns ) , the national science foundation through grants dmr-0605696 ( dns ) and dmr-0529399 ( mpaf ) , and the a.  p.  sloan foundation ( oim ) . some of our numerical simulations were based on the alps libraries @xcite . a faithful formulation of the physical system in the slave particle approach eq . ( [ d1d2 ] ) is a compact u(1 ) lattice gauge theory . we set this up on the 2-leg ladder as follows . denote the vector potential components on the links of the two chains as @xmath256 and @xmath257 , and on the rungs as @xmath258 . in the euclidean path integral , denote the temporal components associated with the sites on the two chains as @xmath259 and @xmath260 . the action for the gauge field contains @xmath261 , @xmath262 , and @xmath263 . for each @xmath264 , the last term gives a 1 + 1d gauge action , and if the cosine is interpreted as a villain cosine , we can replace it with @xmath265 and treat @xmath266 as non - compact with no change in the results . choosing the gauge @xmath267 , the combination @xmath268 is massive , while the combination @xmath269 is massless . having described the gauge sector , we now consider the fermions @xmath7 and @xmath8 with the mean field bands as in fig . [ fig : bands ] and take the continuum limit using the bonding @xmath270 and antibonding @xmath271 fields near the corresponding fermi points @xmath272 . these fields couple to @xmath273 as to the usual gauge field , while the massive combination @xmath268 can be integrated out . the continuum hamiltonian density is @xmath274 where @xmath275 and @xmath276 are the gauge charges of the @xmath7 and @xmath8 respectively . in the dbl[2,2 ] case , there are both bonding and anti - bonding fields present for each species @xmath7 and @xmath8 . the fermi momenta satisfy @xmath277 , where @xmath27 is the original boson density per site . in the dbl[2,1 ] case , the @xmath8 fermions have only bonding fields . the allowed interactions contain general density - density terms , @xmath278 ~,\end{aligned}\ ] ] where @xmath279 sum over all bands of all species , and also the following terms , @xmath280 \\ & + & f \left[\sum_p d_{1p}^{(0)\dagger } d_{2-p}^{(0)\dagger }                      d_{1-p}^{(\pi ) } d_{2p}^{(\pi ) } + \hc \right ] \\   & + & g \left[\sum_p d_{1p}^{(0)\dagger } d_{2-p}^{(\pi)\dagger }                      d_{1-p}^{(\pi ) } d_{2p}^{(0 ) } + \hc \right].\end{aligned}\ ] ] in the dbl[2,1 ] case , all terms that contain @xmath281 are absent . the above qed2-like theory can be analyzed perturbatively in the matter - gauge coupling , e.g. , in a systematic @xmath282 expansion , as was done in ref . for the spinon - gauge treatment of the 1d heisenberg spin chain . from such studies , we often use the following `` amperean '' rule of thumb : gauge interactions modify in a singular way processes that involve fermion fields with oppositely oriented group velocities . if the corresponding gauge currents are parallel ( antiparallel ) , such processes are enhanced ( suppressed ) , which originates from the attraction ( repulsion ) of currents in electromagnetism . as an example , @xmath283 has oppositely charged particles moving in opposite directions and producing parallel gauge currents , so this bilinear is expected to be enhanced compared to the mean field . we will see this explicitly in a solution of the 1d gauge theory by bosonization,@xcite which we pursue instead of the perturbative treatment . we also caution here and will see below that the different fermi velocities and the general allowed density - density interactions complicate the analysis significantly and can independently affect power laws , in addition to the above amperean rule for the gauge interaction effects .    to bosonize,@xcite we write @xmath284 . \label{bosonize}\ ] ] here @xmath285 and @xmath286 are the conjugate phase and phonon fields for each band , while @xmath287 are klein factors , which we take to be hermitian operators that commute with the bosonic fields and anticommute among themselves . the kinetic hamiltonian density becomes @xmath288 .\end{aligned}\ ] ] the density - density interactions eq . ( [ h4fermi ] ) lead to generic terms of the form @xmath289 and @xmath290 , which are strictly marginal and affect power laws . the other interactions eq .  ( [ h4fermii ] ) produce cosines and are written separately for the dbl[2,1 ] and dbl[2,2 ] below . we proceed in the euclidean path integral , choose the gauge @xmath291 and integrate out the @xmath292 . this renders the field @xmath293 massive and essentially pins it to zero . this is the suppression of the charge fluctuations in the gauge theory and realizes the microscopic constraints @xmath294 at long wavelengths . we then have two gapless modes left in the dbl[2,1 ] and three modes in the dbl[2,2 ] , if we can assume that the cosines corresponding to @xmath295 are irrelevant ( see below ) .    to characterize the phases of the original hard core bosons , we examine the single boson correlations , the boson density and current correlations , and the pair boson correlations . the microscopic boson field is written as @xmath296 which we then express in terms of the continuum bosonized fields ( see separate treatments of the dbl[2,1 ] and dbl[2,2 ] below ) . for oppositely moving @xmath297 and @xmath298 , the @xmath161 part has a non - zero projection onto @xmath299 , and gapping out the @xmath299 by the gauge fluctuations enhances the corresponding contribution , in agreement with the amperean rule applied to the oppositely charged @xmath7 and @xmath8 . the fermion densities are @xmath300 the two lines in the expansion separate @xmath149 and @xmath154 parts , and all contributions with the same non - zero momentum are grouped together ( assuming generic distinct bonding and antibonding bands ) . we also consider the particle current on the rungs ( the definition below uses the gauge choice @xmath267 made early on ) : @xmath301 \\ & & = i \sum_{p , p ' } \ ! \left (   d_{\alpha p}^{(0)\dagger } d_{\alpha p'}^{(\pi ) }        - d_{\alpha -p'}^{(\pi)\dagger } d_{\alpha -p}^{(0 ) } \right ) \ ! e^{-i(p k_{f\alpha}^{(0 ) } - p ' k_{f\alpha}^{(\pi ) } ) x } . \nonumber\end{aligned}\ ] ] the continuum expansion resembles that for the density at @xmath154 , but with opposite signs between @xmath302 and @xmath303 ; this will help to distinguish different phases in the analysis of instabilities , secs .  [ subapp : dbl21instab ] ,  [ subapp : dbl22instab ] . bosonized expressions for the contributing bilinears will be given separately in the dbl[2,1 ] and dbl[2,2 ] cases . here we only note that when the constituent particle and hole are on the opposite sides , we again have non - zero projection onto @xmath299 , so the corresponding contribution is enhanced by the gauge fluctuations , in agreement with the amperean rule .    finally , we consider the boson pair operator @xmath304 and expand in terms of the continuum fermion fields . assuming @xmath305 and @xmath306 are nearby , we characterize the pair by its center of mass coordinate @xmath307 along the chains ( this becomes the argument of the continuum fields ) and also by the internal pair structure . to this end , we collect different microscopic contributions that give rise to each term @xmath308 , where @xmath309 are combined band and fermi point indices . using short - hands @xmath310 $ ] , etc . , then @xmath304 contains @xmath311 \\ & & \times     \left [   d_{2c}(x_1 , y_1 ) d_{2d}(x_2 , y_2 )           - d_{2c}(x_2 , y_2 ) d_{2d}(x_1 , y_1)\right ] \\ & = & d_{1a}(x ) d_{1b}(x ) d_{2c}(x ) d_{2d}(x ) \\ & & \times e^{i(k_{1ax } + k_{1bx } + k_{2cx } + k_{2dx } ) x } \phi_{1a , 1b , 2c , 2d}(y_1 , y_2 , \xi ) ~,\end{aligned}\ ] ] where @xmath312 is the relative coordinate along the chains . the `` pair wavefunction '' is @xmath313 with @xmath314 and @xmath315 . in the above , we have chosen to characterize the pair by the momentum along the @xmath115-axis , but not by the momentum along the @xmath316-axis . on the 2-leg ladder , it is easier to visualize pairs by keeping track of both @xmath317 and @xmath318 . a few more words about this characterization on the 2-leg ladder . consider first a 2d setting , where we would write a contribution to a boson - pair @xmath319 by say @xmath320 where @xmath321 is some slowly varying operator of the center of mass coordinate @xmath322 , while @xmath323 is the internal pair function in the relative coordinate @xmath324 . specializing to a square lattice , we could then distinguish s - wave and d - wave ( more precisely , @xmath325 ) pairing by looking at the signs of @xmath326 . on the 2-leg ladder , however , separating out the transverse momentum @xmath5 as in 2d mixes things a little . specifically , a rotation of an @xmath327 pair into an @xmath328 pair can be also achieved by a `` translation '' @xmath329 . to avoid this , we instead define the pair - function by writing @xmath330 where @xmath331 and @xmath332 are the @xmath115-components of @xmath322 and @xmath324 respectively . we then characterize the pair wavefunction @xmath333 by the symmetry under interchanging the two legs of the ladder ( @xmath334 ) . when @xmath333 is even , we call it s - wave , since then the amplitudes for a pair sitting on a @xmath87 diagonal and a pair on a @xmath90 diagonal are the same . when @xmath333 is odd , we call it d - wave , since then the amplitudes for the two diagonals are opposite . the names `` s - wave '' and `` d - wave '' are used in anticipation that it is such local pictures of the pairs on the plaquettes that survive when we build towards 2d . returning to the dbl construction eq . ( [ phipair ] ) , we get s - wave or d - wave character depending whether @xmath335 is even or odd multiple of @xmath181 . we now specialize to the dbl[2,1 ] and dbl[2,2 ] cases in turn .      upon bosonization , we start with three modes @xmath336 , @xmath337 , and @xmath338 ( we drop the bonding band label for the @xmath8 fermions ) . to proceed formally , we change to new canonical variables via @xmath339    the microscopic boson field @xmath340 has the following contributions listed by their momentum ( cf . eq .  [ b ] ) : @xmath341 in each equation , the top and bottom signs are for @xmath342 and @xmath181 respectively . we omit the klein factors for simplicity since we will not need them here . also , we do not show combinations that can be obtained from the above by reversing all momenta . the microscopic boson density has contributions from both the @xmath7 and @xmath8 particles via eq . ( [ rhod ] ) , while the rung current has contributions only from the @xmath7 via eq . ( [ jrung ] ) ; we list all involved bilinears by their momentum : @xmath343 here we include the klein factors , but do not show combinations that can be obtained from the above by reversing all momenta .    when we include gauge fluctuations , @xmath299 becomes massive , locking together @xmath344 and @xmath345 . integrating out @xmath346 , we obtain a generic harmonic liquid theory in terms of two coupled modes @xmath347 and @xmath348 . of the cosine interactions eq . ( [ h4fermii ] ) , we have only @xmath349 assuming this is irrelevant , we have a stable phase with two gapless modes . the potential instability due to this term is considered in sec . [ subapp : dbl21instab ] below . we do not write explicitly the full dbl[2,1 ] theory in terms of the @xmath347 and @xmath348 since even without the additional interactions from eq . ( [ h4fermi ] ) we have a general coupled harmonic system . we can still make some observations about the scaling dimensions of the contributions to the boson and the boson density or current operators . in the initial fermionic mean field , the bands @xmath350 are decoupled and the scaling dimensions of all contributions to the operators @xmath351 and @xmath352 are equal to 1 . the gauge fluctuations effectively set @xmath353 everywhere . consider first the contributions to @xmath351 . for arbitrary couplings in the full harmonic theory , we can argue that the scaling dimensions of the operators @xmath354 are at least @xmath355 , i.e.  larger than the mean field value ( these have no explicit amperean enhancement ) . on the other hand , setting @xmath353 , the scaling dimensions of the operators @xmath356 can be lower than the mean field in accord with the amperean rule ; without further information , we can not say much  in general , these scaling dimensions can be as low as zero . consider now the contributions to the boson density @xmath352 . the zero momentum contribution has scaling dimension @xmath355 . we can argue for arbitrary couplings in the theory that the @xmath46 contribution has scaling dimension larger than the mean field value @xmath355 . on the other hand , the @xmath357 , @xmath45 , and @xmath358 contributions can have smaller scaling dimensions upon setting @xmath353 , in accord with the amperean enhancement rule . we do not have general bounds on the scaling dimensions in the coupled two - mode system , although clearly these are not all independent . furthermore , in all of the preceding discussion , even if the gauge interaction acts to enhance some correlation , the short - range interactions in 1d can act to suppress the correlation , so the ultimate fate is not clear .    as an illustration , we list all scaling dimensions in the case when the @xmath347 and @xmath348 modes are decoupled and characterized by the luttinger parameters @xmath359 and @xmath360 respectively [ our convention is that @xmath361 enters as @xmath362 in the action ] : @xmath363   & = & \frac{3}{8 } ( \frac{2}{3 g_\rho } + \frac{3 g_\rho}{2 } )     + \frac{1}{8 } ( \frac{2}{g_{- } } + \frac{g_{-}}{2 } ) \geq 1 ~ ; \\ \delta[b_{(k_{f1}^{(k_y ) } - k_{f2 } , k_y ) } ]   & = & \frac{1}{8 } ( \frac{2}{g_\rho } + \frac{g_\rho}{2 } )     + \frac{1}{8 } ( \frac{2}{g_{- } } + \frac{g_{-}}{2 } ) \geq 1/2 ~ ; \\ % \delta[\rho_{b , ( 2k_{f1}^{(k_y ) } , 0 ) } ]    & = & \frac{1}{4 } g_\rho + \frac{1}{4 } g_{- } ~ ; \\ \delta[\rho_{b , ( k_{f1}^{(0 ) } - k_{f1}^{(\pi ) } , \pi ) } ]    & = & \frac{1}{2 } ( \frac{2}{g_{- } } + \frac{g_{-}}{2 } ) \geq 1 ~ ; \\ \delta[\rho_{b , ( k_{f1}^{(0 ) } + k_{f1}^{(\pi ) } , \pi ) } ]    & = & \frac{1}{g_{- } } + \frac{1}{4 } g_\rho~ ; \\ \delta[\rho_{b , ( 2k_{f2 } , 0 ) } ] & = & g_\rho ~.\end{aligned}\ ] ] this example is loosely motivated by the observation in fig . [ fig : dbl_nq ] that the boson @xmath4 has roughly similar singularities at @xmath364 and @xmath365 . also , increasing the ring term @xmath18 and decreasing the interchain hopping @xmath50 drives the mean field bonding and antibonding bands to be more similar , suggesting such approximate decoupling of the two modes . continuing with this illustration , the boson correlators at @xmath366 are most enhanced when @xmath367 and @xmath368 giving @xmath369 = 1/2 $ ] ; with such couplings we also find mean - field - like @xmath370 = \delta[\rho_{b , ( k_{f1}^{(0 ) } - k_{f1}^{(\pi ) } , \pi ) } ]   = \delta[\rho_{b , ( k_{f1}^{(0 ) } + k_{f1}^{(\pi ) } , \pi ) } ] = 1 $ ] , while @xmath371 = 2 $ ] . these numbers are qualitatively similar to the singularities detectable by eye in the density structure factor in fig . [ fig : dbl_nq ] , where we see some signatures at @xmath372 , @xmath373 , and no visible signature at @xmath358 . a more careful look at the dmrg data shows that the density singularities at @xmath372 are slightly stronger than the mean field , while the ones at @xmath373 are weaker ; this would suggest @xmath360 somewhat smaller than @xmath374 , which would also be consistent with the assumed irrelevance of the cosine interaction eq . ( [ 4ferme1 ] ) . however the dmrg boson correlator singularity has estimated scaling dimension somewhat smaller than the maximum @xmath375 in this illustration , so one should not take the above too literally . nevertheless , this example gives some sense to the strengths of the observed singularities in fig . [ fig : dbl_nq ] .    finally , consider a boson pair operator ; we list contributing four - fermion combinations by their momentum along the chains : @xmath376 as explained earlier , we characterize each contribution by the internal pair wavefunction eq . ( [ phipair ] ) . again we have omitted the klein factors but can easily restore them when needed . for the zero momentum contribution we have @xmath377 this remains unchanged if we simultaneously change the @xmath316 coordinates of both bosons in the pair , e.g. , if we move a pair that lies entirely in chain 1 to chain 2 or if we turn a @xmath87 diagonal pair into a @xmath90 diagonal pair ; we call this s - wave .    for the contribution at momentum @xmath378 we have @xmath379 \nonumber \\ & & = 4 \sin(k_{f2 } \xi )   \biggl ( \delta_{y_1 = y_2 } ( -1)^{y_1 }           \sin\left[\frac{k_{f1}^{(0 ) } - k_{f1}^{(\pi)}}{2 } \xi \right ]   \nonumber \\ & & \quad\quad\quad\quad\quad         + i \delta_{y_1 \neq y_2 } ( -1)^{y_1 }           \cos\left[\frac{k_{f1}^{(0 ) } - k_{f1}^{(\pi)}}{2 } \xi \right ] \biggr ) ~. \nonumber\end{aligned}\ ] ] this changes sign if we simultaneously change both @xmath316 coordinates in the pair . when @xmath21 and @xmath22 approach each other , the pairs straddling the two chains have larger amplitude . since the amplitude changes sign when we turn a @xmath87 diagonal pair into a @xmath90 diagonal pair , the internal structure of the pair has some d - wave character . the contribution at @xmath380 can be similarly characterized , but is omitted here .    in the free fermion mean field all contributions have scaling dimension 2 . beyond the mean field , the d - wave contribution at @xmath381 is potentially enhanced by setting @xmath353 ; still , its scaling dimension is at least @xmath355 , while we can not say much about the scaling dimension of the s - wave contribution at zero momentum . in the above illustration with decoupled @xmath382 and @xmath383 modes , the d - wave pair has scaling dimension @xmath384 if we set @xmath385 , while the s - wave pair has scaling dimension is @xmath386 if we also set @xmath387 . the latter would be larger if we take smaller @xmath360 , in agreement with the observed dominance of the d - wave - like pair correlations in the dmrg fig . [ fig : dbl_pair ] oscillating at wavevector @xmath158 . let us ask what phase is obtained starting from the dbl[2,1 ] if the interaction eq . ( [ 4ferme1 ] ) is relevant [ the amperean rule does not point one way or another , but additional interactions from eq . ( [ h4fermi ] ) can make this happen ] . then @xmath388 is pinned , with two distinct cases depending on the sign of @xmath389 considered below , while @xmath390 fluctuates wildly . all contributions to @xmath351 contain @xmath390 in the exponent , so the single boson correlation function decays exponentially  there is a charge gap to single boson excitations . the @xmath347 mode still remains gapless , and we have a single - mode harmonic fluid described by the parameter @xmath359 . the boson density shows power law correlations at zero momentum with scaling dimension 1 . the remaining particle - hole bilinears with power law correlations are at wavevector @xmath391 ( scaling dimension @xmath392 ) and at @xmath393 ( scaling dimension @xmath359 ) . the @xmath20 and @xmath394 parts do not distinguish qualitatively between the two possibilities @xmath395 and @xmath396 , while the @xmath397 part does distinguish . explicitly , the @xmath7 particle density is @xmath398 where we have omitted the zero momentum component and parts with exponentially vanishing correlations . on the other hand , the @xmath7 particle rung current is @xmath399    consider now two cases :    : it follows that @xmath400 . in this case @xmath401 shows power law correlations at @xmath45 , while @xmath402 correlations are absent . : it follows that @xmath403 . in this case @xmath401 correlations at the above wavevector are absent , while @xmath402 shows power law at wavevector @xmath381 along the chains . let us finally consider the pair - boson correlations , which are also power law . putting together all contributions to @xmath404 , we have @xmath405 \\ + \sum_p \eta_1^{(0 ) } \eta_1^{(\pi ) }   e^{i ( 2 \phi_\rho + p [ \theta_\rho + ( k_{f1}^{(0 ) } + k_{f1}^{(\pi ) } ) x ] ) }   \phi_{1p(0 ) , 1p(\pi ) , 2r , 2l } ~,\end{aligned}\ ] ] where the appropriate @xmath406 are given in eqs . ( [ phiq0 ] ) and ( [ phiq2pirho ] ) . the pairs at zero momentum have dominant correlations with scaling dimension @xmath407 ; these have s - wave character as far as rotating diagonal bonds is concerned , but the amplitude details depend on the sign of @xmath389 . the pairs at @xmath381 have subdominant correlations with scaling dimension @xmath408 ; these have d - wave character for rotating the diagonal bonds and the details do not depend much on the sign of @xmath389 . summarizing , the resulting phase has a gap to single boson excitations , but has power law correlations for particle - hole and particle - particle composites . the dominant particle - hole composite has scaling dimension @xmath392 and either contributes to the density at wavevector @xmath188 if @xmath395 or to the rung current at @xmath158 if @xmath396 . the dominant particle - particle composite has scaling dimension @xmath407 and represents s - wave pairing at zero momentum . loosely speaking , we can describe this phase by saying that the bosons form s - wave - like pairs and these molecules in turn form a luttinger liquid . the current and/or density fluctuations occur at a wavelength equal to the mean inter - pair spacing along the ladder , as expected for a 1d luttinger liquid ( of pairs ) . the form of the particle - hole fluctuations across the rungs presumably reflects the internal pair structure in the two cases .      in the dbl[2,2 ] case , we start with four bands , @xmath409 , and have four modes upon bosonization eq . ( [ bosonize ] ) . it is convenient to introduce the following canonical variables : @xmath410 with the same transformation for the @xmath161 variables . the cosine interactions eq .  ( [ h4fermii ] ) are @xmath411 where @xmath412 .    proceeding with the analysis as in the dbl[2,1 ] case , upon integrating out the gauge field , @xmath299 becomes massive . if by tweaking the strictly marginal interactions @xmath413 we could render the @xmath295 terms irrelevant , we would end up with a phase with three gapless modes . we call this possible phase dbl[2,2 ] , and it can be analyzed similarly to dbl[2,1 ] . for example , the mean field boson correlation in this phase would read @xmath414 \\        - \cos[(k_{f1}^{(k_{y1 } ) } + k_{f2}^{(k_{y2 } ) } ) x ] \bigr ) ~.\end{aligned}\ ] ] the first term in the brackets comes from @xmath415 and is expected to be enhanced by amperean attraction ( in the bosonization , this term contains @xmath299 in the exponent and is potentially enhanced upon setting @xmath416 ) . on the other hand , the second term in the brackets has no amperean enhancement . as emphasized earlier , besides the gauge fluctuations crudely captured by the amperean rule , other interactions can also change the scaling dimensions . interestingly , and this has important consequences explored below , @xmath417 and @xmath418 carry the same momentum @xmath419 , and also @xmath420 and @xmath421 carry the same momentum @xmath422 . in fact , josephson - like coupling between the first two `` bosonic modes '' is allowed and gives rise to the @xmath423 interaction in eq . ( [ h4fermii ] ) , while josephson coupling between the last two modes gives rise to the @xmath424 interaction . if the gauge physics is the dominant determining factor in the dbl[2,2 ] phase , we expect the state to be strongly unstable to these interactions , and we pursue this scenario below . on the other hand , if the dbl[2,2 ] were stabilized against the @xmath423 and @xmath424 by some density - density interactions dominating over the gauge physics , which is possible and interesting , we would have little intuition from the gauge theory perspective , and we do not pursue this possibility further . if the gauge physics dominates , it appears very likely that the @xmath423 and @xmath424 terms are relevant , as can be seen from their bosonized expressions upon setting @xmath353 ( in accord with the amperean rule ) . let us explore the resulting phase when @xmath423 and @xmath424 flow to large values and pin the fields @xmath425 and @xmath388 . four different possibilities depending on the signs of @xmath423 and @xmath424 are discussed below . in such a phase , the conjugate fields @xmath426 and @xmath390 fluctuate strongly and any operator containing these in the exponent will have only short - range correlations . thus , we conclude that the boson correlations decay exponentially , i.e. , there is a gap to single boson excitations . however , there is still one gapless mode @xmath347 , which we characterize with the luttinger parameter @xmath359 . consider fermion bilinears composed of a particle and a hole of the same species that enter density and current , eqs . ( [ rhod],[jrung ] ) . other than the @xmath427 densities @xmath428 , the following bilinears survive when @xmath429 get pinned : @xmath430 the particle density eq . ( [ rhod ] ) is , upon dropping the zero momentum part and setting @xmath416 , @xmath431 the particle current on the rungs eq . ( [ jrung ] ) is @xmath432    we now discuss different phases that can arise depending on the signs of the couplings @xmath423 and @xmath424 . we set @xmath433 , from which it follows @xmath434 ( the physical results are independent of the choice of @xmath435 ) . there are four cases :    : it follows that @xmath436 , @xmath437 . in this case @xmath438 show power law correlations at @xmath188 with the scaling dimension @xmath392 , while @xmath439 correlations are absent . this is a natural phase coming from the microscopic gauge theory since the fluctuations of the @xmath7 and @xmath8 densities are in sync ; also , when we integrate out the massive gauge field @xmath268 early in the derivation of the continuum theory , we in fact generate such @xmath440 couplings . foretelling the analysis of the pairing correlations below , we propose that this is the s - wave - paired phase observed in the dmrg , sec . [ subsec : s - paired ] . : it follows that @xmath441 , @xmath442 . in this case @xmath443 show power law correlations at @xmath188 , while @xmath439 correlations are absent . the fluctuations of the @xmath7 and @xmath8 densities are out of sync . this phase with spatial modulation of the gauge charge is less natural coming from the lattice theory . it is in principle allowed if the bare couplings are finite , since then such modulation costs only finite energy density that can be offset by some other short - range interactions ; but the cost is large if the microscopic theory is at strong coupling as is usually the case in the slave particle treatments . for example , only the phase with @xmath438 can be realized by our wavefunctions . : it follows that @xmath444 . in this case @xmath445 correlations are absent , while @xmath446 show power low correlations at wavevector @xmath158 along the chains with scaling dimension @xmath392 . the @xmath7 and @xmath8 particle currents are in sync , so this phase is natural in the gauge theory ; for the original bosons , it would have enhanced current - current correlations but not density correlations . however , we have not observed such a phase in the ring model . : it follows that @xmath447 . in this case @xmath445 correlations are absent , while @xmath448 show power low correlations at wavevector @xmath158 along the chains . since the currents are out of sync , this phase is less natural coming from the microscopic gauge theory . let us consider boson pairs . first , for each @xmath449 there are four @xmath450 continuum field combinations that survive when the @xmath429 are pinned : @xmath451 } ~,\end{aligned}\ ] ] where in the last line the top / bottom sign corresponds to @xmath452 . constructing boson pairs via @xmath453 , the field @xmath454 disappears , and we also set @xmath353 . of the sixteen combinations , six carry zero momentum and all have the largest scaling dimension @xmath407 : @xmath455 the second and third lines will likely have smaller amplitudes , since @xmath425 and @xmath388 still fluctuate a little about the pinned values . we now examine the internal pairing structures . corresponding to the first line above , we have : @xmath456 the @xmath457 piece describes a pair formed by bosons in the same chain , while @xmath458 is a pair straddling the two chains . when say @xmath459 , which becomes more accurate as we decrease @xmath50 , the amplitudes for the latter are significantly stronger , i.e. , pairs are predominantly straddling the two chains .    corresponding to the other zero - momentum combinations , we have @xmath460 if @xmath461 , this is independent of whether the two sites are on the same or different chains . on the other hand , if @xmath462 , this has opposite signs if the two sites are on the same versus different chains . all of the above pair - functions are even if both particles are moved perpendicular to the chains . in particular , @xmath87 and @xmath90 diagonal pairs have the same amplitudes . this is what we call s - wave pairing . putting together all zero momentum contributions to @xmath304 , we have @xmath463 ~,\end{aligned}\ ] ] where @xmath464 and eventually drops out since both @xmath465 and @xmath466 are determined uniquely by the signs of the @xmath423 and @xmath424 couplings : @xmath467 , @xmath468 . note also that the contributions in the last four lines will have smaller amplitudes since @xmath469 and @xmath383 fluctuate a little around the pinned values . so the numerically largest contribution is given by the first line , while the smaller contributions from the last four lines will add to it with relative signs that depend on the signs of the @xmath423 and @xmath424 . summarizing , the dbl[2,2 ] is unstable towards a boson - paired phase with s - wave pairs carrying zero momentum . the resulting phase is roughly similar to the s - wave paired phase discussed as a possible instability of the dbl[2,1 ] in the preceding section , but some details are different . for example , coming out of the dbl[2,1 ] the amplitude for rung pairs vanishes , while out of the dbl[2,2 ] the rung pairs have large amplitudes comparable to those of the diagonal pairs . the latter is more similar to what the dmrg finds in the s - wave paired phase in the ring model , sec . [ subsec : s - paired ] . the prediction of dominant density correlations at wavevector @xmath188 also agrees with the dmrg . one additional prediction from the theory is that the pairing and density power laws have exponents that are inverse of each other . while the dmrg estimates for the present system sizes do not satisfy this exactly , there are likely strong finite size effects , and we would like to revisit this with larger systems . in the dbl construction eq . ( [ psidbl ] ) , we can view one determinant as affecting a generalized flux attachment or jordan - wigner ( jw ) transformation and view the other determinant as a ground state of the new jw fermions . in the dbl[2,1 ] case , the @xmath8 determinant is composed entirely from the bonding orbitals and becomes eq . ( [ detjw ] ) ; in particular , it prevents two bosons from being on the same rung . if we consider a hard - core boson model that prohibits double occupancy of rungs , then @xmath470 is the conventional 1d chain jw transformation on this ladder . as can be seen from the bottom panel in fig . [ fig : ed_results ] , our ring model in the dbl[2,1 ] phase appears , by its own dynamics , to strongly suppress double rung occupancy . the optimized power of @xmath153 in the trial wavefunction is also relatively small , making the @xmath471 factor look more like the jordan - wigner . finally , the @xmath7 determinant is composed of the bonding and antibonding orbitals , and its optimized power is relatively large though smaller than @xmath355 , which suggests that the jw fermions are not far from being free despite the large ring exchanges and the restricted rung occupancy . on the other hand , the dbl[2,2 ] would have some double occupancy of rungs and would be an example of a non - trivial jw , but unfortunately this phase appears to be unstable as described in sec . [ subapp : dbl22instab ] . the instability happens because of special kinematic conditions satisfied here , which allow direct josephson coupling and locking between enhanced boson modes . however , we expect that on @xmath0-leg ladders with @xmath225 , dbl phases will exist which can not be described by a conventional jordan - wigner approach . the pure ring model on the 2-leg ladder can be solved exactly . in the absence of boson hopping terms ( @xmath472 ) , the number of bosons on each rung , @xmath473 , is separately conserved . the number of bosons in each chain , @xmath474 or @xmath475 , is also conserved . any rung with @xmath476 or @xmath374 effectively breaks the system into decoupled pieces . consider an isolated segment of @xmath220 rungs with each @xmath477 . the ring hamiltonian is mapped to an xy spin chain by identifying boson configuration @xmath478 with spin up and @xmath479 with spin down : @xmath480 = k \sum_{x=1}^{l-1 } [ \tau^+(x ) \tau^{-}(x+1 ) + \hc ] ~,\ ] ] where @xmath481 are the usual spin-1/2 operators . this is readily solved by free fermions , and the ground - state energy is @xmath482 & = & -2 k \sum_{n=1}^{n_{\rm max } } \cos\frac{\pi n}{l+1 } \\ & = & -k   \left ( \frac { \sin\frac{\pi ( 2 n_{\rm max } + 1)}{2(l+1 ) } }              { \sin\frac{\pi}{2(l+1 ) } } - 1 \right ) ~,\end{aligned}\ ] ] where @xmath483 if @xmath220 is even ( in this case @xmath484 ) , while @xmath485 if @xmath220 is odd ( @xmath486 ) . in either case , the ground - state energy per boson is minimized by making @xmath220 large , with the asymptotic behavior @xmath487}{l } = -\frac{2 k}{\pi } + \frac{k(1 - 2/\pi)}{l } + o(\frac{1}{l^2 } ) ~.\end{aligned}\ ] ]    going back to the @xmath18-only model on an infinitely long ladder , we conclude that for arbitrary boson density @xmath488 , it is advantageous to phase - separate into an empty region and a half - filled region , since this minimizes the energy per particle ( @xmath489 can be treated by particle - hole transformation ) . the phase separation arises because the ring exchanges provide an effective attraction between particles  see also the next section . note that the half - filled region itself is a highly correlated state of bosons . consider two bosons on an otherwise empty ladder . in the absence of the hoppings , @xmath490 , the states @xmath491 are degenerate ground states with energy @xmath492 . each can be viewed as a d - wave pair sitting on a plaquette . small boson hopping can be included perturbatively , and the pair starts to hop along the ladder with amplitude @xmath493 . in fact , inspired by similar considerations in the 1d @xmath66 model in ref .  , the zero momentum bound state can be written exactly for the general @xmath494 model , and is created by @xmath495 ~.\end{aligned}\ ] ] this is defined for @xmath496 and has energy @xmath497 the pair `` size '' is @xmath498 $ ] . note that the bound state does not utilize any rung hopping . also , it competes with a free state of the two bosons where each moves independently and where the total energy is @xmath499 assuming @xmath500 throughout . as an example , for @xmath501 the bound state has lower energy than the state with two free bosons when @xmath502 = 2j$ ] , while @xmath503 = 2.74 j$ ] and @xmath504 increases with increasing @xmath50 .    for a very small density of bosons , it is then natural to propose that for @xmath505 the bosons will form such tightly bound pairs which in turn form a luttinger liquid . this is the picture of the d - wave paired phase . as discussed in the preceding section , in the absence of the hopping ( or for very large @xmath18 ) , the system at any density phase - separates into empty and half - filled regions . this maximizes the ring energy per particle which attains @xmath506 and , in particular , wins over the asymptotic @xmath507 energy per particle in the paired state . bosons in the phase - separated half - filled region form a non - trivial quantum state and can hardly gain any kinetic energy once the hopping is included . in the boson - paired liquid , on the other hand , the pairs can move and gain some kinetic energy , which can then stabilize the system against phase separation . as a crude estimate , we ignore the kinetic energy in the phase - separated case and compare @xmath506 with eq . ( [ ebound ] ) divided per boson and find that the boson - paired liquid is more stable for @xmath508 . thus , there is a sizable window over which the boson - paired liquid is energetically preferable over the free boson liquid and is stable against phase separation .    at very low density we also need to consider more - boson bound states that successively improve their ring energy but lose kinetic energy , and the corresponding liquids could intervene between the boson - paired one and the phase separation . however , it seems likely that a sizable region of the boson - paired phase will remain ; of course , it also competes with other phases such as dbl , as happens in fig . [ fig : phased2 ] for the filling @xmath62 studied by dmrg . let us sketch how one can construct trial wavefunctions for boson - paired liquids such as the d - wave paired phase at low densities described in the last section and in sec . [ subsec : d - paired ] and the s - wave paired phase described in sec .  [ subsec : s - paired ] . let us denote the boson pair wavefunction as @xmath509 . for a given pairing - up of particles we could write @xmath510 times some jastrow - type factor for the center of mass coordinates @xmath511 , and then symmetrize this to obtain a bosonic wavefunction : @xmath512.\end{aligned}\ ] ] even if we take a simpler jastrow which is already a symmetric function of @xmath513 , the symmetrization of the pairing part gives a permanent which is prohibitive in numerical calculations .    an alternative way to construct a symmetric bosonic wavefunction is @xmath514 \\ & \times & { \rm det}[r_1 , r_2 , r_3 , \dots ]   \times { \rm jastrow}[r_1 , r_2 , r_3 , \dots ] ~.\end{aligned}\ ] ] here pairing is realized with the help of a pfaffian , which is a bcs wavefunction for spinless fermions and is specified by some pair function @xmath515 , while the total bosonic symmetry is recovered by multiplying by a second fermionic wavefunction taken for simplicity to be a slater determinant of some orbitals . note that for tightly bound pairs when @xmath516 and @xmath515 are very short ranged the two constructions are essentially similar for statistically significant configurations . the pfaffian and determinant are simple to work with in vmc , which allows one to construct trial wavefunctions also beyond the tightly bound limit . of course , both the spinless fermion pair function @xmath515 in the pfaffian and the orbitals in the determinant need to be judiciously chosen to realize the particular boson pair function @xmath516 . we should also note that such @xmath517 wavefunction does not necessarily represent a paired state of bosons , and one needs to measure the long - distance properties first . for example , we found that for a 1d chain a product of a slater determinant and a @xmath213-wave bcs wavefunction has qlro in the boson correlator when the @xmath213-pairing problem is at weak coupling , while the boson correlator becomes short - range  as is needed for the boson - paired state  only when the @xmath213-pairing problem is at strong coupling . consider for illustration low - density bosons in the @xmath501 limit . we can construct a d - wave paired state by taking @xmath518 in the pfaffian part , while combining the determinant and jastrow into one factor @xmath519 the power of the determinant is understood as @xmath520 . the fermion pair function is @xmath212-type and here requires that the two particles are on different chains ( which is why it is strictly appropriate only in the @xmath501 limit ) . for an isolated pair of bosons , this indeed gives the exact d - wave pair of sec . [ subapp : boundstate ] when @xmath521 . we use this wavefunction with two variational parameters @xmath332 and @xmath68 to connect with the dmrg results at @xmath62 in sec . [ subsec : d - paired ] . one could construct a similar boson - paired state by taking the pfaffian factor to be a solution of a general spinless bcs problem with @xmath212 pairing on the ladder ; this would allow more freedom ( e.g. , if we want to study @xmath522 ) , but has not been explored . we have also tried the @xmath517 construction for the s - wave paired phase of sec . [ subsec : s - paired ] . here we take the determinant to contain both bonding and antibonding orbitals , which resembles the @xmath523 in the dbl[2,1 ] or dbl[2,2 ] containing @xmath7 fermions with dominant hopping along the chains ( in particular , the @xmath7 bonding and antibonding bands become equally populated in the @xmath501 limit ) . we take the pfaffian to be that for a @xmath213-wave spinless bcs problem on the ladder , which resembles @xmath8 fermions trying to enter the dbl[2,2 ] regime but becoming paired . there are many parameters here such as the @xmath8 intra- and inter - chain hoppings and @xmath213 pairings . we have found that the boson s - paired phase from fig . [ fig : phased1 ] and sec .  [ subsec : s - paired ] is reproduced roughly when the @xmath8 bcs problem is in the weak - pairing regime for both the bonding and antibonding bands and has predominant inter - chain pairing . note that the above @xmath524 motivation of the @xmath517 content is very loose and is hardly inspired by any gauge theory thinking but rather found by trial and error . we also mention that there are likely other ways to construct such paired boson phases . for example , coming from 2d where paired bosons can form a true condensate , one natural construction to explore could be @xmath525 . nevertheless , the @xmath517 that we studied gives us complementary confidence in the dmrg identification of the s - wave paired phase .
Abstract:

Answer: developing a theoretical framework to access quantum phases of itinerant bosons or fermions in two dimensions ( 2d ) that exhibit singular structure along surfaces in momentum space but have no quasi - particle description remains as a central challenge in the field of strongly correlated physics . in this paper we propose that distinctive signatures of such 2d strongly correlated phases will be manifest in quasi - one - dimensional @xmath0-leg ladder " systems . characteristic of each parent 2d quantum liquid would be a precise pattern of 1d gapless modes on the @xmath0-leg ladder . these signatures could be potentially exploited to approach the 2d phases from controlled numerical and analytical studies in quasi-1d . as a first step we explore itinerant boson models with a frustrating ring exchange interaction on the 2-leg ladder , searching for signatures of the recently proposed two - dimensional d - wave correlated bose liquid ( dbl ) phase . a combination of exact diagonalization , density matrix renormalization group , variational monte carlo , and bosonization analysis of a quasi-1d gauge theory , all provide compelling evidence for the existence of a new strong - coupling phase of bosons on the 2-leg ladder which can be understood as a descendant of the two - dimensional dbl . we suggest several generalizations to quantum spin and electron hamiltonians on ladders which could likewise reveal fingerprints of such 2d non - fermi liquid phases .