over the past decades scanning tunneling microscopy ( stm ) and spectroscopy techniques have been established as an important experimental tool to study strongly correlated electron systems with a spatial resolution down to the atomic scale . in these experiments the tunneling current between the sample and the stm tip is measured as a function of its position and the applied voltage .
the tunneling current is directly related  @xcite to the local density of states ( ldos ) in the sample ; thus stm experiments provide a tool to investigate the spatial dependence of the ldos , which may arise for example due to impurities .
these spatial modulations can be analyzed in terms of the fourier transform of the tunneling conductance which is directly proportional to the fourier transform of the ldos . using this method of analyzing stm data
one can infer informations about the bulk state of matter like the properties of its quasiparticle excitations .
for example , this has been used to study high - temperature superconductors  @xcite as well as carbon nanotubes  @xcite .
theoretical studies of the ldos and stm in the presence of boundaries have been carried out in particular for one - dimensional systems including luttinger liquids  @xcite , open hubbard chains  @xcite , and charge - density wave ( cdw ) states  @xcite . in these systems
the coupling to impurities is relevant  @xcite , which leads , at sufficiently low energies , to an effective cutting of the system into disconnected pieces . in this sense
the effect of an impurity can be modeled by a boundary condition .
previous studies of the fourier transform of the ldos in both gapless  @xcite and gapped  @xcite systems revealed a pinning of the cdw at the impurity as well as individually propagating spin and charge degrees of freedom .    here
we study a quarter - filled one - dimensional mott insulator in the presence of a strong impurity potential .
it has been established  @xcite that for sufficiently strong interactions double umklapp processes are relevant and generate a gap in the charge sector . as a microscopic realization one may think of a quarter - filled hubbard model extended by including nearest - neighbor interactions , which possesses  @xcite an insulating phase for sufficiently strong repulsions .
experimentally these systems are relevant for the description of organic quasi - one - dimensional conductors , i.e. the bechgaard and fabre salts  @xcite .    in this article
we study the low - energy behavior of the ldos by employing the field theoretical description of the mott insulator .
the spectrum consists of massless spin excitations ( spinons ) carrying spin @xmath0 but no charge and massive charge excitations ( solitons and antisolitons ) carrying charge @xmath1 but no spin . in this framework
the impurity is modeled by a boundary for the collective excitations .
in particular , we will focus on signatures of dispersing quasiparticles in the spatial fourier transform of the ldos , as they have been clearly observed in similar studies of luttinger liquids  @xcite and half - filled mott insulators as well as cdw states  @xcite .
the spectral function of the translationally invariant , quarter - filled mott insulator has been previously derived  @xcite by essler and tsvelik . from the quantum numbers stated above
it is clear that an electron will fractionalize into at least one spinon and two antisolitons .
thus the spectral function exhibits a featureless scattering continuum of at least three excitations , which is also consistent with experimental results  @xcite from angle - resolved photoelectron spectroscopy on bechgaard salts .
this article is organized as follows . in sec .
[ sec : model ] we will briefly discuss the field theoretical description of the quarter - filled mott insulator . in sec .  [
sec : gf ] we present our results for the single - particle green function .
this will then be used to derive an analytic expression for the ldos , which we will discuss in detail in sec .
[ sec : ldos ] .
our main result is presented in eq .   as well as figs .
[ fig : plot1 ] and  [ fig : plot2 ] , which show the fourier transform of the ldos in the presence of an impurity .
surprisingly , we find two dispersing features which follow the dispersion relations and respectively .
these features can be interpreted as arising from an antisoliton or a spinon - antisoliton pair propagating between the position of the stm tip and the impurity .
in addition , we observe a non - dispersing singularity at momentum @xmath2 which is indicative of a pinning of a spin - density wave ( sdw ) at the boundary .
a detailed analysis further reveals dispersing features originating in propagating spinons and antisoliton pairs .
the suppression of these features is attributed to the `` quasi - fermionic '' behavior of antisolitons with equal momenta .
finally we explain the different findings in the spectral function and the ldos . in sec .
[ sec : bbs ] we discuss the effect of possible boundary bound states on the ldos and in sec .
[ sec : hfmi ] we compare our results to the half - filled mott insulator studied  @xcite previously . finally , in sec .
[ sec : discussion ] we set our results in the context of other one - dimensional systems and discuss implications for stm experiments on quasi - one - dimensional materials .
readers mainly interested in the qualitative features of the ldos may start with this section and study the detailed results afterwards .
technical details have been moved to the appendix .
as the underlying microscopic model for the quarter - filled mott insulator one may start with the extended hubbard model  @xcite @xmath3\\    & + u\sum_j n_{j,\uparrow}\,n_{j,\downarrow}+v\sum_j n_{j}\,n_{j+1 } ,    \label{eq : exthm }    \end{split}\ ] ] where @xmath4 and @xmath5 are electron creation and annihilation operators at site @xmath6 with spin @xmath7 , @xmath8 , and @xmath9 . at quarter filling and for sufficiently large on - site and nearest - neighbor repulsions @xmath10 this model possesses a mott insulating phase  @xcite which can be thought of as a microscopic realization of the field - theoretical system we will study below .
the extended hubbard model can be used as an effective model for the description of bechgaard salts  @xcite ; its spectral and transport properties have been investigated in detail using various techniques  @xcite .    instead of working with the microscopic model we will study here the corresponding low - energy field theory .
the continuum description is obtained by focusing on the degrees of freedom around the fermi points @xmath11 .
we introduce slowly varying right- and left - moving fermi fields as @xmath12 and similarly for the electron creation operators . here
@xmath13 denotes the lattice spacing , @xmath14 , and the fermi momentum is given by @xmath15 at quarter filling .
the impurity potential is assumed to be strong such that we can model it by a boundary condition on the continuum electron field @xmath16 using standard bosonization of the right- and left - moving fermi fields it was shown  @xcite that double umklapp processes result in a cosine term in the charge sector , which is relevant at low energies and generates a gap . thus the effective low - energy hamiltonian is given by  @xcite @xmath17\nonumber\\ * & & -\frac{g_\mathrm{c}}{(2\pi)^2}\int_{-\infty}^0 dx\ , \cos\bigl(\beta\phi_\mathrm{c } ) , \label{eq : chargehamiltonian}\\ h_\mathrm{s}&=&\frac{v_\mathrm{s}}{16\pi}\int_{-\infty}^0 dx   \bigl[\bigr(\partial_x\phi_\mathrm{s}\bigr)^2 + \bigr(\partial_x\theta_\mathrm{s}\bigr)^2\bigr ] , \label{eq : spinhamiltonian}\end{aligned}\ ] ] where , as a consequence of , the canonical bose fields @xmath18 are subject to the hard - wall boundary conditions @xmath19 the fields @xmath20 are dual to @xmath18 . the charge and spin velocities @xmath21 , the luttinger parameter @xmath22 , and the coupling constant @xmath23 are functions of the hopping and interactions in an underlying microscopic model . within the field theory
they can be viewed as phenomenological parameters .
experimental estimates  @xcite for the luttinger parameter in the bechgaard salts yield @xmath24 .    in the bosonization procedure leading to the effective low - energy hamiltonian
the explicit relation between the fermi and the bose fields is given by  @xcite @xmath25\nonumber\\ * & & \qquad\qquad\times\exp\biggl[\frac{{\text{i}}}{4}f_\sigma\bigl(\phi_\mathrm{s } + \theta_\mathrm{s}\bigr)\biggr],\label{eq : bosonizationr}\\ * l^\dagger_\sigma&=&\frac{\eta_\sigma}{\sqrt{2\pi}}\ , \exp\biggl[-\frac{{\text{i}}}{4}\biggl(\frac{\beta}{2}\phi_\mathrm{c } -\frac{2}{\beta}\theta_\mathrm{c}\biggr)\biggr]\nonumber\\ * & & \qquad\qquad\times\exp\biggl[-\frac{{\text{i}}}{4}f_\sigma\bigl(\phi_\mathrm{s } -\theta_\mathrm{s}\bigr)\biggr],\label{eq : bosonizationl}\end{aligned}\ ] ] where the klein factors @xmath26 satisfy the anticommutation relations @xmath27 and @xmath28 .
the charge excitations are described by the sine - gordon model on the half - line ( [ eq : chargehamiltonian ] ) . for @xmath29
the cosine term is relevant and generates a gap .
the excitations are massive solitons and antisolitons which carry charges @xmath30 respectively .
they possess a relativistic dispersion relation ; we parametrize their energy and momentum in the usual way using the rapidity @xmath31 by @xmath32 the soliton mass @xmath33 is a function  @xcite of the bare parameters @xmath23 and @xmath22 but will be viewed here as a phenomenological parameter replacing the coupling constant @xmath23 . in the regime
@xmath34 solitons and antisolitons are the only excitations in the charge sector ; for @xmath35 propagating breather ( soliton - antisoliton ) bound states exist as well . at the luther - emery point ( lep )
@xmath36 the charge sector is equivalent to a free massive dirac theory  @xcite . the sine - gordon model ( [ eq : chargehamiltonian ] ) with the boundary condition is integrable@xcite , which we will use below to calculate the necessary correlation functions of right- and left - moving fermi fields . the spin sector ( [ eq : spinhamiltonian ] )
describes massless relativistic spinon excitations which propagate with velocity @xmath37 and carry spin @xmath0 .
we have already assumed spin rotational invariance , i.e. a possible luttinger parameter in the spin sector is set to one .
we note that the half - filled mott insulator studied previously  @xcite possesses the same effective low - energy hamiltonian . the difference to the quarter - filled case studied in this article is given  @xcite by the bosonization relations and for the right- and left - moving fermi fields .
this leads to differing green functions as well as ldos .
we will proceed with the presentation of the results for the quarter - filled mott insulator and present a detailed comparison to the half - filled system in sec .
[ sec : hfmi ] .
in order to derive the ldos below we calculate the time - ordered green function in euclidean space , @xmath38 where @xmath39 is the ground state of in the presence of the boundary and @xmath40 denotes imaginary time . at low energies we can linearize around the fermi points @xmath41 , which results in @xmath42 where e.g. @xmath43 . for the calculation of the ldos we have to set @xmath44 below .
in particular , we will focus on the spatial fourier transform of the ldos as physical properties can be more easily identified .
the four terms in the decomposition then contribute in different regions in momentum space , i.e. @xmath45 and @xmath46 contribute for @xmath47 while @xmath48 and @xmath49 contribute for @xmath50 . in the translationally invariant system right- and left - moving fields in the spin sector decouple which results in @xmath51 .
the presence of a boundary couples right and left sectors and concomitantly the fourier transform of the green function ( [ eq : gflowenergy ] ) acquires a non - zero component at @xmath47 .
we will focus on the @xmath52-part of the green function and the ldos in the following , since this provides a particularly clean way of investigating boundary effects .
we note that the other components of the green function can be calculated in the same way .
due to the spin - charge separated hamiltonian the green function @xmath53 factorizes into a product of correlation functions in the spin and charge sectors .
the correlation functions in the spin sector can be straightforwardly derived using standard methods like boundary conformal field theory  @xcite . on the other hand
, the integrability of the sine - gordon model on the half - line enables us to calculate correlation functions in the charge sector using the boundary state formalism  @xcite together with a form - factor expansion  @xcite .
this yields the following result for the @xmath52-component of the green function @xmath54 the infinite series originates from the applied form - factor expansion , which constitutes an expansion of correlation functions in the number of contributing solitons and antisolitons . in app .
[ app : gf ] we derive explicit expressions for the first three terms @xmath55 which correspond to two - particle processes in the charge sector .
we note that is valid at zero temperature .
the effect of small temperatures @xmath56 on the gapless spin sector can be incorporated using conformal field theory  @xcite as discussed for the half - filled mott insulator in ref .  .
in the next section we use the result to derive the fourier transform of the ldos .
the ldos can now be calculated directly from the green function . in order to analyze the physical properties it is useful to consider the spatial fourier transform of the ldos as features like dispersing quasiparticles
can be more easily identified .
this technique was previously applied to the ldos of luttinger liquids  @xcite as well as cdw states and half - filled mott insulators  @xcite .
the fourier transform of the ldos for positive energies @xmath57 is given by  @xcite @xmath58 here the green function has been analytically continued to real times and we have taken the limit @xmath59 .
we note that due to the assumed spin rotational invariance the ldos does not depend on @xmath60 .
as mentioned before , we will concentrate on the @xmath52-component as it vanishes in the absence of the boundary and hence offers a particularly clean way of investigating boundary effects .
for @xmath61 only @xmath62 contributes and starting from we arrive at our main result ( @xmath63 ) @xmath64 where @xmath65 as well as @xmath66 , @xmath67 , and @xmath68 .
explicit integral representations for the scattering matrix @xmath69 , the boundary reflection matrix @xmath70 , and the function @xmath71 are given in app .
[ app : ff ] .
physically @xmath69 describes the scattering of two antisolitons with relative momentum corresponding to the rapidity @xmath31 , @xmath70 is the reflection amplitude of an antisoliton with rapidity @xmath31 off the boundary [ and thus incorporates the information on the boundary condition ] , and @xmath71 appears in the matrix element [ i.e. form factor ] of the charge part of and with two - antisoliton states [ see ] .
the parameter @xmath72 is defined as @xmath73 .
the normalization constant@xcite @xmath74 has dimension @xmath75 ; thus in order to obtain dimensionless quantities we multiply the ldos by @xmath76 in all figures .
the singularity of the integrands in is smeared out by taking @xmath77 small but finite ; we choose @xmath78 unless stated otherwise .
this mimics the broadening of singularities in experiments due to the instrumental resolution and the finite temperature .
the ldos vanishes for @xmath79 as the fractionalization of an electron creates at least two antisolitons .
the three terms constitute the leading contribution from the form - factor expansion . as can be seen from the double integrals they all originate in two - particle processes in the charge sector .
the number of boundary reflection matrices @xmath70 indicates that @xmath80 involves @xmath81 interactions with the boundary , i.e. the terms contain no , one , or two reflections of charge excitations .
the sub - leading contributions involve at least three particles in the charge sector which either come from a higher number of particles in the intermediate state or from higher - order processes due to the boundary .
similar terms were thoroughly analyzed for the ldos in a cdw state at the lep  @xcite as well as the spectral function in the ising model with a boundary magnetic field  @xcite ; they were found to be negligible in both cases .
we further note that the suppression of sub - leading terms in the form - factor expansion for bulk two - point functions is a well - known feature of massive integrable field theories  @xcite .
we thus expect the higher - order corrections to to be negligible in the low - energy regime we consider here .
: @xmath82 and @xmath83 for @xmath84 and @xmath85 ( upper panel ) and @xmath86 ( lower panel ) .
both plots are shown on the same scale .
we observe a peak at @xmath87 and dispersing features at @xmath88 as well as @xmath89 . for @xmath90 the dispersing features are strongly suppressed . ]
the constant energy scan of the fourier transform of the ldos is shown in fig .  [
fig : plot1 ] .
we observe a singularity at @xmath87 , i.e. @xmath2 , which is indicative of the pinning of the @xmath52-sdw at the boundary ( see sec .  [
sec : ldos0 ] ) .
furthermore there is a peak at @xmath91 .
the constant momentum scans shown in fig .
[ fig : plot2 ] reveal that this peak shows dispersing behavior and splits above a critical momentum @xmath92 [ @xmath93 will be defined in eq .  ] .
more precisely these features follow the dispersion relations and respectively ( see sec .  [
sec : ldos1 ] ) .
the existence of dispersing features is at first surprising since the spectral function  @xcite of the translationally invariant system does not show any dispersing quasiparticles .
the processes contributing to the ldos can be thought of as arising from creating an electron at a position @xmath94 , the fractionalization of this electron into at least one spinon and two antisolitons , the propagation of these elementary excitations through the system , and their subsequent recombination at the position @xmath94 .
this allows a natural interpretation of the leading terms in the form - factor expansion in terms of propagating antisolitons and spinons . in the next sub - sections we will analyze the different terms separately .
afterwards we comment on the different findings in the spectral function and the ldos .     for @xmath85 and @xmath84 .
the curves are constant @xmath95-scans which have been offset along the y - axis by a constant with respect to one another .
the ldos is dominated by the peak at @xmath87 .
we further observe dispersing features at @xmath96 as well as @xmath97 for @xmath92 . ]
let us start with the analysis of the first term in the form - factor expansion , which is shown in fig .
[ fig : plot3 ] .
we first note that this term is dominated by a singularity at @xmath87 , i.e. @xmath2 .
this singularity is due to the pinning of the @xmath52-sdw at the boundary .
as @xmath15 we find a spatial modulation with periodicity @xmath99 as sketched  @xcite in fig .
[ fig : plot3].b . close to @xmath87 the first term in the form - factor expansion and thus the whole ldos behaves as  @xcite @xmath100 we note that the exponent is independent of the luttinger parameter @xmath22 ( see fig .  [
fig : plot3].a ) .
this singularity is similar to the ones observed in the ldos of luttinger liquids  @xcite as well as cdw states  @xcite where , however , the exponents depend on the interaction strength and thus the luttinger parameter .     for @xmath101 .
a ) logarithmic plot close to the singularity at @xmath87 for @xmath102 .
the dotted line is a guide to the eye ; it has gradient @xmath103 .
b ) sketch  @xcite of the @xmath52-sdw pinned at the boundary .
c ) @xmath104 for @xmath85 and different momenta .
we observe a dispersing feature at @xmath105 indicated by small arrows .
d ) sketch of the process giving rise to the propagating feature .
the electron decomposes into two zero - momentum antisolitons ( solid lines ) and one spinon ( dashed line ) reflected at the boundary . ]
furthermore , @xmath106 possesses a weak dispersing feature following ( see  @xcite fig .  [
fig : plot3].c ) @xmath107 as already discussed , the created electron fractionalizes into two antisolitons and ( at least ) one spinon which can propagate through the system . in this way of thinking
the dispersing feature at @xmath108 arises from two antisolitons with zero momentum ( solid lines in fig .
[ fig : plot3].d ) contributing an energy @xmath109 and a massless soliton with momentum @xmath95 ( dashed line in fig .
[ fig : plot3].d ) .
the appearance of @xmath110 in @xmath108 is due to the fact that the spinon has to propagate to the boundary and back , thus covering the distance @xmath111 in time @xmath112 .       for @xmath85 , @xmath101 , and different momenta .
we observe dispersing features at @xmath96 indicated by up - arrows and @xmath97 indicated by down - arrows .
b ) process resulting in @xmath114 .
the electron decomposes into one antisoliton and spinon with zero momentum as well as one antisoliton reflected at the boundary .
c ) process resulting in @xmath115 . ]    the next term in the form - factor expansion contains one boundary reflection matrix @xmath70 .
thus it is natural to interpret it as arising from processes involving the reflection of one antisoliton at the boundary . in fig .
[ fig : plot4].a we show @xmath116 for different momenta .
we clearly observe a dispersing feature ( indicated by up - arrows ) following @xmath117 this feature originates in the process where one antisoliton propagates with momentum @xmath95 while the other antisoliton and the spinon possess zero momentum ( see fig .
[ fig : plot4].b ) .
the first term in @xmath114 is the rest mass of the zero - momentum antisoliton while the second term is the energy of the propagating antisoliton .
furthermore , when @xmath95 exceeds the critical value @xmath93 a second dispersing feature appears ( indicated by down - arrows in fig .  [
fig : plot4].a ) at @xmath118 the critical momentum @xmath93 can be determined from the condition that the spin and charge excitations have the same group velocity @xmath119 which leads to @xmath120 thus this feature can be thought of as arising from an antisoliton with momentum @xmath93 and a spinon carrying momentum @xmath121 , while the second antisoliton stays at rest ( see fig .
[ fig : plot4].c ) .
the fact that the propagating antisoliton and the spinon have the same group velocity results in an increased probability that both simultaneously arrive at position @xmath94 after reflection at the boundary .    the observed splitting of the quasiparticle peak is reminiscent of what is found for the single - particle spectral function in the translationally invariant , half - filled mott insulator  @xcite as well as the ldos of a cdw state in the presence of a boundary in the regime of attractive interactions  @xcite .
the peak splitting is a consequence of the curvature of the antisoliton dispersion ( and thus of the charge gap ) together with the condition @xmath122 .
hence it can not be observed in the luttinger liquid case  @xcite where both sectors are massless .       for @xmath85 , @xmath101 , and different momenta .
we observe dispersing features at @xmath124 ( up - arrows ) and @xmath125 ( down - arrows ) .
processes resulting in ( b ) @xmath126 and in ( c ) @xmath127 . ]
the leading term in second order in the boundary reflection matrix is @xmath128 . as in the previous term
we find two dispersing features .
first , indicated by up - arrows in fig .
[ fig : plot5].a , there is a feature following @xmath129 this feature can be thought of as arising from the propagation of both antisolitons with momenta @xmath130 and thus the same velocity , while the spinon has zero momentum ( see fig .
[ fig : plot5].b ) .
the second feature is found at @xmath131 provided @xmath95 exceeds the critical momentum @xmath132 .
this feature is indicated by down - arrows in fig .
[ fig : plot5].a .
the interpretation is that both antisolitons carry momentum @xmath93 while the excess momentum @xmath133 is carried by the spinon ( see fig .
[ fig : plot5].c ) .
we note that at finite momentum @xmath91 the term @xmath116 dominates the ldos as can be seen by comparing the scales in figs .
[ fig : plot3].c ,  [ fig : plot4].a , and  [ fig : plot5].a , respectively .
more importantly , the dispersing features originating in @xmath116 are much more pronounced than those resulting from @xmath106 and @xmath128 .
we attribute this to the fact that , at least in the naive interpretation of figs .
[ fig : plot4].b and c , the antisolitons involved in processes contributing to @xmath116 possess different momenta as one of them is at rest while the other propagates through the system .
in contrast , the dispersing features in @xmath106 and @xmath128 originate from processes in which the two antisolitons possess equal momenta ( see figs .  [
fig : plot3].d as well as  [ fig : plot5].b and c , respectively ) .
thus their relative rapidity vanishes and the scattering matrix is @xmath134 ( see app .
[ app : ff ] ) .
the resulting `` quasi - fermionic '' behavior leads to a strong suppression of the dispersing features in @xmath106 and @xmath128 via the pauli principle .    finally let us comment on the different findings in the spectral function  @xcite of the translationally invariant system and the ldos in the presence of a boundary .
the spectral function is obtained from the green function @xmath135 via fourier transformation with respect to @xmath136 and @xmath94 .
the contributing processes require the propagation of all excitations , and in particular both antisolitons , from @xmath137 to @xmath138 .
thus one may expect a dispersing feature at @xmath139 from the antisolitons propagating with the same momentum @xmath130 .
in addition , above the critical momentum @xmath140 a linearly dispersing feature following @xmath141 is expected .
indeed these are exactly the conditions  @xcite for the threshold below which the spectral function vanishes . above the threshold
the spectral function is rather featureless , in particular there are no dispersing peaks associated with the propagation of charge or spin degrees of freedom .
the absence of dispersing features may be attributed to the `` quasi - fermionic '' behavior of antisolitons , as all processes contributing to the spectral function will contain antisolitons with equal momenta .
in contrast , the processes contributing to the ldos are given by the fractionalization of an electron at position @xmath94 , the propagation of the two antisolitons and the spinon through the system , and the subsequent recombination at the original position @xmath94 . in particular , processes where the momentum @xmath95 is solely carried by one of the antisolitons or a spinon - antisoliton pair will contribute , resulting in the dispersing features at @xmath114 and @xmath115 , respectively .
the observation of these features thus underlines that the spectral function and the ldos probe rather different physical processes .
while the spectral function requires the propagation of all excitations between two points in the system , the ldos also probes processes where only one or two excitations propagate to the impurity and back .
up to now we have concentrated on the boundary conditions for the electron field or equivalently for the canonical bose fields . in this section
we will consider more general phase shifts for the electron fields which lead to arbitrary dirichlet boundary conditions for the bose fields .
as is well known , such boundary conditions can result in the existence of boundary bound states ( bbs s ) , which are expected  @xcite to manifest themselves in the ldos as features within the spectral gap . here
we will focus on general dirichlet boundary conditions in the charge sector @xmath142 while in the spin sector we still assume @xmath143 . in the boundary state formalism  @xcite applied here the information about the boundary condition
is encoded in the boundary reflection matrix @xmath70 .
in particular , the existence of a bbs will show up as a pole in the physical strip @xmath144 . as was shown in ref
.   this will be the case provided @xmath145 the pole at @xmath146 corresponds to a bbs with energy @xmath147 when calculating the green function the pole of the boundary reflection matrix will lead to two additional terms , @xmath148 and @xmath149 , in the form - factor expansion .
this in turn yields two additional terms in the ldos .     for @xmath150 , @xmath151 , and @xmath84 .
the curves are constant @xmath95-scans which have been offset along the y - axis by a constant with respect to one another . the existence of a bbs results in spectral weight within the two - soliton gap @xmath79 . ]    to be explicit , let us consider the lep @xmath36 in the following .
the three leading terms in the form - factor expansion of the ldos are still given by where the boundary reflection matrix now depends explicitly on @xmath152 @xmath153 this obviously possesses a pole in the physical strip @xmath144 provided @xmath154 corresponding to the existence of a bbs .
the resulting two additional terms in the ldos read ( @xmath155 ) @xmath156 where @xmath157 as well as @xmath158 and @xmath159 .
we note that @xmath160 has the dimension of an inverse length .
the spatial dependence of the terms @xmath148 and @xmath149 ( see app .
[ app : gf ] ) is dominated by @xmath161 ( @xmath162 ) , thus it is natural to interpret @xmath163 as the width of the bbs .
the terms possess a threshold at @xmath164 corresponding to the creation of the boundary bound state and an additional antisoliton .
interestingly , although @xmath165 and @xmath166 separately possess jumps at the threshold , the sum of both terms is continuous .
this is reminiscent of the bbs contributions in the ising model with a boundary magnetic field  @xcite , where a cancellation of singularities yields a smooth spectral function at the threshold .
the ldos in the presence of a bbs is shown in fig .
[ fig : plot6 ] .
we clearly observe a contribution within the gap @xmath79 .
in contrast to bbs s in the cdw state  @xcite there is , however , no singularity at @xmath167 as the underlying processes involve the creation of a spinon , an antisoliton , and the bbs . as before the ldos is dominated by a singularity at @xmath87 and we observe dispersing features following and respectively .
the ldos of a half - filled mott insulator in the presence of a boundary was analyzed in detail  @xcite in ref .  .
as already mentioned , the effective low - energy hamiltonian is also given by
but the expressions for the right- and left - moving fermi fields in terms of the bose fields differ from and . from the physical point of view right- and left - moving fermi fields
create and annihilate individual antisolitons in the half - filled case whereas they create and annihilate antisoliton pairs in the quarter - filled system .
thus the gap in the ldos is given by @xmath33 in the former and @xmath109 in the latter system .
apart from this we find the following similarities and differences .    in both cases
the ldos possesses a singularity at @xmath87 which behaves as @xmath168 independently  @xcite of the luttinger parameter @xmath22 .
this singularity is due to the pinning of the @xmath52-sdw at the boundary , which possesses a wave length @xmath169 or @xmath170 respectively .    beside this singularity the ldos of
the half - filled mott insulator shows dispersing features at @xmath171 , at @xmath172 , and , for @xmath92 , at @xmath173 .
an interpretation is obtained by noting that the electron decomposes into one spinon and one antisoliton , which then propagate through the system . as the leading processes involve only individual antisolitons , the `` quasi - fermionic '' behavior of the antisolitons discussed at the end of sec .
[ sec : ldos2 ] does not affect the ldos .
in contrast , in the quarter - filled system we observe only two dispersing modes following and respectively ( see fig .
[ fig : plot2 ] ) .
all other dispersing features discussed above are strongly suppressed due to the `` quasi - fermionic '' behavior of the antisolitons .
finally , in the presence of a bbs with energy @xmath174 the ldos of the half - filled mott insulator possesses a non - dispersing singularity with @xmath175 which originates from the creation of a spinon and the bbs .
in contrast , in the quarter - filled case the ldos is smooth at the threshold @xmath164 which originates from the simultaneous creation of one antisoliton and the bbs .
in this section we compare qualitative aspects of the ldos of various one - dimensional systems in the presence of a boundary or strong impurity potential . we discuss the spatial fourier transform @xmath176 defined in .
specifically we consider luttinger liquids  @xcite , half - filled mott insulators ( or cdw states )  @xcite , and the quarter - filled mott insulator investigated in sec .
[ sec : ldos ] .
we focus on the @xmath61-component of the ldos as it vanishes in the translationally invariant systems and thus provides a particularly clean way to investigate boundary effects .
( we note that the ldos in the small - momentum regime @xmath50 behaves qualitatively similar . )
the first feature , which all three systems have in common , is a singularity of the ldos @xmath177 at @xmath87 , which originates in the pinning of a charge or spin density wave at the impurity . in the case of a luttinger liquid
the strength of this singularity depends on the luttinger parameter @xmath178 and thus on the interactions between the electrons , whereas in both the half - filled and quarter - filled mott insulator the ldos behaves as  @xcite @xmath179 .
furthermore , as was pointed out in ref .  , the gap of the ldos in the half - filled mott insulator is given by the soliton mass @xmath33 , which equals the thermal activation gap and is half of the gap observed in optical measurements .
in contrast , in the quarter - filled system the ldos possesses a gap of @xmath109 , which is equal to the optical gap but twice the thermal gap .
of the ldos in different one - dimensional systems .
a ) in a luttinger liquid one finds two linearly dispersing modes .
b ) in a half - filled mott insulator ( mi ) there exists a gap @xmath33 .
one observes one linear and one curved dispersion . from the latter a third feature splits off at @xmath180 .
c ) the quarter - filled mott insulator possesses a gap @xmath109 .
one further observes one curved dispersion from which a linear mode splits off at @xmath180 . in both ( b ) and ( c )
the thermal activation gap is given by @xmath33.,width=321 ]    in addition , the ldos possesses various dispersing features ( see fig .  [
fig : plot7 ] ) . in the luttinger liquid case one observes two linearly dispersing modes corresponding to individually propagating spin and charge degrees of freedom respectively ( see fig .  [
fig : plot7].a ) . in the half - filled mott insulator
an electron decomposes into massless spin excitations and one charge excitation with mass @xmath33 .
propagating spin excitations give rise to a linearly dispersing mode with velocity @xmath37 , while the propagating charge excitation results in a curved dispersion relation .
in addition , there exists a critical momentum @xmath180 at which the group velocity of the charge excitation equals @xmath37 . at this momentum a third linearly dispersing feature splits away from the charge mode , which has its origin in the collective propagation of spin and charge degrees of freedom with equal velocities ( see fig .
[ fig : plot7].b ) .
we note that the qualitative behavior discussed for the half - filled mott insulator is also expected in other systems possessing gapless excitations in one sector and gapped ones in the other like , for instance , cdw states .
finally , in the quarter - filled mott insulator an electron decomposes into massless spin excitations and two charge excitations with masses @xmath33 .
this results in a curved dispersion originating in the propagation of one charge excitation while the second one stays at the position of the stm tip .
this dispersion splits at the critical momentum @xmath93 due to the collective propagation of one charge excitation and additional spin degrees of freedom ( see fig .  [
fig : plot7].c ) .
we note that one does not observe individually propagating spin excitations .
the relevant processes would require the two charge excitations to possess equal momenta , which is forbidden by their `` quasi - fermionic '' behavior ( see sec .
[ sec : ldos2 ] ) .    to sum up ,
at sufficiently large momenta @xmath95 we observe two linearly dispersing modes in a luttinger liquid , two linear and one curved dispersion relation in a half - filled mott insulator , and one linear and one curved dispersion in a quarter - filled mott insulator .
this clearly demonstrates that characteristic properties of the bulk state of matter and its electronic excitations can be infered from studying spatial modulations of the ldos in the presence of an impurity .
the ldos discussed in this article is directly related to the tunneling current measured in stm experiments .
if the density of states in the stm tip is approximately independent of energy , the tunneling conductance will be given by the thermally smeared ldos of the sample at the position of the tip  @xcite @xmath181 where @xmath182 denotes the fermi function .
@xmath176 can then easily be obtained via spatial fourier transformation . in this context
the models discussed above apply to quasi - one - dimensional materials at energies above the cross - over scale to three - dimensional behavior .
this situation might be experimentally realized for example in stripe phases of high - temperature superconductors  @xcite , carbon nanotubes  @xcite , bechgaard and fabre salts  @xcite , self - organized atomic gold chains  @xcite , and two - leg ladder materials  @xcite .
as discussed above already the knowledge of qualitative properties such as the presence of a spectral gap and its relation to the thermal activation gap as well as the number and positions of singularities can reveal substantial information about the electronic properties and excitations of quasi - one - dimensional materials . a more detailed analysis of quantitative features like the strength of singularities or the precise positions of dispersing modes can even lead to estimates for the luttinger parameter and the velocities of charge and spin excitations .
there are , of course , experimental limitations .
obviously the experimental resolution and thermal broadenings have to be sufficiently small to clearly distinguish the different singularities and dispersing modes .
moreover , the simplyfied relation does not contain  @xcite the matrix elements for the tunneling from the sample into the tip , which may introduce a non - trivial spatial dependence of the tunneling conductance on the ldos and thus hamper a natural interpretation in terms of propagating excitations .
in conclusion , we have studied a quarter - filled mott insulator in the presence of an impurity , which was modeled by a boundary condition for the electron fields . in this system
we calculated the @xmath52-component of the spatial fourier transform of the ldos using the boundary state formalism together with a form - factor expansion .
the ldos is dominated by a singularity at @xmath2 , which is indicative of a pinning of a sdw at the impurity .
furthermore , we observed several dispersing features above the two - soliton gap @xmath79 .
we identified their physical origin as spin and charge degrees of freedom , which propagate between the stm tip and the impurity , and attributed their relative strength to the `` quasi - fermionic '' behavior of charge excitations with equal momenta .
this also explains the absence of dispersing features in the spectral function  @xcite of the translationally invariant system and underlines that the spectral function and the ldos probe rather different physical processes .
we further showed that a bbs in the charge sector leads to a non - vanishing ldos within the two - soliton gap .
finally , we discussed our results in the context of other one - dimensional systems , i.e. luttinger liquids and half - filled mott insulators , and argued that the measurement of the spatial modulations in the ldos can be used to extract detailed informations about the electronic states of quasi - one - dimensional systems .
i would like to thank fabian essler , markus garst , markus morgenstern , alexei tsvelik , and particularly volker meden for valuable discussions and comments .
this work was supported by the german research foundation ( dfg ) through the emmy noether program .
in this appendix we derive the leading terms @xmath55 in the form - factor expansion of the charge part of the green function , i.e. we determine the correlation function @xmath183 in the sine - gordon model on the half line . here
the operators @xmath184 and @xmath185 are given by the charge parts of and respectively .
we calculate using the boundary state formalism  @xcite together with a form - factor expansion  @xcite .
the same technique was previously applied  @xcite to determine the green function of a cdw state in the presence of an impurity . here
we restrict ourselves to the main steps of the derivation and concentrate on the differences as compared to ref .  ,
to which we refer implicitly for all details and definitions not discussed here .    by performing a rotation in euclidean space we map  @xcite the problem onto the sine - gordon model on the real axis .
the boundary condition is translated into an initial condition encoded in a boundary state . expanding the latter in powers of the boundary reflection matrix and inserting a resolution of the identity between the operators in we find @xmath186 where we have defined the auxiliary functions    @xmath187 & \hspace{-20mm}\times { \left\langle0\right|}o_r(\tau , x_1 ) { \left|\theta_{n},\ldots,\theta_{1}\right\rangle}_{c_n,\ldots , c_1}\ ; ^{c_1,\ldots , c_n}\!{\left\langle\theta_{1},\ldots,\theta_{n}\right| } o_l^\dagger(0,x_2 ) { \left|-\theta'_1,\theta'_1,\ldots,-\theta'_m,\theta'_m\right\rangle } _
{ a_1,b_1,\ldots , a_m , b_m}. \end{split}\ ] ]    here the states @xmath188 ( and so on ) are scattering states of solitons and antisolitons , the indices @xmath189 take the values @xmath190 and label the corresponding u(1 ) charges , and energy and momentum of solitons and antisolitons are given in terms of the rapidities @xmath191 and @xmath192 by .
we note that constitutes an expansion in two parameters : ( i ) the number of particles ( with rapidities @xmath191 ) in the intermediate state and ( ii ) the number of reflections of particles at the boundary [ the number of boundary reflection matrices @xmath193 ) ] .
the scattering of solitons or antisolitons with rapidities @xmath194 and @xmath195 is encoded in the scattering matrix  @xcite @xmath196 .
the boundary reflection matrix  @xcite @xmath197 similarly describes the reflection at the boundary , @xmath198 , where @xmath199 for @xmath200 .
for dirichlet boundary conditions one has  @xcite @xmath201 .
furthermore , @xmath202 , where @xmath203 and @xmath204 are the hamiltonian and the total momentum of the sine - gordon model on the infinite line . the operator @xmath205 ( @xmath185 ) creates ( annihilates ) two antisolitons , i.e. increases ( decreases ) the u(1 ) charge by 2 .
the necessary form factors were derived by lukyanov and zamolodchikov  @xcite . in our conventions
they are given by @xmath206 where the normalization constant @xmath74 and the auxiliary function @xmath71 are stated in app .
[ app : ff ] .
the form factors satisfy the form - factor axioms as stated in ref .
with lorentz spin @xmath207 and semi - locality factor @xmath208 . in the following
we will evaluate the terms @xmath209 , @xmath210 , and @xmath211 which result in @xmath212 , @xmath213 , and @xmath214 , respectively .      the first term to be evaluated
. it does not contain any information about the boundary in the charge sector , i.e. it contains no boundary reflection matrices .
a similar term was calculated in ref .   in the derivation of the spectral function in the translationally invariant system .
after shifting the contour of integration , @xmath216 , we directly obtain @xmath217 where the center - of - mass coordinates are defined by @xmath218 and @xmath219 .
analytical continuation @xmath220 to real times , taking the limit @xmath221 , and fourier transformation with respect to @xmath112 and @xmath94 yield the first term @xmath222 in the form - factor expansion of the ldos .
we start with @xmath210 defined in .
as the operator @xmath205 creates two antisolitons we deduce @xmath224 ; @xmath225 yields @xmath226 .
the matrix element of @xmath185 contains incoming and outgoing particles and thus possesses kinematical poles .
we deal with these terms following smirnov  @xcite and analytically continue the form factor as ( see ref .   for a detailed discussion of this procedure )
@xmath227 here we have already omitted the positive imaginary parts of the rapidities in the form factors in the 2nd to 5th lines as they do not possess kinematical poles .
inserting into @xmath210 yields to different contributions : the first line leads to a term containing three integrations over rapidities and is thus a sub - leading correction ; we will not evaluate it here .
the 2nd to 5th lines yield using the boundary cross - unitarity condition  @xcite @xmath228 as well as unitarity  @xcite @xmath229 @xmath230 in order to proceed let us first assume @xmath231 .
this implies that the boundary reflection matrix does not depend on the u(1 ) charge , i.e. @xmath232 , where @xmath70 is explicitly stated in app .
[ app : ff ] . in particular
, @xmath233 is analytic in the physical strip @xmath234 .
now we shift the contours of integration @xmath235 ( we note that the contribution originating in the first term of does not possess any poles in the physical strip @xmath236 ) and apply the crossing relation in the second form factor , @xmath237 where the phase factor @xmath238 is due to the semi - locality  @xcite of the operator @xmath239 with respect to the fundamental fields creating the excitations .
finally we use to obtain @xmath240 fourier transformation directly yields the second term @xmath241 of the ldos .    on the other hand , if @xmath242 , the boundary reflection matrix @xmath243 possesses a pole at @xmath244 . at the lep
the residue is given by @xmath245 ( see ) .
thus when shifting the contours of integration in we obtain the additional bbs contribution @xmath246 fourier transformation yields @xmath247 .
finally we calculate the leading term involving the reflection of two antisolitons , which originates from @xmath211 .
the analytic continuation of the form factor of @xmath185 reads @xmath249 where the scalar products are evaluated using the faddeev - zamolodchikov algebra  @xcite @xmath250 .
\end{split}\ ] ] therefore the 2nd to 5th line in yield ( we also use the scattering axiom for the first form factor ) @xmath251 let us again first assume @xmath231 . if we shift the contours of integration @xmath252 and use the crossing relation as well as a lorentz transformation , @xmath253 we obtain @xmath254 @xmath255 is now readily obtained via fourier transformation .
starting from the second term @xmath256 related to the bbs is obtained by picking up the poles of @xmath257 and @xmath258 at @xmath259 respectively . at the lep
the result is @xmath260
in this appendix we collect some formulas on the sine - gordon theory used in the present article .
the scattering matrix between solitons and antisolitons possesses the integral representation  @xcite @xmath261,\ ] ] where @xmath73 .
a similar integral representation for the boundary reflection matrix reads  @xcite ( for dirichlet boundary conditions with @xmath262 ) @xmath263,\\    \sigma(\theta)\!\!\!&=&\!\!\!\exp\!\left[2{\text{i}}\int_0^\infty\frac{dt}{t}\ ,      \sin\frac{\theta t}{\pi\xi}\ ,      \frac{\sinh\bigl(\frac{t}{\xi}(1+{\text{i}}\frac{\theta}{\pi})\bigr ) }      { \sinh t\,\cosh\frac{t}{\xi}}\right].\end{aligned}\ ] ] the auxiliary function @xmath71 appearing in the form factors is given by  @xcite @xmath264,\\    \mathcal{c}_1&=&\exp\left[-\int_0^\infty\frac{dt}{t}\ ,      \frac{\sinh^2\frac{t}{2}\,\sinh\bigl(t(\xi-1)\bigr ) }      { \sinh(2t)\,\cosh t\,\sinh(t\xi)}\right],\end{aligned}\ ] ] while for the normalization factor @xmath74 one finds numerically from the integral repesentation derived in ref .
@xmath265 the expressions at the lep are obtained by setting @xmath266 ( @xmath36 ) , i.e. @xmath267 , @xmath268 , @xmath269 , and @xmath270 .
j.  e. hoffman , k. mcelroy , d .- h .
lee , k.  m. lang , h. eisaki , s. uchida , and j.  c. davis , science * 297 * , 1148 ( 2002 ) ; k. mcelroy , r.  w. simmonds , j.  e. hoffman , d .- h .
lee , j. orenstein , h. eisaki , s. uchida , and j.  c. davis , nature * 422 * , 592 ( 2003 ) ; c. howald , h. eisaki , n. kaneko , m. greven , and a. kapitulnik , phys .
b * 67 * , 014533 ( 2003 ) ; m. vershinin , s. misra , s. ono , y. abe , y. ando , and a. yazdani , science * 303 * , 1995 ( 2004 ) .
m. fabrizio and a.  o. gogolin , phys .
b * 51 * , 17827 ( 1995 ) ; s. eggert , h. johannesson , and a. mattsson , phys .
rev . lett . * 76 * , 1505 ( 1996 ) ; s. eggert , _ ibid . _ * 84 * , 4413 ( 2000 ) ; m. guigou , t. martin , and a. crpieux , phys . rev .
b * 80 * , 045420 ( 2009 ) .
f. mila and x. zotos , europhys .
* 24 * , 133 ( 1993 ) ; k. penc and f. mila , phys .
b * 49 * , 9670 ( 1994 ) ; k. sano and y. no , j. phys . soc . jpn . * 63 * , 1250 ( 1994 ) ; phys . rev .
b * 70 * , 155102 ( 2004 ) ; s. ejima , f. gebhard , and s. nishimoto , europhys . lett . * 70 * , 492 ( 2005 ) ; k. sano and y. no , phys . rev .
b * 75 * , 113103 ( 2007 ) .
d. controzzi , f.  h.  l. essler , and a.  m. tsvelik , phys .
86 * , 680 ( 2001 ) ; f.  h.  l. essler , f. gebhard , and e. jeckelmann , phys .
b * 64 * , 125119 ( 2001 ) ; s. andergassen , t. enss , v. meden , w. metzner , u. schollwck , and k. schnhammer , _ ibid .
_ * 73 * , 045125 ( 2006 ) ; t. shirakawa and e. jeckelmann , _ ibid .
_ * 79 * , 195121 ( 2009 )
.      a. schwartz , m. dressel , g. grner , v. vescoli , l. degiorgi , and t. giamarchi , phys . rev .
b * 58 * , 1261 ( 1998 ) ; m. dressel , k. petukhov , b. salameh , p. zornoza , and t. giamarchi , _ ibid . _
* 71 * , 075104 ( 2005 ) .
b. korin - hamzi   , e. tafra , m. basleti   , a. hamzi   , and m. dressel , _ ibid . _
* 73 * , 115102 ( 2006 ) .
f.  h.  l. essler and r.  m. konik , in _ from fields to strings : circumnavigating theoretical physics ( ian kogan memorial collection ) _ , edited by m. shifman , a. vainshtein , and j. wheater ( world scientific , singapore , 2005 ) , vol
.  i.                                j. voit , eur .
j.  b * 5 * , 505 ( 1998 ) ; f.  h.  l. essler and a.  m. tsvelik , phys . rev .
b * 65 * , 115117 ( 2002 ) ; f.  h.  l. essler and a.  m. tsvelik , phys .
* 90 * , 126401 ( 2003 ) ; h. benthien , f. gebhard , and e. jeckelmann , _ ibid . _ * 92 * , 256401 ( 2004 ) ; t. ulbricht and p. schmitteckert , europhys .
* 89 * , 47001 ( 2010 ) .
p. abbamonte , g. blumberg , a. rusydi , a. gozar , p.  g. evans , t. siegrist , l. venema , h. eisaki , e.  d. isaacs , and g.  a. sawatzky , nature * 431 * , 1078 ( 2004 ) ; s. notbohm , p. ribeiro , b. lake , d. a. tennant , k. p. schmidt , g. s. uhrig , c. hess , r. klingeler , g. behr , b. bchner , m. reehuis , r. i. bewley , c. d. frost , p. manuel , and r. s. eccleston , phys .
* 98 * , 027403 ( 2007 ) .