at mesoscopic scale and at low temperature , electrons in metallic samples keep their phase coherence on a length @xmath2 which is bigger than the sample size .
transport and thermodynamic properties of the system are then sensitive to interference between electronic wave functions and present spectacular signatures of this phase coherence . to study these effects
the ring geometry is especially suitable . indeed in the presence of a magnetic flux
@xmath3 through the ring the periodic boundary conditions for electronic wavefunctions acquire a phase factor @xmath4 with @xmath5 the flux quantum @xcite . as a result , the magnetoconductance of a phase coherent ring exhibits quantum oscillations which periodicity corresponds to one flux quantum through the area of the sample @xcite .
the phase of the first harmonics of these oscillations is sample specific so that these harmonics do not survive ensemble average .
in contrast the second harmonics have a contribution which resists this averaging .
this results from the interference between time reversed paths around the ring ( weak localisation contribution ) .
these @xmath6 periodic oscillations were observed in transport measurements on long cylinders or connected arrays of rings @xcite .
their sign corresponds to a positive magnetoconductance in zero field . in the case of singly connected geometries , like full disks ,
the signature of weak localization consists in a single peak of positive magnetoconductance which width corresponds to @xmath7 through the sample @xcite .
magneto - transport experiments on connected systems constitute a very sensitive and powerful probe for the investigation of sample specific signatures of quantum transport .
however , because of strong coupling between the system and the measuring device , quantum corrections represent a small fraction of the conductance ( of the order of @xmath8 where @xmath9 is the dimensionless conductance expressed in @xmath10 units ) which is still dominated by the classical drude response in the diffusive regime @xmath11 .
there exists a number of experiments which can address some of the electronic properties of mesoscopic samples without coupling to macroscopic wires .
this is the case of ac conductance experiments where aharonov bohm rings are coupled to an electromagnetic field . in contrast with the connected case , the response of an isolated system can be dominated by quantum effects .
moreover , the quasi discrete nature of the energy spectrum and the sensitivity to the statistical ensemble ( canonical or grand - canonical ) are new features of isolated mesoscopic systems . in particular it has been emphasized that the average absorbtion of isolated mesoscopic systems is determined by the energy level statistics and its sensitivity to time reversal symmetry breaking by a magnetic field .
the first experiments done in this spirit were performed by coupling an array of disconnected gaas / gaalas rings to a strip - line superconducting resonator @xcite . in such a geometry
the rings experience both ac magnetic and electric field .
the magnetic response of the rings i.e. their orbital magnetism is directly related to persistent currents in the zero frequency limit @xcite . on the other hand the electric response of isolated metallic sample is related to the screening of the electric field inside the metal .
the induced charge displacement is at the origin of the polarisability @xmath12 , defined as the ratio between the induced electric dipole @xmath13 and the applied electric field @xmath14 ( @xmath15 ) .
the polarisability is known to be essentially determined by the geometry of the sample with correction of the order of @xmath16 , with @xmath17 the screening length and @xmath18 the typical size of the system @xcite .
it has been recently predicted that this quantity is sensitive to phase coherence around the ring @xcite and is thus expected to present flux oscillations .
the electric contribution can be in the particular case of gaas rings of the same order of magnitude as the magnetic response @xcite .    to be able to distinguish between the two types of response we have designed a superconducting lc resonator which capacitive part and inductive part are physically separated . in this paper
we present measurements of both magnetic and electric response of aharonov - bohm rings .
note that these experiments are done on the very same array of rings for both types of response , giving us the opportunity to compare them .
preliminary account of measurement of the electric response was given in reference @xcite .
the paper is organised in the following way .
section [ experimentalsetup ] gives a detailed presentation of the sample , an array of aharonov - bohm rings , and the resonating technic used to measure the magnetic and electric response .
results on the non - dissipative part of the flux dependent electric response are presented in section [ magnetopolarisabilty ] .
a comparison with theoretical predictions is given , including frequency dependence .
temperature and electronic density dependence are also investigated .
the next section focuses on the dissipative part of the magnetopolarisability of the rings .
theoretical results for this quantity are derived and compared to the experiment .
the section [ magneticresponse ] is devoted to the measurements of the magnetic response of the same rings .
despite the fact that the signal is then smaller , the magnetic response of the rings is detected and compared to predictions on averaged persistent currents .
we conclude by a comparison between the magnetic and electric response .
[ cols="^,^,^ " , ]         we have studied the electric and magnetic susceptibilities of isolated aharonov - bohm rings . our system is an array of @xmath0 2d rings etched by reactive ion etching in a high mobility algaas / gaas heterojunction .
the characteristics of these rings are given in table [ tab : characteristicsoftherings ] .
they are ballistic in the transverse direction and diffusive longitudinally ( @xmath19 and @xmath20 ) .
it is important to perform a deep etching of the heterojunction ( down to gaas ) in order to minimize high frequency losses , which have been observed to be important in etched algaas . because of etching the electronic density is strongly depressed .
however we are able to recover the nominal density of the heterojunction by illuminating the rings with a infrared diode placed close to the sample in the dilution refrigerator . for each illumination a current of @xmath21 is run through the diode during several minutes .
measurements are done at least one hour after the illumination in order to ensure good stability of the sample .
an upper value of the estimated illumination power coupled to the sample is 600 photons / s with a wavelength of 766  nm . with this setup
we are thus able to perform measurements at different electronic density .
the control on the density is rather qualitative because of the difficulty to calibrate the illumination procedure .
we have checked the effect of illumination on a connected aharonov - bohm ( ab ) ring ( figure [ fig : abillum ] ) . at zero illumination time
the conductance of the ring is zero .
on such a sample we can follow the ab oscillations when the resistance decreases by more than an order of magnitude with illumination . as a consequence a clear effect of illuminating the ring is to increase its conductance .
the fourier transform of the resistance of the ring is shown on figure [ fig : fftabillum ] .
we see for each illumination an oscillation whose periodicity is consistent with a flux quantum @xmath22 in the area of the ring .
however the fourier transform shows that the peak corresponding to this periodicity changes with illumination both in shape and in amplitude .
the amplitude increases with illumination due to the increase of ab oscillations . the fact that the shape , and in particular the width , of the peak changes with illumination is an indication that the width of the ring increases with illumination time . to be more precise
the width of the rings is multiplied by a factor 2 between the first and last curve of figure [ fig : fftabillum ] .
note that the increase of the electronic density has also been shown to induce an increase of the electronic mobility @xcite .
peak is strongly dependent on illumination .
the curves are shifted for clarity.,width=302 ]    in order to study the disorder average we have measured conductance of a single ring and a mesh , representing a two - dimensional square array , etched in the same type of heterojunction than the rings .
the magnetoconductance is shown on figure [ fig : angr ] .
as expected the ab effect disappears under ensemble averaging .
the @xmath6 oscillations on the other hand remains on the mesh . in this case the triangular shape of the magnetoconductance is attributed to weak localisation in the wire of the mesh .
signal disappears with ensemble average , so that in the mesh only the @xmath6 component remains .
note the triangular shape of the magnetoconductance on the mesh .
the curves are shifted for clarity.,width=302 ]          to measure the electric or magnetic response of the rings we couple them to a superconducting micro - resonator and detect the changes in its properties .
this resonator is made by optical lithography .
it consists of a niobium strip - line deposited on a sapphire substrate .
this substrate has been preferred to silicon or gaas because it induces the weakest temperature dependence of the resonance frequency and gives the best quality factor due to the quality of the niobium layer on sapphire .
a schematic drawing is given in figure [ fig : resonator ] .
the width of the wire constituting the resonator is 2 @xmath23 m , the thickness 1 @xmath23 m and the spacing between two adjacent wires is 4 @xmath23 m .
the total length of the capacitance or the inductance is 10 or 20 cm . in this kind of resonator the inductance is physically separated from the capacitance by a distance of 300 @xmath23 m , allowing to submit the sample only to an electric field ( or to a magnetic field ) to measure its electric ( or magnetic ) response .
this separation between magnetic and electric response has been checked by deposition of a paramagnetic system ( dpph ) alternatively on the capacitive and inductive part of the resonator .
a magnetic spin resonance signal was only observed when dpph was on the inductive part .
the resonance frequency of the bare resonator varies between 200 mhz and 400 mhz depending on the geometry .
its quality factor is 10000 at 4.2 k and 200000 at 20 mk .
the resonator can be modelled by an lc circuit of resistance @xmath24 , inductance @xmath25 , capacitance @xmath26 , whose resonance frequency is @xmath27 and quality factor @xmath28 . from the value of the higher resonance frequencies of the resonator
we have estimated that the residual capacitance of the meander line is at least 10 times smaller than @xmath26 .
due to the meissner effect , the dc field just above the resonator is strongly inhomogeneous . in order to minimize this effect ,
we have inserted a thin , @xmath29m - thick , mylar film between the detector and the rings .
this reduces the field inhomogeneity to about @xmath30% , which is of the order of fluctuations in the lithography .      in order to measure their electric response the rings
are placed on top of the capacitance of the resonator .
note that with this procedure the rings are not well aligned with the resonator so that they do not have the same coupling with the capacitance .
this is not a problem as soon as only linear response is investigated .
we note @xmath31 the polarisability averaged over disorder of a ring at the frequency @xmath32 .
the impedance of the capacitance @xmath26 slightly modified by the rings reads  : @xmath33 in this expression @xmath34 is the number of rings coupled to the capacitance , @xmath35 is an averaged coefficient measuring the dielectric coupling between one ring and the capacitance .
the capacitance with the rings is equivalent to a capacitance @xmath36 in series with a resistance @xmath37 .
hence  : @xmath38 the frequency shift due to the rings is then  : @xmath39 the quality factor is determined by  : @xmath40 so that , with @xmath41 at resonance  : @xmath42 provided that @xmath43 .
the electric coupling coefficient is estimated in appendix [ electric_coupling ] . knowing this value and the number of rings coupled to the resonator
, it is possible to evaluate quantitatively the polarisability of the rings by measuring the resonance frequency shift ( equation [ frequencyshift ] ) and the variation of the quality factor ( equation [ qshift ] ) .
when the rings are placed on top of the inductance @xmath25 of the resonator , this inductance is shifted because of their magnetic response @xmath44 according to  : @xmath45 with @xmath34 the number of rings coupled to the resonator , @xmath46 the magnetic coupling coefficient between one ring and the inductance , which has the dimension of the inverse of a volume .
note that , properly defined , the coupling coefficient @xmath46 is of the same order of magnitude than @xmath35 .
more precisely the estimation of @xmath35 and @xmath46 done in appendix leads to @xmath47 , as expected from reference @xcite .
following the same reasoning than for the electric coupling , the properties of the resonator are modified according to  : @xmath48    from previous equations it is in principle possible to measure the absolute value of @xmath12 or @xmath49 .
however when a gaas sample is on the inductive or capacitive part of the resonator the modification of the resonance is dominated by the influence of the substrate . as a consequence
it is very difficult to have an accurate absolute measurement .
nevertheless relative measurements are possible so that the variation of the electric or magnetic response with magnetic field can be detected in a reliable way .
the reflected signal of the resonator is measured with the setup of figure [ fig : reflectedsignal ] and used in a feedback loop to lock the frequency of a rf generator to the resonance frequency .
the setup is summarized in figure [ fig : asservissement ] .
the resonator is coupled capacitively to the external circuit using on - chip capacitances . in order to preserve the quality factor of the resonator we work in a configuration where the resonator is undercoupled .
the rf power injected is sufficiently low ( @xmath50 10 pw ) so as not to heat the sample .
the frequency of the rf generator is modulated at @xmath51 and the signal from the resonator is detected by a lock - in detector at the frequency of the modulation .
the lock - in signal is to first approximation the derivative of the resonance peak : it gives an error signal i.e. this signal is zero at resonance , and changes sign when the frequency of the generator is higher or lower than the resonance frequency . using this signal in a feedback loop the frequency of the rf generator
is locked to the resonance frequency .
this way , by measuring the feedback signal , one has directly access to the shift of the resonance frequency . to enhance the accuracy we modulate the magnetic field by a 1 g ac field oscillating at 30 hz , produced by a small superconducting coil close to the sample , and detect the modulated resonance frequency with a lock - in detector .
at this point we consider that the frequency of the generator is locked to the resonance frequency by the previous setup .
the signal measured is the signal reflected from the resonator . as a consequence it is related to the reflexion coefficient @xmath52 , with @xmath53 the impedance of the resonator and the coupling capacitance and @xmath54 the impedance matched by the external circuit .
we assume that near the resonance frequency the impedance @xmath53 reads  : @xmath55 with @xmath56 the resonance frequency . in the limit @xmath57 , which correspond to a very undercoupled resonator ,
the reflected signal is a linear function of @xmath53 . as a consequence
if the rf signal is frequency modulated at @xmath51 around the resonance frequency @xmath56 the reflected signal at @xmath58 is related to the second derivative of the real part @xmath53 , which is proportional to @xmath59 . this way by measuring the signal at @xmath58 we have access to the quality factor .
however when the frequency modulation is not small compared to the width of the resonance peak or the resonator is not very undercoupled to the external circuit , the relation between the signal at @xmath58 and the quality factor is not straightforward and needs calibration .
in this part we present measurements of the flux dependent polarisability of the rings , which are placed on the capacitive part of the resonator as described in section [ experimentalsetup ] . in this configuration the resonance frequency is decreased by 15 % , due to the dielectric constant of the gaas substrate .
moreover the quality factor drops down to 3000 at 20 mk at zero illumination .
this strong decrease is attributed to dielectric losses in the heterojunction .
the derivative of the resonance frequency of the resonator with the rings is shown on figure [ fig : polarextract ] ( a ) .
this signal is a straight line , on top of which small oscillations are superimposed .
the straight line is due to the field dependence of the penetration length in niobium , which constitutes the resonator .
this behaviour has been verified to be the same with or without the rings .
the oscillating signal is on the other hand attributed to the flux dependent electric response of the rings .
these oscillations are extracted by subtracting the base line ( figure [ fig : polarextract ] ( b ) ) .
note their anharmonicity as well as the existence of an aperiodic signal as illustrated by the fourier transform of the data ( figure [ fig : polarextract ] ( c ) ) showing a well defined peak . in order to focus on this periodic contribution , which is the expected signature of phase coherence , a numerical high pass filter with a cut - off frequency of 0.05 g@xmath60 ( corresponding to the arrow on figure [ fig : polarextract ] ( c ) )
is applied and the signal is then numerically integrated in order to have the frequency shift due to the rings ( figure [ fig : polarextract ] ( d ) ) .
this shift is proportional to the variation of polarisability versus magnetic field according to formula [ frequencyshift ] .
we will return to the aperiodic signal in the section devoted to illumination effect .
the frequency shift due to the rings is periodic with a period of approximately 12.5 g. from the fourier transform ( figure [ fig : polarextract ] ( c ) ) the period of the oscillation is deduced to be consistent with half a flux quantum @xmath61 in a ring with no signature of @xmath22 periodicity , as expected for an aharonov - bohm effect averaged over many rings @xcite .
note the extra broadening ( by more than a factor 2 ) of this @xmath6 peak compared to the measurements on a single connected ring .
we interpret this as resulting from the dispersion in circumfrences in the different rings .
the sign of frequency shift is negative at low magnetic field which means according to formula [ frequencyshift ] that the magneto - polarisability is positive , i.e. @xmath62 at low magnetic field .
the screening is thus better when time reversal symmetry is broken by magnetic field .
the scale of the signal is given by the amplitude of the first oscillation .
from figure [ fig : polarextract ] ( d ) we deduced @xmath63 g@xmath64 . note that this value means detecting a frequency shift of 100hz on a frequency of 350 mhz . with the coupling coefficient estimated in appendix [ electric_coupling ] it leads to the value of the magnetopolarisability @xmath65 , where @xmath66 is the calculated polarisability of a quasi one dimensional ( quasi-1d ) circular ring of radius @xmath67 .       at different value of the parameter @xmath68 . note that the value of magnetopolarisability is zero at zero whatever the level broadening.,width=298 ]
our experiment shows that there is a flux correction to the polarisability of the rings , which is _ positive _ at low field .
let s now compare this result to recent theoretical predictions . since we are using a ring geometry we are going alternatively from a situation where the system presents time reversal symmetry ( at flux values of @xmath69 , with @xmath70 ) to the case where time reversal symmetry is broken by magnetic field .
in the random matrix theory ( rmt ) the first case corresponds to gaussian orthogonal ensemble ( goe ) whereas the second is related to the gaussian unitary ensemble ( gue ) .
so the quantity to be compared with theoretical predictions , which evaluate the variation of a physical variable @xmath71 between goe and gue , is @xmath72 defined as @xmath73 .
note that since the rings are semi - ballistic , the transition with magnetic field may not be exactly from goe to gue .
the polarisability of small metallic grains was studied using rmt first by gorkov and eliashberg @xcite .
the sensitivity of the electrostatic properties of mesoscopic systems to quantum coherence has been emphasized by bttiker for connected geometries @xcite .
the phase coherent correction to the polarisability of isolated systems was recently theoretically investigated .
efetov found that it is possible to relate self consistently this correction to the flux dependence of the screened potential @xcite .
two recent works have calculated this effect in the diffusive regime using linear response formalism ( noat _ et al . _
@xcite ) or supersymmetry techniques ( blanter and mirlin @xcite ) .    in the grand canonical ensemble ( gce )
the chemical potential in each ring is supposed to be constant .
it describes a situation where the rings are connected to a reservoir of particules . _ a priori _ this is not the case in the experiment where the rings are isolated but as the theory is simpler in gce we recall first the result in this statistical ensemble .
no flux dependence for the polarisability is predicted if the rf pulsation @xmath32 is much smaller than the inverse relaxation time @xmath74 .
however when @xmath75 the magnetopolarisability is related to the flux dependence of the diagonal matrix element of the screened potential  : @xmath76 @xmath77 is the disorder averaged square of the diagonal matrix element of the screened potential @xmath78 at energy @xmath23 , the mean chemical potential of the rings .
@xmath79 is the applied electric field .
we note @xmath80 the eigenstates of the unperturbed system .
this matrix element is then given by  : @xmath81 from this expression it appears that the magnetopolarisability is related to the difference of correlation function of the eigenstates with and without time reversal symmetry .
this correlation function has been computed in the diffusive regime within a supersymetric @xmath82-model approach @xcite  : @xmath83 \left [ 1 + 2 \pi_d({\bf r_1},{\bf r_2 } ) \right ] \quad ( goe)\end{gathered}\ ] ] @xmath84 \left [ 1 + \pi_d({\bf r_1},{\bf r_2 } ) \right ] \quad ( gue)\end{gathered}\ ] ] with @xmath85 the volume of the sample , @xmath86 a short range function which decays on the length scale of the mean free path and @xmath87 the diffusion propagator .
the correction due to the short range term @xmath86 has been shown to be negligible @xcite . by considering only the diffusion term
the magnetopolarisability is given by  : @xmath88 note that this derivation of the magnetopolarisability is equivalent to the one used by noat _
et al _ @xcite based on the following rmt argument  : @xmath89 this relation can also be obtained from [ eq : corrgoe ] and [ eq : corrgue ] using the fact that  : @xmath90 due to symmetry properties of the screened potential .
the calculation of the magnetopolarisability using formula [ polardiffuson ] for the case of a quasi-1d ring ( appendix [ quasi1dpolar ] ) leads to  : @xmath91 @xmath92 is a function related to the geometry and the dimension of the sample . using this expression and the value of table [ tab : characteristicsoftherings ] we have @xmath93 .
in our experiment the rings are isolated , the number of electrons in each ring is supposed to be constant , so that the result of the canonical ensemble ( ce ) should apply . at t=0 and zero frequency
the flux dependent correction to polarisability is found to be zero
. however at @xmath94 the gce result is recovered .
the complete frequency dependence of the magnetopolarisability in the ce has been recently derived by blanter and mirlin @xcite .
following their reasoning but taking into account the level broadening @xmath74 we can write  : @xmath95 with @xmath96 a function which depends only on the statistic of energy levels : @xmath97\end{gathered}\ ] ] @xmath98 and @xmath99 are respectively the two and three levels correlation function , known from rmt @xcite . by evaluating this expression versus frequency at different value of level broadening we get the results shown on figure [ fig : polarce ] .
the behaviour at low value of the level broadening is in qualitative agreement with result of reference @xcite .
in particular the magnetopolarisability is found to be zero at zero frequency ( in our calculation the value of @xmath100 is at least 25 time smaller than @xmath101 ) .
the present experiment was performed at @xmath102 , the ce magnetopolarisability is equal at most to 50 % of the gce value ( in the limit of small level broadening ) . as a consequence the expected value for @xmath103 is then @xmath104 which is of the same order of magnitude as the experimental value .
note that the measurement is not sufficiently accurate to give an estimate of the level broadening by comparing the experimental result with the curves of figure [ fig : polarce ] .
a very interesting extension of the experiment would be to study the magnetopolarisability at different frequencies in order to test the theoretical predictions .
this could be done by working with resonators with smaller inductances .
component of the fourier transform of the signal .
the fitting function used is proportional to @xmath105 with two fitting function for @xmath106 .
first we took the phase coherence length measured on connected sample @xmath107 , which exhibits a @xmath108 behaviour .
the other fitting function is @xmath109 .
the best agreement is found with an exponential decay with a temperature scale of 90 mk .
inset : temperature dependence of the frequency shift due to the rings.,width=298 ]    the temperature dependence of the signal is also investigated .
the magnetopolarisability decreases with temperature ( inset of figure [ fig : polartemp ] ) .
theoretically the effect of temperature on magnetopolarisability has not been studied yet .
we will base our analysis of the temperature dependence on the hypothesis that the amplitude of the signal is related to the phase coherence length @xmath2 in the same way as weak - localization . in this case
the amplitude of the @xmath6 component of the signal is proportional to @xmath110 @xcite . in figure
[ fig : polartemp ] the temperature dependence of this component is shown .
we have tried to fit it using two laws for @xmath106 . first using the behaviour deduced from the measurements on connected wires @xcite we have tried the experimental value @xmath107 which exhibits a @xmath108 dependence , as expected for 1d system @xcite .
it leads to a poor agreement with experimental points .
using for the phase coherence time the result of electron - electron interaction in quantum dots ( od system ) @xcite @xmath111 , leading in the diffusive regime to @xmath112 , gives a better agreement . in this case the temperature scale is found to be 90 mk .
we deduced from this value @xmath113 mk at 18 mk , i.e. much smaller than the level spacing .
the phase coherence length deduced from this analysis is 10 times higher than the length measured on connected sample @xcite .
we relate this difference between the non - connected and connected case to the fact that whereas the connected samples are one dimensional with a continuous energy spectrum due to the strong coupling with the reservoirs , the spectrum of the non - connected rings is discrete . as a concluding remark on this temperature dependance , we want to emphasize the need of a deeper theoretical analysis of the behaviour of the magnetopolarisability versus temperature .               component before numerical integration , for different illumination times.,width=298 ]    using the procedure described in [ sec : rings ] we are able to study the influence of electronic density on magnetopolarisability . on figure
[ fig : fftpolarillum ] the fourier transform of the derivative of frequency shift , when the base line due to the resonator is removed , is shown at different illumination time .
as expected the fourier transform exhibits a @xmath6 peak .
note also the low frequency component which corresponds to the aperiodic signal seen on figure [ fig : polarextract ] ( b ) .
the @xmath6 peak depends on electronic density .
its amplitude shows first an increase and then decreases at high illumination .
moreover the width of the @xmath6 peak increases showing that the rings widen with illumination .
the peak becomes asymmetric suggesting that rings having long circumferences initially not populated contribute to the signal at high illumination .
following the procedure described in section [ magnetopolarisabilty ] we measured the amplitude of magnetopolarisability versus illumination time .
it yields figure [ fig : effetillum ] , which shows the change of the amplitude of the magnetopolarisability . at first the signal increases , then for illumination time above 1400 s the amplitude of oscillation decreases .
we interpret this behaviour in the following way . before illumination
the electronic density in the rings obtained after deep etching of the 2deg is strongly depressed compared to the nominal value . as a consequence an important fraction of the rings are likely to be localized and do not contribute to the magnetopolarisability . in this regime
the signal is small .
after illumination the number of rings contributing to the signal increases so that the frequency shift due to the rings increases . at high enough electronic density
when the rings are sufficiently populated so that they contain delocalized electrons the theoretical results obtained in the diffusive regime are expected to be valid leading to a @xmath8 dependence ( formula [ da ] ) , with @xmath114 the dimensionless conductance . this is a possible explanation for the decrease of the magnetopolarisability observed at high illumination level .
note that we can not exclude also a reduction of the screening length @xmath17 due to illumination . from formula [ da ]
we deduce then a decrease of the signal .
however we believe that the screening length changes very weakly with illumination because this length is essentially determined by the density of states at the fermi energy , which is independent of energy for a 2d system ( see table [ tab : characteristicsoftherings ] ) .
we have analysed so far the @xmath6 periodic component of the magnetopolarisability signal obtained after filtering low frequency ( see figure [ fig : polarextract ] ) . on the other hand
the whole integrated signal is depicted in figure [ fig : polarstep ] .
one can clearly see a triangular shape dependence of the signal with magnetic field superimposed to the oscillations , very similar to the weak - localisation conductance of the connected mesh shown in figure [ fig : angr ] . note that this behaviour is only present at low temperature
, it completely disappears for temperature higher than 300 mk .
the amplitude of this extra signal due to the finite width of the rings strongly increases and sharpens with illumination .
we think that it is reasonable to attribute this evolution to the increase of the width of the rings .
note that a similar evolution has been previously observed in the ac magnetoconductance of ballistic gaas squares @xcite .
versus @xmath32 in gce at different value of the parameter @xmath68.,width=298 ]    , width=298 ]    by measuring the quality factor of the resonator versus magnetic field , we have access to the flux dependent electric absorption ( formula [ qshift ] ) , which is related to the conductance , in the case of an electric dipole , through : @xmath115 the contribution due to the rings ( figure [ fig : fig2 ] ) exhibits the same periodicity as the frequency shift , which corresponds to half a quantum flux in a ring . the low field signal decreases .
it corresponds to a negative magnetoconductance , i.e. opposite to weak - localization .
this surprising sign was pointed out in the context of the magnetoconductance of rings submitted to an oscillating magnetic flux in the discrete spectrum limit @xcite .    to explain this result one has to take into account the level spacing distribution in a disordered system @xcite .
the sign and amplitude of the typical variation of electric absorption are understandable using the fact that level repulsion in a disordered system is higher in gue than in goe . following reference @xcite the flux dependent electric absorption in a system described by eigenvalues @xmath116 and the corresponding eigenfunctions @xmath80
could be written in a linear response regime  : @xmath117 with @xmath118 .
we will first consider this expression in the gce where an average over the chemical potential is computed . with this procedure @xmath119 and @xmath120 , where @xmath121 is the range over which the average over the chemical potential is done .
the first term in the right hand side of equation [ absorption ] then reads  : @xmath122 note that in this sum the energies @xmath123 and @xmath124 have to belong to the range @xmath125 $ ] . using this constraint we replace the sum by an integral  : @xmath126 with @xmath98 the two energy level correlation function . in this expression
we neglect the flux dependence of the matrix element and note its average value @xmath127 . with this approximation equation [ nondiag ]
reads  : @xmath128 the debye term of equation [ absorption ] is equal in the gce to  : @xmath129 in the gce in the dynamical regime the flux correction to the polarisability is given by formula [ dagce ] at @xmath130 . using the following correlation function @xcite  : @xmath131 \pi_d(\mathbf{r_1},\mathbf{r_2 } )
\quad ( goe)\end{gathered}\ ] ] we have @xmath132 .
hence : @xmath133 which can be evaluated numerically ( figure [ fig : electricabsorption ] ) .    for isolated rings we have to apply the result of ce .
it is then possible to estimate the correction to electric absorption by using the same treatment as for the real part of polarisability .
it leads to  : @xmath134\end{gathered}\ ] ] numerical estimation of this formula at different value of the level broadening leads to the behaviour shown on figure [ fig : absorpce ] .
the electric absorption is always negative at low frequency and may change sign at low value of @xmath68 .    in order to compare these calculations with our experimental result we will compute the ratio @xmath135 .
it is worth noting that in our experiment this quantity is given according to equations ( [ frequencyshift],[qshift ] ) by @xmath136 and is independent of the number of rings coupled to the resonator or the electric coupling coefficient @xmath35 .
experimentally we find @xmath137 at illumination time zero and -0.23 after 420 s of illumination . theoretically , at @xmath138 , @xmath139 and @xmath140 so that the expected value of @xmath141 is around 1 .
it corresponds to small @xmath68 . for higher value of this parameter
the ratio is of the same amplitude or higher . as a consequence
the predicted behaviour is consistent with the experimental value for the sign , but the theoretical amplitude is too high by more than a factor 2 .
this conclusion is different from our previous statement where the frequency dependence of the real part of magnetopolarisability was not taken into account @xcite .
.,width=298 ]            due to the design of the resonator we have also the opportunity to investigate the magnetic response of the same aharonov - bohm rings used for the measurements of the electric response . in this case the rings are placed on top of the inductive part of the resonator .
note that to do so we have to warm - up , cool down and re - illuminate the rings . as a consequence , strictly speaking , the rings are not the same than for the measurement of the electrical response because the electronic density and the disorder realisation in each ring is not exactly the same from one run to the other .
nevertheless due to the fact that we are measuring an ensemble average quantity the change in disorder realisation of each ring does not modify the result of the experiment .
moreover we have checked ( on the electric response measurement ) that , for the same illumination procedure , the result varies within a 15 % range from one run to the other .
the flux dependent orbital magnetism at a frequency of 350 mhz is then detected . in this configuration
the quality factor of the resonator is only 500 .
we do not understand this strong increase of magnetic losses .
this low quality factor decreases the accuracy of our measurements of the resonance frequency .
moreover it prevents precise measurements of the flux dependence of the dissipative part of magnetic response of the rings . as a consequence in this part we present only the flux dependent non dissipative part of the magnetic response .
note that we can not rule out the fact that the signal measured in this configuration could be partially due to electric response of the rings .
however due to the small value of the residual capacitance of the meander line and the bad electric coupling in this geometry this electric component is estimated to be at least 20 times smaller than when the rings are placed on top of the capacitance .
moreover the very different shape of the electric and magnetic signals is a strong evidence that we are indeed measuring essentially the magnetic response of the rings .
the signal measured at zero illumination , after subtracting the base line due to the resonator , is shown in the lower part of figure [ fig : magneticsignal ] .
inspired by our previous analysis we decompose the measured field dependent part of the signal into an aperiodic and a periodic part which corresponds to @xmath6 in a ring ( figure [ fig : magneticsignal ] ) .
we interpret the @xmath6 component as the contribution of electronic trajectories enclosing the whole ring . on the other hand the triangular shape signal
could be due to the contribution of trajectories confined in the finite width of the ring .
the amplitude of the @xmath6 periodic component of the signal is @xmath142 .
we deduce from formula [ dfmagnetic ] and evaluation of the magnetic coupling done in appendix [ magnetic_coupling ] that the flux dependent magnetic response of the ring is @xmath143 m@xmath144 . in the following
we first assume that the main contribution to this signal is due to the flux derivative of the persistent currents @xcite and then discuss finite frequency effects .
if the flux dependence of persistent currents is @xmath145 , we deduce  : @xmath146 with @xmath147 the surface of the ring .
we find then a _ diamagnetic _ average persistent current , the amplitude of which is : @xmath148 na .
the aperiodic component of the signal corresponds on the other hand to low field paramagnetism .
let s now compare our result for the average persistent currents to other experimental results and to theoretical predictions .
a @xmath6 periodic diamagnetic persistent currents has been also observed in arrays of metallic rings @xcite .
the expected amplitude of the averaged current due to repulsive electron - electron interactions from first order hartree - fock calculation @xcite is @xmath149 na , this value is expected to be decreased by higher order terms .
considering on the other hand theoretical predictions for non interacting electrons @xcite the expected value is between @xmath150 na and @xmath151 na . in both cases the currents are predicted to be paramagnetic .
the rather small difference between interacting and non interacting electrons is very specific to the gaas rings where the number of electrons is small .
the measured signal is consistent for the amplitude but not for the sign ( unless assuming attractive interactions ) with theoretical predictions .
it may also be important to take into account an effect of frequency for the flux dependent orbital magnetism .
in fact by applying a formalism very similar to the one used for magnetopolarisability the variation of the real part ( non - dissipative ) of the susceptibility of a ring submitted to an oscillating magnetic flux in ce without interactions is given by @xcite  : @xmath152 with @xmath153 the matrix element of the current operator .
it is then possible to apply the same reasoning as for the real part of polarisability , and to use the fact that @xmath154 , so that  : @xmath155          \right)\end{gathered}\ ] ] the evaluation of this expression is easily deduced from the evaluation of magnetopolarisability and leads to figure [ fig : magneticce ] .
frequency induces a strong decrease of the magnetic signal for frequencies of the order of @xmath156 but does not seem to induce a sign change of @xmath157 .
note that in strong contrast with the electric response the magnetic susceptibility is maximum at zero frequency , which corresponds to persistent currents .
it would be important to investigate the effect of finite frequency on the contribution due to electron - electron interactions on persistent currents .
recently it has been suggested that the measured currents may be due to a rectifying behaviour of the rings : a high frequency noise leads then to a dc current @xcite .
noise also induces dephasing .
a recent paper by kravtsov and altshuler @xcite predicts that those two quantities , average persistent current and dephasing measured by the saturation of the phase coherence time , are related in a simple way @xmath158 .
@xmath159 is a constant giving the sign of the persistent current and @xmath160 the dephasing time at zero temperature . using the value of @xmath161 ns at 20 mk , deduced from measurements on connected sample , and considering the orthogonal case ( absences of spin orbit , then @xmath162 ) we deduce an expected value for persistent currents of -0.14 na .
the predicted persistent current is then _
diamagnetic_. the sign and amplitude are then consistent with our experimental findings . on the other hand if we take the value deduced from the temperature dependance of the magnetopolarisability of non - connected rings , which is not the case considered theoretically , we deduce @xmath163 na , smaller by an order of magnitude than the experimental value
.      the influence of electronic density on the magnetic response of the rings has been investigated by illuminating them .
different illumination times are shown in figure [ fig : magneticillum ] .
we observed that the triangular envelope of the signal changes sign and increases with illumination . for each illumination
the fourier transform of the signal exhibits a component which is consistent with half a flux quantum periodicity ( figure [ fig : magneticillumfft ] ) .
one sees however that with illumination the shape of the @xmath6 peak in the fourier transform is modified .
the peak broadens with electronic density indicating that the width of the rings increases .
note that this width is always consistent with the one deduced from etching and depletion effects .
the power of the fourier transform integrated in the @xmath6 zone is constant within 10 % .
so the amplitude of the @xmath6 signal is constant but its shape is modified .
it indicates that the amplitude of the persistent currents does not depend much on electronic density .
the sign change of the low frequency part of the flux dependent magnetic response of the rings , going from low field paramagnetism to diamagnetism , is not understood .
we have presented measurements of electric and magnetic responses of aharonov - bohm rings etched in a 2deg .
they present a flux dependent correction to screening .
this correction is positive in low field which means that screening is enhanced when time reversal symmetry is broken by a magnetic field .
the sign of the effect is consistent with theory for isolated rings at finite frequency .
the value of magnetopolarisability is @xmath164 , with @xmath165 the calculated polarisability of a quasi-1d circular ring of radius r. the temperature dependence of magnetopolarisability is consistent with @xmath166 .
the behaviour versus electronic density is compatible with a @xmath8 dependence of magnetopolarisability .
the magnetic response has been measured on the very same array of rings .
the rings exhibits a signal consistent with _
diamagnetic _ average persistent currents of amplitude @xmath167 na .    because the measurements are done on the same rings it is possible to compare the electric and the magnetic signal .
the experimental ratio between the frequency shift due to the electric or magnetic response is around 10 , a value consistent with theoretical expectations taking into account the electric and magnetic coupling coefficient ( appendix [ magnetic_coupling ] and [ electric_coupling ] ) and the ratio between the typical matrix element of the screened potential and the current operator which leads to @xcite : @xmath168 with @xmath169 the vacuum impedance , @xmath170 the drude conductance and @xmath171 the fine structure constant .
we have thus shown that the mesoscopic electromagnetic response of gaas rings is dominated by the flux dependent polarisability instead of orbital magnetism .
this would not be the case in metallic rings , where , due to the very short screening length , the mesoscopic electric response is negligible .
the low field diamagnetic sign of the orbital magnetism needs further investigation both on the experimental and theoretical sides .
we thank b. etienne for the fabrication of the heterojunction . we acknowledge fruitful discussions with l.p .
lvy , ya.m .
blanter , s. guron and g. montambaux and technical help of m. nardone and p. demianozuck .
in this appendix we evaluate the magnetic coupling of one square ring with the resonator in the configuration of the experiment . the inductance @xmath25 is modelled by two cylindrical wires separated by a distance of @xmath172 ( see figure [ fig : res_anneaux ] ( a ) ) .
a ring is submitted to the magnetic field of those wires .
let s first evaluate the mutual inductance @xmath173 between the ring and the resonator . using ampere theorem the magnetic field generated by a current @xmath174 in the the inductance is easily calculated .
the flux of this magnetic field through a ring of size @xmath175 located at the ( 0,0 ) point is then  : @xmath176 if the ring is located at a point ( x , y ) @xmath173 reads  : @xmath177 the ring submitted to a magnetic field @xmath178 acts as a magnetic dipole @xmath179 .
this dipole is equivalent to the ring with a current @xmath180 , so that the flux in the inductance is now @xmath181 .
we deduce from these results that @xmath182 with @xmath183 the inductance of the meander - line . from the resonance frequency and the calculation of the capacitance we deduce @xmath184h , .
the rings are not perfectly well coupled to the inductance so that they are not all located at @xmath185 .
moreover because of the mylar sheet inserted between the rings and the resonator the rings are not in the plane of the resonator . to take this into account the inductance is averaged over the x position of the rings and @xmath186m@xmath187 m . within this
approximations : @xmath188 m@xmath144 .
is determined by the polarisability of the rings and the electric field generated by the capacitance at the ring position .
, width=302 ]    in this appendix we evaluate the electric coupling coefficient @xmath189 of one ring with the capacitance @xmath26 of the resonator .
the capacitance is modelled by two cylindrical wires of radius @xmath24 and separated by a distance @xmath172 , one wire with a linear charge of @xmath190 , the other one @xmath191 .
the electric field outside the wires is the one generated by two lines of lineic charge @xmath190 and @xmath191 separated by a distance @xmath192 determined by @xmath193 @xcite . in our case @xmath194 pf . using gauss theorem
we can easily calculate the electric field in the plane of the rings in every point @xmath195 outside the wires  : @xmath196 a ring submitted to this field generates an electric dipole @xmath197 with @xmath12 the polarisability of one ring , so that the rings submitted to electric field act as an ensemble of dipole .
we model them by two line of lineic charge @xmath82 and @xmath198 separated by a distance @xmath175 , and such that @xmath199 ( see figure [ fig : res_anneaux ] ) .
note that to do so the electric field has to be constant on the scale of the rings : this is roughly the case .
evaluating the potential @xmath200 created by these two wires between each side of the capacitance , and using the relation @xmath201 we have for rings located at ( x , y )  : @xmath202 to have the capacitance shift induced by one ring we have to divide the previous result by the number of rings @xmath203 with @xmath204 the length of the capacitance .
moreover as the ring are imbedded in gaas - algaas we have to divide our result by the dielectric constant of the substrate @xmath205 .
we can now evaluate the electric coupling coefficient @xmath35 defined by @xmath206 by averaging over the x - position of the rings and considering that the rings are located between 1.5 @xmath23 m and 2.5 @xmath23 m in the y - direction from the resonator . within these approximations
@xmath207 m@xmath144 .
note that the previous result is very close to the value of the magnetic coupling coefficient @xmath46 .
in this part we evaluate the magnetopolarisability given by formula [ polardiffuson ] for a quasi-1d ring .
the diffusion propagator at frequency @xmath32 is given by  : @xmath208 @xmath147 and @xmath156 are respectively the surface of the ring and the mean level spacing . @xmath209 and @xmath210 are the eigenvalues and eigenvectors of the diffusion equation  : @xmath211 we consider a 2d ring of perimeter @xmath18 , radius @xmath67 and width @xmath212 , with @xmath213 . in this case
the solutions of the diffusion equation are  : @xmath214 with @xmath215 and @xmath216 .
the modes corresponding to @xmath217 are given by  : @xmath218 @xmath219 is the coordinate along the ring , and @xmath220 the radial coordinate . in our description
the ring corresponds to @xmath221 $ ] .
we consider reflecting border in the @xmath220 direction.the corresponding eigenvalue is  : @xmath222\ ] ] @xmath223 is the thouless energy .
the mean charge density ( average over the width of the ring ) in the ring submitted to an electric field @xmath79 is  : @xmath224 note that using this density we recover the classical result for a quasi-1d ring @xmath165 . in the thomas - fermi approximation
we deduce the mean screened potential  : @xmath225 using this relation and the formula for the diffuson one can do the calculation analytically . because of the form of @xmath78 only the mode ( m=0,n=1 ) remains and leads to  : @xmath226 with @xmath227 .


the energy loss of hard partons propagating through the soft medium created in heavy - ion collisions has long been regarded as a promising tool to gain information on the medium properties @xcite .
radiative energy loss , i.e. the idea that medium - induced radiation predominantly carries away energy from a hard parent parton has been rather successful in describing not only the hadronic nuclear suppression factor for central collisions but also the effects of changing medium geometry @xcite .
calculations within dynamical evolution models in various formalisms @xcite have improved on this result and show also agreement with measured hard back - to - back correlations @xcite and the measured suppression of protons @xcite .    nevertheless , there are indications that a radiative energy loss picture fails to describe the suppression of heavy quarks as seen in the measurements of the single - electron spectra @xcite .
energy loss due to elastic collisions with the medium @xcite has been suggested as a possible solution to this problem .
such calculations indicate a large component of elastic energy loss also for light quarks and gluons .
however , the pathlength dependence of elastic energy loss does not agree with the observed response of the suppression to a large pathlength bias as probed in back - to - back hadron correlations @xcite .    in the light of these inconclusive findings
, it would be desirable to have a baseline calculation of elastic energy loss which treats the relevant physics in as detailed manner as possible .
in particular , many analytical calculations of elastic energy loss assume eikonal propagation of a hard parton , however in a more realistic description even hard partons , @xmath2  gev , can be deflected from their path by interactions with the medium . similarly , in processes like @xmath3 there is not only a change of energy between incoming and outgoing parton but also a change in parton species which influences the subsequent interaction probability of the outgoing parton .
finally , in scattering processes in which the energy is distributed almost evenly among the outgoing partons , it is not justified to use the picture of energy from a hard parton being lost into the medium . clearly , one also has to account for situations in which both partons after a scattering process are hard or both become part of the medium , or in which the identity of the leading parton changes .    in the present manuscript , we present work towards such a baseline computation . for this purpose
, we develop a monte carlo ( mc ) simulation for the hard parton s interaction with the medium , as this is an appropriate framework to account for all the above mentioned effects .
in addition , we also try to to give estimates on the uncertainties of such a computation .
similar , in some aspects even more ambitious approaches are the mc models jewel ( jet evolution with energy loss ) @xcite and martini ( modular algorithm for relativistic treatment of heavy ion interactions ) @xcite , both of which consider pqcd interactions of a parton shower with a thermal medium and include both elastic and radiative energy - loss components . in jewel ,
one models shower evolution and medium - induced radiation in a static ( no flow , @xmath4 constant ) medium , and in martini the vacuum showering is assumed to cease when the hydrodynamic evolution starts , while our approach focuses on the leading partons . like martini , also we utilize a realistic well - tested 3-dimensional hydrodynamical description of the medium ( in our case azimuthal symmetry in central collisions , and longitudinal boost symmetry but 3-d flow ) instead of a static medium applied in jewel .
in the following , we make the assumption that the medium can be characterized as a thermal gas of quasi - free quarks and gluons and that its interactions with a hard parton can be treated as incoherent @xmath5 processes in leading - order pqcd . we stress that these assumptions need not be realized for the medium produced in ultrarelativistic heavy - ion collisions  any discrepancy between the calculation presented here and the data will reflect this .
we model the medium by means of a hydrodynamical calculation @xcite with the initial state obtained from the saturation of the phase space with minijets @xcite .
this provides the medium temperature @xmath4 at each spacetime point . for the moment , we assume that only the scatterings with a partonic medium are relevant and ignore the possibility that a hadronic medium might interact with a hard parton
. thus the mean free path of hard partons becomes infinite once the medium temperature drops below the critical temperature @xmath6 .
given the local @xmath4 , the distribution of scattering partners is given by the fermi - dirac ( bose - einstein ) distribution for quarks ( gluons ) respectively .
we also account for the 3-d hydrodynamical flow of the fluid cell .
our approach , most specifically the assumption that the relevant interactions can be described in terms of partonic pqcd interaction , shares many features with parton cascade models ( e.g. @xcite ) .
however , unlike in parton cascade models where the partonic pqcd interactions are assumed to describe the whole off - equilibrium bulk medium , we assume that the bulk is essentially thermalized and use explicit pqcd interactions only to account for the non - equilibrium dynamics between hard parton and medium .    in the following , we first give a more detailed account of the hydrodynamical simulation , the parton - medium scattering rates and details of the mc implementation .
then we present our results for the partonic and hadronic nuclear modification factor @xmath7 and for qcd - matter tomography , and also discuss the relevant uncertainties . some technical details regarding the scattering rate calculations
can be found in the appendix .      in the present exploratory stage
we consider only the cleanest possible case , ( nearly ) central @xmath8 collisions in the mid - rapidity region . the produced qcd matter , which evolves and flows 3-dimensionally , can to a first realistic approximation then be described by using the nonviscous , azimuthally symmetric , longitudinally boost - invariant hydrodynamic framework .
the 1 + 1 d hydrodynamical framework , which we apply here for the qcd medium evolution , is presented in detail in ref .
the initial energy density @xmath9 for the hydrodynamics is here obtained from pqcd - minijet production @xcite and final - state saturation as computed in the ekrt - model @xcite .
the transverse profile of @xmath10 is assumed to scale with the number of binary collisions , and thermalization is assumed to take place at formation ( i.e. at saturation ) .
for central au+au collisions at @xmath11  gev , we thus initialize the hydrodynamics at ( see table 1 in @xcite ) @xmath12  fm/@xmath13 , where the longitudinal proper time is defined as @xmath14 .
the transverse flow is assumed to be zero at @xmath15 .
the equation of state ( eos ) in our hydrodynamic background contains a high-@xmath4 qgp phase , a low-@xmath4 hadron resonance gas ( hrg ) phase , and a mixed phase between them .
the qgp here is an ideal gas of gluons and three flavours of massless quarks and antiquarks obeying a bag - eos , while the hrg is an ideal gas of all hadrons and hadron resonances up to @xmath16  gev .
the maxwell construction gives the mixed phase , and the bag constant @xmath17  mev fixes the phase transition temperature to @xmath18  mev .
we terminate the hydrodynamic evolution at a single decoupling temperature @xmath19  mev . as discussed in ref .
@xcite , the measured pion and kaon spectra and proton multiplicity in central au+au collisions at rhic can be reproduced very nicely with such hydrodynamical set - up , once the thermal hadron spectra have been computed via the cooper - frye decoupling procedure and once the strong and electromagnetic hadronic decays have been accounted for . for more details , as well as for plots of the initial energy density and of the phase boundary locations in the @xmath20 plane , we refer the reader to ref .
@xcite .
our simulation of energy losses of high - energy partons in the produced qcd matter rests largely on one particular physical quantity , namely the scattering rate for a high - energy parton of a type @xmath21 , @xmath22 where all possible elastic and inelastic partonic processes @xmath23 are accounted for by summing over all types of collision partners , @xmath24 in the initial state , and over all possible parton type pairs @xmath25 in the final state . in general , the scattering rate depends on the frame , and in particular on the high - energy parton s 4-momentum @xmath26 , on the flow 4-velocity @xmath27 and on the temperature @xmath28 of the fluid at each space - time location @xmath29 . in the notation
below , we leave these dependences implicit but in the simulation they are fully accounted for . ignoring pauli blocking and bose enhancement
, we express the scattering rate as follows ( see also ref .
@xcite )    @xmath30    where @xmath31 is the spin- and colour - summed / averaged scattering amplitude which depends on the standard mandelstam variables @xmath32 , @xmath33 and @xmath34 .
the distribution function for a thermal particle of momentum @xmath35 is denoted by @xmath36 , which is the bose - einstein distribution @xmath37 for gluons and the fermi - dirac distribution @xmath38 for quarks .
since the initial state spins and colours have been averaged over in the squared matrix elements , the spin- and color - degeneracy factors in @xmath39 and @xmath40 are the standard ones , @xmath41 for gluons , and @xmath42 for quarks and antiquarks , while the possible flavour - degeneracy needs to be determined separately in each process . in our simulation
, we evaluate the scattering rate always in the local rest - frame of the fluid , so that @xmath43 and @xmath44 is the energy of the high - energy parton @xmath21 in this frame .
the singularities arising in forward and backward scatterings of massless partons  appearing in the scattering amplitudes at @xmath45  need to be regularized in some physically meaningful way . for the moment , an exact field - theoretical solution to this non - equilibrium problem is not known .
one possible treatment is to apply resummed thermal ( equilibrium ) propagators separately for the longitudinal and transverse components , as is done in the amy approach @xcite .
since the problem is , however , still phenomenological , we choose here a more transparent approach ( tested also in jewel before ) and introduce a thermal - mass - like overall cut - off scale @xmath46 for the momentum phase space , such that @xmath47 . here
@xmath48 is the strong coupling constant and @xmath49 is a parameter of the order of one but which we can vary .
as we can not control the running of @xmath50 , we keep the strong coupling constant fixed with momentum scale but vary both @xmath48 and @xmath49 in different runs of the mc code to test the model sensitivity to these parameters . since @xmath51 ,
the cut - off conditions are equivalent to @xmath52 with the requirement @xmath53 .
these are encoded in the function @xmath54 above , which is defined as @xmath55    for the evaluation of the lorentz - invariant part in eq .  , it is useful to recall that the scattering rate can be expressed in terms of the scattering cross sections . for the process @xmath23 , where the lorentz - invariant cut - off conditions are imposed , and the identities of the final state particles are accounted for , the cross section is @xmath56 and we can rewrite eq .   as @xmath57 where @xmath58 is the thermal particle density in the rest frame of the qcd fluid , and where we have expressed the mandelstam variable as @xmath59 , i.e. in terms of the angle @xmath60 between the hard parton and the thermal parton , and their energies in the fluid rest frame .
changing the integration variable @xmath61 to @xmath32 , it is straightforward to arrive at the following result @xmath62 where @xmath63.\ ] ] the analytical results for the relevant elements of @xmath64 are listed in the appendix .
the remaining @xmath65 integral is easiest to do numerically .
the obtained scattering rates @xmath66 for quarks , antiquarks and gluons , along with the contributions @xmath67 from different processes , are shown in fig .
[ fig : rates_vs_t ] as a function of temperature @xmath4 and in fig .
[ fig : rates_vs_e1 ] as a function of the hard parton energy @xmath44 .     of gluon ( above ) and
quark ( below ) as a function of temperature @xmath4 for a qcd plasma at rest .
flavour- and quark - antiquark -summed contributions from different processes and the analytical estimates for the total rates discussed in the text are also shown ( solid lines).,title="fig:",width=340 ]   of gluon ( above ) and quark ( below ) as a function of temperature @xmath4 for a qcd plasma at rest .
flavour- and quark - antiquark -summed contributions from different processes and the analytical estimates for the total rates discussed in the text are also shown ( solid lines).,title="fig:",width=340 ]     of gluon ( above ) and quark ( below ) as a function of the parton energy @xmath44 , for a qcd plasma at rest at @xmath68 gev .
flavour - summed contributions from different processes are shown .
the analytical estimates for @xmath66 ( see text ) are shown by the solid lines . in the upper panel ,
the analytical estimate for the @xmath69 process ( dashed line ) is the same as for the quark - antiquark summed @xmath70 process . in the lower panel ,
the dashed lines show the contributions from the @xmath71 process and an individual @xmath72 process .
, width=340 ]     of gluon ( above ) and quark ( below ) as a function of the parton energy @xmath44 , for a qcd plasma at rest at @xmath68 gev .
flavour - summed contributions from different processes are shown .
the analytical estimates for @xmath66 ( see text ) are shown by the solid lines . in the upper panel ,
the analytical estimate for the @xmath69 process ( dashed line ) is the same as for the quark - antiquark summed @xmath70 process . in the lower panel ,
the dashed lines show the contributions from the @xmath71 process and an individual @xmath72 process .
, width=340 ]    to get a feeling of the general characteristics of the system , we consider also the following analytical estimates . in the process @xmath23
, the average of a quantity @xmath73 in the plasma at a temperature @xmath4 can be computed in the fluid rest frame as @xmath74    at the high - energy limit , @xmath75 , the relevant pqcd cross sections all behave as @xmath76 , and @xmath77 ( see appendix , take the limit @xmath78 , and @xmath79 ) . in the boltzmann approximation
we thus get @xmath80 from which we see that , in this limit , @xmath81 ) .
the last result in eq .
indicates that the collisions of the high - energy parton with the medium partons happen on the average in an angle @xmath82 .    using the same approximation ,
we can easily obtain analytical estimates also for the scattering rates , @xmath83 \nonumber \\ & \approx & ( g_g + 2\cdot3\cdot g_q\cdot{\textstyle \frac{4}{9 } } ) \frac{9}{2}\frac{\alpha_s^2t^3}{\pi m^2 } = 144\frac{\alpha_s^2t^3}{\pi m^2}\\ \label{analytic_gammaq } \gamma_{q_i } & \approx &   \gamma_{{q_i}g\rightarrow { q_i}g } +    \sum_{j}[\gamma_{{q_i}{q_j}\rightarrow { q_i}{q_j } } + \gamma_{{q_i}{\bar q_j}\rightarrow { q_i}{\bar q_j}}]\nonumber\\ & \approx & ( g_g\cdot{\textstyle \frac{4}{9 } } + 2\cdot3 \cdot g_q\cdot({\textstyle \frac{4}{9}})^2)\frac{9}{2}\frac{\alpha_s^2t^3}{\pi m^2 } =   64 \frac{\alpha_s^2t^3}{\pi m^2}\end{aligned}\ ] ] where the factors @xmath84 account for quarks and antiquarks times the number of quark flavours .
these estimates are shown in figs .
[ fig : rates_vs_t ] and fig .
[ fig : rates_vs_e1 ] , together with the actual rates applied in our simulation ( once the boost to the local rest frame has been taken into account ) .
we observe that the analytical rates describe the numerical ones ( both the total rate and the ones from the main processes ) quite well in the kinematic and temperature range of rhic .    in a static homogeneous qgp ,
the average number of collisions experienced by the high - energy parton of an average path length @xmath85 , would , in the high - energy limit , be @xmath86 .
taking @xmath87  fm and using eqs .
- with @xmath88 and @xmath89 , we arrive at the range 145 collisions for gluons and 62 collisions for quarks at @xmath90  mev .
this gives us a simple order - of - magnitude estimate of the expected number of elastic scatterings for a hard parton .      to initiate the hard massless parton of a type @xmath21 in each event
, we sample the lo pqcd single - jet production spectrum ( for more details , see @xcite ) at @xmath91 and at @xmath92 , @xmath93 \frac{1}{1+\delta_{kl } } , \end{split}\ ] ] where @xmath94 with @xmath95 containing the squared matrix elements listed in eqs .
( a-1)-(a-8 ) .
for the parton distribution functions ( pdfs ) , we use the cteq6l1 set @xcite .
for simplicity in this exploratory work we , however , neglect the nuclear effects to the pdfs @xcite , since these are small in comparison with the ones arising from the final state interactions with the medium . the initial rapidity @xmath96 is randomly generated in the range @xmath97 $ ] from a flat distribution .
( this is consistent on the level of the longitudinal boost - symmetry approximation made for the hydrodynamical part . )
this fixes the hard - parton energy @xmath98 and polar angle @xmath99 of its momentum vector .
the azimuth angle @xmath100 is evenly distributed between @xmath101 $ ] .    in order to keep track of the hard partons interactions with the medium in the coordinate space , the initial spatial position of the hard parton has to be fixed as well .
we start to follow the hard parton at the initial longitudinal proper time @xmath15 of our hydrodynamical model . since in the c.m .
frame all hard partons are produced in the lorentz - contracted overlap region at @xmath102 , the longitudinal position at later times ( before the first collision at @xmath103 ) is assumed to be determined by the longitudinal momentum only .
the initial time and longitudinal coordinates for the hard parton are thus @xmath104 and @xmath105 .
the initial position in the transverse plane is sampled from the nuclear overlap function @xmath106 where @xmath107 is the impact parameter . as we focus only on the central collisions , in the following @xmath108 .
the nuclear thickness function @xmath109 is defined as usual , @xmath110 where for the nuclear density @xmath111 we use the standard woods - saxon distribution @xmath112 with @xmath113 , @xmath114 , @xmath115  fm and @xmath116  @xmath117 .
the hard parton propagates through the plasma in small time steps @xmath118 , during which we propagate the parton in position space .
since we consider only incoherent scatterings here , the probability for not colliding in the time interval @xmath118 is assumed to be given by the poisson distribution @xmath119 where @xmath66 is the total scattering rate for the hard parton of the type @xmath21 .
hence the probability to collide at least once during the time @xmath118 is @xmath120 . for small enough @xmath118
we can assume that there will be at most one collision .
however , as explained above , we calculate the scattering rates in the local rest frame of the quark - gluon plasma fluid element and we must boost also the time step @xmath118 to the same frame , @xmath121 where the lorentz factor @xmath122 appears , @xmath123 being the 3-dimensional flow velocity vector of the hydrodynamical plasma in the c.m .
frame of the colliding nuclei .
should a scattering happen , the probability @xmath124 for a given type of scattering process is determined by the partial scattering rates as @xmath125    we produce the 4-momentum @xmath35 of the scattering partner from the thermal medium according the energy and angular distributions defined by the differential scattering rate @xmath126 ( see eq .  ) of the process .
we then boost the 4-momenta of the hard parton @xmath26 and the thermal particle @xmath35 to the c.m .
frame of the partonic collision , rotating the frame so that @xmath26 is on the @xmath127-axis .
it is now easy to determine the scattering angle @xmath128 by sampling the distribution @xmath129 using @xmath130 , and sampling the azimuthal angle @xmath131 from a flat distribution .    after determining the scattering angles @xmath128 and @xmath131
, we can construct the final state 4-momenta @xmath132 and @xmath133 and lorentz transform these back to the original c.m .
frame of the colliding nuclei .
the final state parton with highest energy is then chosen as the new hard parton to be propagated further , for which we repeat the procedure outlined above with the next timestep .
note that the type of the hard parton may change here .
( we also note that in our code , we have the possibility to follow the recoil partons as well but at sufficiently high parton momenta ( @xmath2  gev ) , we do not see a significant effect on the observables studied here . )
this procedure continues until the surrounding matter has cooled down to the critical temperature @xmath6 . at @xmath134
the system is in this work in a mixed phase where energy density @xmath135 keeps decreasing while the temperature stays constant .
it is _ a priori _ not clear how the interaction with the medium in the mixed phase should be modelled properly in pqcd .
however , we will demonstrate _ a posteriori _
( see fig .
[ raapartonqgp ] ahead ) that within the assumption of an effective temperature , the effect of the mixed phase is small because the medium density is small as compared to the initial density . thus , unless there is an anomalous enhancement of parton - medium scattering in the mixed phase by non - perturbative processes , a detailed treatment of the mixed phase is not important , and we can also ignore the interactions of the hard parton in the hadron gas . in the following ,
we treat the interactions in the mixed phase by using an effective temperature @xmath136 where @xmath137 is the quark - gluon plasma degrees of freedom and @xmath138 is the bag constant .
our simulation ends when @xmath139 and the system enters the hadron gas phase , where we assume no significant interaction between parton and medium .    the outcome of the procedure described above is a medium - modified distribution of high - energy partons , @xmath140 . in order to compare with observables
, we have to hadronize this distribution . in the current leading - particle case , this can be done by convoluting the obtained partonic distribution with the fragmentation function @xmath141 : @xmath142 where @xmath143 is the fraction of the final parton momentum @xmath144 available to the hadron with momentum @xmath1 , and @xmath145 is the fragmentation scale . in this work
we have used the kniehl - kramer - ptter ( kkp ) fragmentation functions @xcite .
in the following , we choose a lower bound of the parton transverse momentum , @xmath146 , as 5.0 gev to account for the observation that at low @xmath1 the main hadron production mechanism is not independent fragmentation .
partons which fall below @xmath147 in our simulation are assumed to be thermalized and become part of the medium , and we do not follow these partons anymore .
we average all observables across the rapidity window @xmath148 $ ] which corresponds to the phenix acceptance . to make sure we account properly for all the possible partons falling into this final rapidity window , we choose the initial @xmath149 .
we refer to all partons above both the @xmath146 cut and within the final rapidity acceptance as punch - through partons.     is shown by the dashed lines.,width=340 ]    to verify the magnitude of the rapidity shifts experienced by the high-@xmath98 partons , we plot in fig .
[ fig : rapidity ] the initial rapidity distribution of the punch - through partons in our default set - up with @xmath88 and @xmath150 .
we can see that about 13 % of the punch - through partons originate from outside of the final rapidity interval .
on the one hand , this suggests that the eikonal approximation  where all partons would remain within the interval @xmath151  is a fairly good first approximation . on the other hand ,
if one aims at a precision analysis of the single - hadron distributions or , even more importantly , at correlations of two high-@xmath1 hadrons , the non - eikonal dynamics should be obviously taken into account . an attempt to quantitatively understand the rapidity shift seen in fig .
[ fig : rapidity ] can be found in appendix  [ sec : yshifts ] .
also the leading - parton flavour - conversions are usually neglected in the eikonal energy - loss simulations .
interestingly , this effect is not symmetric between quarks and gluons : due the collision rates , there are more flavour conversions related to gluons than to quarks , and a quark from a @xmath152 leading - parton conversion is then more likely to become a punch - through parton than a gluon from a @xmath153 conversion .
our simulation shows that for a simplified homogeneous static plasma cylinder of , say , radius @xmath87  fm and @xmath154  mev , as much as 40 % of those original high-@xmath144 gluons which will become punch - through partons , will turn into quarks . in the full simulation with 3-d hydrodynamical background medium , however , only about 4% of all punch - through partons have undergone this effect .
we thus conclude that in simulations for single - hadron distributions at rhic , where a realistically evolving background medium is considered , the leading - parton flavour - conversions can be neglected .
nuclear modification factor both at hadron and parton level with strong coupling constant @xmath88 .
phenix data is from @xcite.,width=340 ]    in fig .
[ raahadron ] we present the nuclear modification factor @xmath155 , @xmath156 for @xmath88 , @xmath150 , both on the hadronic level and on the partonic level , and compare it with the @xmath157 data from phenix @xcite .
it is evident that for this default choice of @xmath158 and @xmath49 elastic collisions are not enough to account for the data but that they should not be neglected , either .
this result is in agreement with the calculation of purely elastic energy loss in the arnold - moore - yaffe ( amy ) formalism averaged over hydrodynamics @xcite where a similar value of @xmath159 is found for @xmath160 .
note that a comparison with models not using a constrained fluid - dynamics model for the bulk evolution or the same @xmath158 is not meaningful .
it can further be deduced from the figure that since the suppression as a function of @xmath1 ( or partonic @xmath144 ) is flat , the convolution of the medium - modified parton spectrum with the fragmentation function does not cause a significant effect in @xmath155 . we will thus in the following focus on the partonic suppression pattern , with the implication that the partonic suppression provides a reasonable estimate for the suppression also on the level of hadrons . as a technical point
, we note that in order to collect enough statistics also at high-@xmath144 , the partonic results in fig .
[ raahadron ] consist of three dedicated runs , corresponding to @xmath161  gev .     and
@xmath150 , including only the scatterings in the qgp phase ( squares ) and in the qgp and mixed phase but not the in hadron gas ( circles).,width=340 ]    in fig .
[ raapartonqgp ] we demonstrate , as discussed above , that the mixed phase is in this set - up indeed causing only a small additional suppression and , consequently , that it is justified to neglect the hadron gas effects in the present study .
the results shown in this and in the following figures are from runs with @xmath162  gev .
for the regulator - mass parameter value @xmath163 and strong coupling constant values @xmath164.,width=340 ]    in fig .
[ raapartonalfa ] we show the partonic nuclear modification for three different choices of constant @xmath158 .
as expected , we observe a rather strong dependence on the value of @xmath158 .
the behavior is quite straightforward ; doubling the coupling approximately doubles the relative suppression . with @xmath165
we find a suppression due to elastic energy loss already below the data . on the other hand , with @xmath166
the effect is rather weak .
we thus focus on @xmath88 in the following when examining other model parameters .     for mass parameter values @xmath167 and strong coupling constant @xmath168.,width=340 ]
as is obvious the other significant model parameter in our simulation is @xmath49 , controlling the momentum cutoff scale @xmath169 which regulates the scattering cross sections at smallest and largest angles . like @xmath158
, it has a direct effect on the scattering cross sections of a hard parton with medium partons .
however , it could be expected to influence energy loss in a different way .
while the total cross section can become very large if small - angle forward scattering processes are taken into account , such processes do not transfer a large amount of energy and are hence not efficient for energy loss .
thus , a smaller cutoff mass is expected to increase the total cross section without affecting large - angle scattering whereas a change in @xmath158 affects both forward and large - angle scattering in equal measure .    a comparison of different values of the cutoff mass in fig
.  [ raapartonmass ] shows , however , that this parameter has an unexpectedly large effect on the nuclear modification factor . the suppression increases with decreasing mass ; this indicates that the cutoff mass is not yet in the region where forward scattering , which is inefficient for energy loss , takes place .
( above ) @xmath168 ( middle ) and @xmath170 ( below ) .
all particles are moving to the positive @xmath29-direction.,width=321 ]     ( above ) @xmath168 ( middle ) and @xmath170 ( below ) .
all particles are moving to the positive @xmath29-direction.,width=321 ]     ( above ) @xmath168 ( middle ) and @xmath170 ( below ) .
all particles are moving to the positive @xmath29-direction.,width=321 ]    in fig .
[ xyinit ] we show the geometry underlying the suppression in terms of the initial spatial distribution of punch - through partons for different values of @xmath158 .
this figure can be considered as a tomography image , where we have rotated the position of each particle in transverse @xmath171-plane in such a way that they all move to the same direction , to the direction of positive @xmath29-axis . with @xmath172 ( upper panel ) , the quark - gluon plasma is quite transparent in transverse direction as far as elastic collisions are concerned , thus the initial distribution of punch - through partons is centered around the origin .
doubling the value of @xmath158 to 0.3 ( middle panel ) makes the distribution move off - center .
thus , the medium introduces a surface bias which means that hard partons from the center can not punch through the matter anymore .
overall , the plasma becomes more opaque , as is expected from fig .
[ raapartonalfa ] .
the opaqueness is emphasized even more in the bottom panel , which corresponds to @xmath165 , which reproduces , or slightly overestimates the suppression measured at rhic .
we have also checked the tomography images for different mass parameters ( @xmath49 ) which reproduce the same suppression pattern in @xmath155 as when @xmath158 is varied ( see fig .
[ raapartonmass ] ) .
our conclusion here is that these do not significantly differ from the ones we have shown in fig .
[ xyinit ] .
the results above indicate on the one hand that for our default choice of parameters @xmath168 and @xmath173 , elastic energy loss is not sufficient to account for the observed value of @xmath155 . on the other hand , we see that at least in the framework which is based on individual parton - parton scatterings instead of , e.g. , static scattering centers , the elastic energy losses give a non - negligible contribution to @xmath155 . this has also been observed by a comparable calculation within the amy formalism @xcite in which a combination of radiative and elastic energy loss was needed to account for the data .
however , our results show an unfortunately strong dependence on cutoff parameters which are in essence supposed to effectively account for the physics of a dynamically running @xmath158 and thermal field - theoretical processes generating a screening mass .
thus for a reasonable variation of these parameters there are pairs @xmath174 for which the @xmath155 data can be described by elastic processes alone .
it is thus evident from these findings that a comparison based on the magnitude of the single hadron suppression alone can not resolve the question how significant the contribution of @xmath0 processes to energy loss actually is . to improve upon this problem ,
the next step is to include @xmath175 processes and coherence effects into the simulation .
however , a comparison of the tomographic pictures in this paper ( fig .
[ xyinit ] ) with the results of a radiative energy loss calculation @xcite indicates that the way hard partons probe the medium geometry through energy loss is quite different for radiative and elastic energy loss .
thus , as outlined in a schematical way in @xcite it can be expected that pathlength - dependent observables , like the suppression as a function of the angle of the observed hadron with the reaction plane , or the suppression of back - to - back correlations , will show pronounced differences between elastic and radiative energy loss .
the mc code introduced here provides an excellent tool to explore this idea , and this will be addressed in a future publication .
the spin- and color summed / averaged squared parton scattering amplitudes @xmath176 for gluons and quarks are ( see e.g. @xcite ) @xmath177 inserting the expression for the cross section into the definition of the @xmath178-function we have @xmath179 .
\end{split}\ ] ] several amplitudes contain identical terms .
we list the integrals for each of these terms separately , as it is then easy to combine the appropriate terms to form the final expressions for @xmath180 .
@xmath181 \nonumber \\
\hline \begin{array}{c } 1   \end{array } & 4e_1 e_2 -2m^2
\log{\frac{2e_1 e_2}{m^2}}-2m^2 \quad \\ \begin{array}{c } \frac{ut}{s^2 } \end{array } & \frac{2}{3 } e_1 e_2-\frac{3}{4}m^2+\frac{m^4}{4e_1 e_2}-\frac{m^6}{48e_1 ^ 2 e_2 ^ 2 } \quad \\ \begin{array}{c } \frac{su}{t^2},\frac{st}{u^2 } \end{array } &   -\frac{8e_1 ^ 2 e_2 ^ 2}{m^2 } + 4e_1 e_2 \log{\left(\frac{4e_1 e_2}{m^2}-1 \right ) } + 2m^2 \quad \\
\begin{array}{c } \frac{t^2+u^2}{s^2 } \end{array } &   \frac{8}{3 } e_1 e_2 -2m^2 \log{\frac{2e_1 e_2}{m^2 } }   -\frac{1}{2}m^2-\frac{m^4}{2e_1 e_2}+\frac{m^6}{24e_1
^ 2 e_2 ^ 2 } \quad \\ \begin{array}{c } \frac{t^2+u^2}{tu } \end{array } &   2(4e_1 e_2-m^2 ) \log{\left(\frac{4e_1 e_2}{m^2 } -1 \right ) } - 16e_1 e_2        + 4m^2\log{\frac{2e_1 e_2}{m^2}}+8m^2 \quad \\ \begin{array}{c } \frac{s^2+u^2}{t^2},\frac{s^2+t^2}{u^2 } \end{array } &   \frac{16e_1 ^ 2e_2 ^ 2}{m^2}-8e_1 e_2 \log{\left(\frac{4e_1 e_2}{m^2}-1 \right ) } + 4e_1 e_2 -2m^2 \log{\frac{2e_1 e_2}{m^2}}-6m^2 \quad \\ \begin{array}{c } \frac{s^2+u^2}{su } \end{array } &    -(4e_1 e_2-m^2 ) \log{\left(\frac{4e_1 e_2}{m^2}-1\right ) } + 2e_1 e_2   + m^2 \log{\frac{2e_1 e_2}{m^2}}-m^2 \quad \\ \begin{array}{c } \frac{s^2}{tu } \end{array } & 2(4e_1 e_2-m^2 ) \log{\left(\frac{4e_1 e_2}{m^2}-1\right ) } -8e_1 e_2   + 4m^2 \quad \\
\begin{array}{c } \frac{u^2}{st } \end{array } & -(4e_1 e_2-m^2 ) \log{\left(\frac{4e_1 e_2}{m^2}-1\right ) } + 10e_1 e_2 -3m^2\log{\frac{2e_1 e_2}{m^2}}-5m^2    \quad \end{aligned}\ ] ]
here we try to understand quantitatively the magnitude of the rapidity shifts of punch - through partons shown in fig .
[ fig : rapidity ] . at the high - energy limit @xmath182 , a straightforward calculation in the static plasma approximation ( see sec .
[ subsec : rates ] ) gives the average scattering angle in the cms frame as @xmath183 we set the initial high-@xmath98 parton at zero rapidity and its momentum @xmath184 to point into the direction of the @xmath29-axis , while the beam direction is , as usual , the @xmath127-axis .
having the maximum rapidity shift in mind , we take the medium - parton s momentum @xmath185 to reside in the @xmath186-plane , its energy to be @xmath187 and the collision ( polar ) angle relative to @xmath184 such that @xmath188 , which correspond to the averages in eq .
[ averages ] .
thus @xmath189    in a reverse boost from the cms frame of the parton - parton scattering , the energy and momentum of the scattered high - e parton become @xmath190 where the velocity of the cms frame is @xmath191 , @xmath192 is the lorentz gamma factor and the asterisks denote the cms - frame quantities .
the momentum of the scattered high-@xmath98 parton in the cms - frame can be computed to be @xmath193 \nonumber\\ & + & { \bf \hat e_y}e^*\sin\theta_3^*\cos\phi_3^ * \\ & + & { \bf \hat e_z}[p_{1x}^*\sin\theta_3^*\sin\phi_3^ * + p_{1z}^*\cos\theta_3^ * ] ,   \nonumber\end{aligned}\ ] ] where @xmath194 is the scattering ( polar ) angle relative to @xmath195 and @xmath196 is the azimuthal angle around @xmath195 , measured from the @xmath197-axis , and where the cms - frame quantities are obtained in a boost to @xmath198 as @xmath199    the point here is that although in the cms frame ( where @xmath200 ) the scattering cross section does not depend on the angle @xmath201 , a @xmath201 dependence is generated in the boost back to the plasma rest frame .
the rapidity shift @xmath202 is now easy to obtain numerically from the components @xmath203 and @xmath204 and we find that to a good approximation @xmath205 for a hard parton of , say , @xmath206  gev in a @xmath154  mev plasma , we obtain @xmath207 .
the maximum shifts @xmath208 occur at @xmath209 and @xmath210 , i.e. when the scattering remains in the ( @xmath211)-plane .
similarly , we find that @xmath212 and @xmath213 .    to estimate the maximum possible total rapidity shift @xmath214 induced by @xmath215 scatterings
, we may assume that at the limit of small scattering angles ( small rapidity shifts ) the maximum - shift distribution @xmath216 holds in each scattering .
the rapidity shifts are additive and we can write the probability distribution for @xmath217 as follows : @xmath218 where the probability distribution for a rapidity shift in one collision , @xmath219 where @xmath220 , is obtained from the corresponding flat azimuthal angle distribution @xmath221 by a change of variables according to eq .
[ deltay ] .
now it is easy to obtain @xmath222    for an order of magnitude estimate , we take @xmath223 ( from our simulation , since for this static high-@xmath4 plasma test case the leading - parton flavour - coversions are in fact a non - negligible effect ) , and arrive at the average maximum total rapidity shift @xmath224 .
figure  [ fig : yshifts ] shows our simulation test results for the rapidity spread of punch - through partons which originate from 10 gev partons sent off from the origin into the @xmath29-direction .
the width of the distribution in the static 500 mev plasma case ( squares ) is seen to be about 0.3 , which nicely confirms the order of magnitude of our upper limit estimate @xmath214 .
the smaller width in the simulation follows from the fact that the simulation accounts for scatterings ( corresponding to a continuum of azimuthal incoming angles @xmath225 around @xmath184 ) , where the rapidity shift effects are smaller than in our estimate above .
changing the background into a hydrodynamically evolving more realistic medium ( circles ) reduces the width of the rapidity shift by about a factor 3 .
this is caused by the fact that collision rates , scattering angles @xmath194 and rapidity kicks all decrease with decreasing temperature . finally , in a full simulation , where the high-@xmath98 partons can be produced anywhere in the transverse plane , and where also the surface bias effect is present ( see fig .  [ xyinit ] )
, the rapidity shifts become still somewhat smaller . from these considerations , and realizing that the rapidity spread distributions should be applied to each final - rapidity bin of fig .
[ fig : rapidity ] , we conclude that the degree of non - eikonality ( @xmath197-shifts ) is understood , and also that if we wish to fully account for punch - through partons in the finite rapidity interval @xmath151 , we have to consider an initial window of @xmath226 .
we thank harri niemi for providing us with the hydrodynamic simulations .
j.a . gratefully acknowledges the grant from the jenny and antti wihuri foundation .
this work was also financially supported by the academy of finland : k.j.e.s projects 115262 and 133005 , and t.r.s academy research fellowship with the associated project 130472 .
k.  j.  eskola , k.  kajantie and j.  lindfors , nucl .
b * 323 * , 37 ( 1989 ) .
k.  j.  eskola , k.  kajantie , p.  v.  ruuskanen and k.  tuominen , nucl .
b * 570 * , 379 ( 2000 ) .
t.  matsui , b.  svetitsky and l.  d.  mclerran , phys .
d * 34 * , 783 ( 1986 ) [ erratum - ibid .
d * 37 * , 844 ( 1988 ) ] .