the emerging field of optomechanics concerns with the study of the mechanical effects of light on mesoscopic and macroscopic mechanical oscillators .
this phenomenon has been realized in the optomechanical systems consisting of an optical cavity with a movable end - mirror or with a membrane - in - the middle .
the radiation pressure exerted by the light inside the optical cavity couples the moving mirror or the membrane which acts as a mechanical oscillator to the optical field .
this optomechanical coupling has been employed for a wide range of applications such as the cavity cooling of microlevers and nanomechanical resonators to their quantum mechanical ground state @xcite , producing high precision detectors for measuring weak forces and small displacements and also providing a good approach for fundamental studies of the transition between the quantum and the classical world @xcite .
on the other hand an alternative path to the studies of cavity optomechanics has been provided experimentally by systems consisting of ultracold atomic ensembles trapped in optical cavities @xcite where the excitation of a collective mode of the cold gas plays the role of the vibrational mode of the mirror .
the standing electromagnetic wave inside the cavity forms a periodic potential , the so - called optical lattice , in which the cold atoms exhibit phenomena typical of solid state physics like the formation of energy bands and bloch oscillations @xcite .    in such systems with high finesse cavities
the atom - light interaction is enhanced because the atoms are collectively coupled to the same optical mode .
besides , in the dispersive regime where the laser pump is far detuned from the atomic resonance the excited electronic state of the atoms can be adiabatically eliminated and consequently the only degrees of freedom of atoms will be their mechanical motions @xcite . for low photon numbers when the optical grating produced by the intracavity optical field is very shallow , one can approximately restrict the dynamics to the first two motional modes @xcite . in more recent theoretical investigations @xcite
it has been shown that the simple two - mode model of a bose - einstein condensate ( bec ) can be improved by considering higher motional modes .    in spite of similarities between the two kinds of optomechanical systems ( with a moving mirror and with a bec ) ,
there are some essential differences between them .
firstly , in contrast to the moving mirror of the optomechanical sytems which is placed in the harmonic potential of its spring , the excitation modes of the bec are not based on the presence of such an external harmonic potential @xcite .
secondly , their parameters are realized in different regimes ; the frequency of the excitation mode of the bec , i.e. , recoil frequency , is well below the dipole coupling strength while the oscillation frequency of a moving mirror is of the same order of magnitude as the coupling strength @xcite . on the other hand in addition to pumping the cavity from one of the end mirrors , in the cavities equipped with atomic gas it is also possible to pump the atoms from the side of the cavity@xcite .
one of the most important characteristics of the many - body systems is the two - body collision which can affect the properties of the system . so in order to study the dynamics of a bec gas in a realistic experimental situation
, it is necessary to take it into account . in the optomechanical systems containing the atomic gas
there are two kinds of atom - atom interactions which are the origins of the system nonlinearities .
firstly , due to the atom - photon interaction the potential acting on the condensate depends in a highly nonlocal and nonlinear way on the condensate itself which leads to the long range atomic interaction mediated by the cavity field @xcite . on the other hand
there is an intrinsic nonlinearity due to the s - wave scattering which can take place at arbitrary momentum values and causes a broadening of the atomic momentum distribution @xcite due to the intraband transitions and also can scatter atoms to the other bands ( interband transitions)@xcite .
furthermore , there is another kind of nonlinearitiy due to the finiteness of the particle number of the system which can be neglected in the thermodynamic limit where the total number of atoms is very large @xcite .    in this work , we are going to extend the two - mode model considering the effects of atom - atom interaction for a one - dimensional bec in an optical cavity .
we consider the system in the low photon regime but will not restrict our treatment to the weak atom - atom interaction .
in fact the atomic collisions can populate several nonzero quasimomenta of the energy bands ( intraband transition ) and can also cause the atoms to be scattered to the other bands ( interband transitions ) . in order to have a simplified optomechanical model we take into account just the lowest nozero quasimomenta in the first brillouin zone .
we show that in the bogoliubov approximation the atom - atom interaction not only shifts the energies of the excited modes but also play the role of an oprtical parametric amplifier ( opa ) for the bogoliubov side mode which can affect its normal - mode splitting behaviour @xcite . on the other hand due to the atomic collisions the resonance frequency of the cavity
is shifted which leads to the decrease of the number of cavity photons and the depletion of the bogoliubov mode .
these results are in good agreement with those obtained by the numerical solutions of the full description of gross - pitaevskii equation ( gpe ) @xcite .
the paper is structured as follows . in sec . ii we will first give a thorough theoretical description of the many - body system under consideration and then derive a simplified optomechanical model . in sec .
iii the quantum langevin equations ( qles ) are derived and linearized around the semiclassical steady state . in sec .
iv we study the mean - field solutions and fluctuations of the system and in sec .
v we derive the spectrum of the bogoliubov mode and investigate the normal mode splitting ( nms ) .
finally , our conclusions are summarised in sec .
and the decay rate is @xmath0.,width=3 ]    we are going to study a gas of @xmath1 ultracold bosonic two - level atoms with mass @xmath2 and transition frequency @xmath3 inside the optical lattice of a single - mode , high - finesse fabry - perot cavity with length @xmath4 .
the cavity is driven at rate @xmath5 through one of its mirrors by a laser with frequency @xmath6 , and wavenumber @xmath7 .
we assume the bec is confined in a cylindrically symmetric trap with a transverse trapping frequency @xmath8 and negligible longitudinal confinement along the @xmath9 direction ( fig.[fig : fig1 ] ) . in this way we can describe the dynamics within an effective one - dimensional model by quantizing the atomic motional degree of freedom along the @xmath9 axis only .      in the dispersive regime where the laser pump is far detuned from the atomic resonance ( @xmath10 exceeds the atomic linewidth @xmath11 by orders of magnitude ) ,
the excited electronic state of the atoms can be adiabatically eliminated and spontaneous emission can be neglected @xcite . in the frame rotating at the pump frequency , the many - body hamiltonian reads @xmath12 where @xmath13 is the annihilation operator for a cavity photon and @xmath14 is the cavity - pump detuning .
@xmath15 is the single - particle hamiltonian of an atom inside the optical lattice of the cavity and @xmath16 is the atom - atom interaction that are respectively given by    [ hamil1 ] @xmath17    here @xmath18 is the optical lattice barrier height per photon and represents the atomic backaction on the field , @xmath19 is the vacuum rabi frequency , @xmath20 and @xmath21 is the two - body s - wave scattering length @xcite . the second term in the hamiltonian of eq.([hamil1]a ) is a periodic potential of period @xmath22 .
it is well - known that the eigenfunctions of a particle inside a periodic potential are the bloch functions @xmath23 with eigenvalues @xmath24 where @xmath25 is the bloch band index and @xmath26 $ ] is the quasimomentum of the particle @xcite . if there are @xmath27 periods inside the cavity , then @xmath28 .
using the born - von karman periodic boundary condition @xmath29 and based on the bloch theorem , @xmath30 it is deduced that @xmath31 , where @xmath32 . in this way the quasimomentum @xmath33 in the first brillouin zone can be written in terms of optical wave number , i.e. , @xmath34 where @xmath35 . now using the bloch theorem ,
the eigenfunctions @xmath36 can be written as @xmath37 here @xmath38 is a periodic function with the optical lattice period @xmath39 and reads @xmath40 by substituting eq.([u ] ) into eq.([bloch2 ] ) , the eigenfunctions of the single - particle hamiltonian @xmath15 is obtained as @xmath41    the second quantized atomic wave field @xmath42 can be expanded as a fourier series in terms of the bloch functions , i.e. , @xmath43 here @xmath44 ( @xmath45 ) is the annihilation ( creation ) operator for the atomic field that annihilates ( creates ) a particle in a state determined with the bloch band index @xmath25 and quasimomentum @xmath33 . by using this expansion
the hamitlonian of the system is diagonalized and one can investigate the bloch band structure of the condensate @xcite .
however , the fourier coefficients @xmath46 can only be determined numerically and there is no closed form for the bloch functions . instead , by substituting eq.([blochf ] ) into eq.([psi ] ) one can obtain a fourier expansion of the atomic field in therms of plane waves @xmath47 in the derivation of eq.([psi2 ] ) we have done a bogoliubov transformation @xmath48 where @xmath49 ( @xmath50 ) is the atomic field operator that annihilates ( creates ) a particle in a state determined with the band index @xmath51 and quasimomentum @xmath33 .    by substituting eq.([psi2 ] ) into eq.([h1 ] ) the hamiltonian of the systems can be written as @xmath52 in these sums @xmath53 , @xmath35 , and @xmath54 are the energy eigenvalues in which @xmath55 is the recoil energy .
as seen , the field expansion of eq.([psi2 ] ) can no longer diagonalize the hamiltonian of the system .
in fact that part of hamiltonian which corresponds to the atom - photon interaction remains nondiagonal . due to the interaction with
the optical filed , the atoms are scattered to the nearest bands ( @xmath56 ) while their quasimomentum ( @xmath57 ) remains unchanged .      since the optical potential is symmetric with respect to the origin , the hamiltonian of the system has the parity symmetry .
so if the sytem starts from a homogeneous bec , only the cosine parts of the exponential functions in eq.([psi2 ] ) are excited because of the parity conservation . in the limit of weak photon - atom coupling ,
when either the photon number inside the cavity or @xmath58 is small , the lowest bands @xmath59 can be excited by fluctuations resulting from the atom - light interaction @xcite .
on the other hand , the s - wave scattering populates fluctuations with arbitrary @xmath51 and @xmath57 @xcite . in a very simplified optomechanical model
one can consider the scattering from @xmath60 to @xmath61 in the lowest band @xmath62 . in this way the atomic field operator ( eq.([psi2 ] ) ) can be truncated as @xmath63 where we have done the following bogoliubov transformations    @xmath64    in the case that the system does not have parity symmetry , for example when the bec is inside a ring cavity , one should also consider sine modes with annihilation operators @xmath65 which in our model have been set aside @xcite . by substituting the atomic field operator , eq.([opaf ] ) , into the hamiltonian of eq.([h1 ] ) , one can arrive at the following hamiltonian @xmath66 where    [ ac aa ] @xmath67    the first two terms in eq.([opa ] ) denote the cavity and the pump hamiltonians .
the third term is the energy of the side - modes @xmath68 and @xmath69 which is just the third term in eq.([h2 ] ) with energy eigenvalues @xmath70 given by eq.([epsilon ] ) .
the hamiltonian @xmath71 denotes the atom - photon interaction and @xmath16 is the atom - atom interaction .
furthermore , @xmath72 is the s - wave scattering frequency and @xmath73 is the waist of the optical potential .
this hamiltonian can be further simplified since for weak optical fields and large @xmath74 the depletion of the initial condensate remains weak .
so we can use the bogoliubov approximation and treat the zero - momentum mode classically , i.e. , @xmath75 @xcite . in this way
eqs.([opa ] ) reduce to the following forms    [ hb ] @xmath76    here @xmath77 is the effective stark - shifted detuning . as seen from the hamiltonian of eq.([hb]b ) there are two kinds of optomechanical coupling .
the first term of eq.([hb]b ) is the linear radiation pressure which couples the side mode @xmath68 to the optical field with the optomechanical coupling constant @xmath78 while the second term is the nonlinear optomechanical coupling of the two side modes @xmath68 and @xmath69 with the optical field . in the atom - atom interaction hamiltonian @xmath16
we have disregarded all terms proportional to @xmath79 .
the influence of s - wave scattering has partly appeared as a shift in the side modes energies ( the third and fourth tems in eq.([subhb ] ) ) . since @xmath80 , where @xmath81 is the density of atoms , it is apparent that the energy shits obtained here is the same as that of bogoliubov theory @xcite . on the other hand
the role of @xmath16 in eq.([hb ] ) for the bogoliubov side modes is very similar to that of an optical parametric amplifier ( opa ) in an optomechanical system with the nonlinear gain parameter @xmath82 @xcite . in sec .
v we will show how this parameter can affect the nms behaviour of the coupled bogoliubov mode and the cavity field .
the dynamics of the optomechanical system described by eqs.([hb ] ) is fully characterized by the following set of nonlinear qles , written in the frame rotating at the input laser frequency ,    [ nql ] @xmath83-\eta+\xi,\\ \dot{c}_{10}&=&-(i\omega_{10}+\gamma)c_{10}-i\frac{\sqrt{2n}}{4}u_{0}a^{\dagger}a-\frac{iu_{0}}{2}a^{\dagger}a c_{10}\nonumber\\ & & -\frac{i}{2}\omega_{sw}c^{\dagger}_{10}+f_{10},\\ \dot{c}_{01}&=&-(i\omega_{01}+\gamma)c_{01}-\frac{iu_{0}}{2}a^{\dagger}a c_{01}-\frac{i}{2}\omega_{sw}c^{\dagger}_{01}+f_{01},\nonumber\\\end{aligned}\ ] ]    where @xmath84 and @xmath85 . here
@xmath0 and @xmath11 characterize the dissipation of the cavity field and collective density excitations of the bec , respectively .
the cavity - field quantum vacuum fluctuation @xmath86 satisfies the markovian correlation functions , i.e. , @xmath87 , @xmath88 with the average thermal photon number @xmath89 which is nearly zero at optical frequencies @xcite . besides , @xmath90 and @xmath91 are the thermal noise inputs for the two side modes of bec which also satisfy the same markovian correlation functions as those of the optical noise .
the noise sources are assumed uncorrelated for the different modes of both the matter and light fields .
now we are going to study the weak excitations of the condensate from its ground state .
such excitations include small deviations of both the atomic wave function and the optical field from their respective stationary states .
so we decompose each operator in eqs.([nql ] ) as the sum of its steady - state value and a small fluctuation . by substituting @xmath92 , @xmath93 and @xmath94 into eqs.([nql ] ) one can obtain a set of nonlinear algebraic equations for the steady - state values ,    [ ss ] @xmath95    where we have assumed @xmath96 and @xmath97 are real numbers @xcite and @xmath98 is the effective detuning , @xmath99 and @xmath100 .
eq.([ss]c ) shows that the mean value of the mode @xmath69 is zero .
it is the consequence of the fact that the ground state of the translationally invariant system is also invariant under discrete translation @xcite . on the other hand , the linearized qles for the fluctuating operators take the followinf forms    [ flc ] @xmath101    where @xmath102 is the effective atom - photon coupling parameter . by defining the quadrature fluctuations    [ flc ] @xmath103    with @xmath104 ,
the qles can be written in the compact matrix form @xmath105 where @xmath106^{t}$ ] is the vector of continuous variable fluctuation operators and @xmath107^{t}.\ ] ] is the corresponding vector of noises .
the @xmath108 matrix @xmath109 is the drift matrix given by @xmath110 as is seen from the drift matrix , the side mode @xmath111 has been decoupled from the optical mode .
it is due to the fact that in the qles , eqs.([nql]c ) , there is no standard radiation pressure coupling between the side mode @xmath69 and the optical field .
instead , it is coupled to the radiation field via the nonlinear term , @xmath112 which is as small as @xmath113 of the radiation pressure term , and is deleted during the linearization process ( because the mean field value of this mode , @xmath114 , is zero ) . in this way the side mode @xmath69 acts as a medium level in the dynamics .      in order to study the stationary properties of the system
it is enough to focus our attention on the subspace spanned by the optical mode and the side mode @xmath68 .
it means that we can consider only the upper block of the drift matrix of eq.([a ] ) .
the system is stable only if the real part of all the eigenvalues of the matrix @xmath109 are negative .
these stability conditions can be obtained by using the routh - hurwitz criterion @xcite . due to the linearized dynamics of the fluctuations and since all noises are gaussian the steady state is a zero - mean gaussian state which is fully characterized by the @xmath115 stationary correlation matrix ( cm ) @xmath116 , with components @xmath117 . when the system is stable such a cm is given by@xcite @xmath118     versus the normalized cavity detuning @xmath119 for two values of @xmath120 ( thin line ) and @xmath121 ( thick line ) .
the parameters are @xmath122 , @xmath123 , @xmath124 , @xmath125 , @xmath126 and @xmath127.,title="fig:",width=288 ]   versus the normalized cavity detuning @xmath119 for two values of @xmath120 ( thin line ) and @xmath121 ( thick line ) .
the parameters are @xmath122 , @xmath123 , @xmath124 , @xmath125 , @xmath126 and @xmath127.,title="fig:",width=288 ]     and ( b ) the normalized effective atom - photon coupling @xmath128 ( thick line ) and the mean photon number magnified by @xmath129 ( thin line ) versus the normalized s - wave scattering frequency @xmath130 .
the parameters are the same as those of fig.[fig : fig2 ]
. the cavity detuning has been set at @xmath131.,title="fig:",width=288 ]   and ( b ) the normalized effective atom - photon coupling @xmath128 ( thick line ) and the mean photon number magnified by @xmath129 ( thin line ) versus the normalized s - wave scattering frequency @xmath130 .
the parameters are the same as those of fig.[fig : fig2 ]
. the cavity detuning has been set at @xmath131.,title="fig:",width=288 ]    where @xmath132 and @xmath133 is the matrix of stationary noise correlation functions . for the noise diffusion matrix we have @xmath134 , where @xmath135 $ ] is the diffusion matrix , with @xmath136^{-1}$ ] as the mean number of thermal excitations of the bogoliubov side mode @xmath137 whose frequency of oscillations is given by @xmath138 therefore , eq.([vij ] ) becomes @xmath139    when the stability conditions are satisfied @xmath140 and eq .
( v2 ) will be equivalent to the following lyapunov equation for the steady - state cm @xmath141 equation([lyap ] ) is linear in @xmath116 and can be straightforwardly solved .
however , the explicit form of @xmath142 is complicate and is not reported here .
in this section we first discuss our results based on the numerical solutions of eqs.([ss ] ) for the mean fields and then solve eq.([lyap ] ) to obtain fluctuations in the number of atoms and photons and their entanglement .
we will show how the atom - atom interaction affects the mean fields , cavity resonance and atom - photon entanglement .
we analyse our results based on the experimentally feasible parameters given in ref.@xcite .      here
we are going to study the effect of s - wave scattering on the behaviour of the mean - field values of the optical field and the bogoliubov side mode @xmath68 , given by eqs.([ss ] ) , quantitatively .
we consider @xmath143 atoms distributed in the optical cavity of length @xmath4=178 @xmath144 m with bare frequency @xmath145 corresponding to a wavelength of @xmath146=780 nm . the optical mode is coherently driven at rate @xmath5=80.06@xmath147 with the recoil frequency for rubidium atoms @xmath148 3.57 khz and the one - atom light shift @xmath58=0.96 @xmath147 @xcite .    in fig.[fig :
fig2 ] we have plotted the mean number of photons and the fraction of condensate atoms occupying the bogoliubov mode @xmath68 versus the normalized atomic detuning @xmath119 at a fixed temperature of @xmath149 k , and for two values of s - wave scattering frequencies @xmath82 . as is seen the increasing of the s - wave scattering interaction shifts the resonance frequency of the cavity to the lower values . besides , it causes the number of atoms in the bogoliubov mode to decrease i.e. , leads to the depletion of this mode .
this is in good accordance with the results obtained from numerical solution of the gpe ( fig.1 of ref.@xcite ) .    from eq.([ss ]
a ) one can obtain the resonance condition as @xmath150 such a resonance shift has been obtained in ref.@xcite without considering atomic collisions . as was mentioned in that reference the bec acts as a kerr medium that shifts the empty - cavity resonance .
however , due to the atomic collisions the bogoliubov side mode @xmath68 is depleted to other modes i.e. , the mean value @xmath97 decreases .
so a bec with atom - atom interaction has a different resonance in comparison to a non - interacting one . for clarity
, we have shown the effect of atom - atom interaction on the resonance of the cavity in fig.[fig : fig3](a ) .
for a non interacting bec , i.e. , @xmath151 , the resonance occurs at @xmath152=28966 @xmath147 while for very large values of @xmath82 it goes to @xmath153 . in fig.[fig : fig3](b ) the effective atom - photon coupling parameter @xmath154 ( thick line ) , given by eq.([g ] ) , and the mean photon number of the cavity ( thin line ) have been plotted versus the normalized s - wave scattering frequency .
this figure shows that both the effective atom - photon coupling parameter and the mean number of photons decrease with @xmath82 .
all of these reductions are the direct consequence of the resonance shift exerted by atomic collisions . in the next subsection
it will be shown that increasing the strength of atomic interaction causes the number of fluctuating photons to decrease .
so the decrease of the total cavity photon number is a sign of atomic collisions .
these results coincide with those obtained by numerical solution of gpe @xcite .
( b ) and the the entanglement between the bogoliubov mode @xmath68 and the optical field ( c ) versus normalized cavity detuning @xmath119 for two different values of @xmath120 ( thin line ) and @xmath121 ( thick line ) .
all parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ]   ( b ) and the the entanglement between the bogoliubov mode @xmath68 and the optical field ( c ) versus normalized cavity detuning @xmath119 for two different values of @xmath120 ( thin line ) and @xmath121 ( thick line ) .
all parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ]   ( b ) and the the entanglement between the bogoliubov mode @xmath68 and the optical field ( c ) versus normalized cavity detuning @xmath119 for two different values of @xmath120 ( thin line ) and @xmath121 ( thick line ) .
all parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ]    after calculating the mean - field values we can obtain the elements of the drift matrix @xmath109 and solve for the steady solutions of eq.([na ] ) . as explained before , by solving the lyapunov equation ( eq.([lyap ] ) ) we can obtain the correlation matrix @xmath116 which gives us the second - order correlations of the fluctuations . the correlation matrix corresponding to the upper block of @xmath109 in eq.([a ] )
can be written as @xmath155 where @xmath156 and @xmath157 represent the correlations of the photonic and atomic degrees of freedom respectively , and @xmath158 describes the cross correlations . in this way
we can calculate the incoherent excitation number of photons @xmath159 and the incoherent excitation number of atoms in the bogoliubov side mode @xmath160 on the other hand the bipartite entanglement between the atomic and photonic degrees of freedom can also be calculated by using the logarithmic negativity @xcite @xmath161,\ ] ] where @xmath162^{1/2}$ ] is the lowest symplectic eigenvalue of the partial transpose of the @xmath163 cm , @xmath116,with @xmath164 .     and
the s - wave scattering frequency has been considered to be @xmath120 .
the other parameters are the same as those of fig.[fig : fig2 ] . for clarity , the values of @xmath165 and @xmath166
have been multiplied by 100 and 1000 respectively.,width=288 ]    in figs.[fig : fig4 ] ( a ) and ( b ) we have plotted , respectively , the incoherent excitation number of photons and atoms versus the normalized cavity detuning @xmath119 for two values of @xmath82 . as is seen , fluctuations are small far from resonance but they increase rapidly near the resonance . on the other hand
increasing atom - atom interaction decreases the fluctuations in the number of photons and atoms . since the total number of photons is the mean value @xmath167 that obtained in previous subsection plus the incoherent excitation number of photons , we can conclude that the stronger atom - atom interaction , the weaker optical field of the cavity .
in fact the atom - atom interaction makes larger shifts in the cavity resonance frequency and consequently reduces the cavity field intensity .
hence the decrease of the cavity output provides a direct measure of the atom - atom interaction within the condensate .    in fig.[fig
: fig4 ] ( c ) the entanglement between the bogoliubov mode @xmath68 and the optical field has been plotted versus the normalized cavity detuning @xmath119 for two values of @xmath82 .
again , the maximum of entanglement occurs at resonance . besides , an increasing in the s - wave scattering frequency causes the entanglement of the atoms and photons to decrease .
as it was mentioned above , an increase in atom - atom interaction shifts the resonance of the cavity which reduces the number of photons of the cavity and also causes the depletion of the bogoliubov mode.this reduction in the mean values of photons and atoms leads to a reduction in the effective atom - photon coupling ( fig.[fig : fig3](b ) ) which causes the entanglement of atoms and photons to decrease .
the thermal effects on the incoherent excitation numbers of photons and atoms in the bogoliubov side mode as well as the entanglement between photons and atoms have been illustrated in fig.[fig : fig5 ] . as is seen from the figure , increasing
the temperature causes the fluctuations of the number of photons and atoms to increase and on the other hand reduces the atom - photon entanglement .
furthermore , the thermal effects for the range of temperatures below @xmath168 are negligible .
finally , we are going to obtain the power spectrum of the bogoliubov side mode @xmath68 and also drive its effective frequency and damping rate .
we will show that the coupling between the cavity field and the bogoliubov mode which behaves like a mechanical mirror , leads to the splitting of the normal mode into two modes ( nms ) .
a similar theoretical approach has been done in ref.@xcite for opto - mechanical bose - hubbard hamiltonian ( ombh ) by expanding the atomic wave field in terms of the wannier functions which is valid only for weak atom - field nonlinearity @xcite .
instead , here we use the atomic field expansion , eq.([opaf ] ) , in the momentum space . for this purpose
, we solve the linearized qles for the fluctuations in the displacement operator @xmath169 as @xmath170 where @xmath171 is the fourier transformation of the fluctuations in the total force acting on the bogoliubov mode and @xmath172 is its susceptibility which are given by , respectively ,    @xmath173    where @xmath174 is the effective frequency of the bogoliubov mode , @xmath175 and @xmath176 is its effective damping rate , @xmath177    , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath183 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ] , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath183 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ] , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath183 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ]    , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath184 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ] , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath184 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ] , ( b ) normalized effective bogoliubov mechanical damping @xmath178 and ( c ) normalized power spectrum of the displacement operator of the bogoliubov mode versus the normalized frequency @xmath179 for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) and for @xmath182 and @xmath184 .
the other parameters are the same as those of fig.[fig : fig2].,title="fig:",width=288 ]    to calculate the spectrum , we need the correlation functions of the optical noise sources in the frequency domain    @xmath185    where we have assumed @xmath186 . for the atomic noise sources in the frequency domain we have    @xmath187    where we have omitted the indices of @xmath188 and @xmath189 for simplicity . here , @xmath190^{-1}$ ] is the number of thermal excitations for the bogoliubov mode and @xmath191 is the mechanical bogoliubov frequency of oscillation given by eq.([wm ] ) . by using the relation @xmath192
we can calculate the power spectrum of the displacement operator @xmath193 as follows @xmath194.\end{aligned}\ ] ] the modification of the frequency of the bogoliubov excitations of the condensate due to the radiation pressure given by eq.([omegaeff ] ) is equivalent to the optical spring effect in cavity optomechanical systems with movable mirrors @xcite .    in figs .
[ fig : fig6](a ) and [ fig : fig7](a ) we have plotted the normalized effective bogoliubov mechanical frequency @xmath195 versus the normalized frequency @xmath179 for two values of cavity damping rates @xmath183 ( fig.[fig : fig6 ] ( a ) ) and @xmath184 ( fig.[fig : fig7 ] ( a ) ) , and for two different values of @xmath180 ( thin line ) and @xmath181 ( thick line ) . as is seen from both figures , a higher two - body interaction makes the condensate more robust and the bogoliubov frequency of the condensate does not significantly deviate from @xmath191 .
besides , a decrease in the damping rate of cavity causes the appearance of the peaks in the curves .    in figs.[fig : fig6](b ) and [ fig : fig7](b ) the normalized effective bogoliubov mechanical damping @xmath178 has been plotted versus the normalized frequency @xmath179 for two values of cavity damping rates @xmath183 ( fig.[fig : fig6 ] ( b ) ) and @xmath184 ( fig.[fig : fig7 ] ( b ) ) . in the case of fig.[fig : fig6](b ) where the damping rate of the cavity is nearly equal to the effective atom - photon coupling constant @xmath196 , the effective damping rate @xmath176 is much lower compared to the case of fig.[fig : fig7](b ) where @xmath197 .
it means that the stronger atom - photon coupling , the higher atom loss and hence the higher value of the effective damping .
an experimental observation of this phenomenon of light induced back - action heating and consequent loss of atoms has been reported in ref.@xcite .
it was also found that the atom loss rate is increased near the resonance . on the other hand ,
the effect of atomic collisions on the effective damping rate has been demonstrated in figs.[fig : fig6](b ) and [ fig : fig7](b ) which show that a stronger atom - atom interaction causes the effective damping rate to be increased near the resonance points .
therefore increasing the rate of atomic collisions can help us in cooling of the bogoliubov mode of the bec by the radiation pressure .
finally , we are going to investigate the influence of atomic collisions on the phenomenon of nms in the power spectrum of bogoliubov displacement operator .
the nms is associated with the mixing between the fluctuation of the cavity field around the steady state and the fluctuations of the condensate ( bogoliubov mode ) around the mean field .
the origin of the fluctuations of the cavity field is the beat of the pump photons with the photons scattered from the condensate atoms .
the phenomenon of nms is observable whenever the energy exchange between the two interacting modes takes place on a time scale faster than the decoherence of each mode .    in figs.[fig :
fig6](c ) and [ fig : fig7](c ) we have shown the normalized power spectrum of the displacement operator of the bogoliubov mode @xmath198 versus the normalized frequency @xmath179 for two values of cavity damping rate @xmath183 ( fig.[fig : fig6 ] ( c ) ) and @xmath184 ( fig.[fig : fig7 ] ( c ) ) , and for two values of @xmath180 ( thin line ) and @xmath181 ( thick line ) . in the case of fig.[fig : fig6](c ) where the damping rate of the cavity is nearly equal to the effective coupling constant of atoms and photons @xmath196 , there is no splitting in the power spectrum ( neither for @xmath199 nor for @xmath200 ) while in fig.[fig : fig7](c ) where @xmath197 the nms is appeared for @xmath201 so in order to observe nms firstly the system should be in the strong coupling regime where the atom - photon coupling is larger than the decay rates of photons an atoms .
when this condition is fulfilled , the nms can be observable with increasing atom - atom interaction .
normal mode splitting of a system of large number of atoms coupled to a ring cavity has been observed experimentally in the strong cooperative coupling regime @xcite .
the splitting of the normal mode has been observed by increasing the number of atoms which leads to the increase in the s - wave scattering frequency .
in conclusion , we have done a theoretical investigation on the optomechanical properties of a one - dimensional bose - einstein condensate inside a driven optical cavity considering the effects of atomic collisions . due to the dispersive atom - photon interaction the atoms develop a band structure in the optical lattice of the cavity . in the limit of weak photon - atom coupling
the lowest bands @xmath59 can be excited by fluctuations due to the atom - light interaction . on the other hand ,
the atomic collisions scatter the atoms to the higher band and also populate states with non zero quasimomentum .
it has been shown that there is a nonlinear optomechanical coupling between the nonzero quasimomentum states with the optical field proportional to @xmath113 which can be disregarded in the bogoliubov approximation where the number of atoms in the bec mode is very large .
therefore , there is no radiation pressure coupling between the atomic modes with nonzero quasimomentom and the cavity field in the bogoliubov approximation .    in this way we have obtained a simplified optomechanical model where the bogoliubov side mode @xmath68 is coupled to the optical field through the radiation pressure term and the atom - atom interaction hamiltonian behaves like an opa for the bogoliubov mode with the non - linear gain parameter @xmath82 .
it has been shown that in the strong coupling regime where the effective atom - photon coupling is greater than the cavity damping rate , an increase in the s - wave scattering frequency would lead to the nms in the power spectrum of the bogoliubov mode .
furthermore , the atom - atom interaction causes the depletion of the bogoliubov mode and also shifts the cavity resonance frequency which leads to a decrease in the mean number of cavity photons .
besides , it decreases fluctuations in the number of photons and atoms .
hence the decrease in the cavity output provides a direct measure of the atom - atom interaction within the condensate . on the other hand ,
a stronger atomic collision rate causes a decrease in the effective atom - photon coupling which leads to a decrease in the entanglement between the bogoliubv mode and the optical field .
the authors wish to thank the office of graduate studies of the university of isfahan for their support .        c. genes , d. vitali , p. tombesi , s. gigan , and m. aspelmeyer , phys .
a * 77 * , 033804 ( 2008 ) .
j. d. teufel , j. w. harlow , c. a. regal , and k. w. lehnert , phys .
* 101 * , 197203 ( 2008 ) ; j. d. teufel _ et al _ , nature ( london ) * 475 * , 359 ( 2011 ) ; j. d. teufel _ et al _ , nature ( london ) , * 471 * , 204 ( 2011 ) .
s. de liberato , n. lambert , and f. nori , phys . rev .
a * 83 * , 033809 ( 2011 ) .
m. ludwig , c. neuenhahn , c. metzger , a. ortlieb , i. favero , k. karrai , and f. marquardt , phys .
lett . * 101 * , 133903 ( 2008 ) .
barzanjeh , m. h. naderi , and m. soltanolkotabi phys .
a * 84 * , 023803 ( 2011 ) . c. bradaschia _ et al _ , nucl .
methods phys .
res . a * 289 * , 518 ( 1990 ) ; a. abramovici _ et al _ , science * 256 * , 325 ( 1992 ) ; p. fritschel , proc .
spie * 4856 * , 282 ( 2003 ) .
m. d. lahaye , o. buu , b. camarota , and k. c. schwab , science * 304 * , 74 ( 2004 ) . t. j. kippenberg and k. j. vahala , opt .
express * 15 * , 17172 ( 2007 ) .
f. brennecke , t. donner , s. ritter , t. bourdel , m. kohl , t. esslinger , nature ( london ) * 450*,268 ( 2008 ) .
s. gupta , k.l .
moore , k.w .
murch , and d.m .
stamper - kurn , phys .
* 99 * , 213601 ( 2007 ) .
f. brennecke , s. ritter , t. donner , and t. esslinger , science * 322 * , 235 ( 2008 ) .
s. ritter , f. brennecke , k. baumann , t. donner , c. guerlin , and t. esslinger , appl .
b * 95 * , 213 ( 2009 ) .
o. morsch and m. oberthaler , rev .
phys . * 78 * , 179 ( 2006 ) . c. maschler and h. ritsch , opt
. commun . * 243 * , 145 ( 2004 ) p. domokos and h. ritsch , j. opt .
b * 20 * , 1098 ( 2003 ) . c. maschler and h. ritsch , phys .
rev . lett.*95 * 260401 ( 2005 ) . c. kittel , _ quantum theory of solids _
( wiley , new york , 1963 ) .
a. bhattacherjee , opt . commun . * 281 * , 3004 ( 2008 )
. g. szirmai , d. nagy , and p. domokos , phys .
a * 81 * , 043639 ( 2010 ) .
s. k. steinke and p. meystre , phys .
a * 84 * , 023834 ( 2011 ) .
w. chen , d. s. goldbaum , m. bhattacharya , and p. meystre , phys .
a * 81 * , 053833 ( 2010 ) .
g. konya , g. szirmai , and p. domokos , eur . phys .
j. d * 65 * , 33 ( 2011 ) .
r. kanamoto and p. meystre , phys .
* 82 * 038111 ( 2010 ) . c. j. pethick and h. smith ,
_ bose einstein condensation in dilute gases _
( cambridge university press , 2002 ) .
d. nagy , p. domokos , a. vukics , and h. ritsch , eur .
j. d * 55 * , 659 ( 2009 ) . c. maschler ,
mekhov , and h. ritsch , eur .
j. d * 46 * , 545 ( 2008 ) .
asboth , p. domokos , and h. ritsch , phys .
a * 70 * , 013414 ( 2004 ) .
d. s. goldbaum , k. zhang , and p. meystre , arxiv:0911.3234v1 . i. s. gradshteyn and i. m. ryzhik , _ table of integrals , series and products _
( academic press , orlando , 1980 ) ; a. hurwitz , _ selected papers on mathematical trends in control theory _
, edited by r. bellman and r. kabala ( dover , new york , 1964 ) . c. w. gardiner and p. zoller , _ quantum noise _
( springer , berlin , 2000 ) . c. genes , d. vitali , and p. tombesi , new j. phys . * 10 * , 095009(2008 ) .
j. m. zhang , f. c. cui , d. l. zhou , and w. m. liu , phys .
rev . a * 79 * , 033401 ( 2009 ) .
p. horak , s. m. barnett , and h. ritsch , phys .
a * 61 * , 033609 ( 2000 ) .
j. eisert , ph.d .
thesis , university of potsdam , 2001 ; g. vidal and r. f. werner , phys .
a * 65 * , 032314 ( 2002 ) ; g. adesso , a. serafini , and f. illuminati , phys .
a * 70 * , 022318 ( 2004 ) .
a. bhattacherjee , j. phys .
b : at . mol .
. phys . * 43 * 205301 ( 2010 ) .
j. larson _
et al _ , phys .
lett . * 100 * 050401 ( 2008 ) .
k. w. murch , k. l. moore , s. gupta , and d. m. stamper - kurn , nature phys . * 4 * , 561 ( 2008 ) .
j. klinner , m .
lindholdt , b .
nagorny , and a. hemmerich , phys .
* 96 * , 023002 ( 2006 ) . c. w. gardiner , phys .
a * 56 * , 1414 ( 1997 ) .


until a few years ago there was a common feeling in the condensed matter physics community that we had a very good understanding of the universal properties of one - dimensional systems of fermions through our knowledge of luttinger liquid theory and its instabilities @xcite .
one of the key results of luttinger liquid theory @xcite is that interacting electrons in one - dimension have separate dynamics in the spin and charge sectors at low energies .
these low energy collective spin and charge modes are bosonic and generally propagate with different velocities , with the spin mode typically being the slower of the two .
luttinger liquid theory provides a particularly convenient theoretical framework for discussing the response of the electron gas to a time - dependent local potential , even in higher dimensions : the mathematical trick of  bosonization " makes the evaluation of correlation functions relatively simple , and a potential local on the scale of the fermi wavelength results in scattering of primarily @xmath3-wave type which effectively converts the problem to just one dimension ( the radial coordinate ) .
a few years back it was realized that at finite temperatures a  luttinger liquid system " can take on qualitatively new behaviors in transport @xcite and the single - particle green s function @xcite .
some of the results for the single - particle green s function were anticipated in an earlier work @xcite , but the many novel features that can arise due to finite temperatures were not fully appreciated until the more recent works mentioned above which really started proper the study of the spin - incoherent luttinger liquid @xcite .    the spin - incoherent regime of a one - dimensional interacting electron gas is defined in the following way .
let us assume that at zero temperature our system is spin - charge separated with spin velocity @xmath4 and charge velocity @xmath5 such that @xmath6 .
if the interactions are strong ( regardless of the range of interactions ) it can be shown that @xmath4 can be _ exponentially small _ compared to @xmath5 @xcite . defining a characteristic spin energy @xmath7 and a characteristic charge energy @xmath8 for mean particle spacing @xmath9 ,
it is evidently clear that the spin energy can be exponentially suppressed relative to the charge energy : @xmath10 .
an exponentially separated spin and charge energy makes it possible for the thermal energy to lie in the range @xmath11 where @xmath12 is boltzman s constant and @xmath13 is the temperature .
a one dimensional electron system in the energy window @xmath11 is said to be in the spin - incoherent regime .
if the charge sector is gapless it will be described by a zero - temperature luttinger liquid theory and the resulting system is called a spin - incoherent luttinger liquid ( sill ) .
note that if @xmath14 refers to a spin _ gap _ then a sill is still realized when @xmath11 provided the charge sector remains gapless .
for a general discussion of the sill see @xcite . for sill physics in ferromagnetic and
superconducting hybrid structures see @xcite , and in cold atomic gases see @xcite .
in this article we will only discuss a particularly simple form of a time - dependent spatially local potential whereby its value is changed instantaneously from zero to non - zero .
the prototypical physical realization of this situation is when a photon in the x - ray energy range ejects a deep  core " electron from a metal  instantly changing the local charge by @xmath15 .
the electrons of the metal then rush in to screen the positive potential and thereby change their states .
if the electrons interact strongly with each other the processes of the electrons  rushing in " to screen the local potential will obviously be impacted . according to luttinger liquid theory , interactions qualitatively change the nature of the low energy physics in such a way that
electron - like quasi - particles are absent from the theory @xcite leaving only the collective bosonic excitations . on this basis alone , we expect the response of one - dimensional systems to differ in some ways from their higher dimensional counter parts . in fact , this expectation is correct but to fully appreciate the differences , it is important to review some of basic ideas and qualitative expectations for the singular responses of non - interacting electron systems .
the seminal theoretical ideas of the singular response of a non - interacting electron gas due to a localized core hole were worked out by @xcite and @xcite .
a nice textbook - level discussion can be found in @xcite and a more recent review in @xcite .
the most important results to emerge are : ( i ) there is a threshold frequency @xmath1 above which there is photon absorption and below which there is not .
( ii ) the photon absorption at frequencies just above @xmath1 follows a power law with exponent @xmath2 .
the photon absorption intensity is computed from a fermi s golden rule - type treatment in the electron - photon coupling and from the results just cited above takes the form @xmath16 near threshold where @xmath17 is a constant , and @xmath18 for @xmath19 while @xmath20 for @xmath21 .
a schematic of the possible behaviors is shown in fig .
[ fig : schematic_threshold ] .     schematic behavior of photon absorption intensity , eq.([eq : i_fl ] ) , for an infinite mass core hole in three dimensions .
the exponent @xmath2 given by eq.([eq : a ] ) contains positive contributions from excitonic effects and negative contributions from anderson s orthogonality catastrophe .
it is possible for @xmath2 to be either positive ( a ) or negative ( b ) . in the spin - incoherent regime of one dimensional systems
the absorption is not generally given by eq.([eq : i_fl ] ) but instead exhibits logarithmic corrections and additional physics enters the exponent @xmath2 . ]    the important microscopic details of the electron dynamics and the electron - core hole interaction are encoded in the exponent @xmath2 and the power - law form itself .
the power - law arises from summing an infinite series of diagrams containing logarithmic divergences @xcite .
as we will see , in the sill the power law form is not sacred .
it will be modified by logarithmic corrections , and the physics encoded in @xmath2 will depend on whether or not the one - dimensional system is in the spin - incoherent regime . for a non - interacting electron gas
the exponent takes the form @xcite
@xmath22 ^ 2,\ ] ] where @xmath23 is the phase shift induced by the core hole in the @xmath24 channel of the conduction electrons .
often @xmath3-wave scattering ( @xmath25 ) is the dominant channel and we will assume this is the case throughout this review . as the total amount of screening charge around the core hole must be one unit of charge , there is a friedel sum rule : @xmath26 if the scattering is not too strong , the phase shifts @xmath23 can be computed in the born approximation and one finds @xcite @xmath27 where @xmath28 is the density of states at the fermi energy and @xmath29 is the strength of the electron - hole interaction assumed momentum and energy independent .
( the actual core hole - electron interaction is @xmath30 and therefore attractive . ) within the born approximation and assuming only @xmath3-wave scattering , the exponent becomes @xmath31 ^ 2.\ ] ] evidently , eq.([eq : a ] ) and eq.([eq : a_approx ] ) imply there are competing effects in @xmath2 .
the first term proportional to @xmath32 is positive while the second term proportional to @xmath33 ^ 2 $ ] is negative .
going back to eq.([eq : i_fl ] ) one sees that a positive contribution to @xmath2 tends to make the threshold response more singular , while a negative contribution to @xmath2 tends to make the threshold response less singular .
in fact , the origin of these two contributions are very well understood .
the first contribution to @xmath2 in eq.([eq : a ] ) comes from excitonic effects arising from the attractive interaction between the electrons and the core hole .
the second contribution to @xmath2 in eq.([eq : a ] ) arises from anderson s orthogonality catastrophe @xcite and it tends to suppress the absorption by the small ( vanishing as @xmath34 ) overlap between initial and final states .
note that for a non - interacting fermi sea both excitonic effects and orthogonality effects disappear if @xmath35 and the absorption threshold becomes a simple step .
looking ahead for a moment to luttinger liquid physics , it is instructive to point out that if a collection of fermions is interacting , there will be orthogonality effects if an additional fermion of like charge is added to the system because the other fermions must adjust their state to accommodate the newcomer , but no excitonic effects .
this rearrangement of electronic states is what is responsible for the famous power law suppression of the tunneling density of states in a luttinger liquid @xcite .
finally , it is worth emphasizing that depending on the microscopic details involved , it is possible for @xmath2 to be positive or negative as shown in fig .
[ fig : schematic_threshold ] . if the mass of the core hole is finite and one considers only direct transitions like those in fig .
[ fig : direct_transition ] one finds that the threshold is smeared out and disappears @xcite .
we will find this is not the case when we go to one dimension .
schematic of the possible direct transitions for exciting a hole with ( a ) finite mass and ( b ) infinite mass . here
@xmath36 denotes the fermi energy . in ( a ) an electron
is excited from a valence band up to a conduction band leaving behind a single mobile hole in the valence band . in
( b ) an electron is excited from an infinite mass core level .
the photo - absorption of both process ( a ) and ( b ) differ in one dimension compared to their counterparts in three dimensions .
furthermore , the response in one dimension depends on whether the system is in the luttinger liquid or spin - incoherent luttinger liquid regime . ]
before we turn our attention to the spin - incoherent luttinger liquid , it is useful for orientation to discuss the singular response of its closest cousin , the luttinger liquid .
the behavior exhibited by the one - dimensional luttinger liquid ( ll ) will already take us part of the way to understanding the response of the sill . in all of the results quoted below finite system size effects
will not become important until frequencies smaller than @xmath37 are probed , where @xmath38 is the length of the system and @xmath39 is the fermi velocity .
also , finite temperature will cut - off the singular behaviors when relevant frequencies are of order @xmath40 and smaller . in this sense , finite system size and temperature
can act as a non - zero threshold for tunneling into a finite length system or an infinite system at finite temperature , for example .
let us first consider the simplest kind of singular response : abruptly adding a particle to a ll as occurs during tunneling . a nice pedagogical discussion of the essential physics of tunneling into a ll is given by @xcite .
here we adopt the central elements of their discussion .
the tunneling process is conveniently visualized within a wigner crystal - type picture for spinless electrons , as shown in fig .
[ fig : wigner_schematic ] .
schematic of the tunneling process of spinless electrons into a  wigner crystal " type electron arrangement .
electrons are schematically represented by black dots evenly spaced on either side of the tunneling barrier of strength @xmath41 at @xmath42 .
the tunneling occurs in two steps : a fast process where the electron enters the region @xmath21 from the left followed by a slow process of charge relaxation where each electron shifts over to the right by a distance of one lattice constant . at long times
@xmath43 as discussed in the text .
the resulting behavior of the tunneling density of states is @xmath44 the qualitative behavior of the tunneling density of states looks very much like that in fig .
[ fig : schematic_threshold](b ) with @xmath45 . ]
such a picture implies very strong electron interactions , but the physical result and the final mathematical expression that describes it are actually valid at any interaction strength .
we imagine a process where an electron tunnels from the left of the barrier at @xmath42 to the right of the barrier at time @xmath46 . the rate at which the electron tunnels is determined by the debye frequency of the wigner crystal , @xmath47 , where @xmath48 is the fermi wave vector and @xmath49 is the phonon ( charge ) velocity .
when an electron tunnels the potential near the barrier suddenly ( on a time scale @xmath50 ) changes .
the energy cost of this new configuration is large , @xmath51 , where @xmath52 is the interaction parameter for the spinless luttinger liquid .
( here @xmath53 if the electrons are non - interacting . )
this large energy cost is reduced by the relatively slow process of the electrons relaxing to new equilibrium positions in order to minimize the energy . once the relaxation process is complete , all the electrons
will be shifted to the right by one crystalline period . in the thermodynamic limit
this new wavefunction is orthogonal to the original one which implies at zero energy ( @xmath54 ) there is no tunneling .
this is a concrete example of anderson s orthogonality catastrophe that we earlier mentioned plays a role in the exponent @xmath2 of the threshold behavior of photon absorption . for finite energy
the final state may be in an excited state so the relaxation is not complete and the initial and final wavefunction can have finite overlap .
we will now see how this directly leads to a power law tunneling density of states , one of the  singular " responses of interest in this article .
we can compute the tunneling density of states from the single particle green s function @xcite . in the semi - classical approximation
the electron tunneling takes the path of minimal action which is given by the classical trajectory @xmath55,\ ] ] where @xmath56 is the electron annihilation operator at position @xmath57 and time @xmath58 , and @xmath59 is the creation operator .
the initial density deformation @xmath60 will eventually relax to @xmath61 .
as the processes is classically forbidden , the trajectory runs in imaginary time @xmath62 .
the classical equation of motions for a harmonic lattice give @xmath63 the potential energy of the deformation is @xmath64 ^ 2},\ ] ] which gives @xmath65 . by the virial theorem ,
the kinetic energy is equal to the potential energy . therefore the action describing the slow part of the tunneling is @xmath66 whence @xmath67 . taking the fourier
transform one finds the tunneling density of states @xmath68 . since we have assumed @xmath69 , the tunneling density of states @xmath70 as @xmath34 , as we noted above in the context of the orthogonality catastrophe .
the discussion above illustrates some of the most central aspects of singular responses of interacting one - dimensional systems .
however , there is much more to the story .
let us now suppose that instead of an electron tunneling , we consider an incident photon that abruptly creates a particle - hole pair as shown in fig .
[ fig : direct_transition ] . note that due to the presence of the gap between the hole states and conduction band , there will be a minimum frequency @xmath1 required to create the particle - hole pair .
if we first focus on the electrons in the conduction band , we see that the process is similar to tunneling : conduction band electrons see the sudden appearance of a new electron and they must adjust their positions by relaxing to a new state just as if the electron tunneled from outside the system . as we have seen , this relaxation rate is related to the strength of the electron interactions @xmath71 in the luttinger liquid , and based on the power law density of states we obtained , we expect the threshold behavior will also contain at least a power law contribution with an exponent that contains information about the interactions in the system .
it is also important to consider the interaction between the hole created and the electrons in the conduction band . as holes and electrons
generally have different masses and they attract rather than repel each other , the creation of a hole adds new physics to the electron tunneling physics above .
in fact , to compute the photon absorption spectrum near threshold we must _ simultaneously _ account for the interaction of the electrons with each other and with the hole . in general
, this is a very difficult problem .
fortunately , though , the power of luttinger liquid theory and the commonly used technical methods such as bosonization allow us to make some precise and quite general statements about @xmath72 .
pioneering theoretical work in the study of the fermi - edge singularity in luttinger liquids was carried out more than 15 years ago @xcite , and later works soon followed @xcite .
a recent review of the optical response of low dimensional systems ( including some aspects of one - dimension ) exists @xcite , so here we will only highlight the main issues and results to set the stage for our discussion of the spin - incoherent one - dimensional systems .    in one - dimensional electronic systems
the effect of impurities is especially dramatic compared to three dimensions .
even for an unrealistic non - interacting one - dimensional system the electrons localize leading to insulating behavior . however , for a realistic system which contains repulsive interactions even a single local potential of arbitrary strength will localize the electrons to the half - line @xcite .
the central physics of this result is the relevance of electron backscattering as a perturbation : at low energies the strength of the impurity grows to arbitrarily large values and eventually divides the system into two semi - infinite segments .
this physics is important for the fermi - edge singularity involving an infinitely massive core hole , _ i.e. _ one that is localized in space .    on the other hand , when the mass of the core hole is finite the physics is very different .
the properties of a heavy particle ( holes are typically heavier than electrons ) moving in a luttinger liquid have been investigated @xcite .
it was shown that the mobility of the heavy particle diverges at low energies , implying that electron backscattering from it becomes irrelevant .
we therefore see that the fermi - edge singularity in a luttinger liquid falls into two categories : ( i ) infinite hole mass and ( ii ) finite hole mass .
note that in three dimensions it was also important to divide the physics into the same two categories , but the reason for that division was different . in three dimensions finite
hole mass led to excitonic effects that smeared the edge singularity of the direct transition and eliminated the threshold .
this has to do with a large ( compared to one - dimension ) phase space for recoil effects .
a sharp threshold only exists for the infinitely massive hole . on the other hand , in one dimension
the issue is about the relevance / irrelevance of backscattering and that is decided by the interactions ( repulsive or attractive ) in the ll .
another way to see that it is natural to divide the classes into infinite and finite masses cases is that the symmetry is different .
an infinitely massive hole breaks the translational symmetry of the problem , while a hole with finite mass does not .
one of the first issues addressed regarding the fermi - edge singularity in the ll is whether the threshold still exists for the case of finite hole mass .
this was answered in the affirmative @xcite .
the next logical question is how does the threshold exponent @xmath2 depend on the electron - hole and electron - electron interactions ?
initially it was argued that @xmath2 was independent of the hole mass @xcite , but this was later challenged and corrected @xcite . knowing that any finite mass hole is an irrelevant perturbation @xcite we
can immediately see that only forward scattering interactions will contribute to @xmath2 .
indeed this is found as well as a schematic phase diagram for when @xmath2 is positive or negative @xcite .
the result is @xcite @xmath73,\ ] ] where @xmath74 is an asymmetric phase shift and @xmath75 a symmetric phase shift defined below in the discussion of the spin - incoherent case .
the presence of @xmath74 implies that the exponent depends on the mass of the hole .
note that both @xmath74 and @xmath75 correspond to forward scattering processes , and @xmath69 is the luttinger parameter for repulsively interacting spinless electrons .
we can directly compare the result eq.([eq : a_ll_mobile ] ) with the non - interacting case eq.([eq : a ] ) by taking @xmath53 and @xmath76 ( appropriate for an infinitely massive hole ) .
doing so , we note the form is identical to an s - wave only scattering model . from this comparison
we also see that @xmath75 has the meaning of a phase shift up to factors of 2 and @xmath77 .
note that we are allowed to compare the result for a mobile impurity in a ll with an infinitely massive impurity in a free electron gas in three dimensions because in the latter backscattering is not a relevant perturbation .
they are therefore in a similar  class " .
the situation is different when we turn to immobile infinite mass holes . here , as before , a threshold exists , but now @xmath2 is modified by the relevant backscattering and in fact obtains a universal contribution of 1/8 @xcite , @xmath78 in the next section we will discuss how these two results are modified in the spin - incoherent regime .
but before we do that , it is important to comment on the experimental situation .
unfortunately , there are few published results @xcite and the overall level of agreement between theory and experiment is poor . clearly there is great need for further experiments on fermi - edge singularities on clean , single mode one - dimensional systems .
we hope this article will help to stimulate more .
we finally come to the main topic of this review well prepared from the earlier discussions of singular responses of non - interacting electron gases and luttinger liquids .
the question to ask is what are the similarities and differences between the singular responses of a sill and a ll or fermi liquid ? as we will see , the behavior is qualitatively different from both the ll and fermi liquid making it possible to distinguish them in experiment .
perhaps the simplest and most striking difference is the behavior of the tunneling density of states .
we begin there .
the first theoretical study of tunneling into a sill was reported in the infinite @xmath79 limit of the hubbard model in a lengthy bethe ansatz calculation @xcite . immediately afterwards a simpler and more general method based on bosonization and the imaginary time path integral representation
was given @xcite .
extensions of the method to finite length systems including boundary effects and externally applied magnetic fields have also been reported @xcite .
while the imaginary time path integral formulation is well suited to the regime of the sill , it does not provide a conceptually straightforward generalization to the regime @xmath80 . in order to address that regime , we need a representation of the spin degrees of freedom that are valid at arbitrary energy scales .
such a formulation is now in hand @xcite and it provides a good starting point for numerical studies in the regime @xmath80 .      as is typical of interacting one dimensional electron systems
, we assume the hamiltonian in the strongly interacting case is spin - charge separated @xcite , @xmath81 . here , @xmath82 is identical to that of the ll , and is given by @xcite @xmath83 , \label{eq : h_holon}\ ] ] where @xmath84 and @xmath85 are bosonic fields satisfying @xmath86=-i\frac{\pi}{2 } { \rm sgn}(x - x')$ ] , and @xmath87 is the effective luttinger interaction parameter of the electron gas @xcite .
on the other hand , the spin hamiltonian at arbitrary temperatures is to a very good approximation given by a nearest neighbor antiferromagnetic heisenberg spin chain @xcite , @xmath88 where evidently the spin energy is set by @xmath89 : @xmath90 .
the basic idea is to represent the electron operator as a product of operators that describe the holons @xmath91 [ spinless fermions that naturally arise in the context of strongly interacting fermions and the spin - incoherent regime @xcite ] and the spin degrees of freedom @xmath92 .
the holon operators ( denoted by @xmath93 ) by construction satisfy the equation @xmath94 where @xmath95 is the electron annihilation operator for electrons of spin projection @xmath96 , and @xmath97 the corresponding electron creation operator .
the issue of how to bosonize the electron operator for a strongly interacting system earlier arose in the context of the large @xmath79 limit of the one - dimensional hubbard model where @xcite wrote the electron creation operator as @xmath98 where @xmath99 creates a site on the spin chain eq.([eq : s_chain ] ) with spin projection @xmath96 .
the expression in eq.([eq : penc_bos ] ) can be physically motivated as follows . from eq.([eq : psi ] ) it is clear that the creation of an electron is also accompanied by the creation of a holon .
however , electrons also carry spin so there must be a component of the electron operator that also creates spin .
this is accomplished by @xmath99 . in general
, one has @xmath100 as the object that adds a new site to the spin chain between @xmath101 and @xmath102 .
while this appears physically intuitive , the expression suffers from the drawback that it does not naturally account for the variation of electron density with position in a real electron gas @xcite .
the remedy for this issue @xcite is to define the position at which the spin site is added to the chain eq.([eq : s_chain ] ) in terms of the number of holons to the left of the site , @xmath103 in terms of eq.([eq : l_x ] ) the electron creation and annihilation operators are defined as @xmath104 the operators given above explicitly account for the fact that the spins are attached to electrons , and the formulas are valid at all energy scales .
it is perhaps worth noting that even though the hamiltonian is spin - charge separated , the electron operators are not written as a product of a spin piece and a charge piece because the  spin " pieces @xmath105 also depend on the electron density via eq.([eq : l_x ] ) .    in writing eq.([eq : psi+_x ] ) and eq.([eq : psi_x ] ) no assumptions have been made about the energy scale relative to the spin and charge energies .
we now restrict our considerations to energies small compared to @xmath106 , but arbitrary with respect to @xmath14 . in this case
, we are free to bosonize the holon sector . we have already given the form of the charge hamiltonian in eq.([eq : h_holon ] ) .
the spinless fields @xmath84 and @xmath85 can be related to the holon density as @xcite @xmath107,\ ] ] where the holon fermi wave vector is twice the electron fermi wavevector @xcite @xmath108 and the bosonic fields satisfy the commutation relations @xmath109=i\pi\delta(x - y)$ ] .
since we are interested in low energies with respect to the charge energy , the electron operator may be expanded about the two holon fermi points at @xmath110 , @xmath111 where @xmath112 destroys an holon near the right fermi point and @xmath113 destroys an electron near the left fermi point .
the left and right holon operaters are bosonized as @xmath114},\ ] ] where @xmath115 is a short distance cut - off of order the interparticle spacing @xmath9 . combining the results of eq.([eq : l_x ] ) , eq.([eq : psi_x ] ) , eq.([eq : den_rel ] ) , eq.([eq : psi_rl ] ) and eq.([eq : psi_rl_fields ] ) one obtains the bosonized form of the electron annihilation operator for spin @xmath96 @xmath116}+e^{-i[k_f^hx+\theta(x)]}\right ) z_{l,\sigma}\bigg |_{l=\frac{1}{\pi}[k_f^hx+\theta(x)]},\end{aligned}\ ] ] and an analogous expression for the electron creation operator @xmath117 .
expression eq.([eq : psi_bos ] ) , however , it not quite complete as it does not account for the discreteness of the charge of the electron .
this can be accomplished by interpreting @xmath118 } \to \sum_l z_{l,\sigma}\delta\left({1 \over \pi}[k_f^hx+\theta(x)]-l\right),\ ] ] after which the full electron annihilation operator ( include both left and right moving parts ) becomes @xcite @xmath119},\ ] ] where @xmath120 the expression for the electron annihilation operator eq.([eq : psi_final_bos ] ) is the key result of @xcite who also showed that in the limit of small energies compared to @xmath14 the expression correctly recovers the standard ll formulas for the electron annihilation operator . with eq.([eq : psi_final_bos ] ) correlation functions involving electron operators can be expressed in terms of the correlation functions of the holon and spin sectors at arbitrary temperatures with respect to @xmath14 , but small energies compared to @xmath106 .
the tunneling density of states in a sill is qualitatively different from the ll . while the tunneling density of states in a ll qualitatively resembles the behavior shown in fig.[fig : schematic_threshold](b ) with @xmath45 , the tunneling density of a sill qualitatively resembles the behavior in fig.[fig : schematic_threshold](a ) with @xmath45 , at least for a range of frequencies , @xmath121 and for not too strong interactions . for this reason
, the tunneling density of states may provide the simplest and most transparent test for sill physics in experiment .
tunneling into an infinitely long sill has been considered by @xcite .
deep in the spin - incoherent regime , the first - quantized path integral representation of the bosonized electron s green s function is a convenient method @xcite .
here we instead follow @xcite as it provides a better launching point for numerical studies in the regime @xmath80 , and also readily captures some other important features of the tunneling density of states .
the tunneling density of states , @xmath122 $ ] , where @xmath123 is the fourier transform of the retarded green s function .
the retarded green s function generally contains particle and hole contributions , but for illustrative purpose , we will only focus on the particle contributions .
hole contributions are conceptually and technical similar , and in the spin - incoherent regime it can be shown that @xmath124 for @xmath125 @xcite . here
we focus on the particle contribution : @xmath126 . using the general formulas
eq.([eq : psi_final_bos ] ) and eq.([eq : little_z ] ) we have @xmath127-\phi(x , t ) ] }   e^{-i[(1+\frac{q_2}{\pi})\theta(0,0)-\phi(0,0)]}\rangle,\;\;\;\end{aligned}\ ] ] where we have assumed the time evolution of the spin degrees of freedom are slow compared to the charge degrees of freedom so that we can neglect their time dependence @xcite . in a translationally invariant system , one has @xmath128 . making a change of variables @xmath129 , the summation over the sites of the spin chain results in a delta function , @xmath130 , which immediately kills one of the momentum integrals and sets @xmath131 .
the resulting green s function is @xmath132 where @xmath133 and @xmath134 }   e^{-i[(1+\frac{q}{\pi})\theta(0,0)-\phi(0,0)]}\rangle.\ ] ] since the holon sector is described the the gaussian theory eq.([eq : h_holon ] ) we can make use of the identity @xmath135 for operator @xmath136 and evaluate the correlation function eq.([eq : g ] ) at zero temperature with respect to the charge energy .
[ generalization to small but finite temperatures in the charge sector is trivial @xcite . ]
the expression for the hole - like contributions to the green s function has the same structure as eq.([eq : g_final_trans ] ) and is reported in @xcite . since for low temperatures compared to the charge energy the holon correlations in eq.([eq : g ] ) can be analytically evaluated in a straightforward way , the emphasis for arbitrary temperatures with the respect to the spin energy is on the spin correlations described by eq.([eq : c ] ) . for arbitrary temperatures
the evaluation of @xmath137 is non - trivial @xcite and it is precisely at this point that numerical input for both zero and non - zero magnetic fields would be most helpful . keeping for the moment temperature arbitrary with respect to @xmath14 , but effectively zero compared to @xmath106 one finds for the tunneling density of states @xcite , @xmath138 where @xmath139 - 1 $ ] and @xmath87 is the luttinger interaction parameter appearing in eq.([eq : h_holon ] ) . from the exponent @xmath140 one sees that the frequency dependence @xmath141 which is maximum for @xmath142 .
expanding about this point and doing the guassian integration , one finds @xmath143 which agrees with the result earlier obtained for the infinite @xmath79 limit of the hubbard model @xcite and for general interaction parameter @xmath71 in the spin - incoherent regime @xcite .
note , however , that the result eq.([eq : nu_structure ] ) is valid at @xmath144 .
the only assumptions were that the dominant weight comes from the point @xmath142 and that @xmath145 .
if one is deep in the spin - incoherent regime @xmath146 _ independent of the details and form _ of @xmath147 [ a simple argument for this is given in @xcite ) ] , and one has @xmath148 . hence , the assumptions leading to eq.([eq : nu_structure ] ) are valid and the result indeed applies in the spin - incoherent regime as well . on the other hand ,
if the spins are fully polarized ( along the z - axis , say ) then @xmath149 and the tunneling density of states will necessarily have a different form ( in fact , the familiar ll form ) from eq.([eq : nu_structure ] ) @xcite .
also at @xmath144 and zero magnetic field one has @xmath150 so the logarithmic factor appearing in eq.([eq : nu_structure ] ) is absent and the ll form is obtained .    finally , note that in the spin - incoherent regime where eq.([eq : nu_structure ] ) applies the tunneling response of a sill looks qualitatively like the behavior in fig .
[ fig : schematic_threshold](a ) with @xmath45 for @xmath151 , but it _ does not _ follow the analytical form of eq.([eq : i_fl ] ) because of the logarithmic corrections . if the interactions are both strong and long - ranged , it is possible to obtain @xmath152 in which case the tunneling density of states would look qualitatively like the behavior in fig .  [ fig : schematic_threshold](b ) with @xmath45 although it would still not fit the analytic form in eq.([eq : i_fl ] ) .
[ recall that for zero - range interactions the smallest value @xmath153 can take is @xmath154 , obtained in the infinite-@xmath79 limit of the hubbard model @xcite . ] in this sense , the tunneling response of the sill is in a different universality class from both the fermi liquid and luttinger liquid .
we will now see that the same logarithmic corrections can also appear in the frequency dependence of the photon absorption .
as we discussed earlier , the physics of the fermi - edge singularity falls into two classes : ( 1 ) electron backscattering from the hole is an irrelevant perturbation .
( 2 ) electron backscattering from the hole is a relevant perturbation . for a luttinger liquid these two regimes are distinguishes by whether the hole has a finite mass , case ( 1 ) , or an infinite mass , case ( 2 ) . in the sill
there is a subtlety because the relevance or irrelevance of backscattering from a local potential ( for example , an infinite mass core hole ) depends on the interactions in the system through the luttinger parameter @xmath71 . in particular ,
a local impurity is only a relevant perturbation if @xmath155 [ which is the same as @xmath69 ] @xcite . if @xmath157 the interactions of ll are repulsive , but a local impurity is an _ irrelevant _ perturbation .
therefore , in the sill we effectively have 3 possible regimes for the fermi - edge singularity : ( i ) finite mass mobile hole for which backscattering is irrelevant .
( ii ) infinite mass core hole for which backscattering is irrelevant . (
iii ) infinite mass core hole for which backscattering is relevant .
let us first consider the case of a finite mass hole .
in this case the backscattering is always an irrelevant perturbation . for addressing the threshold behavior of the photon absorption we can therefore neglect backscattering terms and focus only on forward scattering interactions between the electron and the hole . in the initial study of the fermi - edge singularity in a ll
@xcite it was incorrectly concluded that the exponent @xmath2 is independent of the dynamics of the hole , _ i.e. _ its mass .
a more careful treatment showed that indeed the mass of the hole enters @xmath2 @xcite .
the work of tsukamoto _ et al .
_ employed a combined bethe ansatz - conformal field theory study of a zero range interaction between particles .
the asymptotic structure of the energy and critical exponents in this model could be obtained exactly . from the formulas
obtained it was shown explicitly that the exponent @xmath2 depends on the mass of the hole .
perhaps the most important result is the structure of the electron - hole interaction : there are actually two forward scattering terms that emerge .
one is the familiar density - density interaction between hole density and electron density .
the other is less familiar and it describes the interaction of the hole with a current of electrons that is flowing by it . in other words , it is natural to transform to frame co - moving with the hole @xcite . in this frame
the hole sees a current of electrons .
one might think that such a picture would break time - reversal symmetry , but in an actual experiment electrons would be excited near both the right and left fermi points shown in fig.[fig : direct_transition](a ) .
thus , as a whole , time - reversal symmetry is preserved . as the contribution to the absorption edge is the same in both cases , one need only focus on one fermi point .
our discussion below will cover both fermi points in one fell swoop , as the electron operator eq.([eq : psi_final_bos ] ) contains both right and left moving pieces .
the fermi - edge singularity in a spin - incoherent luttinger liquid was discussed by @xcite .
the full hamiltonian for our problem in the frame co - moving with the hole is @xmath158 , where @xmath159 is given by eq.([eq : h_holon ] ) and eq.([eq : s_chain ] ) , and @xmath160 with @xmath161 and @xmath162 . here
@xmath163 is the symmetric part of the forward scattering from the hole and @xmath164 is the antisymmetric part of the forward scattering @xcite .
( in our convention @xmath165 represents the density fluctuations and @xmath166 the particle current . ) again , the antisymmetric part appears since in the frame of the hole , it sees a net current of particles scattering from it .
the `` + '' sign is for a right - moving hole and the `` - '' sign is for a left - moving hole .
the parameter @xmath164 depends on the momentum and mass of the hole , and when it is at rest , @xmath167 @xcite .
the operator @xmath168 creates a hole with spin @xmath96 and @xmath169 annihilates a hole with spin @xmath96 .
the energy of such a hole in its rest frame is @xmath170 , which in the presence of an externally applied magnetic field can depend on @xmath96 .    to compute the photon absorption
we treat the electron - photon interaction in second order perturbation theory .
standard manipulations show @xcite @xmath171 to evaluate the correlations appearing in eq.([eq : threshold ] ) we must diagonalize the hamiltonian .
the hamiltonian @xmath172 can be diagonalized with the unitary transformation @xmath173h^\dagger h\}\ ] ] where @xmath174 and @xmath175 applying this transformation we find @xmath176 , where the only change to @xmath177 is a shift in the hole energy @xmath178 , which is unimportant to us here .
our main concern is with computing the threshold exponent @xmath2 and the functional form of the response .
we begin with the evaluation of the correlation function @xmath179 , where @xmath180 is the imaginary time . at finite temperature ,
@xmath181 \\ & = & \frac{1}{z}{\rm tr}\left[e^{-\beta \bar h }   \bar \psi_\sigma(\tau)\bar h_\sigma(\tau)\bar h_\sigma^\dagger(0)\bar \psi_\sigma^\dagger(0)\right],\end{aligned}\ ] ] where in the second line we have introduced the unitary transformation @xmath182 and used the cyclic property of the trace .
we have already determined @xmath183 , so it remains to determine @xmath184 and @xmath185 .
direct evaluation gives @xmath186 ( up to unimportant multiplicative factors ) and @xmath187 } $ ] .
therefore , the correlation function @xmath188 separates as @xmath189 .
one readily finds @xmath190 at zero temperature .
this factor will enter the threshold frequency @xmath1 in eq.([eq : threshold ] ) and will in general also depend on the external magnetic field .
our main interest here is in the evaluation of the part of the correlation function that will give us the frequency dependence just above threshold , @xmath191 }   e^{i[\delta_a\theta(0)+\delta_s\phi(0 ) ] } \psi_\sigma^\dagger(0)\biggr],\end{aligned}\ ] ] where @xmath192 $ ] .
formally , eq.([eq : c_psi ] ) bears a striking resemblance to the single particle green s function evaluated in the spin - incoherent regime above .
it is worth pausing a moment to understand the physics of eq.([eq : c_psi ] ) . recalling that the operator @xmath193 creates a particle at space - time point @xmath194 in the many - body system and the electron number is related to the @xmath84 field via @xmath195 , with @xmath196 the average particle density
, we see immediately that @xmath197 involves adding an electron _ plus _ an additional  background excitation " ( from @xmath198 ) and then removing the same particle and its additional background at a time @xmath180 later .
the factors @xmath199 contribute additional density fluctuations coming from the motion of the finite mass valence hole .
this factor is absent in the infinite mass limit .
it remains to evaluate the trace in eq.([eq : c_psi ] ) .
we again follow the more general bosonization scheme for strongly interacting systems described above which reduces the to previously obtained results deep in the spin - incoherent regime @xcite .
we have @xmath200 where @xmath201 }   e^{-i[(1+\frac{q}{\pi}-\delta_a)\theta(0)-(1-\delta_s)\phi(0)]}\rangle.\ ] ] in the zero temperature limit with respect to the charge degrees of freedom , we can evaluate eq.([eq : f ] ) using eq.([eq : h_holon ] ) as we did for the single - particle green s function .
the result is @xmath202}.\ ] ] the frequency dependence of @xmath72 can now be computed from eq.([eq : threshold ] ) using @xmath189 .
the computation is essentially identical to that of the tunneling density of states , eq.([eq : tdos ] ) , only @xmath203 and @xmath204 - 1 $ ] .
therefore , we have @xmath205 where the @xmath206 factor tells us there is no photon absorption below below a threshold frequency @xmath1 . from fig.([fig :
direct_transition])(a ) it is clear that @xmath207 is of order the band gap plus the fermi energy of the electrons plus the absolute value of the energy of the hole near the fermi points of the electron gas . as with the expression eq.([eq : tdos ] ) for the tunneling density of states , eq.([eq : i_gen ] ) is the general expression for the photon absorption in an interacting one - dimensionsal system . for temperatures low compared to the charge energy , but arbitrary with respect to the spin energy , the heavy lifting is in the evaluation of @xmath137 as we emphasized in the discussion of the single - particle green s function .
following the same arguments we used there , we see that the frequency dependence @xmath208(1+q/\pi-\delta_a)^2}$ ] which is a maximum for @xmath209 . expanding about this point and doing the guassian integration
, one finds @xmath210 note that the final form of the photon absorption eq.([eq : i_mobile_final ] ) is almost identical to that of the tunneling density of states , eq.([eq : nu_structure ] ) , aside from the factor @xmath75 that appears in the threshold exponent .
it is instructive to compare the similarities and differences of the fermi - edge singularity with a mobile hole in the ll and sill . comparing with the general formula eq.([eq : i_fl ] )
we see that the threshold behavior is different in the sill compared to the ll and even the fermi liquid .
the threshold behavior of the sill contains the same logarithmic corrections that we found earlier in the tunneling density of states .
focusing on the power - law part of the threshold behavior we see that @xmath211 so that it is independent of @xmath74 , unlike the ll case in eq.([eq : a_ll_mobile ] ) .
since the asymmetric phase shift does not appear , this implies the threshold exponent is _ independent of the hole dynamics and therefore the mass of the hole _ in the spin - incoherent case . of course , from eq.([eq : i_mobile_final ] ) it is clear that the effects of @xmath74 are still felt , but only as a prefactor to the frequency dependence .
the effects of an externally applied magnetic field can also be easily incorporated @xcite . for a quantum wire
the magnetic field does not affect the orbital part of the wavefunction until the magnetic length is of the order of the wire width .
we therefore only discuss the coupling of the magnetic field to the spin of the electron , _
i.e. _ a zeeman coupling . the magnetic field will influence the spin correlations described by @xmath137 . in the spin - incoherent regime with zero field
one has @xmath212 , independent of @xmath96 . in a finite external field one instead has @xmath213 where @xmath214 and @xmath215 is the zeeman energy @xcite . while the presence of the magnetic field may shift the threshold energy
, it will not affect the exponent eq.([eq : a_sill_mobile ] ) @xcite .
the magnetic field also gives rise to a frequency scale @xmath216 below which sill behavior is observed and above which behavior characteristic of a spinless ll is observed @xcite , namely the form eq.([eq : i_fl ] ) with the exponent eq.([eq : a_ll_mobile ] ) .
we next turn to the case of an infinitely massive impurity for which backscattering is irrelevant . in this case
the effective electron - hole interaction is given by eq.([eq : e - h ] ) only with @xmath76 as this corresponds to the case of a stationary hole with no backscattering .
the photon absorption in this case is thus given by eq.([eq : i_mobile_final ] ) with @xmath76 .
while an infinite mass hole breaks translational symmetry , the irrelevance of backscattering effectively restores translational symmetry at low energies and the threshold behavior maps onto the finite mass case .    on the other hand , when backscattering is relevant the situation is very different . at low energies the hole acts as a boundary to a semi - infinite one dimensional system @xcite . in this case , the sill maps onto a _ spinless _ electron system @xcite which allows us to immediately exploit @xcite the results obtained for that system .
in particular , one obtains a universal contribution of 1/8 to the threshold exponent @xcite , @xmath217 and the frequency dependence follows eq.([eq : i_fl ] ) _ without _ the logarithmic corrections .
in fact , for the infinitely massive hole with relevant backscattering the only signature of spin - incoherent physics is hidden in the mapping @xmath87 @xcite .
while our theoretical understanding of the spin - incoherent regime has advanced rapidly , the difficulty of obtaining the regime @xmath10 has hampered experimental tests of the theory . to date our best experimental evidence comes from momentum resolved tunneling on cleaved - edge overgrowth quantum wires @xcite , transport in split - gate devices @xcite , and somewhat more speculatively in transport measurements on gated single wall carbon nanotubes with low electron density @xcite .
given that the existing theory now encompasses hybrid structures of sills with ferromagnets and superconductors as well @xcite , the number of experimental groups poised to make contributions to this exciting field is greatly enlarged .
for a recent discussion of numerous possible experiments that would provide strong evidence for spin - incoherent luttinger liquid behavior see @xcite .    on the theoretical side , a number of open issues remain . in particular , it is desirable to have a better understanding of the behavior on temperature scales @xmath218 which is most likely very relevant for many experiments that approach but do not quite reach the spin - incoherent regime
. related to this is the need for a better understanding of the crossover between the luttinger liquid and the spin - incoherent luttinger liquid regimes .
both of these issues will likely require a numerical attack as there are no obvious analytical methods available to address them .
finally there is the issue of spin - orbit coupling that has so far received no attention . for very strong spin - orbit coupling is there novel behavior in the regime @xmath219 ?
these issues await further study .
i am grateful to my collaborators ophir auslaender , leon balents , adrian feiguin , thierry giamarchi , bert halperin , anibal iucci , markus kindermann , karyn le hur , jiang qian , hadar steinberg , dagim tilahun , yaroslav tserkovnyak , and amir yacoby for their important contributions to our work on spin - incoherent luttinger liquids .
i am also indebted to matthew fisher , akira furusaki , leonid glazman , and kostya matveev for many enlightening discussions on this topic .
this work was supported by nsf grant no .
phy05 - 51164 and the lee a. dubridge foundation .