shock waves have been observed in a wide range of environments outside the earth , from the solar wind to the hot gas in galaxy clusters .
however , the mechanism whereby the gas in these environments is shocked has been poorly understood . while shock transitions in the earth s atmosphere are mediated by molecular viscosity ( and hence direct particle collisions ) , those in interstellar space and the solar wind are too dilute to form in this way . in non - relativistic shocks , the role of collisions is effectively played by collective interactions of the plasma with the magnetic field .
this results in a multi - scale shock transition having sub - structure at ion kinetic length scales ( larmor radius or inertial length ) and potentially electron kinetic scales ( inertial lengths or whistler mode ) ( e.g. , schwartz et al .
( this volume ) , treumann 2009 ) .
such plasmas are termed collisionless .
the magnetic fields threading through the charged particle plasmas in space endow the plasmas with elastic properties , much like a fluid .
the kinetic energy of the inflowing gas is dissipated within this fluid via collective interactions between the particles and magnetic field , transferring energy from the magnetic field to the particles .
the collective processes are the result of the dc electromagnetic fields present in the shock transition layer , kinematic phase mixing , and also plasma instabilities ; the last give rise to a rich range of plasma waves and turbulent interactions .
it has long been known that these processes may heat the electrons beyond the mass - proportional value predicted by the rankine - hugoniot jump conditions - ample evidence is found in spacecraft studies of solar wind shocks ( schwartz et al .
1988 ) and multi wavelength spectroscopy of supernova remnants ( ghavamian et al .
2001 , 2002 , 2003 , 2007 ; laming et al .
1996 ; rakowski et al .
2008 ) and galaxy cluster gas ( markevitch et al . 2005 ; markevitch & vikhlinin 2007 ; russell et al .
2012 ) understanding how this process depends on such shock parameters as shock speed , preshock magnetic field orientation and plasma beta has been slow .    in collisionless plasmas ,
the downstream state of the plasma can not be uniquely determined from the upstream parameters because the rankine - hugoniot jump conditions only predict the _ total _ pressure downstream , not the individual contributions from the electrons and ions : @xmath2 . at the limit of a strong shock ,
@xmath3 and @xmath4 are each 4 times their preshock values , so the relative values of @xmath0and @xmath5immediately behind the shock are wholly dependent upon the nature of the collisionless heating processes occurring at the shock transition .
although an mhd description can be used to describe the behavior of the gas far upstream and far downstream of the shock , a more detailed kinetic approach is required for understanding how the dissipation at the shock front transfers energy from plasma waves and turbulence to the electrons and ions .
non - relativistic collisionless shocks can be broadly sorted into three categories : slow , intermediate and fast .
the three types are defined according to the angle between the shock velocity and upstream magnetic field , as well as the relative value of the shock speed compared to the upstream sound speed ( @xmath6 ) and alfvn speed ( @xmath7 ) .
most astrophysical shocks are quasi - perpendicular ( i.e. , they propagate at a nearly right angle to the preshock magnetic field ) , allowing only for fast - mode propagation . in that case , the relevant quantity is the magnetosonic mach number , @xmath8 ( @xmath9 ) .
collisionless shocks have also been classified according to whether the flow speed exceeds the sound speed in the downstream plasma ( kennel 1985 ) . above the critical mach number where the flow is subsonic , the dissipation of flow energy into thermal energy
can no longer be maintained by electrical resistivity , and plasma wave turbulence ( the cause of which are instabilities generated when the electron and ion distribution functions are distorted at the shock transition ) is required ( kennel 1985 ) .
shocks above the critical mach number are termed supercritical , while those below are termed subcritical .
note that even for subcritical shocks , observations suggests that kinetic processes other than resistivity and turbulence contribute to the shock dissipation ( greenstadt and mellott 1987).there are also other , higher critical mach numbers related to the formation of subshocks ( kennel et al .
1985 ) and non - steady cyclic shock reformation beyond the whistler critical mach number ( krasnoselskikh et al . 2002 ) .
the critical mach number for quasi - perpendicular shocks is estimated to be rather low @xmath102.8 ( edmiston and kennel 1984 ) .
the mach numbers of most snr shocks are expected to be well in excess of this value , meaning that they are both fast - mode and supercritical .    as one approaches the high mach numbers expected for astrophysical shocks ( @xmath1120 - 100 for an ambient magnetic field of 3 @xmath12 ) , a greater and greater fraction of the incoming ions become reflected back upstream from the shock front , corresponding to an increasingly turbulent and and disordered shock transition .
physical parameters such as @xmath13 , @xmath14 and @xmath15 no longer jump in an ordered manner ( i.e. , the transitions are no longer laminar ) .
in addition to the reflected particles , the hot ions from downstream become hot enough to escape upstream , further enhancing the population of ions in front of the shock .
these ions , which form a precursor , are now believed to play an essential role in the dissipation of high mach number collisionless shocks . aside from providing the seed population for the acceleration of cosmic rays , the precursor ions are likely generate a variety of plasma waves capable of selectively heating the electrons over the ions , thereby providing an important mechanism for raising @xmath16 above the mass - proportional value of @xmath17 .
taken at face value , the rankine - hugoniot jump conditions predict that electrons and ions will be heated in proportion to their masses :    @xmath18    in collisionless shocks , there are three relevant timescales to consider : the time required for coulomb collisions to isotropize a distribution of electrons , @xmath19 , the time required for coulomb collisions to isotropize distribution of ions , @xmath20 , and finally the time required for the electrons and ions to equilibrate to a common temperature , @xmath21 ( spitzer 1962 ) . after a time scale @xmath19 or @xmath20 , the particles in question attain a maxwellian velocity distribution .
the self - collision time , @xmath22 for electrons of density @xmath3 and temperature @xmath0 is ( spitzer 1962 ) :    @xmath23    where ln@xmath24 is the coulomb logarithm , equation [ massprop ] has been used to write the relaxation time in terms of the shock speed and @xmath25 is the shock speed in units of 1000 kms@xmath26 . for a young snr having @xmath271000 kms@xmath26 , postshock density n@xmath281 @xmath29 and ln @xmath3030
the time required to establish a maxwellian distribution at the electron temperature given by equation 1 is @xmath31 0.004 yrs .
the electrons are isotropized by self - collisions first , then the protons , and finally over a longer timescale the electrons and protons equilibrate to a common temperature .
for this reason , coulomb collisions alone are unable to establish equilibration electron and proton distributions at the shock front .
this equilibration is described by the relation    @xmath32    where @xmath33    the temperatures @xmath0 and @xmath5 equilibrate to a common density - weighted average temperature @xmath34 , given by @xmath35 , where @xmath36 is the mean mass per particle ( @xmath37 for cosmic abundances ) . for mass - proportional heating ,
the timescale given by equation [ tep ] is @xmath382000 yrs for @xmath39 @xmath401000 kms@xmath26 , of similar order but longer than the proton - proton isotropization timescale , @xmath41 .
these are longer than the age of the snr , substantially so at higher @xmath39 indicating that for minimal heating ( i.e. , @xmath16@xmath421/1836 ) the electrons and ions will not equilibrate to @xmath34 during the lifetime of the snr ( itoh 1978 ; draine & mckee 1993 ) .    the arguments above indicate that coulomb collisions are ineffective at both isotropizing the heavy ion particle distributions and equilibrating the electron and ion temperatures at the transtions of collisionless shocks .
the emission spectra of non - radiative snrs ( dominated mostly by x - ray and ultraviolet emission ) should therefore be sensitive probes of the collisionless heating processes at the shock transition .
we consider the observational constraints of these processes below .
the most useful shocks for studying collisionless equilibration processes are those exhibiting detectable emission from the immediate postshock gas . to be diagnostically useful
, the emission should arise from the region close to the shock front , where temperature disequilibrium between electrons , protons and heavy ions is substantial enough to affect both the relatives fluxes and relative velocity widths of emission lines .
the shape and extent of the spatial profile for the different emission species behind the shock is also a useful diagnostic , especially for the uv resonance lines of he  ii @xmath431640 , c  iv @xmath441548 , 1550 , n  v @xmath441238,1243 , and o  vi @xmath431032 , 1038 . for a given shock speed , the distance behind the shock where the emission peaks depends strongly on the initial electron temperature , and electron temperature immediately behind the shock .
the ubiquity of non - radiative snrs , as well as their relatively simple geometry and very high shock speeds , makes these objects the most important laboratories for investigating the efficiency and nature of electron - ion and ion - ion equilibration .
other non - radiative shocks available for study are those occurring in stellar wind bubbles ( for example , wolf - rayet bubbles ) ( gosset et al .
2011 ) and in galaxy cluster shocks ( e.g. , the bullet cluster ( markevitch et al . 2005 ) and abell 2146 ( russell et al .
however , these shocks are far less frequently observed than those in snrs , leaving the latter as the most important objects providing both a broad range of shock speeds and corresponding diagnostic information ( such as proper motions and spatially resolved structure ) .
shocks in star forming regions , such as hh objects and their associated bow shocks , tend to be clumpy and complicated , and are usually located in dense environments ( @xmath45100 - 1000 @xmath29 ) with low enough shock speeds ( @xmath39@xmath46100 - 200 kms@xmath26 ) to be radiative . in those cases ,
the optical and uv emission ( where the most valuable shock diagnostic line emission arises ) is dominated by emission from the cooling and recombination zones far downstream from the shock . in those regions
the cooling of the gas below 100,000k and the accompanying compression of the gas result in a collisional plasma with @xmath47 .
this erases any memory of the initial electron - ion equilibration .
furthermore , a significant fraction of the ly @xmath48 continuum produced in the recomination zone is expected to pass upstream and ionize the preshock gas . in strong shocks
( @xmath39@xmath4980 kms@xmath26 ) , this results in complete ionization of hydrogen ( shull & mckee 1979 ; cox & raymond 1985 ; dopita et al .
1993 ) , thus precluding the use of collisionally excited balmer line emission from neutral h ( described in the next section ) as a temperature equilibration diagnostic .      in the late 1970s
it was discovered that very fast ( @xmath102000 kms@xmath26 ) shocks in young snrs could generate detectable optical emission very close to the shock transition ( chevalier & raymond 1978 ; chevalier , kirshner & raymond 1980 ) , providing a valuable diagnostic of physical conditions at the shock before coulomb collisions or cooling could alter them .
this emission is produced by collisional excitation of h  i as it flows into the shock front .
the cold neutral component does not interact directly with the plasma waves and turbulence at the shock , while the ionized component is strongly heated and compressed by a factor of four ( when the shock is strong ) .
some of the cold h is destroyed by collisional ionization ; however , the rest of the cold h undergoes charge exchange with hot ions behind the shock , generating a separate population of hot h. approximately 1 in every 5 collisions results in collisional excitation to the n=3 level of h , producing h@xmath48 and ly @xmath52 emission .
the h@xmath48 line from the cold neutrals is narrow and reflects the preshock temperature ( @xmath5330,000  k ) , while that from the hot neutrals is broad ( typically @xmath54500 kms@xmath26 ) , and reflects the postchock temperature ( and to a large extent , the velocity distribution ) of the protons ( chevalier & raymond 1978 ; chevalier , kirshner & raymond 1980 ; smith et al .
1991 ; ghavamian et al .
2001 ) ( figure [ fig : balmerprofile ] ) . in balmer - dominated shocks , the broad to narrow h@xmath48 flux ratio is proportional to the ratio of the charge exchange rate to the ionization rate , with the latter being highly sensitive to the electron and proton temperatures .
this makes the broad to narrow ratio , @xmath55 , very sensitive to @xmath16 .
there is only weak dependence of @xmath55 on the preshock h  i fraction and preshock temperature , mainly due to differences in the amount of ly @xmath52 converted into h@xmath48 in the narrow component ( ghavamian et al . 2001 ; van adelsberg et al .
2008 ) .
-0.2 in    the first systematic attempt to use balmer - dominated snrs to infer @xmath16 for collisionless shocks was attempted by ghavamian et al .
( 2001 ) . using the broad h@xmath48 line widths
measured from a sample of balmer - dominated shocks , they constrained the range of plausible shock speeds between the limits of minimal equilibration and full equilibration .
they then predicted the broad - to - narrow ratios for a grid of shock models over this range of @xmath39 and @xmath16 , allowing @xmath16 to be constrained .
the range of broad component h@xmath48 widths observed in snrs ranges from @xmath10250 kms@xmath26 for the slowest balmer - dominated shocks ( cgynus loop ) , to @xmath10 500 kms@xmath26 for intermediate - velocity shocks ( rcw 86 ) and finally @xmath102600 kms@xmath26 for the fastest shocks ( snr 0509@xmath5667.5 ) .
this corresponds to a well - sampled range of shock speeds : nearly a factor of 10 . the primary uncertainty in measurement of the broad component width at low shock speeds ( @xmath57200 kms@xmath26 , as seen in the northeastern cygnus loop ; hester et al .
1994 ) is disentangling the broad and narrow components when they are of comparable width . at high shock speeds ( @xmath582000 kms@xmath26 )
the main difficulty is the baseline uncertainty of the surrounding continuum : if the peak of the broad line is low and the width very large , errors in ascertaining where the broad line merges into the background can lead to underestimates of the broad component width .
the range of @xmath55 for this sample of balmer - dominated shocks lies between 0.4 and 1.2 ( kirshner , winkler & chevalier 1987 ; smith et al .
1991 ; ghavamian et al . 2001 ; 2003 , rakowski , ghavamian & laming 2009 ) .
however , it does not not vary with broad component width in a monotonic fashion .
recently there has been substantial improvement in the modeling of balmer - dominated shocks .
earlier calculations of the broad h@xmath48 line profile assumed that it formed from a single charge exchange ( chevalier , kirshner & raymond 1980 ; smith et al .
1991 ; ghavamian et al .
2001 ) , and it treated the hot and cold neutrals as two separate , distinct populations , with a given neutral belonging to either one or the other .
however , in reality an interaction tree is required to track the number of photons emitted by each neutral over multiple excitations and charge exchanges .
these effects were first incorporated in the balmer - dominated shock models of heng & mccray ( 2007 ) , who also found that charge exchange results in a third population of neutrals having velocity widths intermediate between the hot and cold neutrals .
further improvements in modeling of balmer - dominated shocks were included by van adelsberg et al .
( 2008 ) , who included the momentum transferred by charge exchange between the hot neutrals and protons .
this allowed the bulk velocity of the postshock neutrals to be calculated separately from those of the protons .
inclusion of the momentum transfer showed that for shock speeds @xmath571000 kms@xmath26 , charge exchange effectively couples the fast neutral and thermal proton distributions , while for high shock speeds ( @xmath585000 kms@xmath26 ) , it does so far less effectively .
this results in a fast neutral distribution that is skewed relative to the protons in velocity space and an average velocity that is much higher for the fast neutrals than the protons ( van adelsberg et al .
together , the inclusion of all these effects has enhanced the ability of the models to match the observed broad - to - narrow ratios ( and hence predict @xmath16 ) .
in particular , the newer models can now match the low broad - to - narrow ratio ( @xmath55@xmath590.67 ) observed in knot g of tycho s snr , yielding @xmath16@xmath590.05 , @xmath39@xmath591600 kms@xmath26 .
however , even after all the additional physics is included , the balmer - dominated shock models are still unable to reach the low broad - to - narrow ratios measured along the rims of dem l71 ( @xmath55@xmath590.2 - 0.7 ; ghavamian et al .
2003 ; rakowski , ghavamian & laming 2009 ) .
the electron - ion equilibration in dem l 71 was instead determined via comparison of broad h@xmath48 line fwhm with postshock electron temperatures measured from _
chandra _ observations ( rakowski , ghavamian & hughes 2003 ) .
the most promising explanation advanced for the @xmath55 discrepancy has been added narrow component flux from the shock precursor ( raymond et al . 2011 ; morlino et al .
2012 ) , hitherto not included in the earlier shock models ( ghavamian et al .
2003 , rakowski , ghavamian & laming 2009 ) .
these developments are described in more detail in the next section .
the plot of @xmath16 versus @xmath39for the available sample of balmer - dominated shocks shows a declining trend of equilibration with shock speed ( ghavamian et al . 2007 ; heng et al .
2007 ; van adelsberg et al . 2008 ) .
the trend is described by ghavamian et al .
( 2007 ) as full equilibration for shock speeds up to and including 400 kms@xmath26 , and a declining equilibration proportional to the inverse square of the shock speed above 400 kms@xmath26 .
this description can be characterized in the following way : @xmath60 where the functional form of the the @xmath16 relation is designed to asymptotically transition to @xmath16=@xmath17 at very high shock velocities .
the most up to date plot of @xmath16 versus shock speed , reproduced from van adelsberg et al .
( 2008 ) , is shown in figure  [ fig : plotequil ] ( note that the shock models used in producing these plots do not include contribution from the shock precursor ) .
although the inverse correlation between @xmath16 and @xmath39was largely confirmed by van adelsberg et al ( 2008 ) , there may be some evidence of departure from the @xmath16@xmath50 @xmath61 relation at shock speeds exceeding 2000 kms@xmath26 .
van adelsberg et al .
find that when all three measured broad component widths and broad - to - narrow ratios from sn 1006 are included in the plot ( @xmath39@xmath102200 - 2500 kms@xmath26 ) , a slight upturn in the @xmath16-@xmath39 relation appears . the @xmath16 ratios for those cases
are found to be @xmath100.03 , superficially similar to @xmath62 , rather than @xmath63 .
the reason for the discrepancy is not clear , nor whether this indicates a breakdown in the @xmath61 dependence at shock speeds exceeding 2000 kms@xmath26 . some caveats to consider when interpreting the upturn in @xmath16 seen in figure [ fig : plotequil ] are as follows : for shock speeds @xmath582000 kms@xmath26 collisional ionization and excitation of h are primarily caused by proton ( and to a lesser degree , alpha particle ) impact ( laming et al . 1996 ; ghavamian 1999 ; ghavamian et al . 2001 ; tseliakhovich et al .
experimentally measured cross sections for these interactions have still not been available to high precision ( uncertainties @xmath1020%-30% still exist ) , although more sophisticated theoretical calculations are now becoming available ( see , for example , tseliakhovich et al .
in addition , as the broad component width increases , the h@xmath48 profiles are spread out over an increasing number of pixels , resulting in noisier spectra and greater measurement uncertainty in the broad component width . these larger error bars in turn result in a larger uncertainty in @xmath39 , especially at shock speeds of 2000 kms@xmath26 and higher .
the best way to further constrain the equilibration - shock speed relation is to add new data points to the curve shown in figure  [ fig : plotequil ] . to this end
, we present balmer - dominated h@xmath48 profiles for two additional positions in tycho s snr ( different from those of knot g presented by kirshner , winkler & chevalier ( 1987 ) , smith et al .
( 1991 ) and ghavamian et al .
( 2001 ) ) .
these profiles , shown in figure  [ fig : tychocennwspec ] , were acquired with a moderate resolution spectrograph in 1998 ( for details on the observational setup , see ghavamian 1999 and ghavamian et al .
the profile marked ne was obtained from a clump of h@xmath48 emission located along the northeastern edge of tycho s snr , approximately 1@xmath64 northward of knot g. the clump appears behind the main body of the balmer filaments and exhibits a broad component that is substantially doppler shifted to the red ( 10.7@xmath651.2  , or about 490 kms@xmath26 ) .
the doppler shift reflects the bulk velocity of the hot posthock proton distribution , so the significant velocity shift of the broad component centroid indicates that the shock in the ne has a substantial velocity component into the plane of the sky , i.e. , that the ne shock is located on the far side of the blast wave shell .
the broad h@xmath48 width of the ne shock is 1300@xmath6565 kms@xmath26 ( the smallest broad component width measured in tycho s snr so far ) , with @xmath55=0.85@xmath650.04 .
the nw shock , on the other hand , has a broad h@xmath48 width of 2040@xmath6555 kms@xmath26 ( the largest broad component width measured in this snr so far ) , with @xmath55=0.45@xmath650.15 .
neither of these broad - to - narrow ratios is strictly reproduced by the latest models of van adelsberg et al .
( 2008 ) , with the lowest predicted ratios being for case b ( the assumption of optically thick conditions for ly @xmath52 photons in the narrow component ) and for low equilibrations ( @xmath16@xmath570.1 ) . for these low equilibrations , the corresponding shock speeds for the ne and nw shocks in tycho are approximately 1400 kms@xmath26 and 2250 kms@xmath26 ( using figures 5 and 10 of van adelsberg et al .
2008 ) .
in figure  [ fig : plotequil ] we have added data from the ne and nw shocks in tycho s snr to the @xmath16-@xmath39 plot .
the two data points help fill in a portion of the plot where the data are sparse : the region between approximately 1200 kms@xmath26 and 1500 kms@xmath26 , as well as the region beyond 2000 kms@xmath26 , where the existing data are taken entirely from sn 1006 .
the added point near 1400 kms@xmath26 is fully consistent with the @xmath61 relation , while @xmath16 is not well constrained enough for the point at 2250 kms@xmath26 to contradict the appearance of an upturn at the highest shock velocities .
it is clear that proper characterization of the equilibration - shock velocity curve above 2000 kms@xmath26 will require both higher signal - to - noise spectra on existing balmer - dominated shocks , as well as new data points beyond a broad component width of 2500 kms@xmath26 .
although the inverse squared relation between equilibration and shock velocity appears has been the most salient result of the study of balmer - dominated shocks , there have been discrepant results reported in a small subset of cases . recently the broad h@xmath48 component in the lmc snr 0509@xmath5667.5 was detected for the first time ( helder , kosenko & vink 2010 ) .
broad emission was observed along both the northeastern rim ( fwhm 3900@xmath65800 kms@xmath26 ) and southwestern rim ( fwhm 2680@xmath6570 kms@xmath26 ) , with the former being the fastest balmer - dominated shock detected to date having the characteristic broad and narrow component h@xmath48 emission .
interestingly , the broad - to - narrow ratios for both shocks are exceptionally low , @xmath55=0.08@xmath650.02 and 0.29@xmath650.01 along the ne and sw rims , respectively .
as noted by helder , kosenko & vink ( 2010 ) these ratios are nearly twice as large as the smallest ones predicted by the models of van adelsberg et al .
( 2008 ) , precluding measurement of @xmath16 from the balmer - dominated spectra and implicating excess narrow component h@xmath48 emission in a cosmic ray precursor once again . using rgs spectra of snr 0509@xmath5667.5 acquired with _
xmm - newton _
, they obtained a forward shock speed of approximately 5000 kms@xmath26 in the sw .
this implies a broad h@xmath48 width of 3600 kms@xmath26 , substantially smaller than the measured width of 2680 kms@xmath26 .
helder , kosenko & vink ( 2010 ) suggest that this indicates some thermal energy loss ( @xmath1020% ) to cosmic ray acceleration .
this picture is supported by the presence of nonthermal x - ray emission in their fitted rgs spectra of snr 0509@xmath5667.5 .
the result described above relies upon an accurate disentangling of the bulk doppler broadening from the thermal doppler broadening in the x - ray lines .
the disentangling depends on parameters such as the ratio of reverse shock to forward shock velocities , as well as the ratio of the gradients in the reverse and forward shock velocities , which in turn had to be assumed from evolutionary snr models .
the result , while intriguing , is still significantly uncertain . on the other hand , helder , kosenko & vink ( 2010 ) found that the x - ray shock velocities from the ne could only be reconciled with the observed broad h@xmath48 width there if @xmath16@xmath660.2 , which clearly which predicts @xmath16@xmath67 predicted from equation  [ tetpeq ] .
if this result were to be confimed by future observations , it would present a new challenge in understanding how electron - ion equilibration occurs in fast collisionless shocks .
one possibility , given the presence of nonthermal x - ray emission in the spectra of 0509@xmath5667.5 ( warren & hughes 2004 ; helder , kosenko & vink 2010 ) is that the moderate loss of thermal energy to cosmic ray acceleration may have slightly increased the compression and reduced the temperature at the shock front compared to the case with no acceleration ( decourchelle , ellison & ballet 2000 ; ellison et al . 2007 ) .
both of these effects would tend to render the plasma more collisional , possibly explaining the @xmath16@xmath660.2 result .
however , it is also worth noting that the shock velocity used for obtaining @xmath16 in the ne is subject to the same model dependence and uncertainties as the sw measurements , so similar caution is required in its interpretation .
a similar combined optical and x - ray study of rcw 86 was performed by helder et al .
( 2009,2011 ) . there ,
the broad component h@xmath48 widths were supplemented with electron temperatures measured from _ xmm - newton _ rgs spectra from the same projected locations along the rim .
one of the main results of this study was that the slower shocks ( broad h@xmath48 fwhm @xmath10500 - 600 kms@xmath26 ) showed showed @xmath16@xmath101 , agreeing with earlier results from similar shocks observed both in rcw 86 and elsewhere ( ghavamian et al .
( 1999 , 2001 , 2007 ) ) . however , helder et al . ( 2011 ) found while the results were indeed consistent with low equilibration at the shock front for fast shocks ( @xmath16@xmath590.02 for broad fwhm @xmath101100 kms@xmath26 ) ) and higher equilibration for the slower shocks ( @xmath16@xmath591 for broad fwhm @xmath10650 kms@xmath26 ) , their x - ray derived electron temperatures were inconsistent with @xmath680.3 kev at the shock front contradicting the suggestion of ghavamian et al .
( 2007 ) that shocks above 400 kms@xmath26 may all heat electrons to roughly 0.3 kev . however , a major caveat of these results is that the forward shock in rcw 86 is believed to be impacting the walls of a wind - blown bubble ( williams et al .
2011 ) , resulting in substantial localized variations in shock speed around the rim .
these variations occur as different parts of the forward shock impact the cavity wall at different times .
while the broad component h@xmath48 widths closely trace the current position of the shock front , the x - ray emission behind that shock arises over a much more extended spatial scale , and is sensitive to the history of the forward shock interaction with the cavity wall .
furthermore , narrowband h@xmath48 imagery of rcw 86 with the eso very large telescope ( helder et al .
( 2009 , 2011 ) ) shows a complex morphology of filaments , especially along the eastern side of the snr .
the broad h@xmath48 components of these filaments exhibit substantial , localized variations in line width , ranging from @xmath10600 kms@xmath26 fwhm to 1100 kms@xmath26 ghavamian ( 1999 ) , helder et al .
( 2009 , 2011 ) ) .
these variations reflect localized changes in density and viewing geometry along the line of sight . as
such , uniquely mapping the observed balmer filaments to their corresponding x - ray emission in _ xmm - newton _ data ( especially given the somewhat coarse 10@xmath69 spatial resolution of that instrument ) is fraught with uncertainty .
additional corroboration for these results would be desirable .
perhaps the most important and exciting recent development in our understanding of balmer - dominated shocks has come with the development of new kinetic - based shock models .
blasi et al . (
2012 ) have introduced a kinetic model for following the momentum and energy exchange between neutrals and ions , along with the back - reaction of those neutrals when they pass back upstream and form a fast neutral precursor . rather than assume a maxwellian velocity distribution for the neutrals ( as had been done in previous models , despite the lack of thermal contact between neutrals needed to justify such an assumption ) , both the ion and neutral distributions are computed from their appropriate boltzmann equations . building on these models , morlino et al .
( 2012a , b ) have confirmed what had been suspected earlier ( smith et al . 1994 ; hester et al . 1994 ; sollerman et al . 2003 ) , namely that the broadening of the narrow component beyond the expected ism value ( @xmath1025 - 30 kms@xmath26 instead of 10 kms@xmath26 ) is most likey due to heating in a cosmic ray precursor . in particular ,
morlino et al .
( 2012 ) found that the characteristic charge exchange length of the incoming neutrals exceeds that of the neutrals crossing back upstream , so that the narrow component width is impacted not by the neutral return flux , but rather by heating in the cosmic ray precursor .
raymond et al . ( 2011 ) predicted a similar broadening of the narrow component , though they focused mainly on the contribution of collisional excitation in the precursor to the flux in the narrow h@xmath48 component .
spatially resolved line broadening of the narrow h@xmath48 component was detected in ground - based longslit spectroscopy of knot g in tycho ( lee et al .
in addition , a small ramp - up in h@xmath48 emission was observed ahead of knot g in hst imagery tycho s snr ( lee et al . 2010 ) .
the results , taken together , are strong evidence for the presence of cosmic ray precursors in balmer - dominated snrs .
one prediction of the new balmer - dominated shock models is the existence of a third component of the h@xmath48 emission ( morlino et al .
2012a , b ) . when the hot neutrals escape upstream , they undergo charge exchange with the colder preshock protons .
this results in fast protons with cold neutrals , with the former rapidly equilibrating with preshock protons and pre - heating the gas .
the temperature of the equilibrated protons in the precursor lies between the temperature of the far upstream protons ( @xmath105000 k ) and the far downstream protons ( @xmath10 10@xmath70 - 10@xmath71 k ) , typically @xmath1010@xmath72  k. further charge exchange between these warm protons and the preshock neutrals gives rise to a third , warm neutral component ( neither fast nor slow ) having velocity widths of hundreds of kms@xmath26(morlino et al .
interestingly , the presence of a third h@xmath48 component was first observationally reported by smith et al .
( 1994 ) in their high resolution echelle spectroscopy of balmer - dominated snrs in the large magellanic cloud ( an example of one of the spectra from smith et al .
( 1994 ) is reproduced in figure [ fig:0509_3rdhalphacomp ] ) .
a third h@xmath48 component was also reported in high resolution spectra of knot g in tycho s snr ( ghavamian et al . 2000 ) . in their models of balmer - dominated shock emission , morlino et al .
( 2012a ) found that the importance of the third component relative to that of the broad and narrow components depends strongly on the preshock neutral fraction and @xmath16 , in line with earlier theoretical predictions on properties of a fast neutral precursor ( smith et al .
the fact that the third component has been detected in tycho s snr ( width @xmath10150 kms@xmath26 ) is consistent with the high preshock neutral fraction ( @xmath730.9 ) inferred from the broad - to - narrow ratio of knot g by ghavamian et al .
a similar third component may have been detected in high resolution spectra of snr 0509@xmath5667.5 , where measurement of the narrow component width required the inclusion of an additiona component of width 75 kms@xmath26 ( smith et al .
1994 ) .
the new fast neutral precursor models predict that a substantial fraction of the h@xmath48 excitation in balmer - dominated shocks can arise ahead of the shock , where warm neutrals are excited by electron impact .
interestingly , the relative contribution of the preshock h@xmath48 to the total ( upstream + downstream ) is sensitive to @xmath16 behind the shock .
morlino et al .
( 2012a ) found that up to 40% of the total h@xmath48 flux from a balmer - dominated shock can arise from the fast neutral precursor when @xmath39@xmath102500 kms@xmath26 and @xmath16=1 both upstream and downstream . in these models the preshock contribution to the total flux drops substantially for lower downstream equilibrations for @xmath39@xmath581000 kms@xmath26 ( the slowest shocks considered by morlino et al .
this is generally consistent with the fact that in most cases , shock models not including the precursor h@xmath48 emission have been able to match the observed broad - to - narrow ratios .
in other words , if the postshock temperature equilibration were not low for such fast shocks , the agreement between the observed and predicted @xmath55would have been substantially worse for such remnants as tycho s snr and sn 1006 .
sn 1006 is an example of a snr accessible to uv spectroscopy due to its galactic location , 450 pc up from the galactic plane , and therefore with relatively low extinction due to intervening dust and gas .
the hopkins ultraviolet telescope ( hut ) observed the uv resonance lines of he  ii @xmath431640 , c  iv @xmath441548 , 1550 , n  v @xmath441238,1243 , and o  vi @xmath431032 , 1038 ( raymond et  al .
1995 ) emitted from the balmer dominated filament in the nw quadrant .
the line showed doppler broadening consistent with that of the h @xmath48 broad component observed in the optical , indicating insignificant ion - ion equilibration .
laming et  al .
( 1996 ) were also able to infer the degree of electron - ion equilibration at the shock .
-0.2 in    he  ii @xmath431640 , by virtue of its relatively high excitation potential @xmath10 48 ev ) , is excited only by electrons , and its intensity is therefore directly related to the electron temperature .
c  iv @xmath441548 , 1550 , n  v @xmath441238,1243 , and o  vi @xmath431032 , 1038 , by contrast , have much lower excitation potentials of @xmath74 and 12 ev , so although these ionization states are established by electrons , the line emission in these transitions can also be excited by impacts with hot protons and @xmath48 particles , and the intensity ratio of he  ii @xmath431640 to c  iv @xmath441548 , 1550 , n  v @xmath441238,1243 , and o  vi @xmath431032 , 1038 can be sensitive to the post - shock electron - proton temperature equilibration .
the spatial distribution of the uv resonance line emission , when spatially resolved , provides additional constraints on the degree of electron - ion equilibration at the shock front . for shocks slower than 1500 kms@xmath26 , @xmath75 at the shock front results in both a more rapid rise and higher maximum in emissivity of the uv resonance lines with distance behind the shock ( see figure  [ fig : ovispatial ] ) .    laming et  al .
( 1996 ) calculated impact excitation cross sections for protons and @xmath48 particles colliding with li - like ions , using a partial wave expansion with the coulomb - bethe approximation , and applying a unitarization procedure following seaton ( 1964 ) .
they found a degree of equilibration of order @xmath76 or less , which implied for a 2250 km s@xmath26 shock an electron temperature immediately postshock of @xmath77 k , in very good agreement ( and actually predating ) the optical results discussed above in subsection 4.1 .
the ion - ion equilibration in sn 1006 was revisited by korreck et  al.(2004 ) , using higher spectral resolution fuse data comprising o  vi @xmath431032 , 1038 and the broad ly  @xmath781025 emission lines .
they found a slightly broader line profile in ly @xmath52 , implying less than mass - proportional heating and possibly a small degree of ion - ion equilibration .
sn 1987a represents another sn / snr in a region of the sky accessible to uv observations .
hst cos observed the he  ii @xmath431640 , c  iv @xmath441548 , 1550 , n  v @xmath441238,1243 and n  iv @xmath441486 lines emitted from the reverse shock ( france et  al .
when combined with optical spectroscopy of h @xmath48 , the @xmath16 ratio at the shock is determined to be in the range 0.14  0.35 , significantly higher than similar ratios coming from balmer dominated forward shocks .
france et  al .
( 2011 ) argued that a different equilibration mechanism is likely at work .
considering the relative youth of sn 1987a , and the fact that the reverse shock is the origin of the emission , significant populations of cosmic rays and associated magnetic field amplification are unlikely .
in fact , in the expanding ejecta the magnetic field is likely to be very weak , leading to a very high aflvn mach number shock .
as will be discussed below in connection with shocks in galaxy clusters , electron heating in such a case is likely to be due to acceleration in the cross - shock potential .
the cross - shock potential is effective at heating electrons , and so may explain the higher @xmath16 in sn 1987a
.    the forward shock of sn 1987a has also been observed in x - rays with the grating instruments on chandra ( e.g. zhekov et  al .
2009 ) . in general electron heating well below complete equilibration is seen , though precise interpretation is difficult because one observation sees emission from shocks at a variety of different velocities , due to irregularities in the density of the surrounding medium .
the inverse relationship between the temperature equilibration and shock speed is an interesting result from studies of balmer - dominated snrs .
however , the applicability of this result to both fully ionized shocks and shocks undergoing efficient cr acceleration ( @xmath5850% of their energy transferred to crs ) remains unsettled .
recently , vink et al . (
2010 ) used a two - fluid - model for cosmic rays and thermal gas to simulate the effect of cosmic ray acceleration on the temperature and ionization structure of fast , non - relativistic shocks .
they found that if 5% of the shock energy were to be channeled into cosmic rays ( the minimum needed if snrs are the dominant source of cosmic rays ) then approximately 30% of the postshock pressure must reside in cosmic rays ( corresponding to a ratio of cosmic ray to total postshock pressure , w , of 0.3 ) .
for w=0.3 , vink et al .
( 2010 ) predicted a lowering of the average temperature of the postshock gas to @xmath1070% of the value when the cosmic ray contribution is ignored .
this is a significant alteration of the postshock temperature profile , and should result in much more rapid equilibration of electrons and protons close to the shock .
however , do the effects described above actually occur in snr shocks ?
one of the principle lines of evidence cited by vink et al .
( 2010 ) in support of this picture was the result found in rcw 86 by helder et al .
( 2009 ) . in that snr
, the broad h@xmath48 widths of balmer - dominated filaments were found to be nearly 50% smaller than the minimum allowed given their x - ray proper motions .
filaments in the ne of the snr exhibited broad h@xmath48 widths of 1000 kms@xmath26 , but their apparent x - ray counterparts , which exhibited strong x - ray synchrotron ( nonthermal ) emission , exhibited proper motions indicating shock speeds of 3000 - 6000 kms@xmath26 .
this result , along with the theoretical prediction that x - ray synchrotron emission requires shock speeds of at least 2000 kms@xmath26 ( aharonian et al .
1999 ) was taken by helder et al .
( 2009 ) and vink et al .
( 2010 ) as evidence for substantial energy loss ( w@xmath100.5 ) from the balmer - dominated shocks to cosmic rays .
however , this association has now been refuted by subsequent multi - epoch optical imagery of the h@xmath48 filaments , which have failed to show the kind of high proper motions seen in the nonthermal x - ray filaments ( helder et al .
2012 , in preparation ) . instead
, they show proper motions consistent with shock speeds predicted by the broad h@xmath48 widths without energy loss to cosmic rays(@xmath10600 - 1200 kms@xmath26 ) , implying that w@xmath790.2 .
the association between the balmer - dominated shocks studied spectroscopically by helder et al .
( 2009 ) and the x - ray filaments was due either to coincidental spatial alignment , or due to sudden deceleration of the outer shock in rcw 86 during its encounter with the surrounding cavity wall ( williams et al .
2011 ; helder et al .
2012 , in preparation ) .
-1.2 in    the lack of association between the balmer - dominated filaments and the non - thermal x - ray filaments in rcw 86 raises some important questions about the feasibility of using balmer - dominated shocks to study electron - proton equilibration in cases where 50% or more of the thermal energy is diverted to cosmic rays . in snrs such as sn 1006 ( koyama et al .
1996 ; katsuda et al .
2010a ) , tycho s snr ( warren et al . 2005 ; katsuda et al .
2010b ) , cas a ( vink & laming 2003 ) and rx j1713.7@xmath563946 ( koyama et al . 1997
; slane et al . 1999 ; tanaka et al . 2008 ) the presence of strong synchrotron x - ray filaments has been interpreted as evidence for highly efficient cosmic ray acceleration .
the narrowness of the synchrotron filaments most likely reflects the short emitting lifetimes of the ultra high energy electrons ( energies @xmath1010 - 100 tev ) as they spiral in the postshock magnetic field ( vink & laming 2003 ) .
the detection of @xmath80-ray emission from the shells of sn 1006 ( acero et al .
2010 ) and rx j1713.7@xmath563946 ( aharonian et al . 2006 ; abdo et al .
2011 ) has shown that cosmic rays are accelerated to energies as high as 100 tev in these snrs . in all cases thermal x - ray emission has been exceptionally faint due to the very low inferred preshock densities ( @xmath570.1 @xmath29 ) , making it more likely that the overall x - ray emission will be dominated by synchrotron radiation from the most energetic cosmic rays .
snrs expanding into such low density media can propagate at the high shock speeds required for cosmic ray acceleration ( @xmath582000 kms@xmath26 ; aharonian et al .
1999 ) for a longer time , allowing their structure to become modified by the back pressure from the cosmic rays .
investigating the temperature and ionization structure of such shocks with balmer line spectroscopy requires finding balmer - dominated shocks exhibiting x - ray synchrotron radiation .
all of the known snrs exhibiting balmer - dominated shocks have also been imaged at x - ray wavelengths with _
chandra _ or xmm , allowing reasonably detailed searches for shocks emitting both h@xmath48 and synchrotron x - ray emission ( the latter producing hard continuum that is dominant at energies of 2 kev and higher ) . a detailed comparison for all balmer - dominated snrs has not yet been published . however
, even a cursory comparison between the narrowband h@xmath48 and hard x - ray images of these snrs shows a distinct _ anticorrelation _ between shocks emitting in these two bands .
for example , overlaying h@xmath48 and _ chandra _ images of tycho s snr acquired during the same epoch ( 2007 ) ( figure [ fig : tycho_balmer ] ) shows little or no correlation between the prominent balmer - dominated filaments on the eastern and northeastern edges and the non - thermal x - ray filaments ( e@xmath543 kev ) circling the remnant .
the balmer - dominated filaments ( shown in red in figure [ fig : tycho_balmer ] ) on the eastern side of tycho s snr are seen projected 30@xmath69 - 1@xmath64 inside the edge of the nonthermal x - ray filaments ( marked in green ) , an indication that portion of the shell along the line of sight has significantly decelerated .
the balmer filaments are seen at the outermost edge of the thermally emitting x - ray ejecta ( marked in blue ) , but only at locations where little or no nonthermal x - ray emission is present .
the bright optical filament known as knot g ( at the far left edge of figure [ fig : tycho_balmer ] ) is the only location where balmer filaments and x - ray synchrotron emission appear coincident .
however , upon closer inspection the anticorrelation between the balmer line and synchrotron emission can be seen in knot g as well : the upper half of the filament , where balmer line emission is strongest , exhibits minimal synchrotron emission , while the opposite is true in the lower half of the filament .
the enhanced nonthermal emission inside of knot g may be due to the strong recent deceleration of knot g , where the snr is currently propagating into a strong density gradient at the outermost edge of an h  i cloud ( ghavamian et al .
the lack of optical / x - ray synchrotron correlation is especially striking given that the balmer - dominated filaments in tycho s snr have a high enough shock velocity ( @xmath101800 - 2100 kms@xmath26 ) to accelerate particles to tev energies .
the anticorrelation between the optical and nonthermal x - ray emission can be observed in other snrs as well , including sn 1006 , where recently x - ray proper motions have been measured along the entire rim by katsuda et al .
( 2012 ) . as with tycho s snr ,
the locations of the balmer - dominated filaments and the nonthermal x - ray filaments along the nw rim of sn 1006 are mutually exclusive . instead
, the balmer - dominated shocks are closely associated with thermal x - ray filaments having a proper motion consistent with a shock velocity of 3300@xmath65200@xmath65300 kms@xmath26 ( statistical and registrational uncertainties , respectively ) for a distance of 2.2 kpc .
this result is in excellent agreement with the shock velocity of 2890@xmath65100 kms@xmath26 determined from the broad h@xmath48 width and broad - to - narrow ratio of the nw filament by ghavamian et al .
such close agreement is a strong indication that little substantial energy has been lost from the thermal plasma to cosmic ray acceleration , similar to optical proper motion studies from rcw 86 ( helder et al .
, in preparation ) .
from the above discussion it appears that balmer - dominated snrs , while offering powerful diagnostics of @xmath16 and @xmath39 , are not useful for investigating equilibration the extreme cases of strongly cosmic - ray modified shocks .
in fact , the very condition allowing for the detection of the balmer line emission - presence of neutral gas ahead of the shock - is also responsible for limiting the fraction of shock energy lost to cosmic ray acceleration .
quantatitive evaluations of this effect by drury et al .
( 1996 ) and reville et  al .
( 2007 ) show that when the preshock gas is significantly neutral , alfvn waves driven by the cosmic rays ahead of the shock are dissipated by ion - neutral damping .
as long as the charge exchange frequency , @xmath81 ( @xmath82 ) is larger than the alfvn wave frequency , @xmath83 ( @xmath84 ) the ions and neutrals oscillate coherently , and ion - neutral damping is not important . however , when @xmath85 the neutrals are left behind by the ions in the alfvn wave motion , and during the incoherent oscillation between the two , charge exchange exerts a drag on the alfvn waves , damping them . the condition required for alfvn waves to not be strongly damped in the precursor
can be written out as @xmath86 where @xmath87 is the alfvn speed of the ions ahead of the shock and @xmath88 is the preshock neutral density .
given that cosmic rays resonantly scatter off alfvn waves having doppler shifted frequencies comparable to their gyrofrequency , and that the cosmic ray gyrofrequency is related to its energy via @xmath89 , the inequality above can be cast in terms of the energy , @xmath90 , below which a significant fraction of the cosmic ray flux out of the shock is reduced by ion - neutral damping : @xmath91 where we have set @xmath928.4@xmath9310@xmath94 @xmath95 s@xmath26 ( kulsrud & cesarsky 1971 ) , @xmath96 is the preshock magnetic field strength in units of 3 @xmath12 , @xmath15 and @xmath97 are the total preshock density and neutral fraction , and @xmath98 is the preshock temperature in units of 10@xmath99  k. for balmer - dominated snrs , where recent models have required moderate amplification of the preshock magnetic field ( @xmath1003 - 5 ; ghavamian et al .
2007 ) and where the preshock temperature may exceed 20,000  k ( raymond et al . 2011 ) , @xmath1014 tev for the typical case where @xmath1020.5 @xmath29 .
snrs exhibiting nonthermal x - ray emission are believed to contain cosmic rays with energies of tens of tev , so @xmath1034 tev is certainly high enough reduce the effectiveness of balmer - dominated shocks in producing nonthermal x - ray emission .
however , as noted earlier , a modest back pressure from cosmic rays is required to explain the width of the h@xmath48 narrow component line , as well as the low broad - to - narrow ratios seen in some snrs ( rakowski et al .
2009 ; raymond et al .
in fact , one model for electron heating in fast collisionless shocks requires at least some feedback from the cosmic rays in order to explain the moderate heating of electrons in snrs , as well as the inverse squared relationship between @xmath16 and @xmath39 ( ghavamian et al .
2007 , described in the next section ) . furthermore , as pointed out by drury et al .
( 1996 ) , the ion - neutral damping of alfvn waves in the precursor is unimportant for cosmic rays which have already exceeded @xmath90 . since the acceleration time for cosmic rays shortens considerably with shock speed ( @xmath104 ; malkov & drury 2001 ) , the fastest balmer - dominated shocks are more likely to have accelerated particles beyond @xmath90 and hence will begin to exhibit nonthermal x - ray emission and cosmic - ray modified shock structure .
a good example is the aforementioned snr 0509@xmath5667.5 , where the shock speeds exceed 5000 kms@xmath26 ( helder , kosenko & vink 2010 ) and nonthermal x - ray emission from cosmic ray accelerated electrons is detected from the forward shock .
the forward shocks in more evolved balmer - dominated snrs ( such as sn 1006 and tycho s snr ) will have swept up more mass and slowed down to speeds @xmath572000 kms@xmath26 , by which point @xmath105 will have lengthened and the shocks will be less cosmic - ray dominated .
given the lack of in situ measurements of the particle distributions in snrs , the electron heating mechanisms in these shocks have been studied primarily via numerical methods . on one hand ,
a number of studies have focused on electron heating in relativistic shocks , with the aim of modeling high energy emission from gamma - ray bursts ( e.g. , gedalin et al . 2008 ; sironi & spitkovsky 2011 ) .
these shocks are in a different area of parameter space than the snr shocks discussed here , and the physics governing the electron heating in relativistic shocks is substantially different . at the very high alfvnic mach numbers characteristic of gamma - ray bursts , the shock transition becomes very thin ( less than an electron gyroradius ) .
electrons in this case may be once again accelerated by the cross shock potential , similar to the very low mach number case . on the other hand ,
a number of other studies consider non - relativistic shocks relevant to snrs ( @xmath570.01c ) , where accelerated particles such as cosmic rays or solar energetic particles ( seps ) may play an important role in establishing their shock structure .
these studies have sought to identify plasma waves capable of boosting electrons to mildly relativistic energies ( e.g. , amano & hoshino 2010 ; riquelme & spitkovsky 2011 ) , with the objective of understanding how electrons are injected into the cosmic ray acceleration process .
this is a different ( though related ) question from what we consider here , namely how electrons are promptly heated to temperatures @xmath105@xmath9310@xmath70  k at the shock front ( ghavamian et al . 2007 ; rakowski et al . 2008 ) .
this limits our consideration of the work done so far to two broad scenarios of electron heating in fast , non - relativistic collisionless shocks .
one scenario is based on lower hybrid wave heating in the cosmic ray precursor ( laming 2000 ; ghavamian et al .
2007 ; rakowski et al . 2008 ) while the other is based on counterstreaming instabilities ahead of the shock ( e.g. , the buneman instability , cargill & papadapoulos ( 1988 ) , matsukiyo 2010 ; dieckmann et al .
we discuss these mechanisms below in turn .
the most significant result of the balmer - dominated shock studies , the inverse squared relation between @xmath16 and @xmath39 , places a useful constraint on the range of plausible equilibration mechanisms at the shock front .
the simplest way to obtain @xmath16@xmath106 is to set @xmath107 at shock speeds of 400 kms@xmath26 and higher , while allowing @xmath1 to rise according to the rankine - hugoniot jump conditions , @xmath108 .
the requirement that @xmath16=1 at @xmath39=400 kms@xmath26 gives @xmath1090.3 kev for @xmath39@xmath54400 kms@xmath26 , independent of shock velocity ( ghavamian et al .
although there may be marginal evidence of a departure from this relation at shock speeds exceeding 2000 kms@xmath26 ( van adelsberg et al . 2008 ) , a velocity - independent heating of electrons in snr shocks is an important clue to the nature of plasma heating processes in fast collisionless shocks .
it suggests that plasma processes ahead of the shock front are an important ( if not dominant ) source of electron heating in snrs ( ghavamian et al .
2007 ; rakowski et al . 2009 ) .    as mentioned earlier ,
strong interstellar shocks are expected to form a precursor where cosmic rays crossing upstream give rise to alfvn waves and turbulence ( blandford & eichler 1987 ; jones & ellison 1991 ) , compressing and pre - heating the gas before it enters the shock .
as long as the shock is strong ( @xmath1101/4 ) and cosmic ray pressure does not dominate the postshock pressure ( @xmath111b / b does not greatly exceed unity , with @xmath5720% of the postshock pressure provided by cosmic rays ) the thermal heating within the precursor does not depend strongly on shock velocity
. the limited range of narrow component h@xmath48 widths observed in balmer - dominated snrs over a wide range in shock speeds ( sollerman et al . 2003 ; raymond et al .
2011 ) is consistent with the relative insensitivity of the preshock heating to shock speed . since the widening of the h@xmath48 narrow component line
is now believed to arise in a precursor where the gas is heated by the damping of cosmic - ray driven waves ( wagner et al .
2008 ; raymond et al . 2011 ; morlino et al .
2012 ) , it stands to reason that perhaps the physical processes generating a constant electron heating with shock speed ( @xmath1090.3 kev ) also originate within the cosmic ray precursor .
the above argument was used by ghavamian et al .
( 2007 ) and rakowski et al .
( 2008 ) to advocate for a heating model where lower hybrid waves within the cosmic ray precursor preheat electrons to a constant temperature before they enter the shock front .
this model was based on the conception of mcclements et al .
( 1997 ) , who suggested that lower hybrid waves driven by the reflected population of nonthermal ions could generate lower hybrid waves ahead of the shock , pre - heating electrons and injecting them into the cosmic ray acceleration process .
the condition for generating such waves is that the shock be quasi - perpendicular , and that the reflected ions form a beam - like ( gyrotropic ) distibution .
a similar scenario was suggested by ghavamian et al .
( 2007 ) and rakowski et al .
( 2008 ) , but with one crucial difference : the reflected particles considered are ultra - relativistic cosmic rays rather than suprathermal ions .
the lower hybrid waves are electrostatic ion waves which propagate perpendicular to the magnetic field and whose frequency is the geometric mean of the electron and ion geofrequencies , @xmath112 . the group velocity of these waves is directed primarily along the magnetic field lines ( @xmath113 ; laming 2001 ) and the waves can simultaneously resonate with ions moving across the field lines and electrons moving along the field lines .
although the growth rate of lower hybrid waves is generally small ( rakowski et al . 2008 ) , their group velocity perpendicular to the magnetic field ( and hence the shock front ) can be on the order of the shock velocity ( @xmath114@xmath39 ) .
this allows the lower hybrid waves to remain in contact with the shock for long periods of time , attaining high intensities capable of effectively heating the electrons ( mcclements et al .
1997 ; ghavamian et al .
2007 ) .    in the case of cosmic rays ,
the time spent by the electrons in the precursor is @xmath115 .
the kinetic energy acquired by the electrons in the precursor is @xmath116 , where @xmath117 is the momentum diffusion coefficient of electrons ( ghavamian et al .
2007 ) . for lower hybrid wave turbulence , @xmath118 ( karney 1978 ; ghavamian et al .
2007 ; rakowski et al . 2008 ) , so that @xmath119 , as needed to account for the inverse squared relationship between equilibration and shock speed .
note that under the assumption that nonlinear amplification of the preshock magnetic field is not too strong ( @xmath120 ) , @xmath121 is that of bohm diffusion , which scales as 1/b , so that @xmath122 is also approximately independent of b.    during the past decade more refined models of cosmic ray acceleration have shown that a non - resonant mode of alfvn waves , having a higher growth rate than the previously considered resonant mode ( skilling 1975 ) , can be excited by cosmic rays in the precursor ( bell & lucek 2001 , bell 2004 , 2005 ) . unlike for the resonant case , the non - resonant amplification allows for @xmath123 , driving preshock magnetic fields to values as high as 1 mg ( vink & laming 2003 ; berezhko et al .
2003 ; bamba et al . 2005
; ballet 2006 ) .
such magnetic fields are hundreds of times stronger than the canonical preshock magnetic field of 3 @xmath12 and high enough to account for the observed narrowness of x - ray synchrotron - emitting rims in such snrs as sn 1006 ( assuming the narrowness is due to rapid cooling of high energy electrons behind the shock ; see ballet 2006 and morlino et al .
additional studies have suggested that non - resonant amplification may dominate early in the life of the snr , while resonant amplification may take over during the sedov - taylor stage of evolution ( amato & blasi 2009 ; schure et al . 2012 ) , though in either case , @xmath12410 is readily attained .
such a strong magnetic field effectively reduces the acceleration time for particles , and is very well suited for explaining how cosmic rays can reach the knee in the cosmic ray spectrum near 10@xmath125 ev ( bell & lucek 2001 ; eriksen et al .
2011 ) .
an important factor influencing the effectiveness of lower hybrid wave heating of electrons is the orientation of the preshock magnetic field relative to the shock front .
lower hybrid wave heating is only effective in perpendicular shocks ( vink & laming 2003 , ghavamian et al .
2007 ; rakowski et al .
given their spherical global geometry , snr blast waves generally propagate at a range of angles to the interstellar magnetic field .
x - ray observations and models of such snrs as sn 1006 ( orlando et al . 2007
; petruk et al .
2008 ) have indicated that perpendicular shocks are far more effective at accelerating cosmic rays than parallel shocks .
although the detailed implications of such differences have not yet been worked out for the lower hybrid wave heating model , rakowski et al .
( 2008 ) argue that even for quasi - parallel shock geometries the cosmic ray current driving the nonresonant alfvn waves will generate a significant perpendicular magnetic field ahead of the shock ( such a possibility has also been inferred from numerical simulations ; riquelme & spitkovsky 2011 ) .
this would allow lower hybrid wave growth to overtake modified alfvn wave growth for arbitrary orientations of the far upstream magnetic field , and allow for a more ubiquitous role for lower hybrid wave heating of electrons .
the amplification of the preshock magnetic field well beyond its far upstream value introduces an interesting possibility : effective lowering of the alfvnic ( and hence magnetosonic ) mach number of the shock . for the balmer - dominated shocks , where analysis of the optical spectra has shown that at best only a moderate fraction ( @xmath5720% ) of the shock energy
has likely been channeled into cosmic rays , the widening of the h@xmath48 narrow component has been interpreted as nonthermal broadening caused by the lowest frequency waves in the precursor ( ghavamian et al . 2007 ; though see raymond et al ( 2011 ) for a thermal intepretation ) . to explain the 30 - 50 kms@xmath26 widths of the h@xmath48 narrow component , the preshock magnetic field must be enhanced by a factor of a few . for the non - resonant alfvn waves in the bell ( 2004 ) mechanism
, the magnetic field energy density immediately behind the shock is given by ( schure et al .
2012 ) @xmath126 where @xmath127 is the fraction of the shock ram pressure channeled into cosmic rays . solving this expression for b
gives @xmath128 , where @xmath129 is the shock speed in units of 1000 kms@xmath26 . for @xmath1300.1 - 0.2 , @xmath451 @xmath29 , a postshock compression factor of 4 and balmer - dominated shock speeds @xmath102000 kms@xmath26
, this gives @xmath1314 - 10 ahead of the shock .
correspondingly , @xmath132 can increase by nearly the same factor , so that @xmath133 can be reduced by as much as an order of magnitude .
treumann & jaroschek ( 2008 ) describe the physical picture in this case as that of the shock having to prevent an increasing number of ions from crossing the shock jump by deflecting an increasing number of them at higher and higher mach numbers .
this deflection is necessary so that the ability of the shock to dissipate the inflowing energy is not overwhelmed . by deflecting these ions back upstream into a precursor ,
the net inflow of momentum and energy density into the shock is reduced , reducing the net difference in velocity between the inflowing and outflowing ions .
this effectively reduces the mach number in the frame of the upstream medium .
note that the ions in the precursor are only mildly compressed ( wagner et al .
2008 ; morlino et al .
2012 ) , which only weakly counteracts the rise in b. in addition , @xmath132 only scales as @xmath134 , but scales directly as b. the result is that _ given the compelling evidence for enhanced preshock magnetic fields in snrs shocks , the mach numbers of these shocks may be overestimated by as much as an order of magnitude . _ as we describe in section [ sec:5 ] , a unified description of solar wind and snr shocks , where the physics of electron - ion temperature equilibration occurs over a similar range in mach numbers and involves a similar range of physical processes , may be possible .
similar to the lower hybrid wave model , the buneman instability - driven wave model focuses on the region immediately ahead of a quasi - perpendicular shock . however , unlike the lower hybrid wave model , the buneman instability models consider the reflected nonthermal ion distribution , rather than ultrarelativistic cosmic rays . in the latter model , @xmath1020% of the ions
are reflected backstream against the incoming electron and ion plasma ( papadapoulos 1988 ; cargill & papadapoulos 1988 ) . here
the upstream plasma is not electrically neutral due to the positive charge of the reflected ion distribution . in such cases ,
a drift is induced between the electrons and ions .
the microinstabilities excited by this configuration depend upon the relative size of the electron thermal speed relative to the electron - ion drift velocity .
the buneman instability occurs when the drift velocity of the reflected ions relative to the upstream electrons exceeds the thermal speed of the upstream electrons ( @xmath135)(cargill & papadapoulos 1988 ) , a condition which occurs for very high mach number ( @xmath13650 ) shocks .
if the reflected ion current upstream is strong enough , the electron current generated to counteract it may produce a large enough drift between the preshock ions and electrons to cause a secondary buneman instability when the ion speed exceeds the electron thermal speed ( dieckmann et al .
the buneman instability generates electrostatic plasma waves which damp by rapidly heating the preshock electrons to @xmath137 kev , where @xmath138 is the shock speed in units of 1000 kms@xmath26 . until their thermal speed matches the electron - ion drift speed , at which point the instability saturates .
the rapid heating of the electrons perpendicular to the magnetic field results in @xmath1391 and makes it possible for an ion acoustic instability to occur between the preshock electrons and either the reflected ions or preshock ions ( cargill & papdapoulos 1988 ) . the waves generated by the ion acoustic instability can then transfer a substantial fraction of the shock energy ( tens of percent ) into electron thermal energy .
this process occurs over a length scale of @xmath140 ( as opposed to @xmath141 for the cosmic ray precursor ) , resulting in a @xmath16@xmath590.2 , independent of shock speed .
_ this is in strong disagreement with the equilibrations obtained for the balmer - dominated shocks . _    a number of other electron heating mechanisms , such as the modified two - stream instability and electron - cyclotron drift instability ,
have been proposed for collisionless shocks based on results from particle in cell ( pic ) simulations ( umeda et al .
2012a , 2012b ) .
a unified picture proposed by matsukiyo ( 2010 ) suggests that electrons can be strongly energized at low mach numbers ( @xmath14210 ) via a modified two - stream instability , where the velocities of the reflected / incoming ions and the electrons are lower than the thermal speed of the electrons and the electrons are able to damp out the buneman instability . in this case , obliquely propagating whistler mode waves are excited , having frequencies between the electron cyclotron frequency and the lower hybrid wave frequency .
when the electron - ion drift velocity and electron thermal speed become nearly equal , the electron - cyclotron drift instability becomes important ( umeda et al .
2012a ) , exciting waves with frequencies that are multiples of the electron cyclotron frequency . at higher mach numbers the electron thermal speed is lower than that of the ions , and the buneman instability / ion acoustic wave process descrbied earlier
is predicted to take over .    the amount of electron heating predicted by the buneman instability / ion acoustic wave model scales as @xmath143 ( cargill & papadapoulos 1988 ; matsukiyo 2010 ) , so that for shocks in the 2000 kms@xmath26 - 10,000 kms@xmath26 range , @xmath1442 - 50 kev .
this is clearly at odds with @xmath1450.3 kev observed between 400 kms@xmath26 and 2000 kms@xmath26 for balmer - dominated shocks .
one explanation for this discrepancy is that growth of the buneman - like and two - stream instabilities described above requires that the reflected ions form a distribution function with a positive gradient at some velocity ( laming 2000 ) .
this distribution forms when specularly reflected ions have a mostly monoenergetic , beamlike configuration . at the low mach numbers in the solar wind ( @xmath5710 ) , where the shock structure is laminar
, the reflected ions closely resemble a monoenergetic beam . however , at the higher mach numbers , where the shock front is more turbulent and disordered , the reflected ions are likely to have a greater spread in energy and are probably less beamlike ( laming 2000 ) .
this would lead to suppression of buneman - like instabilities .
however , this line of reasoning is still speculative , and the real explanation for the lack of agreement between the observed @xmath16 and those predicted by models in this section remains to be explored .
cosmic - ray driven processes may ultimately provide a better explanation for electron heating at snr shocks than those involving reflected suprathermal ions .
the cross shock potential arises from the charge separation produced by the different gyroradii for ions and electrons as they cross the shock transition .
it accelerates electrons into the post shock layer , and can be a means of electron heating at subcritical shocks with an approximately laminar structure . at supercritical shocks ,
the time dependence and non - locality introduced reduces the degree of electron heating .
however at sufficiently alfvn mach number ( where the shock transition becomes thin ( on the order of the electron convective gyroradius or inertial length ) , significant electron heating may again occur . in the absence of magnetic field amplification by cosmic rays , this might be expected to happen at snr shocks .
however it is much more likely in environments where the plasma beta is low , such as galaxy clusters .
it may also occur in cases where a significant population of cosmic rays is unlikely due to the low age of the shock , as in gamma - ray bursts .
from the beginning , a detailed study of the fastest collisionless shocks has been hampered by one inherent limitation : they occur in objects which are too remote for in situ study .
although some collisionless shocks in our solar system reach alfvnic mach numbers as high as 30 ( such as those around saturn ; achilleos et al .
2006 ; masters et al .
2011 ) , there are no physical phenomena in our solar system energetic enough to produce the type of shocks seen in snrs ( mach numbers @xmath10100 or more if no preshock magnetic field enhancement in the cr is assumed ) .
in addition , the range of plasma betas attainable in the solar system is larger than the range of plasma betas attained in the interstellar medium .
another fundamental difference between solar wind and snr shocks is the fact that the former are short - lived phenomena confined to small spatial scales ( millions of km ) in a curved ( bow shock ) geometry , whereas the latter are sustained for thousands of years , on spatial scales of parsecs , often well - described by a planar geometry .
this results in an irreducible difference between the two types of shocks : particles crossing back and forth between upstream and downstream can remain in contact with snr shocks for long periods of time , allowing accelerated crs to acquire much more energy in snr shocks than in solar wind shocks .
this potentially allows the crs to create shock precursors with properties needed to heat electrons and influence @xmath16 .
the heating of electrons at the earth s and other planetary bow shocks has been the subject of much theoretical and observational work .
typical features of the electron temperature change , @xmath146 , observed at the bow shock as observed by @xcite include ( a ) an approximate relationship between heating and the incident solar wind energy @xmath147 , where @xmath148 is the component of the solar wind s velocity incident upon the shock , and ( b ) a relationship between the change in temperature normalized by the incident energy and the fast magnetosonic mach number , @xmath149 .
a similar approximate relationship holds between the normalized electron temperature change and the alfvn mach number @xmath133 , especially for shocks with a low plasma @xmath52 , which is the ratio between thermal and magnetic pressures .
recent work by @xcite shows that this relationship with @xmath133 holds well at saturn s bow shock .
this is particularly interesting as @xmath133 at saturn is often much larger than at earth .    in figure
[ fig : combined ] we plot the ratio of electron and ion temperatures downstream of the shock against the magnetosonic , and alfvn mach numbers , as well as the upstream flow velocity relative to the shock .
the data in these figures are taken from bow shock crossings of the isee spacecraft , and is a subset of those listed in @xcite , consisting of 61 crossings for which all the necessary data is available .
it is well known that quantities other than those displayed here may be more appropriate , specifically the change in electron temperature over the change in total temperature @xmath150 or even @xmath151 are better correlated with inverse mach numbers than @xmath152 @xcite .
nevertheless , we use the latter quantity here as it enables a comparison with data from extra - solar system and outer planetary shocks at where less data are available . furthermore , the approximate inverse dependence upon @xmath8 , @xmath133 , and @xmath39 is still quite apparent in this data .
it is interesting to note that the relationship with the mach numbers is much more favourable than the dependence upon @xmath39 , indicating that the mach number is the more relevant quantity for organizing the relationship between @xmath16 and shock strength .
is plotted versus shock speed ( left ) , alfvnic mach number ( center ) and magnetosonic mach number ( right ) .
green symbols show data from crossings of earth s bow shock ( @xcite ) , while the black symbols show data from crossings of saturn s bow shock ( @xcite ) .
shock speeds for the saturnian bow shock are based on a solar wind model and an assumed shock speed with respect to the spacecraft of 100 kms@xmath26 , and ion temperatures are based on electron distribution measurements and the application of the rankine - hugoniot conditions ( see masters et al .
( 2011 ) for a full discussion of shock parameter derivations at saturn s bow shock ) .
red symbols show data acquired from balmer - dominated snr shocks ( van adelsberg et al .
2008 ) , and assume @xmath132=9 kms@xmath26 , @xmath153=11 kms@xmath26 .
, title="fig:",width=480 ] -.5 in    many mechanisms may be involved in the heating of electrons in solar system shocks .
proposed mechanisms involve acceleration of electrons by a cross - shock potential @xcite , wave turbulence @xcite , microinstabilities @xcite , and electron trajectory scattering @xcite .
the existence of a cross - shock potential may be deduced from the generalized ohm s law in which a gradient in electron thermal pressure gives rise to an electric field . examining the energetics of electrons crossing the shock may simplified by working in the de hoffmann - teller frame of reference @xcite , defined as the frame in which the shock is at rest , and in which the magnetic field and plasma flow velocity are ( anti-)parallel . in this case
the the electric field is dominated by that generated by the electron pressure gradient , and the work done on electrons crossing the shock is determined by the cross - shock potential .
additional mechanisms are required to scatter electrons to pitch angles that are more perpendicular , and to flatten the distribution so that empty regions of phase space are filled .
this results in a distribution whose temperature is controlled to a large extent by the de hoffmann - teller frame cross - shock potential and downstream density .
in addition to direct measurement of the electric fields within the shock , @xcite , comparison of upstream and downstream electron phase space distributions have shown that these are consistent with electron acceleration by a cross - shock potential in the de hoffmann - teller frame @xcite .
figure [ fig : combined ] may indicate that similar mechanism(s ) heat the electrons in solar wind and in snr shocks .
_ this is especially appealing when we remember that @xmath133 in snrs may be overestimated due to preshock amplification of magnetic field by cosmic rays . _ in their study of the terrestrial bow shock and interplanetary shocks , schwartz et al ( 1988 ) found that the electron - ion temperature equilibration organizes best by @xmath154 . given the difficulty in determining the mach numbers of snr shocks , the equilibration dependence on shock strength has been characterized via the shock speed instead , and found to obey @xmath155 1/@xmath39@xmath156 ( @xmath39 is much more accurately known than the mach numbers ) .
subsection 4.1 outlined a model for snr electron heating , where the cosmic ray diffusion coefficient @xmath157 is assumed independent of @xmath39 . from quasi - linear theory ( blandford and eichler 1987 ) , @xmath158 where @xmath148 is the energy density of turbulence ( @xmath159 ) and @xmath160 is the cosmic ray velocity . for resonant amplification ,
we evaluate @xmath148 at @xmath161 , where @xmath162 is the gyroradius and @xmath163 the parallel component of the cosmic ray velocity . for relativistic cosmic rays , where @xmath164 , this results in @xmath165 , assumed constant with @xmath39
. however , for nonrelativistic suprathermal particles , @xmath160 will most likely be proportional to the shock velocity @xmath39 , which with the same assumptions leads to @xmath1661/@xmath39 for solar wind shocks ( as opposed to @xmath167 in snrs ) .
this argument is admittedly loose , and should not be be viewed as much more than a hypothesis to motivate further work .    in our arguments above we have made considerable assumptions about @xmath121 .
the most obvious one is that @xmath157 as written above applies to parallel shocks , whereas we are most likely dealing with quasi - perpendicular cases .
this may reduce the difference anticipated between solar wind and smr shocks , depending on the turbulence spectrum ( summarized in appendix a of rakowski et al .
2008 ) .
the degree of cosmic ray magnetic field amplification at snr shocks required to bring snr data points in figure  [ fig : combined ] into alignment with solar wind data points is approximately an order of magnitude or less at the highest velocities considered .
given the degree of magnetic field amplification expected from cosmic ray acceleration , this appears to be highly plausible . in the case of saturated nonresonant instability ( bell 2004 , 2005 ) by resonant scattering ( luo & melrose 2009 )
, @xmath168 is expected . in the case of nonresonant saturation ,
higher values , but strongly dependent on @xmath169 , are predicted .
saturation by electron heating ( i.e. the @xmath133 where growth of lower hybrid waves becomes greater than the growth rate of magnetic field ) leads to similar magnetic field enhancement , with @xmath170 ( rakowski et al .
such magnetic field amplification is less likely at solar wind shocks .
the suprathermal particle densities are lower in solar wind shocks , and the ambient magnetic fields are higher , much closer to where the nonresonant instability would saturate ( if not already beyond it ) .
collisionless shocks occur over a vast range of length scales , with those in galaxy clusters being among the largest . while the shock speeds in the galaxy cluster shocks are similar to those in supernova remnants ( up to 4000 kms@xmath26 ; markevitch 2005 ; markevitch & vikhlinin 2007 ; russell et al .
2012 ) , they occur in environments that are substantially different from both the solar wind and the ism .
these differences can be encapsulated via the plasma beta , defined as the ratio of the thermal pressure to the magnetic pressure of a plasma ( @xmath171 ) . utilizing solar wind parameters listed by bruno & carbone ( 2005 )
, this ratio ranges from around unity at 1 au under fast solar wind conditions ( wind velocity @xmath10900 kms@xmath26 ) to around 20 for the quiescent wind ( wind velocity @xmath10300 kms@xmath26 ) .
the plasma @xmath52 of the ism is close to that of the fast solar wind , @xmath1721 - 4 ( assuming a b=3@xmath36 g , n= 1 @xmath29 and t=10@xmath99  k ) . on the other hand
, the electron temperature of the intracluster medium ( icm ) ranges from  10@xmath173k - 10@xmath71k ( 1 - 10 kev ) , with number densities steeply declining from @xmath10 10@xmath174 @xmath175 near the cluster centers to @xmath1010@xmath176 @xmath175 at the outer edges .
the corresponding sound speed is close to 1000 kms@xmath26 , nearly two orders of magnitude higher than in the general ism .
the magnetic fields measured in galaxy clusters are actually close that of the ism , typically on the order of a microgauss ( carilli & taylor 2002 ) .
therefore , @xmath17750 kms@xmath26 in galaxy clusters , so that @xmath1781 and the magnetic field pressure has negligible contribution to the dynamics of shocks in galaxy clusters .
this puts galaxy cluster shocks in a different region of parameter space than solar wind and snr shocks .
clusters such as 1e0657@xmath5656 ( the bullet cluster ) and a520 show strongly enhanced x - ray emission from collisionless shocks , formed during major mergers when gas from one cluster plunges through gas from the other ( markevitch et al .
2005 ; markevitch & vikhlinin 2007 ; russell et al .
shocks moving mostly along the plane of the sky have a favorable viewing geometry and appear as giant curved structures hundreds of kiloparsecs in length .
the large clumps of infalling gas drive bow shocks into the cluster gas , which has already been heated to at least 1 kev , and produces thermal bremsstrahlung emission peaking close to that energy .
this is another major difference between collisionless shocks in the ism and those in the icm .
while the alfvnic and magnetosonic mach numbers of snr shocks are very difficult to determine due to the lack of available observational constraints on magnetohydrodynamic quantities upstream ( such as @xmath14 , @xmath13 and @xmath15 ) , those in galaxy clusters can readily be measured by spectral analysis of x - ray emission upstream .
comparison of this emission to that of the enhanced postshock region gives the density contrast between the downstream and upstream gas ( i.e. , @xmath179 ) . this density contrast yields the sonic mach number , @xmath180 , via solving the rankine - hugoniot jump conditions ( russell et al . 2012 ) : @xmath181 where @xmath80 is the ratio of specific heats of the cluster gas .
measurement of these jumps from x - ray observations have yielded @xmath1821.5 - 3 for the bullet cluster ( markevitch et al . 2005 ) , abell 520 and
abell 2146 ( russell et al .
2012 ) . using these estimated mach numbers ,
the shock velocity itself can be calculated via @xmath39=@xmath183 , where @xmath153 is the upstream sound speed as inferred from the x - ray spectra .
this yields shock speeds ranging between 2500 kms@xmath26 and 4000 kms@xmath26 , similar to the fastest known balmer - dominated shocks .
however , fits to the x - ray spectra behind these shocks show prompt electron - ion equilibration at the shock front , consistent with @xmath16=1 ( markevitch et al .
2005 ; markevitch & vikhlinin 2007 ) , despite the extremely high shock speeds involved .
this result can only be reconciled with equilibration measurements from the solar wind and snrs if the equilibration depends on the mach number , rather than @xmath39 .
furthermore , given the low mach numbers found in the galaxy clusters , it is plausible that the shock transitions in these cases are laminar ( as opposed to turbulent like the snr and fastest solar wind shocks ) , with electron heating occurring efficiently at the shock front via the same type of cross - shock potential as seen in the slowest solar wind and slowest snr ( @xmath39@xmath53 400 kms@xmath26 ) shocks .
this is of course a speculation ; further insight into collisionless cluster shocks may be obtained via numerical simulations for the appropriate conditions .
there have been exciting advances in the study of electron - ion temperature equilibration in collisionless shocks during the past few years .
perhaps most notable has been the growing realization that temperature equilibration and cosmic ray acceleration may be intertwined processes .
optical studies of collisionless shocks in partially neutral gas ( balmer - dominated shocks ) have shown that the electron - ion temperature equilibration is a declining function of shock speed , well characterized as @xmath16@xmath184 .
this relationship most likely arises due to electron heating ahead of the shock that is nearly independent of shock speed above 400 kms@xmath26 .
cosmic ray precursors , with moderately amplified preshock magnetic field and density , are the most logical sites for electron heating in snr collisionless shocks .
the transition to fully equilibrated snr shocks at speeds below 400 kms@xmath26 may be due to a less turbulent and more laminar structure at low shock speeds and mach numbers .
this allows the electrons to experience a more uniform cross - shock potential , and hence higher energization compared to higher shock speeds and mach numbers , where the shock jump is more turbulent and disordered .
the magnetosonic mach numbers of snr shocks may match those in solar wind shocks if there is a approximately an order of magnitude increase in the alfvn speed of the preshock gas in snrs compared to the average ism value .
this is possible if there is a moderately amplified preshock magnetic field ( @xmath1010 ) , and is readily provided by compression and heating in a cosmic ray precursor . in solar wind shocks , the precursor is due to suprathermal , non - relativistic ions , resulting in a @xmath16@xmath185 seen in the solar wind .
while balmer - dominated shocks have allowed us to elucidate some of the physics of electron - ion temperature equilibration , ion - neutral damping limits most of those observed to cases where the shock structure has not become strongly modified by cosmic ray acceleration .
given this limitation , _ electron - ion equilibration studies of fast , collisionless shocks in fully pre - ionized gas would be highly desirable_. such a sample would allow the equilibration to be studied over a range of speeds where shocks are increasingly affected by feedback from the accelerated cosmic rays .
in such circumstances it is unclear what would happen to the @xmath16 versus @xmath39  relation .
if , as predicted by amato & blasi ( 2009 ) , bell s non - resonant cosmic ray instability takes over from the resonant instability at the highest shock speeds , then additional electron heating may occur in the fastest shocks ( @xmath39@xmath585000 kms@xmath26 ) , resulting in substantial deviation from the @xmath16@xmath184 relation .
such deviations may already have been seen in the fastest ( @xmath39@xmath582000 kms@xmath26 ) shocks , where there is evidence that @xmath16 does not settle down to @xmath17 , but rather @xmath100.03 .
other deviations may have been detected in snr 0509@xmath5667.5 , where @xmath16 for a 5000 kms@xmath26 shock has been estimated to be @xmath100.2 , substantially higher than predicted by the inverse squared relation .
however , the study of such shocks would be challenging . without an h@xmath48 broad component to constrain the range of @xmath39
, shock speeds would have to be determined by other means , such as proper motion studies . that would require x - ray and or uv imagery of snrs with well - constrained distances ( such as those in the lmc or smc ) , over multiple epochs .
it would also require spectroscopy of the forward shocks in these snrs , in order to constrain both the electron temperature ( via x - ray continuum fits and and uv emission line ratios ) and the ion temperature ( e.g. , via he  ii , c  iv , n  v and o  vi resonance lines ) .    an important test of our ideas concerning electron heating by cosmic ray generated waves in a shock precursor would be to measure electron temperatures at snr shocks exhibiting strong cosmic ray modification and substantial magnetic field amplification ( @xmath186100 ) . several snrs show x - ray filaments produced by synchrotron radiation from cosmic ray electrons , and are generally distinct from those shocks with strong balmer emission .
the absence of neutral material ahead of these shock means that optical and uv emission is weak , and electron temperatures will have to be measured from x - ray spectra .
difficulties arise in distinguishing thermal electron bremsstrahlung from cosmic ray electron synchrotron emission , requiring data of high signal to noise .
further complications in some snrs ( e.g. cas a ) stem from scattering of bright x - ray emission from the ejecta such that it coincides spatially with emission from the forward shock .
such scattering may either be local , due to snr dust , or instrumental , due to telescope imperfections .
this leaves sn 1006 as the best likely target for such an observation , since due to its evolutionary state , only the outer layers of ejecta have encountered the reverse shock .
the development of missions like solar orbiter and solar probe plus will allow _ in situ _ measurements of shocks in the solar wind , most likely associated with coronal mass ejections , much closer to the sun .
these will probe a different parameter regime , where the magnetic field pressure dominates over the gas pressure ( low @xmath52 , similar to ism shocks ) .
as such , measurements here might yield insights into the properties of similar plasma in the precursors of snr shocks where the magnetic field has been amplified by cosmic rays .
g. acknowledges support by hst grant hst - go-11184.07-a to towson university .
jml acknowledges support by grant nnh10a009i from the nasa astrophysics data analysis program , and by basic research funds of the office of naval research .
k.  a.  eriksen , j.  p.  hughes , c.  badenes , r.  fesen , p.  ghavamian , d.  moffett , p.  p.  plucinsky , c.  e.  rakowski , e.  m.  reynoso , and p.  slane . .
_ astrophysical journal _ , 728:0 l28 , february 2011 .
doi : 10.1088/2041 - 8205/728/2/l28 .
k. france , r. mccray , s.  v. penton , r.  p. kirshner , p. challis , j.  m. laming , p. bouchet , r. chevalier , p.  m. garnavich , c. fransson , k. heng , j. larsson , s. lawrence , p. lundqvist , n. panagia , c.  s.  j. pun , n. smith , j. sollerman , g. sonneborn , b. sugerman , and j.  c. wheeler .
_ astrophysical journal _
, 743:0 186 , december 2011 .
doi : 10.1088/0004 - 637x/743/2/186 .
e.  a.  helder , j.  vink , c.  gh .  bassa , a.  bamba , j.  a.  m.  bleeker , s.  funk , p.  ghavamian , k.  j.  van der heyden , f.  verbunt , and r.  yamazaki . .
_ science _
, 325:0 719 , august 2009 .
doi : 10.1126/science.1173383 .
a.  masters , s.  j. schwartz , e.  m. henley , m.  f. thomsen , b.  zieger , a.  j. coates , n.  achilleos , j.  mitchell , k.  c. hansen , and m.  k. dougherty . .
_ journal of geophysical research ( space physics ) _ , 116:0 a10107 , october 2011 .
doi : 10.1029/2011ja016941 .
h.  r. russell , b.  r. mcnamara , j.  s. sanders , a.  c. fabian , p.  e.  j. nulsen , r.  e.  a. canning , s.  a. baum , m. donahue , a. edge , l.  j. king and c.  p. odea . .
_ monthly notices of the royal astronomical society _ , 423:0 236255 , june 2012 .
doi : 10.1111/j.1365 - 2966.2012.20808.x .
t.  tanaka , y.  uchiyama , f.  a.  aharonian , t.  takahashi , a.  bamba , j.  s.  hiraka , j.  kataoka , t.  kishishita , m.  kokubun , k.  mori , k.  nakazawa , r.  petre , h.  tajima , and s.  watanabe . .
_ astrophysical journal _ , 685:0 9881004 , october 2008 .
doi : 10.1086/591020 .
j.  s.  warren , j.  p.  hughes , c.  badenes , p.  ghavamian , c.  f.  mckee , d.  moffett , p.  plucinsky , c.  e.  rakowski , e.  reynoso , and p.  slane . .
_ astrophysical journal _ , 634:0 376389 , november 2005 .
doi : 10.1086/496941 .
b.  j.  williams , w.  p.  blair , j.  m.  blondin , k.  j.  borkowski , p.  ghavamian , k.  s.  long , j.  c.  raymond , s.  p.  reynolds , j.  rho , and p.  f.  winkler . .
_ astrophysical journal _
, 741:0 96 , november 2011 .
doi : 10.1088/0004 - 637x/741/2/96 .
c.  s. wu , d.  winske , m.  tanaka , k.  papadopoulos , k.  akimoto , c.  c. goodrich , y.  m. zhou , s.  t. tsai , p.  rodriguez , and c.  s. lin . .
_ space science reviews _
, 37:0 63109 , january 1984 .
doi : 10.1007/bf00213958 .