an exotic and distinctive feature of interacting one - dimensional systems is the phenomenon of spin - charge separation where the elementary excitations of the fermionic system are decoupled charge and spin bosonic modes that propagate at different velocities.@xcite at low energies , luttinger liquid theory@xcite has been quite successful at describing the properties of one - dimensional systems such as quantum wires@xcite and nanotubes.@xcite one of the elegant features of the theory is the ease with which one can obtain exact results , thanks to the quadratic form of the hamiltonian describing the low energy properties , @xmath0 , \label{1dhamltn}\ ] ] where @xmath1 stands for the charge and spin sectors and @xmath2 and @xmath3 are bosonic fields satisfying @xmath4=-i\frac{\pi}{2 } \delta_{ij}{\rm sgn}(x - x')$ ] , and representing charge / spin density fluctuations and charge / spin current density fluctuations , respectively , and @xmath5 are the collective mode velocities .
the parameters @xmath6 describe the interactions in the 1-d system . for repulsive interactions @xmath7 , and for su(2 ) symmetry in the spin sector @xmath8 .
generally , the theory is applied under the assumption that both the spin and charge modes operate at low energy .
but at low densities , @xmath9 , where @xmath10 ( with @xmath11 the bohr radius for the material ) , strong interactions open up a wide window between the characteristic energies of the charge and the spin sectors by suppressing the spin exchange energy @xmath12 while enhancing the charge energy @xmath13 .
when the temperature @xmath14 is raised into this window of energy , @xmath15 ( where @xmath16 boltzmann s constant ) , the spin sector consists of thermally randomized spins that are no longer described by ( [ 1dhamltn ] ) while the charge sector effectively becomes a spinless luttinger liquid ( ll ) .
a one - dimensional system in this regime is known as a spin - incoherent luttinger liquid ( sill).@xcite    several properties of the sill have already been established theoretically,@xcite but the experimental effort has been slowed due to the difficulty of reaching the aforementioned window of energy
. nevertheless , there are experimental indications of this regime in high quality quantum wires grown with the cleaved - edge technique,@xcite in recent split - gate devices,@xcite and possibly even in low - density carbon nanotubes.@xcite given the recent progress in manufacturing hybrid devices consisting of superconductors ( sc ) and ferromagnets ( fm ) , and one - dimensional systems such as nanotubes , we propose here that the properties of the sill be experimentally investigated by studying junctions of superconductors and ferromagnets with a sill . as we show below , these sill hybrid systems lead to distinctive behaviors ( relative to a ll or non - interacting 1-d system ) that will allow the sill to be identified in various types of transport measurements , and in the dynamics of a coupled ferromagnet .    in early work on transport in the spin - incoherent regime
, matveev found that a sill adiabatically connected to non - interacting leads resulted in a universal reduction of the conductance of a single - channel quantum wire from @xmath17 to @xmath18.@xcite this effect has been ascribed to the deactivation of the spin channel , whereby the spin excitations are reflected while the charge ones propagate freely . in a recent related work
, the authors studied the problem of a sill adiabatically connected to two superconductors .
it was found that the critical josephson current through such a sc - sill - sc systems suffers the same fate : the critical value of the josephson current is halved relative to the ll case.@xcite these universal results follow from the assumption of adiabaticity of the contacts between the sill and the leads / sc , which may not be satisfied in some experimental situations .    in this work
we study the transport properties of a sill in the opposite limit where it is contacted via tunnel junctions to a superconductor and/or a ferromagnet .
we will discuss several scenarios that may be implemented experimentally with existing technologies .
the central aim of this work is to contribute to our fundamental understanding of the sill by addressing its properties in hybrid structures , and to add to the arsenal of experiments that can be carried out to probe it and infer its existence .
our main results are the following .
for the case of a fm tunnel coupled to a sill , we compute the spin current pumped into the sill due to a time - dependent magnetization in the ferromagnet
. the algebraic form of the expression for the pumped spin current is the same as that for non - interacting fermions , or a luttinger liquid , but the coefficients appearing have a temperature dependence characteristic of the sill .
we argue the spin transport is diffusive and under this assumption compute the renormalization ( due to the pumped spin current and its back - flow ) of the gilbert damping @xmath19 of the magnetization motion in the fm .
this also has a characteristic temperature dependence in the spin - incoherent regime . for the case of a sc tunnel
coupled to a sill , we evaluate the ac and dc josephson current in a wire that can be either side - coupled ( bulk - coupled ) or end - coupled to the sc .
we find the dc josephson current suffers an exponential decay in real space due to the incoherence of the spin sector ( in contrast to the inverse length suppression that was found for adiabatically contacted scs@xcite ) .
finally , we show that the ac josephson effect also contains important information about spin - incoherent effects in its voltage and temperature dependence .
this paper is organized as follows . in sec .
[ sec : bosonization ] we review bosonization for a strongly interacting electron system for energy scales small compared to the charge energy , but arbitrary compare to the spin energy . in sec .
[ sec : fm ] we study spin pumping and magnetization dynamics in the fm - sill system . in sec .
[ sec : sc - sill ] we study ac and dc josephson effects in the sc - sill - sc system . finally in sec .
[ sec : discussion ] we discuss prospects for experimental realizations of the physics we discussed and directions for future work .
some technical details and results appear in various appendicies .
in this paper we will be studying a strongly interacting one dimensional electron system coupled to either a ferromagnet or a superconductor .
strong interactions imply@xcite that @xmath20 and therefore it is possible to be in a regime of temperatures that may not be small compared to the spin energy @xmath12 , but are still small compared to the charge energy @xmath13 .
this means that we are no longer free to use standard bosonization procedures@xcite for the electron operator because these formulas rely on the assumption that the relevant energy scales on which the system is probed are small compared to both spin and charge energies .    for a strongly interacting system
we require a formalism that applies to all energy scales relative to the spin energy ( and in particular to @xmath21 ) , but only low energy scales to relative to the charge energy .
such a formalism was recently developed by matveev , furusaki , and glazman,@xcite and applied to the evaluation of the spectral function and tunneling into a strongly interacting electron system at arbitrary temperatures relative to the spin energy.@xcite in this section we review the essential elements of their work and show how to represent the electron operator in a strongly interacting one dimensional electron system .
the general expressions for electron correlation functions obtained in this way provide a launching point for numerical studies in the intermediate ( with respect to spin ) temperature regime @xmath22 , which is not well handled by existing analytical methods .      as is typical of interacting one dimensional electron systems
, we assume the hamiltonian in the strongly interacting case is spin - charge separated,@xcite @xmath23 . here
, @xmath24 is identical to that of the ll , and is given by eq .  .
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 @xmath25 where evidently the spin energy is set by @xmath26 : @xmath27 .
the basic idea is to represent the electron operator as a product of operators that describe the holons @xmath28 ( spinless fermions that naturally arise in the context of strongly interacting fermions and the spin - incoherent regime@xcite ) and the spin degrees of freedom @xmath29 .
the holon operators ( denoted by @xmath30 ) by construction satisfy the equation @xmath31 where @xmath32 is the electron annihilation operator for electrons of spin projection @xmath33 , and @xmath34 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 @xmath35 limit of the one dimensional hubbard model .
et al._@xcite wrote the electron creation operator as @xmath36 where @xmath37 creates a site on the spin chain with spin projection @xmath33 .
the expression can be physically motivated as follows . from 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 @xmath37 . in general
, one has @xmath38 as the object that adds a new site to the spin chain between @xmath39 and @xmath40 .
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 matveev , furusaki , and glazman showed@xcite that a remedy for this issue is to define the position at which the spin site is added to the chain in terms of the number of holons to the left of the site , @xmath41 in terms of the electron creation and annihilation operators are defined as @xmath42 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 @xmath43 also depend on the electron density via .      in writing eqs .
, 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 @xmath13 , but arbitrary with respect to @xmath12 . in this case
, we are free to bosonize the holon sector .
the hamiltonian for the holons must then necessarily take the low energy form@xcite @xmath44 , \label{eq : h_holon}\ ] ] where the interaction parameter @xmath45 of the holon is related to the interaction parameter of the charge sector as@xcite @xmath46 , and the spinless fields @xmath47 and @xmath48 can be related to the holon density as@xcite @xmath49,\ ] ] where the holon fermi wave vector is twice the electron fermi wavevector@xcite @xmath50 .
the bosonic fields satisfy the commutation relations @xmath51=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 @xmath52 , @xmath53 where @xmath54 destroys an holon near the right fermi point and @xmath55 destroys an electron near the left fermi point .
the left and right holon operaters are bosonized as @xmath56},\ ] ] where @xmath57 is a short distance cut - off of order the interparticle spacing @xmath58 . combining the results of , , , and one obtains the bosonized form of the electron annihilation operator for spin @xmath33 @xmath59}+e^{-i[k_f^hx+\theta(x)]}\right)\nonumber \\
\times z_{l , s}\bigg |_{l=\frac{1}{\pi}[k_f^hx+\theta(x)]},\end{aligned}\ ] ] and an analogous expression for the electron creation operator @xmath60 .
expression , 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 @xmath61 } \to \sum_l z_{l , s}\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 @xmath62},\ ] ] where @xmath63 the expression for the electron annihilation operator is the key result of the this section , obtained earlier by matveev , furusaki , and glazman , who also showed that in the limit of small energies compared to @xmath12 the expression correctly recovers the standard ll formulas for the electron annihilation operator.@xcite with 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 @xmath12 , but small energies compared to @xmath13 .
see appendix  [ app : spin_current ] and appendix  [ app : andreev_tdos ] for examples of the evaluation of the single particle green s function near different types of boundaries using the bosonization formula and the boundary conditions to be discussed next .      in this work we will be primarily interested in two types of boundary conditions on the electron operators : ( 1 ) open or  hard wall " boundary conditions appropriate for tunnel junctions and ( 2 ) andreev boundary conditions appropriate for adiabatically contacted ( no electron backscattering ) superconductors .
other , more general ,  intermediate " boundary conditions are possible though less generic.@xcite to help orient the reader , in each of the two cases above we will briefly summarize the result appropriate for the ll regime for the purpose of drawing contrast with the strongly interacting , finite temperature sill regime .
our discussion of the boundary conditions extends the results of ref .  [ ] .
open boundary conditions , or  hard wall " boundary conditions result in electron waves that are perfectly reflected at the boundary .
this implies there is no charge current or spin current through the boundary . for concreteness ,
let us assume that our boundary is located at @xmath64 with the interacting one dimensional system living on @xmath65 .
then an electron traveling to the left with spin @xmath33 will be reflected to a right moving electron with the same spin @xmath33 @xmath66 with a phase shift @xmath67 that depends on details of the boundary scattering potential.@xcite from the standard ll bozonization formulas ( in our convention ) @xmath68},\ ] ] we see that implies that @xmath69 , and therefore also that @xmath70 and @xmath71 . in computing correlation functions at the boundary
, this means that we must take @xmath72 and @xmath73 to be non - fluctuating quantities , allowing only @xmath74 and @xmath75 to fluctuate .    in the strongly interacting case
, the reflection condition implies that @xmath76 , where we have used . since the spin site @xmath40 is related to the @xmath47 field via and we have @xmath77 at @xmath64 .
this implies any correlation function involving electron operators evaluated at the boundary will depend only on @xmath78 or its hermitian conjugate , and the fluctuations of the field @xmath48 at @xmath64 .
we will apply these boundary conditions in a calculation of the single particle green s function in appendix  [ app : spin_current ] .    for completeness , below we give the expansions of the bosonic fields for a system of finite length @xmath79 with open boundary conditions at @xmath64 and @xmath80:@xcite @xmath81 where @xmath82 , @xmath83=\delta_{nm}$ ] , and @xmath84=i$ ] . for a semi - infinite system , we take @xmath85 and the discrete sums over @xmath86 become integrals over momentum @xmath87 .      in a certain sense
the andreev limit is the opposite limit of  hard wall " scattering from a sc - m , or sc - ll , or sc - sill interface . whereas open boundary conditions imply an incident electron
is perfectly reflected , andreev boundary conditions imply that an electron is perfectly absorbed by the sc ( with a concomitant reflected hole of the opposite spin).@xcite however , compared to the open boundary conditions the andreev boundary conditions are more subtle as they involve an energy scale , the superconducting gap @xmath88 , that under many circumstances can not be taken to be infinitely large relative to the energy scales of interest ( such as @xmath89 , or an applied voltage ) in the interacting one - dimensional system .
one of the most important consequences of a finite @xmath88 is an energy - dependent reflection coefficient , which ultimately leads to the proximity effect in the normal material.@xcite in the context of interacting one dimensional systems , andreev boundary conditions have been discussed by a number of authors.@xcite the conclusion of these works is that when a left moving electron with spin @xmath33 is reflected as a right moving hole with the opposite spin there is an energy dependent phase shift @xmath90 ( proportional to the momentum difference @xmath91 with respect to the fermi point)@xcite that multiplies a factor@xcite @xmath92 that encodes the phase @xmath93 of the superconducting order parameter ( assumed non - zero for @xmath94),@xcite @xmath95 where @xmath96 is the superconducting coherence length .
the function @xmath97 for @xmath98 and @xmath99 for @xmath100 .
the boundary conditions are valid only at energy scales much smaller than @xmath88 , which we will assume throughout this work .
applying the boundary conditions to the ll case ( with bosonized electron operator below ) gives @xmath101 and @xmath102 .
therefore , we find that much like the situation of perfect reflection there is no spin - current through the interface , but there is is a net charge current .
moreover , analysis of singlet superconductivity and spin density wave correlation functions using in the ll regime shows that there are suppressed spin fluctuations near ( distances less than @xmath103 ) the interface.@xcite it is perhaps worth noting that if the sc - ll interface has a very weak electron backscattering , the interactions in the ll tend to renormalize the interface scattering.@xcite    for strongly interacting electrons the boundary condition implies @xmath104 and @xmath105 at the interface , where we have again used and also .
we show in appendix  [ app : andreev_tdos ] that we recover our earlier results@xcite for sc - sill junctions in the andreev limit using this formalism .    finally , for completeness and for unification of our notation , below we give the expansions of the bosonic fields@xcite for andreev boundary conditions on a system of length @xmath79 ( see fig .
[ fig : sc_tunnel ] ) with superconducting order parameter phase difference @xmath106 : @xmath107(b_m + b^\dagger_m)+\theta^0 , \nonumber   \\
\phi(x)&=&i\sum_{m=1}^\infty \sqrt{\frac{1}{2k_\rho m}}\sin\left[\frac{m \pi } { l}\left(x+{\xi \over 2}\right)\right](b_m -b^\dagger_m)\nonumber \\ & & \hspace{5cm}+\phi(x),\end{aligned}\ ] ] where @xmath108 , @xmath83=\delta_{nm}$ ] , and @xmath109=i$ ] . the topological number @xmath110 and the total spin of the system @xmath111 must satisfy the constraint @xmath112even.@xcite we have also included the proximity effects@xcite via the length @xmath103 in the field expansions that appeared earlier in the general expression for the boundary conditions appropriate to andreev reflection.@xcite    note that compared to the field expansions for open boundary conditions , the expansions for andreev boundary conditions have @xmath47 and @xmath48  switched "
. this can be understood in simple physical terms : @xmath47 and @xmath48 are conjugate fields so if one is a constant , the other is strongly fluctuating .
thus , the  switching " of the fields comes from  opposite " nature of the two boundary conditions in the charge sector . for open boundary conditions
there is no charge current through the boundary , while for andreev boundary conditions the charge current on either side of the boundary is unchanged by its presence because there is no electron backscattering there .
had we been concerned with the bosonization of the spin sector , we would have found that due to the absence of spin current through the boundary that sector would have expansions appropriate for open boundary conditions similar to .@xcite    having spent the time to develop the formalism used in our calculations , we now turn our attention squarely to the physics of hybrid junctions involving spin - incoherent luttinger liquids .
we first discuss a junction consisting of a ferromagnet tunnel coupled to a spin - incoherent luttinger liquid .
we consider the set - up shown schematically in fig.[fig : fm ] whereby a ferromagnet is coupled via a tunnel juntion to a spin - incoherent luttinger liquid .
a similar situation was considered by bena and balents in the context of a fm - ll junction.@xcite they found that when the magnetization @xmath113 of the fm acquires a time dependence ( perhaps by the application of external fields ) , spin current is pumped into the adjoined ll at the rate @xmath114 where @xmath115 and @xmath116 are temperature dependent parameters depending on the interface tunneling and interactions in the ll . below
we show that an identical expression is obtained for a sill and we derive explicit expressions for @xmath115 and @xmath116 in this case .
we schematically indicate how the temperature dependence crosses over from the ll to the sill case .
we note that was earlier derived for a non - interacting metal attached to a fm . in this case ,
@xmath115 and @xmath116 are temperature independent to lowest order.@xcite thus , the temperature dependence of @xmath115 and @xmath116 encodes information about the interactions and may also indicate whether the attached 1-d system is in the ll or sill regime .
we now turn to a derivation of these results .
we assume that the ferromagnet is itinerant , such as fe or co , and can be described by a stoner - type model in which there is a different density of states for spin - up and spin - down electrons .
the different density of states can be absorbed into distinct effective tunneling matrix elements , @xmath117 , for spin - up and spin - down electrons.@xcite this also implies that the local action@xcite describing the spin up and spin down electrons are identical to that of non - interacting electrons @xmath118 where @xmath119 is the inverse temperature and @xmath120 are the bosonic fields describing local fluctuations at the tunneling point , @xmath64 , of spin up and spin down electrons in the fm . for the sill we have the following spin - charge separated form of the local action @xmath121 where @xmath122 is the luttinger parameter of the charge sector of the sill , @xmath123 is the bosonic charge field , and @xmath124 is the action for the spin sector which will actually play no role in the evaluation of correlation functions in the spin - incoherent regime as it effectively drops out leading to  super universal spin physics".@xcite    for an arbitrary spin quantization axis relative to the magnetization @xmath125 it is useful to define @xmath126 which projects the spin quantization axis onto the magnetization direction @xmath127 . here @xmath128 is a @xmath129 matrix and @xmath130 are the pauli spin matricies .
the tunneling hamiltonian then takes the form@xcite @xmath131 where @xmath132 ( @xmath133 ) annihilates / creates an electron in the fm and @xmath134 ( @xmath135 annihilates / creates an electron in the sill .
the tunneling is assumed to occur at @xmath64 as shown in fig.[fig : fm ] .
the @xmath136 tunneling matrix @xmath137 is then @xmath138 the tunneling hamiltonian can be expressed more explicitly as @xmath139 where @xmath140 is the time - dependent magnetization of the ferromagnet , and @xmath141 and @xmath142 .
( from here onwards we will assume @xmath143 . )
the spin current operator is obtained from the relation @xmath144 $ ] , where @xmath145 is the spin density in the sill at the boundary and summation over repeated spin indices is understood .
it follows then that the spin current operator is given by @xmath146 the spin current , @xmath147 \rangle$ ] , is obtained@xcite from second - order time dependent perturbation theory which gives @xmath148-\int \frac{d \omega}{h } c(\omega)e^{-i\omega t}\times\nonumber\\\left(\text{im } ( u_2^*u_1)\vec{m}(\omega)+\frac{|u_2|^2}{2}\vec{m}(\omega)\times\vec{m}(t)\right),\end{aligned}\ ] ] where @xmath149 is the fourier transform of the retarded greens function @xmath150\rangle$ ] .
if we assume the typical frequencies of the magnetization precession are in the ghz range@xcite ( which corresponds to energies of roughly 100 mk ) or smaller , then this is a small energy scale in the problem and we may safely expand @xmath149 for small @xmath151 . it can be easily checked that the zero frequency terms cancel exactly , leaving only the contributions linear in @xmath151 , provided we drop terms proportional to @xmath152 and all higher powers . upon integration over frequency ,
the linear frequency terms are converted to time derivatives yielding expression where the coefficients @xmath115 and @xmath116 are given by @xmath153 since only the imaginary part of @xmath149 is odd with respect to frequency , only this piece will contribute to @xmath154 and those parameters will therefore be real quantities . the final step is to compute the temperature dependence of @xmath155 in the spin - incoherent regime .
since we are working within linear response , the commutator in @xmath156\rangle$ ] is evaluated in the state where there is no tunneling between the ferromagnet and the sill , and a hard - wall at @xmath64 . by
writing @xmath157 $ ] , and noting that the long time behavior of both terms has the same functional dependence , we can easily extract the long time behavior@xcite of @xmath158 using and .
( a more careful calculation that yields the same result is given in appendix  [ app : spin_current ] . ) at finite temperatures ( small with respect to the fermi energy of the ferromagnet and charge energy of the sill , but large compared to @xmath159 ) , we have @xmath160 and @xmath161^{1/(2k_\rho)}$ ] .
the crucial difference with the luttinger liquid is that the exponent for the correlations @xmath162 have changed : in the case of a ll @xmath163 is replaced by @xmath164 .
the remainder of the computations carry through exactly as they would for a ll and we find for @xmath165 @xmath166 \propto \hbar \omega \left(k_b t\right)^{\delta^{sill}},\ ] ] where @xmath167 .
this implies @xmath168 quite generally so that @xmath169 where @xmath170 interpolates between the ll and sill regimes as the temperature is swept through the spin energy @xmath171 .
the temperature dependence of the exponent @xmath170 is shown schematically in fig .
[ fig : delta ] .
it is interesting to note that @xmath172 for @xmath122 in the luttinger liquid regime implying that less spin current is pumped as the temperature decreases . on the other hand , so long as the system remains in the spin - incoherent regime the opposite behavior may be obtained if @xmath173 : since @xmath174 _ more _ spin current is pumped as the temperature is lowered ( for @xmath175 ) .
this is related to the diverging density of states at the boundary@xcite when @xmath173 . in gated cleaved - edge overgrowth
quantum wires it appears possible to lower @xmath176 down to values of order @xmath177 and so it may be possible in the spin - incoherent regime to tune between the positive and negative exponent regimes.@xcite this qualitative difference should be easily seen in experiment .
having computed the spin current pumped into the sill by a time - dependent magnetization vector , it is important to ask how the spin accumulation in the sill in turn affects the magnetization motion .
we address this question by computing the renormalization of the gilbert damping constant due to the spin flow into the sill .
we again closely follow the notation of bena and balents@xcite to clearly establish a connection to the ll case .
the landau - lifshitz - gilbert equation for magnetization @xmath178 precessing around effective magnetic field @xmath179 is @xmath180 where the spin current @xmath181 flows ( @xmath182 is given by for @xmath183 and is non - zero only for time dependent @xmath178 ) into the sill and @xmath184 is the saturation magnetization of the ferromagnet . here @xmath185 is the gyromagnetic ratio which is typically equal to its free electron value,@xcite @xmath186 , in transition metal ferromagnets , and @xmath187 is the dimensionless gilbert damping parameter in the absence of the spin current .
its value is typically of order @xmath188 .
the spin backflow due to spin accumulation in the interacting 1-d system can be described by boundary conditions on left and right moving spin currents and is expressed as@xcite @xmath189 where @xmath190 is the spin chemical potential in the wire related to the magnetization by @xmath191 , with @xmath192 the ( generally temperature dependent ) spin susceptibility , and @xmath193 describes the effective exchange coupling between electrons in the sill and the fm.@xcite the current @xmath194 arises from the electron tunneling contribution to the spin current.@xcite note that the back - scattered spin current contains two terms : ( i ) a term arising from electron transfer from the fm to the interacting 1-d system and ( ii ) a second term arising from exchange between the local spin density and the precessing magnetization vector @xmath178 . since @xmath181 this implies the spin current itself has a contribution due purely to exchange effects , even in the absence of electron transfer.@xcite    our goal in this section is to express as @xmath195 and determine the renormalized parameters @xmath196 and @xmath19 . from it
is evident that we must find @xmath197 .
we have already computed @xmath182 , eq .
, in the previous subsection .
now we must determine @xmath198 , , which depends on @xmath190 and @xmath194 .
we start with @xmath190 which is a function of the spin dynamics in the interacting 1-d system .
it has earlier been shown that su(2 ) invariant electron backscattering leads to diffusive spin behavior in the weakly interacting regime.@xcite in the strongly interacting regime where the spin sector may effectively be described by a heisenberg spin chain coupled to phonon distortions@xcite the spin dynamics has also been shown to be diffusive.@xcite in both limits , the mean free path @xmath199 . following the arguments@xcite that the spin - incoherent regime can be approached from  below " ( @xmath200 ) , we expect the diffusion length @xmath40 to saturate at the interparticle spacing , @xmath58 , for @xmath175 .
since the diffusion constant @xmath201 , this implies the diffusion constant in the sill becomes independent of temperature and takes the value @xmath202 .
note that in the limit @xmath203 , there is no spin diffusion , as the spin excitations can not propagate in the system .
spin will simply pile up at the boundary of the fm without moving further into the interacting 1-d system . unlike our earlier results for the temperature dependence of pumped spin into the sill , eq . with @xmath204 , which were valid for vanishing spin velocity , here
our results are qualitatively dependent on keeping @xmath205 finite , although still small enough to be in the spin - incoherent regime .
we also assume the finite length sill is characterized by a spin - flip time @xmath206 ( due to impurities , spin - orbit effects , etc ) .
the diffusion equation for spin in the sill is then @xmath207 with the boundary conditions @xmath208 at @xmath64 and vanishing spin current @xmath209 at @xmath80 .
see fig .
[ fig : fm ] for the set - up .
the solution to this equation is simple to obtain and is given by @xmath210}{\kappa \sinh[\kappa l]}\left(\frac{1}{d_s
\chi(t)}\right ) \vec i(x=0),\ ] ] where @xmath211 . if the precession frequency @xmath151 is small compared to the inverse spin relaxation time , then to a good approximation @xmath212 . at the boundary @xmath64
we have @xmath213 where @xmath214 \frac{1}{\kappa d_s \chi(t)}$ ] .
note that @xmath215 so the spin current decays exponentially with distance into the sill on a length scale set by the product of the diffusion constant and the spin - relaxation time . for a long spin - relaxation time
, this length scale can be large compared to the inter - particle spacing .
we have now determined all parts of the spin back flow , except the tunneling current @xmath194 , which we now do .
the tunneling current is proportional to the imaginary part of the fourier transformed correlation function@xcite @xmath149 defined below eq .
, @xmath216\nonumber \\
\propto |u_1|^2(k_b t)^{\delta(t)+1 } \sinh\left(\frac{\mu_s}{4 k_b t}\right ) \nonumber \\
\times \left|\gamma\left(1+\frac{\delta(t)}{2}+i\frac{\mu_s}{4\pi k_b t}\right)\right|^2 \nonumber \\ \propto \mu_s |u_1|^2(k_b t)^{\delta(t)}\left|\gamma\left(1+\frac{\delta(t)}{2}\right)\right|^2,\end{aligned}\ ] ] where we have used the results of appendix [ app : spin_current ] and taken the limit @xmath217
. evidently , the main effect of the spin - incoherent physics is to change the exponent @xmath170 to the spin - incoherent value , @xmath204 , so that we have @xmath218 in the spin - incoherent regime .
let us define @xmath219 which is temperature dependent and @xmath220 independent .
at this point the determination of @xmath196 and @xmath19 is identical to the ll case.@xcite we find @xmath221 and @xmath222 where @xmath223 $ ] , @xmath224 $ ] , and @xmath154 are given in eqs.([eq : a_gen],[eq : a_temp ] ) .
the temperature dependence in @xmath225 is entirely contained in @xmath226 and @xmath103 .
see table  [ tb : ll_sill ] below for a comparison of the ll and sill regimes .
[ h ]    .temperature dependence of key quantities in ll - fm and sill - fm hybrid systems .
here @xmath227 is a constant that depends on @xmath228 , and @xmath229 . [ cols="^,^,^,^,^,^",options="header " , ]     as with the case of a ll , for a sill we expect there to be little renormalization of @xmath196 relative to @xmath185 , but the smallness of the gilbert damping @xmath187 means this may obtain a significant temperature dependent correction depending on @xmath230 .
one could easily imagine other scenarios such as a fm - sill - fm junction , a fm - sill - m junction , or even a fm - sill - sc junction .
however , because the spin transport is generally diffusive , for junctions whose length @xmath79 is long compared to @xmath231 , the two contacts to the sill will more or less behave independently from the point of view of spin transport .
to the extent that the two leads are coupled , it is evident from the discussion of the previous section that most of the physics ( and corresponding general formulas ) present for a ll system also apply to the sill system only with some modifications in the temperature dependence of parameters appearing in the theory .
while this may seem like a somewhat trivial result , it is not .
the temperature dependence serves as a means to determine whether spin - incoherent physics is present in the system .
in particular , we note quite generally that if @xmath232 the temperature dependence of many quantities ( e.g. , @xmath233 and those derived from them ) change _ qualitatively _ relative to the expectations for a ll .
such qualitative differences should be observable in experiment .
it is worthwhile to step back and emphasize some of the essential differences between the sill and ll cases .
first , we note that the different temperature dependence can be traced to three quantities : ( 1 ) the correlation function @xmath158 defined below eq .  .
( 2 ) the temperature dependence of the spin susceptibility @xmath93 . ( 3 ) the temperature dependence of the spin diffusion constant @xmath229 .
so long as the sill is tunnel contacted to the fm , the correlation function @xmath158 will always appear at lowest order in perturbation theory and carry along with it the characteristic temperature dependence of the sill .
this quantity appears in both the description of the pumped spin current and the magnetization dynamics . on the other hand , when we are specifically interested in how the spin propagates in the sill [ as we saw for the backscattered spin current ] the diffusion constant @xmath229 enters , and also the spin susceptibility via the einstein relation for the spin conductivity , @xmath234 .
we remark that the spin transport is generically diffusive in the spin - incoherent regime , while it may be either ballistic or diffusive in the ll regime.@xcite    having now flushed out what we feel are the central considerations and results for fm - sill hybrid structures , we now turn our attention to sc - sill hybrid structures .
throughout this work , we assume that all energy scales are small compared to the superconducting gap , @xmath88 , unless explicitly stated otherwise , such as in the limit of a short junction ( defined below ) .
we are primarily interested in the case where the superconductor is contacted via a tunnel junction with a spin - incoherent luttinger liquid .
let us begin our discussion of sc - sill hybrid systems by considering the simplest case : where the fm in fig .
[ fig : fm ] is replaced by a sc .
earlier we studied the properties of such a junction in the andreev limit.@xcite we found a number of remarkable properties , including a completely universal ( independent of the charge and spin hamiltonians ) tunneling density of states in the spin - incoherent regime.@xcite in the opposite limit of a tunnel contact between the sc and sill , no proximity effects of the sc are felt other than those perturbative in the tunneling.@xcite ( in contrast to the andreev limit where the behavior of the pair correlations and tunneling density of states follows directly from the form of the field expansions imposed by the boundary conditions.@xcite ) as the proximity effects are already weak in the andreev limit ( they decay exponentially fast with distance , unlike the power law decay in the ll@xcite ) we do not pursue the very weak proximity effects in the spin - incoherent regime in the limit of a weak tunneling sc - sill junction .    on the other hand ,
it is worthwhile to briefly discuss the behavior expected in opposite limit of nearly perfect andreev reflection with weak ( perturbative ) electron backscattering at the sc - sill interface . in the ll context , this was addressed by vishveshwara _
et al._@xcite in a sc - ll - m system who found dips in the conductance of the junction as a function of gate voltage .
these dips correspond to multiple electron reflections off the sc interface and occur at voltages corresponding to integer values , @xmath235 , of the inverse electron traversal time , @xmath236 , where @xmath9 is an integer counting the number of traversals in the ll.@xcite for a sc - sill - m system we would expect similar qualitative effects , although with an overall suppression in the conductance.@xcite finally , we briefly remark that the conductance per channel of an adiabatic sc - sill junction should drop to 1/2 the value expected for a sc - ll junction:@xcite @xmath237 .
recall that the conductance of the adiabatic sc - ll interface @xmath238 generalizes the result for the adiabatic sc - m interface@xcite @xmath239 realized for @xmath240 .
we now turn our attention to the final topic of this work , sc - sill - sc junctions of the types shown in fig .
[ fig : sc_tunnel ]
. our primary focus will be on the behavior of the ac and dc josephson effects in such a hybrid structure . in the context of sc - ll - sc structures there
have already been a number of theoretical works@xcite and some experimental results.@xcite    earlier we studied the dc josephson effect in a junction like that shown in fig .
[ fig : sc_tunnel](a ) only with andreev boundary conditions at the interfaces of the sill with the scs.@xcite we found that in spite of the fact that the pair correlations decay exponentially with distance from the boundary of the sc , the josephson critical current scaled inversely with the length of the sill and was reduced by a factor of 2 relative to either the sc - ll - sc or sc - m - sc case ( which have identical critical currents that also scale inversely with the length).@xcite we remarked that this  robustness " is essentially related to the fact that the superconductor phase difference only couples to the charge degrees of freedom in the sill , which remain completely coherent .
this effect can also be taken to support the notion@xcite that in some sense superconductors are spin - charge separated.@xcite yet another perspective on this result can be obtained by noting that the andreev boundary conditions on a system of finite length can be mapped onto a system with _
periodic _ boundary conditions of twice the length.@xcite in this way , the josephson response becomes equivalent to the problem of persistent currents for _ spinless _ electrons in a ring threaded by magnetic flux , with the superconductor phase difference being simply related to the flux .
therefore , the dc josephson current for andreev boundary conditions is essentially given by the physics of _ spinless _ electrons which is why spin - incoherent physics does not manifest in a dramatic way .
the factor of 2 reduction in the critical current and josephson period are a direct consequence of the mapping of the sill onto spinless electrons .
we now turn our attention to the case of open boundary conditions at the interfaces between the scs and sill . in this case phase information
must be carried between the two superconductors by cooper pair tunneling between them .
because the _ pairs _ must hop from one sc onto the sill , then onto the other sc , the josephson current will occur at 4th order in this tunneling process .
moreover , because the cooper pair tunneling process creates a pair of electrons locally in the sill near the barrier , many fourier modes of the electron operator will come into play .
this is in contrast to the case of andreev boundary conditions where only the zero mode enters in the evaluation of josephson current.@xcite this crucial difference will turn out to dramatically suppress the josephson response relative to what one would have for a ll when the sill is tunnel coupled to scs . in short ,
we find the critical current decays as an exponential of the junction length ( compared to a power of the inverse junction length for a ll at zero temperature ) with the characteristic decay length the inter - particle spacing .
the sc - sill - sc tunnel junctions shown in fig .  [ fig : sc_tunnel ] are modeled by the following hamiltonian , @xmath241 where @xmath242 are @xmath33-wave bcs hamiltonians for the superconductors and @xmath243 is the sum of eqs . and .
the tunneling hamiltonian @xmath244 is given by @xmath245 we assume the scs support a phase difference of @xmath246 and have the same superconducting gap @xmath88 .
the josephson current is obtained from the relation @xmath247 where @xmath248 is the electronic charge and @xmath249 is the partition function .
we employ imaginary time perturbation theory , where terms dependent on @xmath93 appear in fourth order,@xcite @xmath250 where @xmath251 is a two - particle ( cooperon ) propagator in the one - dimensional interacting system .
the  similar terms " account for time ordered permutations and spin projections ( there are 23 of them , @xmath252 all together ) .
the propagation of cooper pairs in the superconductors is described by the time - ordered anomalous green function @xmath253 , @xmath254 where the expectation value is taken with respect to @xmath255 and the electron operators in are evaluted at @xmath64 for @xmath256 and @xmath80 for @xmath257 .    of the 24 terms
( plus their hermitian conjugates ) that appear in the partition function , eq .  , physical considerations aid us in choosing the most relevant ones depending on the length of the sc separation distance @xmath79 .
there are two important length scales in the junctions we study that give rise to corresponding time scales : the superconducting coherence length @xmath103 and the junction length @xmath79.@xcite we define  long " to mean @xmath258 and  short " to mean @xmath259 .
we now turn to a discussion of the general expression for @xmath251 in the limit of long and short junctions .    for long junctions the tunneling into the sill
is  fast " and the propagation of the cooper pairs is  slow " .
the relevant two - body propagator is @xmath260 where @xmath32 is the bosonized electron annihilation operator given by eq . and the averaging is taken with respect to @xmath243 .
we have @xmath261 , @xmath262 , and @xmath263 .    in the opposite limit of a short junction @xmath259 ,
the propagation through the sill is  fast " compared to the  slow " tunneling of the cooper pairs , so the two - particle propagator approximately separates into a product of two single - particle green s functions@xcite @xmath264 where we have @xmath261 , @xmath262 , and @xmath265 ( note the re - ordering of electron operators relative to ) .    _
end contacted sill when @xmath266 _ : we first consider the case of an  end - contacted " sill , shown in fig.[fig : sc_tunnel]a .
for this geometry the expansions of the holon field operators are the appropriate ones .
our assumption that all energy scales are small compared to @xmath88 ( specifically @xmath267 , and @xmath268 ) enables us to approximate the anomalous green s function with a delta function in time , @xmath269 where @xmath270 is the density of states of the sc @xmath271 at the fermi energy in its normal state .
this approximation forces the times the electrons tunnel into or out of the sill to coincide , @xmath272 and @xmath273 . without loss of generality ,
we let @xmath274 and define @xmath275 . using the expression for the electron operator , eq . , the propagator reads    @xmath276    the calculation of the cooper pair propagator is carried out in the same manner as that of the green s function presented in appendix [ app : spin_current ] .
we first evaluate the expectation value of the spin chain site creation / annihilation operators , @xmath277 in the spin - incoherent regime .
to this end it is convenient to fully exploit the symmetries and boundary conditions of the problem .
we first note that the open boundary conditions at @xmath278 force @xmath279const at both these points .
as before , we take this constant to be zero so that all @xmath47 fields effectively drop out of eq .. the momentum integrals in can then be done trivially to give delta functions on the sites @xmath280 which sets @xmath281 and @xmath282 . therefore , we have @xmath283 which in the spin - incoherent regime can be evaluated quite simply by appealing to the physics of the spins . here
@xmath284 measures the amplitude for two spins of opposite orientation introduced at @xmath80 to arrive at @xmath64 .
deep in the sill regime where all the spins are randomized and spin exchange is highly suppressed , this is just the probability of finding neighboring sites of the spin chain housing oppositely aligned spins , @xmath285 from which it follows that @xmath286 the remaining correlations over the charge degrees of freedom can be computed using the expansions and the identity @xmath287/2}$ ] valid when @xmath288 $ ] commutes with both @xmath289 and @xmath290 .
the final result is @xmath291 where @xmath292 is a short distance cut - off that ensures ultra violet convergence of the integrals over the charge fluctuations in the wire . up to unimportant phase factors and other overall multiplicative constants ,
the josephson current in the wire is then @xmath293 where the integral has been cut off at times @xmath294 , the time it takes for charge to propagate between the two ends of the sill .
compared to the ll case , the most important difference is that the critical current scales as an exponential of the length , @xmath295 rather than a power law @xmath296 in the ll regime.@xcite a measurement of the length dependence of the critical current then serves as clear signature of spin - incoherent physics ( or lack thereof ) , provided one has some knowledge of the density to infer @xmath297 .
we also note that the superconducting phase difference @xmath93 appears as an argument to the @xmath298 function , rather than the saw - tooth form we found for andreev contacts.@xcite the @xmath298 form follows from the assumption of tunneling contacts.@xcite    _ bulk contacted sill when @xmath258 _ : we now consider the case of a  bulk - contacted " sill , shown in fig.[fig : sc_tunnel]b . for this case , we assume that the perturbations induced by the contacts are irrelevant , which is so for @xmath299.@xcite otherwise , if the perturbations induced by the contacts are relevant this case reduces to an  end - contacted " sill discussed above .    starting again with we note that translational symmetry implies that the spin correlations satisfy @xmath300 which allows us to shift the sums : @xmath301 in .
the summation over @xmath302 leads to a delta function setting @xmath303 and the summation over @xmath304 leads to a delta function setting @xmath305 . making the change of variables @xmath306 and @xmath307 one can do the integration over @xmath308 which sets @xmath309 .
the result is @xmath310-[\phi(0,\tau)-\phi(l,0)]\right\ } } \rangle .
\label{eq : pi_bulk_gen}\end{aligned}\ ] ] the expression is valid quite generally under the assumption that the contacts are irrelevant perturbations and do not lead to the  end contacted " result described earlier .
if we now specialize to the spin - incoherent case we have @xmath311 substituting into and integrating over momentum , @xmath312\nonumber \\
\times e^{i2\left\{\left [ \theta(0,\tau)-\theta(l,0)\right]-[\phi(0,\tau)-\phi(l,0)]\right\}}\rangle .
\end{aligned}\ ] ] finally recalling that for @xmath313 we can replace the sum over @xmath40 by an integral,@xcite we find @xmath314 \ln2 }   e^{i2\left\{\left [ \theta(0,\tau)-\theta(l,0)\right]-[\phi(0,\tau)-\phi(l,0)]\right\ } } \rangle , \nonumber \\\end{aligned}\ ] ] which can be evaluated ( up to phase factors ) to give @xmath315 where @xmath316.$ ] the josephson current in the wire is then @xmath317 where we have again cut off the integral at times @xmath318 .
evidently , the critical current again scales as an exponential of the length , @xmath319 rather than a power law@xcite @xmath296 in the ll regime .
this luttinger liquid result corrects the result originally obtained by fazio _
et al._@xcite by taking into account the proximity effect .
comparing the critical currents of the end contacted case with the bulk contacted case , one sees that the length dependence is not materially different : they only differ in the inconsequential power law that multiplies the exponential .
_ end contacted sill when @xmath259 _ : to compute the josephson current we must evaluate the single - particle green s functions that appear in . since @xmath103 is a property of the sc which we assume is tunnel contacted to the sill , it is independent of the properties of the sill .
we will further assume @xmath320 .
the opposite limit of the wire length being shorter than the inter - particle spacing is not well motivated physically .
the single - particle green s functions appearing in can readily be evaluated following the method of appendix  [ app : spin_current ] .
this gives @xmath321}\rangle , \nonumber \\ & = & \frac{1}{2\pi \alpha_c } e^{-\frac{k_f^h l}{\pi}\ln 2 }    \left(\frac{\eta}{1+e^{-\eta}e^{-v_\rho \pi \tau / l}}\right)^{\frac{1}{4k_\rho}},\end{aligned}\ ] ] which immediately leads to the josephson current @xmath322 ^ 2 , \label{eq : endjoscurntshort}\end{aligned}\ ] ] with a critical current @xmath323 .
_ bulk contacted sill when @xmath259 _ :
for a bulk contacted sc in the short wire limit we again apply the formula where the green s functions are those appropriate for an infinite system , computed earlier in the literature in ref .  [ ] .
the result is @xmath324 ^ 2 , \label{eq : bulkjoscurntshort}\end{aligned}\ ] ] which evidently has the critical current @xmath325 , where @xmath326 $ ] .    comparing the results for bulk and end contacted , as well as short verses long wires
, we find the most important result is how the josephson critical current scales with the length of the wire .
suppressing the relatively unimportant power law decay that multiplies the dominant exponential decay we find @xmath327 so that the decay of critical current with the length of the sill is twice as fast as for a short wire .
this originates in the fact that for a short wire electrons propagate independently so the spin incoherence affects each electron independently , rather than a single coherent pair as in the case of the long wire .
it is also worth emphasizing that the assumption of tunneling contacts always results in a josephson current that is proportional to the product of the bare conductances of the two contacts with a sinusoidal dependence on the phase difference of the superconducting order parameter , @xmath328
. therefore it is the length dependence of the critical current that provides a smoking gun signature of spin - incoherence in the dc josephson effect with tunnel contacts . for andreev ( adiabatic ) contacts , the length dependence is identical to those of a ll or non - interacting one dimensional system.@xcite instead , the sill physics is manifest clearly in the flux dependence of the josephson current which takes on a saw - tooth form of half the usual period.@xcite      the ac josephson effect occurs when there is a finite voltage @xmath329 across the sc - sill - sc system .
the josephson phase acquires a time dependence @xmath330 leading to a josephson current that oscillates in time .
a sub - gap dissipative current is also induced , but at small voltages this can be estimated to be small.@xcite as in the case of the dc josephson effect , the qualitative features of the ac josephson effect also depend on whether the wire is short or long , as defined in the dc case .
_ long wire case _ :
we have seen that the qualitative features of the dc josephson current in the spin - incoherent regime for bulk and end contacts are not very different .
that is also true of the ac josephson effect .
the ac josephson current can be computed from@xcite @xmath331,\end{aligned}\ ] ] where @xmath332 is the cooper pair propagator , eq .
, evaluated at real times .
the ac josephson current can be decomposed into sinusoidal and cosinusoidal components , @xmath333,\ ] ] where @xmath334,\;\;\;\;\\ \label{eq : j_c } j_{c}&=&{\rm re}\left[\int_0^\infty dt ' \left(e^{-i evt ' } \pi(t ' ) - e^{+i evt ' } \pi(t')\right)\right].\end{aligned}\ ] ] in the luttinger liquid regime @xmath335 and @xmath336 oscillate with voltage across the junction , as does the amplitude @xmath337.@xcite the frequency of the oscillations of @xmath338 with voltage depends on the relative spin and charge velocities .
one period occurs when the spin and charge parts differ by @xmath339 , that is when @xmath340.@xcite in the spin - incoherent regime , we have @xmath203 , implying vanishing voltages will lead to oscillations and they may cease to be observed . the lack of amplitude oscillations will persist into the spin - incoherent regime , which has effectively only one velocity , the charge velocity . on the other hand , if @xmath205 is not too different from @xmath341 ( say , @xmath342 ) , then one can expect to find a temperature dependence of the voltage oscillations that reveals spin - incoherent physics in a way analogous to coulomb drag@xcite or charge fluctuation noise:@xcite at temperatures below the spin energy , there will be amplitude oscillations in the josephson current as a function of voltage , while for temperatures above the spin energy there will be no such oscillations because the spin mode effectively does not propagate .
this  washing out " of the high frequency ( because of the ratio @xmath343 ) oscillations with temperature is the signature of spin - incoherent physics in the ac josephson effect .
_ short wire case _ : for a short wire , the voltage dependence is independent of the properties of the 1-d system,@xcite be it a ll or a sill . in this case
the ac josephson current is @xmath344 where @xmath345 is the zero voltage critical current for the short wire given in sec .
[ sec : dc - josephson ] , and @xmath346 is an elliptical integral . recall that for a short wire , @xmath347 .
thus , the ac josephson effect is only effective at revealing spin - charge separation in the long wire limit , and there are no new spin - incoherent features that appear relative to the ll aside from the length dependence of the dc critical current that enters in .
in this paper we have touched on what we believe are some of the most easily observed consequences of spin - incoherent behavior in ferromagnet / spin - incoherent luttinger liquid and superconductor / spin - incoherent luttinger liquid junctions . for the case of fm - sill junctions , we computed the spin current pumped into the sill as a result of magnetization dynamics and the effect of spin accumulation in the sill on the parameters governing the magnetization dynamics .
we found that for interactions in the sill with @xmath348 the temperature dependence of the spin current and magnetization dynamics is _ qualitatively different _ from the case of a ll and should thus be observable experimentally .
some of the key differences between fm - sill and fm - ll systems are summarized in table  [ tb : ll_sill ] .
the crossover from fm - ll to fm - sill in the exponent @xmath349 governing the temperature dependence of several key quantities is shown in fig .
[ fig : delta ] .    in the case of sc - sill junctions
our results greatly extend those obtained earlier by us.@xcite in that earlier work we were concerned only with the case of adiabatic ( andreev ) contact of the spin - incoherent luttinger liquid to the superconductor . here
we have developed those results further and also discussed the opposite limit of tunnel contacts to the sc . in the tunneling limit
we have computed the ac and dc josephson response in the geometries shown in fig  [ fig : sc_tunnel ] .
we find that in contrast to the case of adiabatic contacts , the tunnel contacts lead to a josephson critical current that is exponentially suppressed with the length of the sill region .
this difference arises from the fact that in the andreev limit the dc josephson effect is determined solely by the zero modes of the hamiltonian , while for tunnel contacts the non - zero modes generally dominate the response .
these non - zero modes enter because in the tunneling process an electron is created locally near the tunnel barrier and thus requires many wavevectors to build its wavepacket .    with the aim of providing a general discussion of junctions of ferromagnets or superconductors with a strongly interacting one dimensional system ,
we have couched many of our calculations in the recently developed scheme for bosonizing strongly correlated electron systems.@xcite this formalism is valid for arbitrary temperatures with respect to the spin energy @xmath12 , but requires that the temperatures remain small compared to the charge energy @xmath13 . in order to adapt that formalism to the systems we discussed here , we extended those results to express open and andreev boundary conditions in a language valid for arbitrary temperatures with respect to the spin sector . through several examples , we showed how various correlations functions could be evaluated and verified that in the spin - incoherent regime the results properly reduce to the results obtained using the world - line picture.@xcite    one of our primary motivations for using the language of ref .
[ ] is to help provide a clear starting point for numerical studies of the spin - incoherent regime and the many interesting ( and probably experimentally relevant ) crossovers that occur between it and the luttinger liquid regime . within this formalism
the charge physics can be computed analytically via a standard bosonization scheme , but the spin sector must be addressed numerically for temperatures of order the spin energy @xmath12 .
the types of correlation functions that must be computed numerically are those that add and remove a site ( or multiple sites ) from a spin chain , such as @xmath350 , that appear in the evaluation of a single - particle ( or multi - particle ) greens function .
( see eq .  , for example .
) we would like to emphasize that numerical studies of strongly interacting one - dimensional electron systems with appreciable temperature compared to the spin energy is an entirely untouched area and is now ripe for investigation .
on the experimental side , the lack of a clear experimental  smoking gun " observation of the sill remains a key issue to be addressed .
however , there are mounting experimental indications we are not far away.@xcite our best numerical estimates@xcite suggest that many low density quantum wires sit right on the edge of the spin - incoherent regime and perhaps all that is needed is a focused experimental effort to search for its signatures , rather than any key technical breakthough .
one of our aims in this work it to highlight certain classes of hybrid structures where spin - incoherent physics should be observable .
finally , we would like to close with what we regard as some of the outstanding theoretical issues surrounding the spin - incoherent luttinger liquid .
perhaps the main one is the behavior on temperature scales @xmath22 that we already alluded to above .
related to this is a better understanding of the crossover between the luttinger liquid and the spin - incoherent luttinger liquid regimes .
both of these will likely require a numerical attack as there are no obvious analytical methods available to address them .
there is also 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 @xmath351 ?
the subject of noise in hybrid structures involving a sill will be discussed in a forthcoming work.@xcite    we thank leon balents for discussions on spin transport in one dimensional systems and oleg starykh for discussions on pair correlations in superconductors .
this work was supported by nsf grants phy05 - 51164 , dmr-0606489 , the lee a. dubridge foundation , and the welch foundation .
in this appendix we compute the correlator @xmath158 and its fourier transform @xmath149 using the general formalism for bosonization of strongly correlated electrons in one dimension developed by matveev , furusaki , and glazman who also applied it to the evaluation of the single particle green s function for an infinite system and its fourier transform.@xcite we summarized the main results of the bosonization scheme in sec .
[ sec : bosonization ] .
the calculation here is for a semi - infinite , or finite system with open boundary conditions , so is different in detail from what has been discussed in refs .
[ ] , but the basic elements of the bosonization are the same .
this appendix is meant to illustrate clearly in a specific example which parts of correlation functions can be computed analytically at finite temperatures and which pieces in general require methods yet to be developed , or a numerical attack . as was already discussed in sec .
[ sec : bosonization ] the chief difficulty lies in computing the correlations at arbitrary temperatures in the spin sector .
it is hoped that the details given here will provide a good starting point for those skilled in numerics to enter the study of strongly interacting one dimensional systems where there is currently _ no _ quantitative understanding of the regime @xmath352 .
given the typical values of @xmath353 present@xcite in quantum wires , this
intermediate " temperature regime may turn out to be the most relevant experimentally .    in the evaluation of the spin current ( pumped from a ferromagnet into an adjoined one dimensional system coupled via a tunnel junction ) , the correlation function @xmath354\rangle\ ] ]
arises in lowest ( second ) order perturbation theory .
it contains the difference of the product of single particle greens functions for both the fm and the interacting one dimensional system @xmath355 ,   \end{aligned}\ ] ] where the brackets @xmath356 denote a thermal average computed with the open boundary conditions described in sec .
[ subsec : open ] .
the correlators @xmath357 and @xmath358 can be computed by standard bosonization methods.@xcite the result is @xmath359 where @xmath360 is a high energy cut off .
our real objects of interest here are the boundary greens functions @xmath361 and @xmath362 evaluated at arbitrary temperature with respect to @xmath12 , but small compared to @xmath13 . to evaluate these we make use of expression for the electron operator .
straight forward substitution gives @xmath363 } z_{l_1,s } z^\dagger_{l_2,s } e^{-i[(1+\frac{q_2}{\pi})\theta(0)-\phi(0)]}\rangle.\end{aligned}\ ] ] open boundary conditions at @xmath64 prevent fluctuations of the @xmath47 field which sets it to a constant , which we take to be 0 . with @xmath47 set to a constant
, the integrals over momentum can be done trivially resulting in @xmath349-functions for @xmath364 which kill those sums and selects @xmath281 .
the final result is remarkably simple @xmath365}\rangle.\ ] ] similar manipulations yield @xmath366}\rangle,\ ] ] where we have made explicit the time dependence in the spin correlators involving @xmath37 and @xmath367 .
the general theoretical challenge is then to compute those correlation functions which involve adding and removing a site from the spin chain , which is an interacting many - body problem analogous to fermi edge physics@xcite only it involves knowledge of the hamiltonian on potentially many energy scales .
however , for @xmath20 the time evolution of the spin degrees of freedom are very slow compared to the charge degrees of freedom , regardless of temperature , and can be neglected for times @xmath368 .
note that in the limit @xmath369 , the spin dynamics can always be neglected .
furthermore , if @xmath370 , then the spins become essentially non - interacting , and therefore non - dynamical . then the correlation functions simplify considerably . in the spin - incoherent regime
one has@xcite @xmath371 and @xmath372 independent of time which can be straightforwardly generalized to include an applied external magnetic field.@xcite the result is @xmath373 where @xmath374 is the probability of having spin projection @xmath33 .
it takes the values @xmath375 , where @xmath376 is the zeeman energy . note that the correlation functions @xmath377 and @xmath378 in the spin - incoherent regime are identical to their values in the infinite system.@xcite this is because in the spin - incoherent regime the spins are non - dynamical , so the open boundary conditions which tend to suppress fluctuations have essentially no effect on the spins which are rendered non - dynamical by @xmath370 .
finally , we are left to evaluate @xmath379}\rangle$ ] , which is identical to @xmath357 and @xmath358 only with the exponent changed @xmath380}\rangle = \frac{(-i)^{1/2k_\rho}}{2\pi \alpha_c}\left[\frac{\left(\pi k_b t/\epsilon'_c\right)}{\sinh\left ( \frac{\pi k_b t t}{\hbar}\right)}\right]^{1/2k_\rho},\ ] ] where @xmath381 is a high energy cut off for the charge sector in the sill .    combining the results , , , and
, we find @xmath382^{\delta^{sill}+2},\ ] ] where @xmath383 is given in table  [ tb : ll_sill ] , and @xmath384 is an effective high energy cut off given by @xmath385^{2k_\rho/(2k_\rho+1)}$ ] . for making the following formulas more compact , we define @xmath386 .
we are now ready to compute the fourier transform of , @xmath387 . making the change of variables
@xmath388 , we have @xmath389^{\tilde \delta}}.\ ] ] the integral is standard @xmath390^{\tilde \delta } } = 2^{\tilde \delta-1}\gamma(1-\tilde \delta)\frac{\gamma\left ( { \tilde \delta \over 2}-i\frac{\hbar \omega}{2\pi k_b t}\right)}{\gamma\left(1-{\tilde \delta \over 2}-i\frac{\hbar \omega}{2\pi k_b t}\right)}.\ ] ] the result of the integral can be transformed to a more convenient form using the following identities for complex number @xmath391 : @xmath392 and @xmath393 .
this gives @xmath394 { 1 \over \pi},\end{aligned}\ ] ] where the sin can be expanded to pick out the real and imaginary parts .
selecting the imaginary part of @xmath149 we find @xmath395\propto \frac{1}{(2\pi \alpha_c)^2}\left(\frac{\hbar}{\pi k_b t}\right ) \left(\frac{\pi
k_b t}{\tilde \epsilon_c}\right)^{\tilde\delta}\nonumber \\
\times \sinh\left(\frac{\hbar \omega } { 2\pi k_b t}\right)\bigg|\gamma\left({\tilde \delta \over 2}+i\frac{\hbar \omega}{2\pi k_b t}\right)\bigg|^2,\end{aligned}\ ] ] which for @xmath165 gives @xmath396\propto \hbar \omega ( k_b t)^{\delta^{sill}},\ ] ] in agreement with our , and the general analytical form earlier obtained by bena and balents@xcite up to the exponent describing the temperature dependence which is changed in the spin - incoherent regime , as illustrated in fig  [ fig : delta ] .    the calculation above is valid provided the temperature is larger than the level spacing of the finite length wire : @xmath397 . for @xmath398
the finite length results of ref .
[ ] can readily be generalized to finite temperatures .
no simple power laws emerge , but rather a more complicated dependence involving hyperbolic trig functions .
in this appendix we show that the bosonization scheme introduced in sec .
[ sec : bosonization ] recovers our earlier results for the single - particle greens function , the tunneling density of states , and the pair correlations.@xcite      we are interested in computing the fourier transform of the single - particle green s function @xmath399 where here @xmath400 and @xmath401 measure the distance from the boundary with the sc . taking @xmath402 and using the general formula we have @xmath403-\phi(x , t)]}\nonumber \\
\times   e^{-i[(1+\frac{q_2}{\pi})[k_f^hx+\theta(x,0)]-\phi(x,0)]}\rangle,\;\;\;\end{aligned}\ ] ] where the @xmath47 and @xmath48 fields have the expansion given in .
we take the limit @xmath85 as we did in ref .
[ ] so that the zero modes play no role . since the holon sector is described the the gaussian theory we can make use of the identity @xmath404 for operator @xmath289 and evaluate the correlation functions at zero temperature with respect to the charge energy but infinite temperature with respect to the spin energy : @xmath405 , that is we take both the spin energy and the temperature to zero but alway maintain @xmath406.@xcite in the spin - incoherent regime , the spin correlations are translationally invariant so we have @xmath407 . making a change of variables @xmath408 , the summation over the sites of the spin chain results in a delta function , @xmath409 , which immediately kills one of the momentum integrals and sets @xmath410 .
we re - label the remaining momentum variable @xmath91 , @xmath411 }   e^{-i[(1+\frac{q}{\pi})\theta(x,0)-\phi(x,0)]}\rangle.\;\;\;\end{aligned}\ ] ] next , we define a variable @xmath412 and re - express the integration in terms of this variable , @xmath413 }   e^{-i[y\theta(x,0)-\phi(x,0)]}\rangle.\end{aligned}\ ] ] finally , noting that @xmath414)=\int _ { -\infty}^{\infty } \frac{dy}{2\pi } e^{-i[\pi l - ( \theta(x , t)-\theta(x,0))]y}$ ] we see that eq .   correctly reproduces eq.(4 ) of ref .
[ ] after we recall the relations @xmath415 and @xmath416.@xcite the identical result therefore follows for the greens function and the tunneling density of states derived from it .
here we compute the pair correlation , @xmath417 , a distance @xmath400 from the superconductor boundary . by direct substitution of
and we have @xmath418-\phi(x)]}\nonumber \\ \times e^{-i\phi(x ) } z_{l_2,\downarrow } e^{i[(1+\frac{q_2}{\pi})[k_f^hx+\theta(x)]-\phi(x)]}\rangle.\end{aligned}\ ] ] at the boundary , we showed in sec .
[ subsec : andreev ] that andreev reflection implies that @xmath419 . since @xmath420 at the boundary we may fluctuate , we expect also that the result @xmath419 approximately holds near ( within @xmath103 ) of the boundary
. hence , @xmath421}\nonumber \\ \times e^{-i\phi(x ) } z^\dagger_{l_2,\uparrow } e^{i[(1+\frac{q_2}{\pi})[k_f^hx+\theta(x)]}\rangle.\end{aligned}\ ] ] we now assume that we are in the spin - incoherent regime where the spin correlations are translationally invariant : @xmath422 . then repeating the steps used above to compute the single particle greens function ( changing the summation variable @xmath408 , killing the momentum integral with @xmath349-function ) , we find @xmath423 } e^{-i2\phi(x ) } \rangle.\end{aligned}\ ] ] we then make use of the integral representation of the @xmath349-function , @xmath424\right ) } = \delta \left(l-\frac{2}{\pi}[k_f^hx+\theta(x)]\right)$ ] to obtain @xmath425\right)e^{-i2\phi(x ) } \rangle,\ ] ] which for @xmath426 can be evaluated by taking the discrete sum to an integral : @xmath427 . doing so we recover the result ( 9 ) of ref .
[ ] after again recalling the relations @xmath415 and @xmath428.@xcite