recently there was much interest in rectification in nanoscale systems .
@xcite nonlinear mesoscopic transport exhibits much interesting physics such as a peculiar magnetic field dependence of the current @xcite and negative differential resistance for the rectification current at low voltages.@xcite another motivation for the investigation of mesoscopic diode or ratchet @xcite effect comes from possible practical applications in nanoelectronics and energy conversion .
following the pioneering paper by christen and bttiker @xcite most attention has focused on a simpler case of a fermi - liquid conductor . at the same time , the diode effect requires a combination of spacial asymmetry and strong electron interactions in the conductor .
hence , one may expect a stronger ratchet current in strongly interacting luttinger liquid systems .
this expectation has been confirmed by a recent study of transport asymmetries in one - channel quantum wires.@xcite refs .
have focused on a one - channel luttinger liquid in a linear conductor in the presence of a single asymmetric scatterer .
this is the conceptually simplest situation giving rise to rectification . at the same time , changing geometry may increase asymmetry and hence the rectification current . in this paper
we consider an asymmetric setup based on a y - junction of three quantum wires with three impurities .
we show that a stronger diode effect can be achieved in such system than in a linear luttinger liquid and rectification is possible even in the case of symmetric point scatterers .
we focus on the simplest one - channel y - junctions .
more complicated luttinger liquid junctions , such as y - junctions of single - wall carbon nanotubes , are also of interest . in particular , it might be easier to make such junctions in a reproducible way .
y - junctions are among basic elements of electric circuits , however , a theoretical investigation of luttinger liquid y - junctions@xcite has begun only during the last decade . by now , there is a good understanding of linear conductance near various fixed points as well as tunneling density of states.@xcite in this paper we extend the previous research to the problem of transport asymmetries .
specifically , we consider a setup of the type shown in fig .
we assume that one of the three terminals is kept at zero voltage .
ac voltages with amplitudes @xmath2 and @xmath3 , @xmath4 , are applied to the remaining two terminals . in a general case ,
a dc - current is generated in each of the three wires in the junction , fig .
1 . the three currents are different but can be computed in a similar way .
we calculate such rectification dc currents for various types of y - junctions .
we focus on the limit of a low - frequency ac voltage bias . in order to determine the amplitude of the rectification effect in that limit
, it is sufficient to find the difference of the dc currents at the opposite dc bias voltages , i.e. , compare the current when the potentials at the terminals are time - independent and equal 0 , @xmath2 and @xmath3 with the current when the potentials are 0 , @xmath5 and @xmath6 ( cf .
refs .  ) .
this corresponds to a dc current generated by low - frequency square voltage waves .
we find that in some classes of y - junctions the rectification current is higher than in a linear wire with the same strength of the repulsive electron interaction at the same voltage bias . in a certain interval of low voltages
the rectification current exhibits a power - law dependence on the bias : @xmath7 .
the exponent @xmath8 can be negative .
clearly , such dependence with a negative exponent can not extend all the way to zero voltage as @xmath9 at @xmath10 .
hence , the rectification current reaches a maximum at a certain voltage .
we demonstrate that in junctions without time - reversal symmetry , the maximal rectification current can be comparable with the total current @xmath0 even for low voltages and modest interaction strength in the wires .
in particular , such rectification current can be achieved in the `` island setup '' , illustrated in fig .  1 , which can be experimentally realized in quantum hall systems @xcite .
for comparison , in a linear wire , the rectification current @xcite is always much lower than @xmath11 .
the diode effect in y - junctions is not as strong in the presence of time - reversal symmetry as in its absence . still , the maximal rectification current is higher than in a linear conductor .
the paper is organized as follows : first , after a brief qualitative discussion , we describe our setup and formulate a model .
we then review the properties of y - junctions and derive a general expression for the rectification current .
next , we apply the general formalism to the model of three weakly connected wires .
then we determine the leading contributions to the rectification current at different interaction strengths and matching conditions at the junction .
we discover several regimes with different voltage dependences of the current .
finally , we discuss how to build a junction with the maximal rectification current .
is applied to the first wire .
wires 2 and 3 are connected to the ground .
we calculate the dc current in wire 1 . one can view ab , cd and ef as edges of an integer quantum hall liquid .
this corresponds to an `` island junction''.@xcite , width=192 ]    we consider a junction of @xmath12 one - channel wires . in the subsequent sections we will specialize to the case of @xmath13 and assume that the system is spin - polarized or , equivalently , that charge carriers are spinless fermions . in all cases
we assume that long range coulomb interaction in the conductor is screened by the gates .
in such situation the applied voltage bias affects not only the current but also the charge density in the system .
rectification is possible , if a change of the voltage bias sign results in a change of not only the sign but also absolute value of the current .
the latter obviously requires left - right asymmetry in the system .
asymmetry may have two origins : asymmetry caused by the way how the bias is applied and geometric asymmetry due to , e.g. , an asymmetric scatterer.@xcite in particular , one can consider a situation in which several terminals are kept at zero electric potential and the potential of the other terminals changes between @xmath14 and @xmath5 .
the charge density is different for the opposite voltage signs .
this in turn affects the current and leads to rectification .
such density - driven rectification is possible in a broad range of situations including systems without electron interaction in the presence of time - reversal symmetry .
it has been investigated in the context of a carbon nanotube junction in the noninteracting - electrons approximation in ref .  .
at a high voltage bias a strong rectification effect was found . in this paper ,
on the other hand , we focus on the low - voltage regime and , in particular , the universal behavior at low bias . in that limit the density - driven rectification effect
is of little importance as it results in a small rectification current @xmath15 in a noninteracting system .
besides , the density - driven rectification requires a specific way to apply bias . if we are interested , for example , in the transformation of incoming electromagnetic radiation into a dc current then clearly the potential oscillates in all terminals .
thus , in this paper we focus on the rectification mechanism due to geometric asymmetry .
y - junctions have a `` built - in '' geometric asymmetry .
indeed , let us focus on the current @xmath16 in one of the three wires connected by the junction .
the current @xmath17 equals the sum of the currents in the other two wires .
thus , the same current enters the junction through two wires and leaves through one only .
this means that the two sides of the junction are not equivalent since they correspond to 1 and 2 wires .
thus , we may expect that geometric mechanism of rectification applies to any y - junction .
however , as we demonstrate below , the geometric mechanism only works , if electron interaction is present or time - reversal symmetry is broken .
indeed , let us consider a time - reversal invariant system of non - interacting spinless fermions in an @xmath12-terminal mesoscopic junction .
we calculate the current between terminals @xmath18 and @xmath19 .
we compare the current in the situation when terminals @xmath18 are kept at the voltage @xmath2 and terminals @xmath20 are kept at zero bias ( case 1 ) with the current in the situation when terminals @xmath21 are kept at zero bias and terminals @xmath19 are kept at the voltage @xmath2 ( case 2 ) . for non - interacting fermions
the current reduces to the sum of the single particle contributions corresponding to each energy in the window @xmath22 .
we thus compare such contributions for two opposite voltage biases .
the wave function of an electron incoming from terminal @xmath23 with the momentum @xmath24 is @xmath25 , where @xmath26 is the coordinate in wire number @xmath27 and @xmath28 is a unitary scattering matrix .
the above notations for the wave function imply that the probability to find an electron with the wave function @xmath29 in wire number @xmath30 is determined by the outgoing wave and is proportional to @xmath31 . for @xmath32 ,
the probability is determined by both the incoming and outgoing waves .
the time - reversal symmetry implies that @xmath33 is also a solution of the schrdinger equation .
the solution is made of @xmath12 incoming waves and one outgoing wave .
it can be represented as a linear combination of waves @xmath34 with different @xmath27 .
hence , the outgoing wave satisfies the equation @xmath35 .
the absence of outgoing waves in the channels with numbers @xmath36 implies that @xmath37 .
hence , @xmath38 , i.e. , @xmath39 is a symmetric matrix .
now we can compare the contributions to the current from particles with the energy @xmath40 for two opposite signs of the bias .
electrons , incoming from different terminals , are not coherent . hence , in case 1 , @xmath41 . in case 2 , @xmath42 .
the currents are equal from the symmetry of the scattering matrix and hence the diode effect is absent .        in the presence of electron interaction , the rectification effect becomes possible .
this can be understood already from a simple model with two point scatterers of unequal strength in a linear wire ( fig .
similar rectification mechanisms operate in more complex junctions of luttinger liquids .
we will assume that long - range coulomb interactions are screened by the gates and thus the charge density depends on the voltage . in fig .
2a we consider the situation with the incoming current from the left .
backscattering off two impurities results in a `` staircase '' charge density profile ( we have averaged over friedel oscillations ) .
the `` staircase '' goes up as one moves from the right to the left . for the opposite voltage sign ,
the incoming current arrives from the right , fig .
we again have a `` staircase '' charge density profile but now the `` staircase '' goes down as one moves from the right to the left . in the absence of electron interactions
, the transmission coefficients must be the same in both cases and rectification is absent .
let us now consider the effect of electron interactions in a simple mean - field hartree picture .
since long - range coulomb interactions are screened , the relation between the charge density @xmath43 and the electric potential @xmath44 is local . assuming small charge density variations
, we thus find @xmath45 .
incoming electrons are scattered by a combined potential of the two impurities and the hartree potential @xmath44 .
the latter depends on the incoming charge density and hence the applied voltage .
hence , for opposite voltage signs , electrons feel different backscattering potentials .
this results in different transmission coefficients and thus rectification .
another mechanism of transport asymmetries comes from time - reversal symmetry breaking .
the latter is possible in the presence of a magnetic field .
we illustrate asymmetric transport in the absence of the time - reversal symmetry with a model system depicted in fig .
while it is not a realistic model of any mesoscopic conductor , it is very simple and exhibits strong transport asymmetries .
we assume no electron interactions in the model .
we consider a box with four holes in its sides .
each hole corresponds to a terminal in a more realistic description of a junction .
we assume that all charge carriers , entering the box through the holes , have exactly the same speed , perpendicular to the wall in which the hole is made . the speed @xmath46 , where @xmath2 plays a role , similar to the electrostatic potential of a terminal in a realistic junction . at @xmath10 the speed @xmath47 and
hence no current is injected . in a y - junction ,
charge density , injected from each terminal , depends on its electrostatic potential .
a comparison between opposite voltage signs in a realistic junction corresponds in our model to the comparison between the situations with the charge injected through the right and left holes ( figs .
3a , c ) and through the upper and lower holes ( figs .
3b , d ) . in the absence of a magnetic field ,
the current through each hole is zero in both situations : currents of the particles , injected through the opposite holes , cancel as shown in figs .
let us now turn on a magnetic field in shaded areas ( figs .
if the charge carriers are injected from the upper and lower holes then the current remains zero even in the presence of a magnetic field ( fig .
3d ) . indeed
, the carriers never enter the region with the field .
however , if the carriers arrive from the left and right holes then they are deflected by the field ( fig .
3c ) and the electric currents are nonzero in all four holes .
this illustrates how magnetic field can result in transport asymmetry .
the above arguments make it plausible that the rectification effect is the strongest if both strong electron interaction is present and time - reversal symmetry is broken .
this is confirmed by our calculations below for y - junctions of quantum wires with and without time - reversal symmetry .
in this section we formulate our model and review basic properties of y - junctions .
we consider a y - junction with the action    @xmath48,\ ] ]    where @xmath49 are the actions of three uniform wires and @xmath50 describe three impurities in the wires close to the junction point @xmath51 .
the above action must be supplied by matching conditions for three wires .
they will be discussed below .
the action of a uniform wire is given by the equation    @xmath52,\end{aligned}\ ] ]    where @xmath53 are the operators of incoming and outgoing chiral electron fields , @xmath54 is the effective fermi velocity , @xmath55 defines the interaction strength .
we set @xmath56 and the electron charge @xmath57 in most of the following text .
we assume that the long range part of the coulomb force is screened by the gates .
hence , the electron density depends on the voltage bias .
the impurity hamiltonians @xmath58,\ ] ] where @xmath59 is the potential of the impurity in wire number @xmath24 , the impurity being located close to the junction at @xmath51 .
we only take the backscattering part of the impurity hamiltonian into account as the forward scattering terms do not affect our results .
we assume the same interaction strength and fermi velocity in each wire .    in order to treat the case of strong interaction
it is convenient to bosonize @xcite the action in terms of the chiral fields @xmath60 such that @xmath61 , where the opposite signs should be chosen for the in- and out - fields , the commutator @xmath62=\mp2\pi i\delta(y - x)\delta_{kl}\ ] ] and @xmath63 plays the role of the effective fermi momentum and determines the average charge density in the wires .
@xmath64 are klein factors , necessary to ensure the proper fermi commutation relations .
the local densities of in- and out - moving particles in point @xmath65 are @xmath66 .
the action now assumes the form @xmath67;\end{aligned}\ ] ] @xmath68)+h.c.,\ ] ] where @xmath69 , with real @xmath70 and @xmath71 , are of the order of the fourier components of the asymmetric potential , @xmath72 .
note that @xmath71 can be nonzero even for a symmetric potential @xmath73 in contrast to the situation considered in refs . .
for example , for @xmath74 , @xmath75 .
the above expression for the backscattering operators @xmath50 includes multi - particle backscattering processes.@xcite such multi - particle contributions to the action are inevitably generated under the action of the renormalization group by the interplay of a short - range coulomb interaction and impurity potential ( see ref .   for a discussion ) .
all prefactors @xmath70 have the same order of magnitude and are proportional to @xmath73 .
we assume that backscattering amplitudes @xmath70 have the same order of magnitude in all three wires .
the use of the fields @xmath76 at @xmath51 in the backscattering operators is justified , if the distance from the impurities to the junction is lower than the scale @xmath77 set by the voltage bias .
we assume that @xmath70 are sufficiently small so that a perturbative expansion in powers of @xmath70 can be developed for the calculation of the current .
the conditions on @xmath70 will be formulated below .    following the standard notation conventions , the hamiltonian , corresponding to eq .
( [ bos ] ) , can be written as    @xmath78,\ ] ]    where the dimensionless interaction strength @xmath79 and @xmath80 .
the above action alone is not enough to describe the system and matching conditions at the junction are necessary . as we will see , matching conditions are subject to several restrictions .
the most general matching condition has the form @xmath81 , where @xmath82 is an arbitrary function .
below we will focus on the situation in which the system is close to a fixed point . then
only most relevant operators should be kept in the matching conditions and hence all derivatives of the bose - fields @xmath83 can be neglected .
next , we note that the shift of any field @xmath84 by a constant @xmath85 is a gauge transformation that does not change the physics of the system .
hence , the matching conditions must be invariant with respect to such gauge transformations .
this singles out the boundary conditions of the form @xcite @xmath86 where @xmath87 is a matrix with real matrix elements . with this matching condition ,
it is easy to solve the equations of motion for the chiral fields @xmath88 in the absence of the interaction ( @xmath89 ) and impurities , for arbitrary initial conditions for fields @xmath84 . substituting the solution into the commutation relations ( [ comm ] ) for in- and outgoing fields
, one finds that they are compatible , provided that the matrix @xmath87 is orthogonal .
the same condition for the matching matrix must hold for a general case with an arbitrary interaction strength @xcite since the problem can be reduced to a model of non - interacting bosons by diagonalizing eq .
( [ ham ] ) in each wire ( see the next section ) .
thus , we established that @xmath87 is a real orthogonal matrix .
we next demonstrate that the sum of the elements of each of its rows and columns equals one.@xcite it is again convenient to consider the situation without electron interaction and tunneling . in that case
the currents of in- and outgoing electrons in each wire are @xmath90 .
charge conservation implies that the sum of the currents in all wires is zero at @xmath51 .
taking into account eq .
( [ cond ] ) one finds that @xmath91 .
multiplying eq .
( [ cond ] ) by @xmath92 and repeating the same argument one also finds that @xmath93 .
again , electron interactions do not affect the above result @xcite .
the origin of the matching conditions can be understood if one notes that they can be imposed by adding to the hamiltonian a term of the form    @xmath94),\ ] ]    where @xmath95 is a large constant .
each cosine tends to keep its argument at zero .
this can be achieved simultaneously for each cosine only if their arguments commute .
this happens for a unitary @xmath87 .
the cosines describe tunneling between different wires .
thus , charge conservation implies @xmath93 .
this provides an alternative derivation of the matching conditions .    since we plan to separately investigate time - reversal invariant and non - invariant systems , we next need to determine what matching conditions satisfy the time - reversal symmetry
this is easy as the time - reversal transformation corresponds to the change of variables @xmath96 . from eq .
( [ cond ] ) we then see that in time - reversal - invariant systems @xmath97 . taking into account that @xmath87 is orthogonal
we conclude that in time - reversal invariant systems it is also symmetric.@xcite the action ( [ l ] ) is always time - reversal invariant .
thus , the behavior of the whole junction with respect to time reversal is determined solely by the symmetry of the matching matrix @xmath87 .    at this point
we are in the position to give a full classification of fixed - point matching matrices.@xcite for time - reversal invariant systems , it is convenient to use a parametrization from ref .  .
we discover three possibilities with time reversal symmetry : @xmath98    @xmath99    @xmath100    where @xmath101 , @xmath102 and @xmath103
. eq .  ( [ identity ] ) corresponds to three disconnected wires .
the cases of @xmath104 correspond to a junction of two wires and a detached third wire .
the physics is the same as for one linear wire and will not be discussed below .    in the absence of time - reversal symmetry we use the parametrization from ref .  :
@xmath105 where @xmath106 , @xmath107 and @xmath108 .
the case of @xmath109 , @xmath110 corresponds to the island setup , depicted in fig .
1 . note that the matching matrix ( [ 122 ] ) has the form ( [ no - time ] ) with @xmath111 so we will consider it with the case of no time - reversal symmetry .    in general , matching matrices contain negative matrix elements .
this means that the incoming current in one wire may suppress outgoing current in that or other wires .
such situation is possible in the presence of andreev scattering , if , for example , a part of the junction is superconducting .
negative matrix elements were also predicted in y - shaped beam splitters for cold bosonic atoms @xcite .
note that for noninteracting electrons no negative matrix elements are allowed @xcite since in the absence of interaction the elements of the matrix @xmath87 reduce to the squares of the absolute values of the elements of the scattering matrix ( cf .
in order to calculate the current we need to include a voltage bias .
we model fermi - liquid leads by assuming that electron interaction in the wires is zero at large distances @xmath65 from the junction.@xcite we will assume that leads 2 and 3 are kept at a zero voltage .
the bias voltage @xmath112 is applied to lead 1 .
the results do not change in a more general model with the potential @xmath112 applied to lead 1 , the potential @xmath113 , @xmath4 , applied to lead 2 , and a zero potential applied to lead 3 .
we will use the language of zero @xmath114 since it is simpler .
our approach can be easily extended to a general @xmath114 .
the calculations will be based on a renormalization group approach with the voltage @xmath2 playing the role of the infrared cutoff .
thus , we will assume that the temperature @xmath115 .
we would like to emphasize that the role of the voltage does not reduce solely to that of a cut - off ; otherwise the density - driven rectification would be lost ( see a discussion after eq .
( [ int - repr ] ) ) .
it is well known that in luttinger liquids a voltage bias can play a more prominent role than just a cut - off ( see , e.g. , refs . ) . at the same time , the leading contribution to the total current at zero temperature can be estimated by simply setting the renormalization group cutoff to the value of the voltage @xcite .
below we find that the rectification current is comparable with the total current in certain regimes in the absence of the time - reversal symmetry .
this certainly means that in those regimes the rectification current can be computed by assuming that the voltage plays the role of a cutoff only ( we however do not make such an assumption ) .    for a general @xmath114
, we do not expect the current to depend significantly on the temperature at @xmath115 .
such dependence can emerge at particular values of @xmath114 .
in particular , in the island setup ( fig .
1 ) , this happens at @xmath116 .
this can be seen from the keldysh perturbation theory @xcite .
the infrared cutoff in some integrals in the perturbation expansion is set by @xmath117 and in some others by @xmath118v)$ ] .
hence , at @xmath116 our approach applies only for @xmath119 in the island setup .
we want to find the current in the first wire . at low frequencies the current conserves
and hence it is sufficient to find the current in lead 1 .
it is given by the sum of the chiral incoming and outgoing currents @xmath120 .
the incoming current can be found from the landauer formula and is linear in voltage .
it does not contribute to the rectification current .
we thus focus on the outgoing current in lead 1 .
it can be found with a generalization of the approach of ref .  .
our approach is also related to that of ref . .
in what follows we set the temperature @xmath121 to simplify notations .
our method can be easily generalized to finite temperatures .
let us introduce an auxiliary field @xmath122 : @xmath123 the field @xmath124 satisfies a simple matching condition @xmath125 .
thus , the auxiliary field is a chiral field propagating through the junction . in the stationary regime
the average time - derivative of any operator is zero .
let us now consider the operator @xmath126 , where the integration extends between points taken in the noninteracting leads .
its meaning is the charge carried by the mode @xmath124 .
the equation @xmath127\rangle$ ] , where the hamiltonian @xmath128 describes the whole system including the leads , wires and scatterers , reduces to the following relation : @xmath129\rangle.\ ] ] the left hand side is the outgoing current we want to find .
the first term in the right hand side is linear in the voltage bias and can not contribute to the rectification current .
thus , we have to calculate only the second contribution to the right hand side .
in other words , our problem reduces to the calculation of the average backscattering current whose operator equals @xmath130 , where @xmath131 \nonumber\\ = & i[u_l(f^{k\dagger}_of^k_i)^n\exp(i\alpha_l+in\{\sum_p ( m^{-1})_{kp}\phi^p_o(0 ) -\phi^k_o(0)\ } ) \nonumber \\ & + h.c.,\hat o ] \nonumber \\
= & i n u_l ( f^{k\dagger}_of^k_i)^n \exp(i\alpha_l+in[\sum_p ( m^{-1})_{kp}\phi^p_o(0)-\phi^k_o(0 ) ] ) \nonumber\\ & \times [ ( m^{-1})_{k1}-\delta_{k1}]+ h.c . \ , .\end{aligned}\ ] ] note that for any @xmath23 the current operator expresses via _ all _ backscattering operators in all three wires .
in order to find the average current @xmath132 , we will assume that the backscattering operators @xmath50 are absent at @xmath133 and are then gradually turned on . without the operators
@xmath50 the system can be viewed as an equilibrium one in the ground state of an appropriate effective hamiltonian . in order to find it we introduce another auxiliary field with the structure , similar to @xmath124 : @xmath134 in the absence of the backscattering operators @xmath50 ,
the operator @xmath135 commutes with the hamiltonian .
it can be understood as the charge of a chiral mode propagating through the junction . in other words
, it is an additive integral of motion .
hence , the system can be described by a gibbs distribution with an appropriate thermodynamic potential conjugated to @xmath136 .
the physical meaning of that thermodynamic potential is the applied voltage bias @xmath2 . at zero temperature
, one finds that the system is in the ground state of the effective hamiltonian @xmath137 , where @xmath138 is the actual hamiltonian of the y - junction .
the current , i.e. , the average of the sum of the operators @xmath139 $ ] , can be now calculated with the keldysh technique @xcite .
it is convenient to apply the interaction representation @xmath140 .
the interaction representation introduces time - dependence into the operators @xmath50 and @xmath141 according to the rule : @xmath142)\rightarrow   \nonumber\\ & \exp(iv\hat at)\exp(in[\phi^k_i(0)-\sum_p m_{kp}\phi^p_i(0)])\exp(-iv\hat at ) \nonumber\\ & = \exp(in[\sum_p ( m^{-1})_{kp}\phi^p_o(0)-\phi^k_o(0)+vt(m_{k1}-\delta_{k1})]).\end{aligned}\ ] ] after such time - dependence is added into all backscattering operators , the contribution of the form @xmath143 can be removed from the action by a linear shift of all fields @xmath144 .
this does not mean , however , that the voltage would only enter the action through the time - dependence of @xmath50 . indeed ,
since the charge density depends on the voltage , the amplitudes @xmath145 of the operators @xmath50 may get corrections proportional to the small voltage bias .
this effect is discussed in ref . .
the current can now be found with a perturbative expansion in powers of @xmath145 from the standard expression @xmath146 where @xmath147 is the evolution operator from time @xmath148 to time @xmath149 and @xmath150 is the ground state of the effective hamiltonian @xmath151 . to complete the calculation
, we need the green functions determined by the quadratic part of the action ( [ lk ] ) and the matching conditions ( [ cond ] ) .
they have been found in ref .  .
a bogoliubov transformation of the form @xmath152 $ ] in each wire allows one to obtain free chiral fields @xmath153 with correct luttinger liquid commutation relations.@xcite the correlation function of the incoming fields at @xmath51 is given by the free particle relation @xmath154/\tau_c+\delta),\ ] ] where @xmath155 is the ultraviolet cutoff time of the order of the inverse bandwidth , @xmath156 is infinitesimal .
the matching conditions for the new fields have the form @xcite @xmath157 @xmath158^{-1}[(1+g)m-(1-g)].\ ] ] in the time - reversal invariant case @xmath159 . in the absence of the time - reversal symmetry , @xmath160 has the same general structure ( [ no - time ] ) as @xmath87 and satisfies the same set of constraints but the matrix elements are different : @xmath161 @xmath162    at this point
one can write an expression for the current as an expansion in powers of backscattering amplitudes @xmath163 . evaluation of the terms of that expansion is technically difficult and not very informative . indeed ,
neither the amplitudes nor the ultra - violet cutoff are known exactly . as a result
, it is only possible to estimate the order of magnitude of each contribution . at the same time
, such estimation can be performed without explicitly calculating integrals over the keldysh contour and will be sufficient to find the leading power dependence of the rectification current on the bias voltage at small @xmath163 and low voltages near each of the fixed points
. such estimation will be the focus of the remaining sections .
it is convenient to use a renormalization group point of view for the calculation of the current.@xcite in the renormalization group procedure , the coefficient @xmath163 in the hamiltonian and the operators @xmath141 scale as @xmath164 , where @xmath165 is the energy scale and @xmath166 the scaling dimension . at the energy scale @xmath167
the renormalization group procedure stops .
different perturbative contributions to the current can be estimated from the scattering theory.@xcite the current can be expressed as an infinite sum of the contributions , proportional to the products of different combinations of @xmath163 : @xmath168 strictly speaking , it is not enough to include only backscattering operators from the bare tunneling hamiltonian ( [ tun ] ) in the expansion ( [ current - scale ] ) .
all operators generated by the renormalization group must also be included .
their general structure is @xmath169 ) \nonumber\\ = & w_l\exp(i\alpha_l+i\sqrt{g}\sum_k\sum_p n_k[\tilde m_{kp}-\delta_{kp}]\tilde\phi^p_i).\end{aligned}\ ] ] for simplicity we use the same notation @xmath170 for the phases of @xmath171 and @xmath145 .
note that it is always sufficient to keep only two different values of @xmath24 in the sum in eq .
( [ gen ] ) , since @xmath172=0 $ ] .
this means that each possible operator @xmath173 can be generated from a product of two operators @xmath50 ( [ tun ] ) and hence @xmath174 .
since klein factors do not change the scaling dimensions of @xmath171 , we omit them in eq .
( [ gen ] ) .    to complete the calculation of the current we need to find @xmath166 . a straightforward calculation based on eq .
( [ correl ] ) yields for the operator ( [ gen ] ) with an arbitrary choice of @xmath175 : @xmath176 in the absence of the time - reversal symmetry the above expression greatly simplifies : @xmath177.\ ] ] in the low - voltage regime , the main contribution comes from the operators with the smallest @xmath166 .
one easily sees that it equals @xmath178 and is achieved , if two of the coefficients @xmath179 with the remaining one being @xmath180 .    in the time - reversal invariant system ,
the expression for the scaling dimensions is more complicated .
it can be simplified by setting @xmath181 and @xmath182 . using the matching matrix ( [ time - reverse ] ) , one finds : @xmath183 ^ 2}{b+c}=\frac{g}{1+\alpha+\alpha^2}(k\alpha - m)^2.\ ] ] we will assume for simplicity that @xmath184 is irrational .
we will thus avoid situations with @xmath185 at some @xmath24 and @xmath27 .
certainly , such situations correspond to the operators @xmath171 with no @xmath76-dependence .
such operators can not affect transport .
the renormalization group procedure only applies if all @xmath163 are small at the initial energy scale @xmath186 and remain small up to the scale @xmath2 . in the case without time - reversal symmetry , this condition can be easily expressed in terms of @xmath187 : @xmath188 the above equation assumes that @xmath189 , i.e. , @xmath163 is relevant . if @xmath190 then the only restriction on @xmath163 is @xmath191 .
what happens in the presence of the time - reversal symmetry will be addressed in subsequent sections .
the above discussion applies to all contributions to the current , both even and odd in the voltage bias .
we are interested in the even contribution , i.e. , the rectification current @xmath192/2 $ ] .
what terms of the perturbation expansion contribute to the rectification current ?
the answer to this question can be obtained from symmetry considerations .    in what follows
we will need the commutation relations for the klein factors in ( [ tun ] ) . using the commutation relations @xmath193 and
@xmath194=-i\pi { \rm sign}(x - y ) $ ] , an expression for @xmath195 in terms of @xmath196 and the baker - hausdorff formula , one finds @xmath197)f_i^k f_o^{k\dagger}.\ ] ]    next , some terms of the perturbative expansion ( [ current - scale ] ) are zero identically as only certain combinations of the vertex operators in @xmath50 and @xmath141 produce non - zero results after averaging with respect to the quadratic part of the action ( [ bos ] ) .
the condition is well known : @xmath198    let us now apply the above results to possible contributions of different orders of the perturbation theory to the rectification current .
strong limitations emerge for second order contributions .
indeed , eq .  ( [ nonzero ] ) implies that any second order contribution , proportional to an operator @xmath163 , must also contain its hermitian conjugated operator @xmath199 .
hence , the phase factors @xmath200 cancel each other and drop out from the expression for the current .
thus , we can just set @xmath201 .
let us now compare the currents at the bias voltages @xmath2 and @xmath5 . in the case with the bias @xmath5
, we make the change of variables @xmath202
. this transformation does not affect the form of the quadratic part of the action ( [ bos ] ) and the linear matching conditions .
let us also change the order of the klein factors in the tunneling and current operators using the commutation relation ( [ klein - comm ] ) and redefine @xmath203 .
since @xmath163 pairs up with @xmath204 in the nonzero second - order contribution to the current , one can easily see that the phase factors @xmath205)$ ] in the commutation relation ( [ klein - comm ] ) drop out . at the same time , if we omit the phase factors @xmath200 and @xmath205)$ ] from the action and current operators , we discover that the action assumes precisely the same form as for the voltage bias @xmath2 . on the other hand ,
the current operator @xmath141 ( [ backscattering ] ) changes its sign .
this means that the average of @xmath141 also changes its sign and hence the second order term of the order @xmath206 does not contribute to the rectification current .
the above argument assumes that the backscattering amplitudes @xmath163 are independent of the voltage bias . since the injected charge density depends on the voltage in our setup , @xmath163 can exhibit a weak linear dependence on the bias voltage .
this means that second order terms can in fact lead to a `` density - driven '' rectification effect but its amplitude is suppressed by an additional factor @xmath207 , where the ultra - violet cutoff @xmath155 is of the order of the inverse band width .
thus , the second order contribution to the ratchet current scales as @xmath208    beyond the second order , no additional general restrictions on possible contributions to the rectification current can be derived .
particular contributions to the ratchet current can disappear at particular values of the phases @xmath170 and interaction strength @xmath209 .
is applied to wire 1 , and bias @xmath210 to wire 2 .
wire 3 is connected to the ground .
we calculate the dc current in wire 1.,width=192 ]    so far we discussed general features of y - junctions .
we now determine the rectification current in different setups . as a warm - up exercise
we investigate the simplest situation of three almost disconnected wires , fig .
this situation corresponds to the matching matrix ( [ identity ] ) .
the rectification current is much weaker in this limit than for other matching conditions . at the same time
, the qualitative picture is quite similar .
there is no current at all in our original model with the action ( [ l ] ) and backscattering operators ( [ tk ] ) as it describes three disconnected conductors .
we thus modify the model in this section . instead of the backscattering operators ( [ tk ] ) we consider weak tunneling between the wires .
the tunneling operators have the form @xmath211 where @xmath212 , @xmath213 are klein factors and we use the convention @xmath214 .
the commutation relations for the klein factors are @xmath215 .
we assume that wires 3 is kept at zero voltage .
the chemical potential of the first wire is @xmath112 .
the second wire experiences the bias @xmath113 .
note a different meaning of the amplitudes @xmath145 and phases @xmath216 from the model with three weak scatterers discussed above . in eq .
( [ 3-tun ] ) , the phases @xmath216 describe the aharonov - bohm effect due to the magnetic field through the junction . on the other hand , in the model with weak scatterers , all information about the magnetic flux is contained in the matching conditions .
the voltage bias can be included in the action in the form of the time - dependence of the tunneling operators , similar to the discussion above .
we will not need explicit expressions below .
the operator of the current , tunneling between wire 1 and wires 2 and 3 , is @xmath217    the current can be estimated from a renormalization group procedure that stops at the scale @xmath167 . at that scale @xmath218 , where @xmath219 characterizes the interaction strength .
the second order contribution to the rectification current can now be easily found and scales as @xmath220 where @xmath221 .
the @xmath222-th order contribution to the current can not exceed @xmath223 , @xmath224 .
thus , the leading higher - order contribution corresponds to @xmath225 , @xmath226 interestingly , at low voltages and @xmath227 , the third order contribution exceeds the second order contribution .
the above discussion ignored the issue of the time - reversal symmetry .
depending on the presence or absence of the magnetic field through the junction , the system can have or have no time - reversal symmetry.@xcite it is instructive to investigate the effect of the symmetry breaking . in the presence of the time - reversal symmetry
the magnetic flux is zero and all @xmath228 .
indeed , the time - reversal transformation can be represented as @xmath229 . since @xmath230
, we conclude that @xmath231 at zero magnetic field . on the other hand , in the presence of the magnetic field @xmath232 .    the second order contribution to the rectification current
does not depend on the phases @xmath216 .
let us compare the third order contributions for different values of @xmath216 .
nonzero third order contributions to the current originate only from the product of all three operators @xmath50 or all three operators @xmath233 in the perturbative expansion of the keldysh expression for the current ( [ keldysh ] ) . the first contribution , @xmath234 , is proportional to @xmath235 , @xmath236 , and the second contribution is proportional to @xmath237 , @xmath238 . let us change the voltage sign .
simultaneously , we change @xmath239 , change the order of klein factors in all operators and redefine @xmath240 .
this transformation is equivalent to changing the sign of @xmath216 in all operators and simultaneously changing the overall sign of the tunneling hamiltonian ( [ 3-tun ] ) but not the sign of the current operator .
hence , the third order current @xmath241 .
we see that the current is an even function of the voltage in the presence of the time - reversal symmetry , i.e. , at @xmath242 .
the third order current is an odd function of the voltage and does not contribute to the rectification effect , if @xmath243 .
thus , for three weakly connected wires , magnetic field suppresses rectification .
we will see below that typically it has an opposite effect and rectification is stronger in the absence of the time - reversal symmetry .
note that in the above example , the ratchet current returns to its maximal value at @xmath244 .
let us finally compare our results with the case of a linear wire with a high asymmetric potential barrier or , equivalently , two weakly connected wires . in that case
the rectification current is given @xcite by eq .
( [ 3-i-2 ] ) .
thus , at @xmath227 the rectification effect is stronger in a y - junction than in a linear wire with the same interaction strength .
here we address matching matrices of the form ( [ no - time ] ) .
in particular , the results of this section are relevant for the island setup , fig .
our analysis also applies to the time - reversal - invariant situation with the matching matrix ( [ 122 ] ) .
below we find that at sufficiently low voltages the rectification current scales as eq .
( [ 3-no - time ] ) with @xmath245 being the impurity potential strength , @xmath2 the voltage , @xmath219 the luttinger liquid parameter , characterizing electrostatic interactions , @xmath155 the ultraviolet cutoff of the order of the inverse band width , and @xmath246 is determined by the matching conditions according to eq .
( [ tildea ] ) .
similar to the previous section we need to compare the leading second and third order contributions to the rectification current .
they are determined by the operators of the form @xmath247)$ ] with the minimal scaling dimension ( [ z - min - no - time ] ) .
the leading second order contribution scales as @xmath248 , where @xmath246 is given by eq .
( [ tildea ] ) and satisfies @xmath249 .
the leading third order contribution comes from the terms in the perturbation expansion proportional to the product of all three @xmath50 ( or all three @xmath233 ) .
it scales as    @xmath250    at small @xmath245 the second order contribution always exceeds higher - order contributions .
interestingly , however , at greater @xmath245 which are still within the region of validity of the perturbation expansion , the third order contribution is leading .
let us show that the third order contribution dominates .
we will focus on the largest @xmath163 accessible with the perturbation expansion : @xmath251 at @xmath252 when @xmath163 is irrelevant ; and @xmath253 at @xmath254 when @xmath163 is relevant . in all cases
the electron interaction is repulsive , i.e. , @xmath219 .
let us first consider the case of @xmath252 . according to eq .
( [ tildea ] ) , @xmath255 .
. the ratio of the second order contribution to the third order contribution scales as @xmath257 .
thus , the third order contribution dominates indeed .
the case of @xmath254 is also easy . in the limit of @xmath253
we find @xmath258 and @xmath259 .
interestingly , the third order current may become comparable to the total ac current @xmath0 through the junction at @xmath254 .
note also that the exponent @xmath260 in the voltage dependence of the dc current ( [ 3-no - time ] ) is negative at @xmath261 .
this corresponds to the dc current _ increase _ as the ac voltage decreases .
the above calculation applies for the voltage interval @xmath262 at @xmath252 and @xmath263}/\tau_c$ ] at @xmath254 .
the left inequality is dictated by the applicability of the luttinger liquid model .
the right inequality is determined by the validity of the perturbation theory .
one can not calculate the rectification current at lower voltages with the above perturbative approach .
however , it is obvious that @xmath9 at @xmath10 . from this
we conclude that there is a bump on the voltage dependence of the rectification current at @xmath261 : it grows as a function of the voltage at low voltages and achieves its maximal value of the order of @xmath11 at @xmath264}$ ] . in the last equation we omitted a power of @xmath155 as this
does not lead to a confusion .
finally , let us compare the result with the case of a linear wire with the same interaction strength @xmath209 in the presence of an asymmetric impurity . in that case
the maximal rectification current , calculated in the perturbative regime in ref . , scales as @xmath265 at @xmath266 ( we omit powers of @xmath155 ) .
the rectification current does not exceed @xmath267 at @xmath268 . in a y - junction , at @xmath269
the rectification current @xmath270 at @xmath271 is always greater than in a linear wire . on the other hand
, a negative @xmath272 implies @xmath273 .
the maximal rectification current corresponds then to @xmath274 and scales as @xmath275 .
thus , the maximal ratio of the rectification and total currents is always higher in a y - junction than in a linear wire in the luttinger liquid regime @xmath276 .
this case is more complicated than the situation without the symmetry .
the summary of the results is given in section ix.c .
we will need to take into account many different backscattering operators .
thus , it is important to classify them .
the classification relies on the matching conditions at the junction .
the matching matrix is given by equation ( [ time - reverse ] ) where the expressions for @xmath277 in the parametrization @xcite are written under the equation .
it is always possible to redefine @xmath184 in such a way that @xmath278 , where @xmath279 is an arbitrary transposition of @xmath280 . indeed ,
the change of the variable @xmath281 exchanges @xmath282 and @xmath283 ; @xmath284 exchanges @xmath282 and @xmath285 ; @xmath286 exchanges @xmath283 and @xmath285 . any other transposition is a superposition of the above three .
since @xmath87 is an orthogonal matrix , at least one of the elements must be negative or zero . otherwise , different rows can not be orthogonal . without loss of generality
we may assume that @xmath287 and hence @xmath288 .
this can always be achieved by renumbering the wires and redefining @xmath184 .
note that @xmath289 .
similarly , without loss of generality we may assume that @xmath290 .
hence , @xmath291 .
note that @xmath292 .
thus , the ranges of the parameters that we consider are : @xmath293 , @xmath294 , @xmath295 and @xmath296 .
as has already been discussed , we need to deal with two types of backscattering operators in the action and the current operator :    @xmath297 @xmath298)$ ] ;    @xmath299 @xmath300+in_3[\phi^3_o(x=0)-\phi^3_i(x=0)])$ ] .
note that we do not need to include a contribution of the form @xmath301 $ ] in the exponent in @xmath171 .
we will first concentrate on the operators of the first type and determine the leading contribution to the ratchet current in the absence of the operators of the form @xmath171 .
next , we will check what changes after @xmath171 are taken into account .
we summarize our findings in sec .
ix.c .
the scaling dimensions of the operators @xmath163 are simply    @xmath302    the lowest scaling dimension is @xmath303 with @xmath304 .
it sets the scale of the maximal @xmath245 in the perturbation theory , @xmath305 . strictly speaking we need to know the scaling behavior of the operators @xmath171 to determine the maximal allowed value of @xmath245 in the perturbative regime .
we will see in the next subsection that the operators @xmath171 do not change the above expression for the maximal @xmath245 .
we thus immediately find the leading second - order contribution to the rectification current @xmath306    what about the third and higher orders ?
it turns out that the leading contribution is third order and comes from the operators @xmath307)$ ] , @xmath308)$ ] and @xmath308)$ ] .
it scales as @xmath309 - 2}\ ] ]    to see that ( [ third - order - main ] ) is the leading contribution , we first note that @xmath310 , eq .
( [ third - order - main ] ) , exceeds @xmath311 at the maximal allowed @xmath305 . indeed , at such @xmath245 , @xmath312 while @xmath313}>(v\tau_c)^{1 + 3g\times 1/3}>(v\tau_c)^2 $ ] .
next , let us compare the contribution ( [ third - order - main ] ) with other higher order contributions to the current , @xmath314^{z_l - g[1-a]}\right).\ ] ] if an operator with the scaling dimension @xmath315 $ ] enters the above expression then ( [ third ] ) is smaller than ( [ third - order - main ] ) . indeed , any contribution to the rectification current that contains such operators satisfies the inequality @xmath316}=(v\tau_c)^{1+g[1 + 2\alpha]/[1+\alpha+\alpha^2]}<(v\tau_c)^{1 + 3g[1-a]}$ ] .
thus , only contributions with the operators @xmath317)$ ] and @xmath318)$ ] have to be considered .
if one of the contributing operators has @xmath319 or @xmath320 then it is obvious that the contribution can not exceed ( [ third - order - main ] ) .
we thus consider the case with all @xmath321 .
a nonzero contribution to the current requires @xmath322+(\sum_l n_{3l})[\phi^3_o(x=0)-\phi^3_i(x=0)]=0 $ ] . for a general transcendental @xmath184
this equality is satisfied only if @xmath323 .
thus , we consider contributions in which the operator @xmath324)$ ] enters the same number of times as the operator @xmath325 and the operator @xmath326)$ ] enters the same number of times as the operator @xmath327 .
hence , the phases @xmath328 drop out from the final answer .
let us now change the voltage sign , perform the transformation @xmath329 , change the order of the klein factors in each term and redefine @xmath330 . changing the order of the klein factors introduces complex conjugate phase factors into @xmath145 and @xmath331 .
hence , we can ignore both those phase factors and @xmath216 since neither affects the final result .
on the other hand , if we ignore the phase factors , we discover that the backscattering part of the hamiltonian does not change under our transformation while the current operators @xmath141 change their signs .
this means , in turn , that the contribution to the total current we are calculating changes its sign when the bias voltage changes its sign .
hence , it does not contribute to the _ rectification _
thus , eq .  ( [ third - order - main ] ) describes the main contribution to the rectification current .
we now find the leading contribution to the rectification current that contains operators @xmath171 .
first of all , let us check that renormalized operators @xmath332 remain small at @xmath333}/\tau_c$ ] .
indeed , @xmath334 and @xmath335 . hence , @xmath336}<(v\tau_c)^{1 - 2g/3}<(v\tau_c)^{1/3}\ll 1 $ ] .
thus , all renormalized operators remain small and the maximal @xmath245 for which the perturbation theory can be used was found correctly in the previous subsection . in what follows we focus on the case with @xmath245 of the order of its maximal allowed value .
the opposite limit @xmath337 is trivial since in that limit the second order contribution in @xmath245 always dominates the ratchet current .
next , let us check how the operators @xmath171 affect second order contributions to the rectification current .
the corresponding second order contribution @xmath338 ^ 2v^2/\tau_c^3\le ( u\tau_c)^4/\tau_c\sim ( v\tau_c)^{4 - 4g[1-a]}/\tau_c$ ] .
we need to compare it with the current ( [ third - order - main ] ) , @xmath339}/\tau_c$ ] .
one easily sees that @xmath340>1 + 3g[1-a]$ ] since @xmath341<1/3 $ ] .
thus , the second order contribution to the rectification current can be neglected in comparison with the third order contribution even after @xmath171 are taken into account .    at the same time , in a certain region of parameters the dominant higher - order contribution to the rectification current contains an operator @xmath171 . as we will see , such contribution contains exactly one operator @xmath171 .
we thus begin our analysis by excluding contributions which contain three or more @xmath171 s . indeed , any such contribution @xmath342}/\tau_c$ ] .
we need to compare this estimate with @xmath339}/\tau_c$ ] .
one easily sees that @xmath343>[1 + 3g(1-a)]$ ] and hence we do not need to take into account contributions with three or more operators @xmath171 .
thus , it remains to consider higher - order contributions with one or two operators @xmath171 .
we first consider contributions with exactly one operator @xmath344)$ ] , where @xmath345 can be set to zero as discussed above .
the contribution also contains at least two operators @xmath163 .
it is sufficient to consider the case with all @xmath346 in the definition of the operators @xmath163 @xmath297 .
indeed , a contribution with an operator with @xmath347 can be increased by substituting it with @xmath348 operators @xmath349)$ ] .
moreover , we can assume that for each @xmath24 only operators @xmath350 or only operators @xmath351 enter .
let us denote the number of the operators of the form @xmath352)$ ] in the expression for the contribution under consideration as @xmath353 .
then we can estimate the contribution to the rectification current as @xmath354^{l_3}[uv^{g[1-b]-1}]^{l_1}[uv^{g(1-c)-1}]^{l_2},\ ] ] where @xmath355 ^ 2/(1+\alpha+\alpha^2)$ ] . the same argument as in the previous subsection shows that the maximal contribution corresponds to @xmath356 .
then for a general @xmath184 , the product of the vertex operators @xmath357 gives a nonzero result after averaging with respect to the quadratic part of the action only if @xmath358 and @xmath359 .
note that we can assume that both @xmath360 and @xmath361 are nonzero and @xmath362 .
otherwise the operator @xmath171 would have the same from as one of the operators @xmath163 and the analysis from the previous subsection would apply .    we want to compare the contributions ( [ third - order - main ] ) and ( [ est - one - w ] ) . at @xmath363}/\tau_c$ ] we find @xmath364 = \gamma(s_1,s_3)\log ( v\tau_c ) ; \nonumber \\ & \gamma = 1+\frac{g}{1+\alpha+\alpha^2}[-5\alpha^2+|s_1|(1-\alpha^2)+(s_1+s_3\alpha)^2].\end{aligned}\ ] ] if @xmath114 is positive for every choice of @xmath365 then @xmath310 is the main contribution to the rectification current . if @xmath366 for a certain choice of @xmath367 then the leading contribution comes from @xmath368 .
we thus want to investigate at what conditions @xmath366 and find what choice of @xmath360 and @xmath345 minimizes @xmath114 .
that choice determines the power dependence of the rectification current ( [ est - one - w ] ) on the voltage .
first , let us prove that @xmath114 is minimal , if @xmath369 .
indeed , let us compare @xmath370 with @xmath371 , where @xmath372 and @xmath373 is arbitrary .
one finds @xmath374 \nonumber\\ & \ge\frac{g}{1+\alpha+\alpha^2}(1-\alpha)[(|p_1|-1)(1+\alpha)-(1-\alpha)]>0.\end{aligned}\ ] ] thus , we can focus on @xmath375 ( the case of @xmath376 is completely analogous ) .
next , we prove that a negative @xmath377 is minimal at @xmath378 .
indeed , a negative @xmath377 implies that @xmath379/[1+\alpha+\alpha^2]<0 $ ] .
hence , @xmath380 .
taking into account that @xmath381 , we find that @xmath382 is minimal at @xmath378 .
this allows one to establish that @xmath114 is minimal at @xmath378 .
the remaining task is simple .
we just determine at what conditions @xmath383 is negative , i.e. , we need to investigate the inequality @xmath384}{1+\alpha+\alpha^2}<0.\ ] ] one easily sees that @xmath114 can only be negative , if @xmath385 , i.e. , for relatively weak repulsive electron interaction .
@xmath114 is negative in the largest interval of @xmath184 at @xmath386 . in that case ,
@xmath366 for @xmath387 . the contributions with @xmath171 matter only if @xmath388 .
they give rise to the rectification current of the form ( [ answer2 ] ) at negative @xmath114 .
the last question we must address in this subsection concerns the role of the contributions to the rectification current with two operators of the form @xmath299 .
an estimation of such contributions is similar to ( [ est - one - w ] ) : @xmath389^{l_3}[uv^{g[1-b]-1}]^{l_1},\ ] ] where @xmath390 ^ 2/(1+\alpha+\alpha^2)$ ] are the scaling dimensions of the operators @xmath391+is^{x , y}_3[\phi^3_o(x=0)-\phi^3_i(x=0)])$ ] .
as above , we first divide @xmath392 by @xmath310 and ask when the ratio is greater than 1 : @xmath393=\mu(s^x_1,s^x_3,s_1^y , s_3^y)\log ( v\tau_c ) ;   \nonumber \\ & \mu=2+\frac{g}{1+\alpha+\alpha^2 } ( -7\alpha^2+l_1(1-\alpha^2)+\sum_{r = x , y}(s_1^r+s_3^r\alpha)^2).\end{aligned}\ ] ] if @xmath394 then @xmath395 .
obviously , a negative @xmath396 implies that @xmath397 and hence @xmath398 .
let us now compare @xmath396 and @xmath383 at @xmath399 .
one finds    @xmath400    one easily checks that the above difference is positive at @xmath399 .
hence , contributions with one or no operators @xmath171 are always more important than contributions with two such operators .
we found that the leading contribution to the rectification current depends on the sign of @xmath383 , eq .
( [ inv - ineq ] ) . for strong repulsive interaction ( @xmath401 ) , @xmath383 is always positive .
if @xmath402 , at sufficiently low voltages @xmath403}/\tau_c$ ] the current scales as    @xmath404    at a negative @xmath114 and low voltages @xmath271
@xmath405    the exponents in both voltage dependences are negative .
this is related to the fact that our calculations are only valid in an interval of low voltages , @xmath406}/\tau_c$ ] .
similar to section viii the @xmath1-@xmath2 curve for the rectification current exhibits a bump at @xmath271 .
the maximal rectification current at @xmath402    @xmath407    at @xmath366    @xmath408    one can easily verify that the maximal current ( [ answer3],[answer4 ] ) exceeds the maximal possible rectification current at the same voltage @xmath2 in a linear wire @xcite @xmath409^{3g},[v\tau_c])$ ] .
we have found the rectification current in the absence ( [ 3-no - time ] ) and in the presence ( [ answer1],[answer2 ] ) of the time - reversal symmetry in y - junctions . in all cases
the maximal rectification current is greater than in a linear wire with the same interaction strength and bias voltage . in the absence of the time - reversal symmetry
the rectification current can be comparable with the total ac current @xmath0 for sufficiently strong interaction strength , i.e. , it achieves its maximal possible order of magnitude .
this reflects the fact that both electron interaction and time - reversal symmetry breaking facilitate rectification . for most values of parameters
the rectification current is a nonmonotonous function of the bias voltage .
our calculations are valid in the vicinity of various fixed points in the low voltage regime . in a general case ,
a junction is controlled by a stable fixed point at low voltages . for repulsive interaction of spin - polarized particles
there is only one stable fixed point : three disconnected wires.@xcite we found a stronger rectification effect near that fixed point than for two weakly connected wires .
however , the current is low for weakly connected wires .
the diode effect is much stronger in the vicinity of unstable fixed points .
thus , it is important to understand how to tune the system close to those fixed points . some of them may be tricky to realize experimentally . indeed ,
as discussed above , negative elements in the matching matrix @xmath87 imply andreev reflection and could be obtained in a hybrid normal - superconductor structure or in cold atom systems @xcite . at the same time , it is straightforward to make an `` island junction '' @xcite , fig .  1 , with positive matching matrix elements @xmath410 for example , one can use three line junctions between three quantum hall systems .
the role of impurities is played by three constrictions in the junctions .
alternatively , one can use a single quantum hall island confined between edges ab , cd and ef ( fig .  1 ) .
the rectification current scales as @xmath411 at @xmath412 .
time - reversal symmetry is broken in such setup .
thus , the rectification current can be made comparable to the total ac current even in the low voltage regime .
this is the main result of the paper .
most of the time we ignored phase factors @xmath413 in tunneling operators .
as the example of three weakly connected wires shows , for special values of the phases the rectification effect is suppressed .
it has the same order of magnitude for other values of the phases .
we considered the simplest example of a junction : 3 spin - polarized wires .
it would be interesting to generalize our results to a system with spin .
this may result in a more complicated behavior as the minimal model with spin includes 6 channels : two for each wire .
still , we expect a similar physics .
in particular , a time - reversal - invariant system with spin can be obtained from two copies of the `` island junctions '' with opposite spins and chiralities . in the absence of the interaction between the copies
, the problem reduces to a spin - polarized island junction .
an actual realization , based , e.g. , on a topological insulator , must involve interaction between opposite spins .
we expect that such interaction does not change the qualitative picture .