ion pressure anisotropies are ubiquitous in heliospheric and astrophysical plasmas . in the absence of coulomb collisions , the magnetic moment of ions , @xmath4 ( @xmath5 , where @xmath6 is the ion velocity perpendicular to the local magnetic field @xmath7 and @xmath8 ) is an adiabatic invariant .
thus , if ion collisions are infrequent , the ion velocity distributions parallel and perpendicular to the magnetic field decouple , making @xmath9 ( where @xmath10 and @xmath11 are the pressure components of ions perpendicular and parallel to @xmath7 ) .
examples of systems where ion pressure anisotropies are important are low - luminosity accretion disks around compact objects @xcite , the intracluster medium ( icm ) @xcite , and the heliosphere ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) . in these systems
the ion pressure anisotropy is expected to play a key role in the large scale dynamics of the plasma .
for instance , @xcite pointed out that pressure anisotropies can give rise to an anisotropic stress that can contribute to angular momentum transport and plasma heating in low - luminosity accretion disks .
this physics is not included in standard mhd models of accretion disks .
velocity space instabilities limit the growth of the pressure anisotropy and give rise to small scale fluctuations in the magnetic field .
these fluctuations can , in turn , affect the large scale transport properties of the plasma by modifying the mean free path of particles .
magnetic field amplification in a collisionless system generically drives @xmath12 .
there are two ion - scale instabilities that can be excited in this regime : the mirror and the ion - cyclotron ( ic ) instabilities @xcite .
the mirror instability consists of non - propagating , strongly compressional modes .
their fastest growing wave vectors @xmath13 are oblique to @xmath7 , with the magnetic variations parallel to @xmath7 being much larger than the perpendicular fluctuations , @xmath14 . for a bi - maxwellian distribution of ions , and
assuming cold electrons , the threshold condition for mirror instability growth is given by @xmath15 , where @xmath16 ( @xmath17 ) is the ion temperature perpendicular ( parallel ) to @xmath7 , and @xmath18 @xcite .
the ic instability , on the other hand , consists of transverse electromagnetic waves , with the fastest growing @xmath13 being preferentially parallel to @xmath7 ( e.g. , * ? ? ?
whether the ic or mirror instability sets in first depends on how fast these instabilities grow for a given set of plasma conditions .
@xcite showed that for @xmath19 , the growth rate of the ic instability is @xmath20 , where @xmath21 is the cyclotron frequency of the ions .
in addition , for @xmath22 , the threshold anisotropy depends very weakly on @xmath23 .    in the linear regime
, the dominant instability will be determined by which threshold is reached first as @xmath7 is amplified .
these estimates for the threshold conditions imply that , in the linear regime , the ic instability should dominate for @xmath24 , while the mirror instability should dominate for @xmath25 .
although both the mirror and ic instabilities have been extensively studied in the linear regime , a complete theory of their nonlinear evolution and saturation is still lacking . in this paper , we are particularly interested in the question of how the mirror and ic modes behave after the initial phase of exponential growth .
indeed , in most astrophysical scenarios where the magnetic field is amplified , the growth of the field occurs on time scales much longer than the growth time of the relevant kinetic instabilities .
therefore , most of the evolution of the velocity space instabilities happens in a nonlinear regime , where the conditions typically assumed in linear studies , such as a homogeneous plasma or a bi - maxwellian distribution of particle velocities , are not necessarily satisfied .
moreover , if the field is amplified by order unity or more , a quasi - linear analysis is not applicable .
one of the key questions we are interested in addressing in this case is whether the mirror instability saturates via pitch - angle scattering that violates magnetic moment conservation @xcite or via a nearly @xmath4-conserving rearrangement of the magnetic field , with fluctuations @xmath26 @xcite .
we shall see that both of these saturation processes can in fact be important .
the nonlinear regime of the ic and mirror instabilities has been studied both theoretically ( e.g. , * ? ? ?
* ; * ? ? ?
* ) and numerically @xcite . in most numerical studies , however , the pressure anisotropy is treated as an initial condition , which decreases as the instabilities grow and saturate .
thus , in these approaches the nonlinear saturation can not be followed for a time much longer than that of the initial exponential growth .
moreover , many of the previous numerical studies have focused on one - dimensional simulations .    in this paper
we study the long term , nonlinear evolution of the mirror and ic instabilities using two - dimensional particle - in - cell ( pic ) simulations in driven systems . in order to self
- consistently explore a time much longer than the initial exponential growth , we continuously induce the growth of @xmath27 by amplifying the mean field @xmath28 during the simulation .
we concentrate on amplification by incompressible plasma motions .
thus , our simulations impose a shear velocity on the plasma , which sustains the growth of the magnetic field and maintains an overall positive pressure anisotropy .
this growth is intended to mimic magnetic field fluctuations in a turbulent plasma and mhd - scale instabilities that amplify the magnetic field , like the magnetorotational instability ( mri ; * ? ? ? * ) or convective instabilities driven by anisotropic thermal conduction in dilute plasmas @xcite .
we suspect that our results are also relevant to cases in which perpendicular heating or parallel cooling drives a plasma mirror and/or ic unstable , as in , e.g. , the solar wind .
@xcite conducted an analogous study to ours using hybrid simulations . in their study , expansion of the simulation box decreased the mean magnetic field , so that on average @xmath29 .
this in turn led to the growth and nonlinear saturation of the firehose instability .
also , @xcite conducted a similar hybrid simulation study where the simulation box is expanded along the directions parallel and perpendicular to @xmath7 . with this setup , @xmath10 becomes larger than @xmath11 on average , and both the mirror and ic instabilities can grow and reach saturation .
our work is complementary to these hybrid simulation studies and focuses on the process in which the condition @xmath12 is achieved by field amplification due to shearing plasmas .
we also consider higher values of @xmath30 @xmath31 , relevant to accretion disks , the icm , and the near - earth solar wind .
our calculations are less directly applicable to heliospheric measurements of mirror modes driven unstable by rapid heating through the earth
s bow shock ( e.g. , * ? ? ?
this is because we assume that the pressure anisotropy is generated on a timescale long compared to the ion cyclotron timescale , which is probably not true at collisionless shocks .
our paper is organized as follows . in  [ sec : numsetup ] we describe the simulation set up , emphasizing the key physical and numerical parameters .
[ sec : results ] shows our results and  [ sec : conclu ] presents our conclusions . throughout the paper
we compare some of our simulation results with linear theory predictions appropriate to our pic simulations .
these have artificially low ion to electron mass ratios ( @xmath32 . the linear theory results are based on the linear vlasov solver developed in @xcite .    during the completion of this work ,
@xcite presented calculations of firehose and mirror saturation in shearing plasmas with very similar results to those that we present here .
we use the electromagnetic , relativistic pic code tristan - mp @xcite in two dimensions .
the simulation box consists of a square box in the @xmath33 plane , containing plasma with a homogeneous magnetic field @xmath34 initially pointing along the @xmath35 axis .
since we want to simulate a magnetic field that is being amplified in an incompressible way , we impose a shear motion of the plasma so that the mean particle velocity is @xmath36 , where @xmath37 is a shear parameter with units of frequency and @xmath38 is the distance along @xmath35 . from flux conservation ,
the @xmath39-component of the mean field evolves as @xmath40 , implying a net growth of @xmath41 .
although we present simulations resolving the @xmath33 plane only , we also tried runs where the @xmath42 plane was resolved . in those cases the isotropization efficiency was substantially lower
this is because , if the @xmath42 plane is resolved , the growing component of @xmath7 is perpendicular to the plane of the simulation . this way the angle between the relevant wave vectors , @xmath13 , and @xmath7 can not be 0 .
this over - constrains the wave vectors that are allowed to exist , artificially reducing the isotropization efficiency .
an important physical parameter in our simulations will be the ratio of the initial ion cyclotron frequency to the shear frequency , @xmath43 .
we refer to this as the magnetization . in typical astrophysical environments @xmath44 . due to computational constraints
, however , we will use values of @xmath43 much smaller than expected in reality , although still satisfying @xmath45 .
the dependence of our results on the ratio @xmath43 will be carefully assessed .
ccccccccc beta6mag93 & 1 & 6 & - & 93 & 0.05 & 14 & 30 & 12 + beta6mag670a & 1 & 6 & - & 670 & 0.15 & 14 & 10 & 18 + beta6mag670b & 1 & 6 & - & 670 & 0.15 & 14 & 10 & 24 + beta20mag93a & 1 & 20 & - & 93 & 0.05 & 10 & 60 & 35 + beta20mag93b & 1 & 20 & - & 93 & 0.05 & 10 & 30 & 35 + beta20mag670a & 1 & 20 & - & 670 & 0.05 & 7 & 30 & 36 + beta20mag670b & 1 & 20 & - & 670 & 0.05 & 7 & 30 & 48 + beta20mag2000 & 1 & 20 & - & 2000 & 0.05 & 7 & 50 & 55 + beta80mag670 & 1 & 80 & - & 670 & 0.05 & 5 & 30 & 24 + beta80mag1340 & 1 & 80 & - & 1340 & 0.05 & 5 & 30 & 24 + betai20magi240mass10 & 10 & - & 20 & 240 & 0.05 & 5 & 20 & 25 [ table:1d ]    in standard mhd simulations where shear plasma motions are imposed , shearing periodic boundary conditions would be used along the @xmath38 direction ( see , e.g. , * ? ? ?
in that case , the flow velocities at the @xmath38-boundaries of the box are matched using galilean transformations of the fluid velocity . in the case of relativistic pic simulations
this can not be done in a self - consistent way .
the reason is that , under a relativistic change of reference frame , the current @xmath46 transforms in a way that is inconsistent with simply transforming the particle velocities .
we instead implement _ shearing coordinates _ , in which the grid moves with the shearing velocity @xmath47 . in this new frame
the average plasma velocity vanishes everywhere in the box and simple periodic boundary conditions are allowed both in the @xmath38 and @xmath39 axes . in these shearing coordinates ,
maxwell s equations gain additional terms ( see appendix a of * ? ? ?
* ) , becoming @xmath48 @xmath49 where @xmath50 is the speed of light , and @xmath51 is the electric field .
in addition to the modifications to maxwell s equations , the forces on the particles also acquire an extra term : @xmath52 where @xmath53 , @xmath54 , and @xmath55 are the particle s momentum , velocity , and charge , respectively .
the third term on the right hand side of equation [ eq : ind ] and the fourth term on the right hand side of equation [ eq : amp ] are proportional to time and arise from the motion of the shearing coordinate grid points with respect to those of the non - shearing coordinates .
indeed , as time goes on , the @xmath38 axis of the shearing coordinates is gradually tilted ( as seen from the non - shearing frame ) ; therefore the @xmath56derivatives of the fields ( @xmath57 ) in the shearing coordinates must include a time dependent term that accounts for the evolution of the @xmath38 axis .
the last terms in equations [ eq : ind ] and [ eq : amp ] , which are proportional to the coordinate @xmath39 , will be neglected for the following reasons . in our simulations
the box size will typically be a few times the larmor radius of the ions , @xmath58 .
thus , @xmath59 ( where @xmath60 is the thermal velocity of the ions ) , which is much smaller than unity since we are interested in the regime @xmath61 and @xmath62 . as a result ,
the last term in equation [ eq : ind ] will be much smaller than the term on the left hand side , especially since @xmath63 .
the last term of equation [ eq : amp ] , on the other hand , is not necessarily much smaller than the displacement current ( left hand side term ) , mainly because we expect @xmath64 .
however , if the characteristic timescale for the mirror and ic modes is close to @xmath65 ( below we will check that this is indeed the case ) , the last term in equation [ eq : amp ] will be much smaller than the first term on the right hand side ( the ratio between these terms scales as @xmath66 ) , so we choose to neglect it in the limit @xmath67 .
the fact that the neglected term can still be comparable to the displacement current implies that we may be excluding a relevant contribution to the charge density in gauss law .
however , the physics of interest in this paper is non - relativistic and charge separation does not play a role . as a result , @xmath68 ( equivalent to @xmath69 in the _ non - shearing _ coordinate system ) . neglecting the terms proportional to @xmath39
is also required by the periodic boundary conditions in the y - direction .
the existence of these terms at all is a consequence of the incompatibility of the galilean invariance of the shearing box and the lorentz invariance of the pic calculations . for consistency we will also neglect the term @xmath70 in equation [ eq : amp ] , which is also comparable to the displacement current ( in the case where @xmath51 evolves on a timescale close to @xmath65 )
can be integrated using simple time and space interpolations of @xmath71 .
this way , after this modification is implemented , the numerical algorithm used by tristan - mp continues to be second order accurate in time and space .
] our set of simulations are summarized in table 1 .
the physical setup is defined by the magnetization parameter , @xmath43 ( the ratio of the initial cyclotron frequency to the shear rate ) , the ratio between the initial particle and magnetic pressures , @xmath72 , the mass ratio between ions and electrons , @xmath73 , and the initial alfvn velocity of the plasma @xmath74 , relative to the speed of light @xmath50 , where @xmath75 and @xmath76 is the mass density .
the numerical parameters of our runs are defined by the spatial resolution ( @xmath77 ) , box size relative to the ion larmor radius ( @xmath78 ) , and number of particles per cell ( @xmath79 ) , where @xmath80 is the plasma frequency of the electrons and @xmath81 is the spacial separation of the grid points .
the simulations shown in this paper are described in table 1 .
these simulations are a subset of a higher number of runs , where we used different combinations of numerical and physical parameters , which confirmed numerical convergence .
we want to explore the regime @xmath82 . in order to do so ,
most of our runs use the mass ratio @xmath83 . since for @xmath83
both species will be essentially indistinguishable , we initially give ions and electrons maxwellian energy distributions with the same temperature , and use the parameters @xmath84 and @xmath85 to quantify the pressure of the plasma perpendicular and parallel to @xmath7 .
simulations with @xmath83 of course do not allow us to study the physics of electron isotropization .
however , in
[ sec : compmime ] we use an example of our @xmath86 runs to show that @xmath83 calculations reproduce the ion - related phenomena fairly well , with @xmath87 and @xmath88 playing the same role of @xmath89 and @xmath90 in the @xmath86 runs .
we defer a detailed study of electron scale pressure isotropization to future work .
we divide our runs into cases with three different initial betas : @xmath91 , 20 , and 80 .
our analysis focuses on the nonlinear structure of the ic / mirror generated fluctuations , the evolution of the pressure anisotropy @xmath92 , and the conservation of the ion magnetic moment .       runs .
panels @xmath93 and @xmath94 show the evolution of @xmath95 ( @xmath96 ; solid ) and @xmath97 ( @xmath98 ; dotted ) for runs beta6mag93 and beta6mag670a , respectively ( @xmath99 and @xmath100 correspond to black , red , and green , respectively ) .
for the same runs , panels @xmath50 and @xmath101 show @xmath102 ( @xmath103 ; black - solid ) , compared with the linear mirror and ic thresholds ( @xmath104 and @xmath105 , in red and green , respectively ) for pair plasma and growth rates @xmath106 and @xmath107 , respectively .
panels @xmath108 and @xmath109 show @xmath110 and @xmath111 , where the subscripts @xmath112 and @xmath113 denote the components parallel and perpendicular to @xmath28 , respectively .
after the plasma pressure anisotropy exceeds the mirror and ic thresholds , there is an initial phase of exponential growth followed by a secular phase .
the pressure anisotropy initially grows but saturates at the linear threshold for the mirror / ic instabilities.,width=340 ]     runs with @xmath114 and @xmath115 ( beta6mag93 and beta6mag670a in table [ table:1d ] ) .
we show the evolution of the true average magnetic moment @xmath116 @xmath117 ( black line ) and an effective global magnetic moment , @xmath118 @xmath119 ( red line ) , where the subscript @xmath120 denotes an average over all particles .
the ion magnetic moment @xmath116 is conserved until @xmath121 , at which point it decreases on the same timescale as the mean magnetic field ( see fig .
[ fig : muderivs ] ) .
the close similarity between @xmath116 and @xmath118 shows that there are not strong correlations between the particle velocity @xmath122 and magnetic field @xmath123 ; such correlations are more prominent at higher @xmath124 where the ic instability is subdominant ( see fig . [
fig : mubeta20]).,width=340 ]    , for different @xmath124 , compared to the growth rate of the mean magnetic field @xmath125 ( black line ; this quantity is the same for the different @xmath124 runs ) .
@xmath126 is shown for our simulations with @xmath127 and 80 ( beta6mag670a , beta20mag670a , and beta80mag1340 , respectively ) .
the vertical - dotted lines mark the beginning of the saturated state for each run based on when the growth of the fluctuations becomes roughly secular . in the saturated state ,
the magnetic moment changes on a timescale comparable to that of the mean magnetic field ( to within @xmath128 ) . at earlier times , however , the magnetic moment is reasonably well conserved.,width=283 ]    as discussed above , the expectation is that the ic instability should play an important role in the @xmath91 case , and it should become significantly less important in the @xmath129 and 80 cases . the contributions from the ic and mirror modes can be seen from figure [ fig : fieldsbeta6 ] .
this figure shows 2d images of the field fluctuations @xmath130 ( @xmath131 ; @xmath132 stands for  volume average " ) and the plasma density fluctuations @xmath133 for run beta6mag670@xmath50 (  @xmath134 " stands for the component @xmath38 , @xmath39 , or @xmath100 of @xmath135 ) at two different times : @xmath136 ( upper panels ) and @xmath137 ( lower panels ) .
we can see that at the earlier time ( @xmath136 ) the three components of @xmath135 have about the same amplitude ( see figures [ fig : fieldsbeta6]@xmath93 , [ fig : fieldsbeta6]@xmath94 , and [ fig : fieldsbeta6]@xmath50 ) , but appear to be dominated by different mode orientations , indicating the simultaneous presence of ic and mirror modes .
the mirror modes are expected to have @xmath135 mainly in the plane of the simulation ( the plane of @xmath13 and @xmath28 ; * ? ? ?
on the other hand , the ic modes have the three components of @xmath135 of comparable magnitude , so their presence can be most clearly revealed by @xmath138 ( figure [ fig : fieldsbeta6]@xmath50 ) . indeed , while @xmath139 and @xmath140 are dominated by a combination of oblique ( mirror ) and quasi - parallel ( ic )
modes , @xmath138 appears to be dominated by only quasi - parallel modes .
thus , initially the mirror and ic instabilities contribute comparably to @xmath135 . at the later time ( @xmath137 ) , however , the fluctuations become dominated by an oblique wave vector , with the ic instability playing a subdominant role .
the density fluctuations at all times are @xmath141 , and seem to correlate primarily with the mirror modes .
the relative contribution of the mirror and ic instabilities can also be seen from figure [ fig : beta6 ] , which shows the time evolution of different volume - averaged quantities for runs with @xmath142 ( first column ; run beta6mag93@xmath93 ) and @xmath143 ( second column ; run beta6mag670@xmath93 ) . in figures [ fig : beta6]@xmath93 and [ fig : beta6]@xmath94 we plot @xmath95 for the two magnetizations . in both cases
there is an initially exponential growth of @xmath144 until @xmath145 .
both instabilities have about the same growth rates , which are @xmath146 .
note that the dominant growth rate of the mirror / ic modes is nearly the same in the two simulations with different values of @xmath21 .
the growth rate is thus set by the background shear rate @xmath37 instead of @xmath21 .
this can be understood by noting that the smaller growth rate modes have a smaller anisotropy threshold .
therefore , as @xmath123 and @xmath147 grow , modes with smaller growth rate will begin to grow first .
this implies that the dominant modes will be those that can reach a significant amplitude to stop the growth of @xmath147 on the shear timescale @xmath65 .
after the exponential phase , the growth becomes dominated by the mirror modes and is secular ( as predicted by * ? ? ?
this transition from exponential to secular is clearly seen at @xmath148 for @xmath149 in figure [ fig : beta6 ] . despite the similar exponential growth rate for the mirror and ic instabilities in the linear regime ,
the mirror modes are more robust and dominate in the nonlinear regime .
this dominance of the mirror modes is remarkable because it occurs even for @xmath24 , where the linear analysis predicts an important contribution of the ic instability ( in the nonlinear regime the mean magnetic energy has been amplified by about one order of magnitude , so that @xmath24 ) .
we suggest below ( fig . [
fig : spectrabeta6 ] ) that departure from the bi - maxwellian energy distribution assumed in the linear instability analysis is likely playing an important role in this sub - dominance of the ic instability .
run beta6mag670a at two different times .
_ top : _ energy spectra as a function of the energies perpendicular and parallel to the magnetic field , @xmath150 ( black - dotted line ) and @xmath151 ( black - solid line ) . _
bottom : _ energy spectra as a function of @xmath152 .
bi - maxwellian fits are shown in red .
there is a clear deviation from the bi - maxwellian at the highest energies , which grows in time .
there is , however , little deviation from a maxwellian at @xmath153 , in contrast to some of the 1d runs discussed in the appendix .
comparison with figure [ fig : spectrabeta20 ] suggests that the high energy deviation from maxwellian is due to the ic instability , since it is less prominent at higher @xmath154 when the ic instability is sub - dominant.,width=340 ]    one difference in the evolution of @xmath135 for the two different magnetizations is that the end of the exponential phase occurs at somewhat higher amplitude in the @xmath142 case .
this difference implies that , for smaller values of @xmath43 , the mirror and ic fluctuations are less efficient in suppressing the pressure anisotropy , requiring larger values of @xmath144 . at the end of this section
, we will explain in further detail this dependence of @xmath144 on magnetization by focusing on the mechanism by which the isotropization occurs .
one consequence of the difference in @xmath144 is that the lower magnetization runs take longer to reach the point when @xmath102 ( @xmath103 ) saturates .
thus in the @xmath142 run the anisotropy @xmath102 has more time to grow , reaching larger values before saturating .
this can be seen in the evolution of @xmath102 , shown in black lines in figures [ fig : beta6]@xmath50 and [ fig : beta6]@xmath101 for runs beta6mag93@xmath93 and beta6mag670@xmath93 , respectively . for both magnetizations , there is an initial overshoot in @xmath102 , whose amplitude is larger for the @xmath142 run . after the overshoot
, @xmath102 evolves in a similar way for both magnetizations .    since in the nonlinear regime @xmath135
is dominated by mirror modes , after the saturation one would expect the pressure anisotropies to behave according to the marginal stability conditions of the mirror instability . since in our simulations both the ic and mirror modes grow at @xmath155 , we calculated the ic and mirror thresholds for equivalent growth rates in the pair plasma case with magnetizations @xmath142 and @xmath143 .
thus , figures [ fig : beta6]@xmath50 and [ fig : beta6]@xmath101 show that the threshold for mirror ( red ) and ic instability ( green ) in the cases @xmath156 and @xmath157 , respectively .
we see that the mirror threshold coincides fairly well with the saturated @xmath102 for both magnetizations .
these results imply that the linear marginal stability condition for the mirror instability is fairly accurate in determining @xmath102 in the nonlinear regime . in realistic astrophysical setups , where @xmath21 can be many orders of magnitude larger than @xmath37 ,
the expectation is that @xmath102 will follow the @xmath158 mirror threshold , given by @xmath159 .
this differs from the results shown in figure [ fig : beta6 ] , which are only appropriate for pair plasmas and modest magnetization .
( red ) for @xmath160 and @xmath149 ( run beta6mag670a ) .
the good agreement between the two results implies that the anisotropic stress coupling to the background shear is the primary mechanism for particle heating.,width=283 ]    as discussed in the introduction , it is not known whether the mirror instability saturates via pitch - angle scattering ( see , for instance , * ? ? ?
* ) or whether it cancels the appearance of a pressure anisotropy by substantially modifying the structure of the magnetic field @xcite . in the latter case
@xmath135 should grow secularly until @xmath161 , with the breaking of @xmath162invariance not necessary for the regulation of @xmath163 .
figures [ fig : beta6]@xmath108 and [ fig : beta6]@xmath109 show the evolution of @xmath111 and @xmath110 for @xmath114 and @xmath115 , respectively , where @xmath164 and @xmath165 refer to the magnetic field component perpendicular and parallel to the mean field @xmath28 .
we see that @xmath166 @xmath167 dominates in the nonlinear regime , as expected from the larger amplitude of the oblique mirror modes , and @xmath168 the saturation amplitude is @xmath169 , which favors the scenario where the non - conservation of @xmath116 is the key mechanism for the isotropization of the plasma pressure , at least if the background field has been amplified significantly .    to quantify the time variation of the average ion magnetic moment we define @xmath170 where @xmath171 stands for average over all the particles .
@xmath116 corresponds to the actual average of the magnetic moment .
thus if a nearly @xmath162conserving process dominates the mirror saturation , @xmath116 should remain fairly constant .
@xmath118 , on the other hand , can be thought as an _ effective _ global magnetic moment that ignores the fluctuations in @xmath172 or the correlation between @xmath122 and @xmath123 due to particles collecting in mirrors .
@xmath118 necessarily decreases since the mirror and ic instabilities suppress the growth that @xmath173 would have if the particles were only affected by the mean field @xmath28 .
thus , if @xmath116 is nearly conserved , @xmath116 and @xmath118 should evolve quite differently , with @xmath174 . on the other hand ,
if there is significant pitch angle scattering @xmath175 and both should decrease on the same timescale that @xmath176 increases .
figures [ fig : beta6b]@xmath93 and [ fig : beta6b]@xmath94 show @xmath116 and @xmath118 for our @xmath160 runs with magnetizations @xmath142 and @xmath143 , respectively .
we see that for both magnetizations the difference between @xmath116 and @xmath118 is very small ( only @xmath177 a few @xmath178 difference ) .
this implies that there is little spatial correlation between @xmath122 and @xmath123 .
we shall see below that @xmath116 and @xmath118 differ somewhat more at higher @xmath154 where the mirror instability dominates over the ic instability .
figure [ fig : muderivs ] compares the rate of change of the ion magnetic moment with that of the mean magnetic field for runs with @xmath1796 , 20 , and 80 .
this comparison is important since it is the evolution of the mean field that is driving the velocity space instabilities in our calculations . if the total thermal energy varies on a timescale long compared to the mean magnetic field , maintaining marginal stability to the mirror instability via pitch - angle scattering implies @xmath180 $ ] .
the term @xmath181 is exactly @xmath182 for the canonical high magnetization mirror instability threshold but is somewhat different in our pair plasma , modest magnetization simulations ( hence the use of @xmath183 ) .
this expression for @xmath184 implies that pitch - angle scattering should lead to the magnetic moment varying at a rate slightly less than that of the mean magnetic field .
this is consistent with our numerical results in figure [ fig : muderivs ] in the saturated state where the magnetic field grows secularly in time .    at earlier times , when the magnetic field fluctuations due to the mirror and ic instability are smaller
, @xmath116 is approximately conserved even though the mirror and ic are present and grow exponentially . in particular
, there is a temporal lag of @xmath185 between the onset of the exponential growth of the mirror and ic instabilities and the onset of pitch angle scattering that decreases @xmath116 .
our interpretation of this result is that the mirror / ic fluctuations must reach a sufficient amplitude , roughly @xmath186 , in order for pitch angle scattering to be effective . for the mirror instability ,
this violation of @xmath116 conservation by finite amplitude fluctuations is likely due to the stochasticity of particle orbits that sets in for large amplitude ( albeit low frequency ) fluctuations ( e.g. , @xcite ) . for the ion cyclotron instability ,
the need for finite amplitudes before significant scattering sets in is a consequence of the well - known result that the scattering rate for cyclotron - frequency fluctuations is @xmath187 .
the breaking of local @xmath116-invariance for this calculation ( @xmath87 , @xmath188 ) might be affected by the subdominant ( but still significant ) contribution of the ic modes to the magnetic fluctuations .
however , figure [ fig : muderivs ] shows that the magnetic moment varies at the same rate as the mean magnetic field for all of the @xmath124 we have simulated .
we will discuss the higher @xmath124 results in more detail in the next two subsections .
at @xmath189 for @xmath190 and magnetization @xmath114 ( panel @xmath93 ; run beta20mag93b ) and @xmath149 ( panel @xmath94 ; run beta20mag670b ) .
panels @xmath50 and @xmath101 show the same quantities at @xmath191 .
the black arrows show the direction of @xmath28 .
the higher magnetization run produces mirror modes whose wavevectors are more perpendicular to @xmath28 and with longer wavelengths
. a migration to larger wavelengths in time can also be seen.,width=340 ]    an interesting question is whether the non - linear evolution of the velocity space instabilities makes the energy distribution of the particles substantially different from a bi - maxwellian spectrum .
this question is particularly relevant in terms of understanding the suppression of the ic modes in the nonlinear regime as @xmath87 goes from 6 to @xmath192 .
indeed , it has been shown@xcite that a significant departure from a bi - maxwellian distribution can increase the anisotropy threshold for the growth of the ic modes and that a change of this kind is expected given the resonant character of the ic instability .
figure [ fig : spectrabeta6 ] shows that the particle spectrum does differ somewhat from a bi - maxwellian for our @xmath160 run , mainly due to a significant excess of particles at the highest energies . in future work
, we will study whether or not the spectral variation in figure [ fig : spectrabeta6 ] can completely account for the decrease in the amplitude of the ic modes in the saturated state .
figure [ fig : spectrabeta6 ] also shows that the ion velocity distribution does not deviate significantly from a bi - maxwellian at low @xmath152 , as has been found in some 1d studies of the saturation of the mirror instability ( e.g. , @xcite ) .
we discuss this in more detail in the appendix .
finally , figure [ fig : anistress ] shows the evolution of the total heating rate of the particles as a function of time ( black line ) . for comparison , in the red line
we plot the particle heating rate by the anisotropic stress , @xmath193 @xcite .
we see a reasonable agreement between the two results , implying that the anisotropic stress plays the dominant role in particle heating .
other mechanisms , like wave - particle interaction , must play a secondary role in the total particle energization .
this is true for all of the @xmath124 we have simulated .
we now consider a somewhat weaker magnetic field case with @xmath129 . in this case , the field fluctuations are dominated by the mirror instability at all times . this can be seen in figures [ fig : beta20b2s]@xmath93-@xmath50 which show the evolution of the three components of @xmath194 for @xmath190 runs with magnetizations @xmath43=93 , 670 , and 2000 , respectively ( runs beta20mag93@xmath93 , beta20mag670@xmath93 and beta20mag2000@xmath93 in table 1 ) . in all three cases
@xmath195 is small compared to @xmath194 , showing the subdominant role of the ic instability .
figures [ fig : beta20b2s]@xmath93-@xmath50 also show the characteristic transition between exponential growth ( with @xmath196 ) and secular growth ( with @xmath197 ) for the mirror modes .    as in the @xmath91 case ,
the transition between the exponential and secular regimes occurs at smaller amplitudes and in a smoother way in the more strongly magnetized runs .
figures [ fig : beta20b2s]@xmath101-@xmath109 also show that , as in the case of @xmath91 , after the overshoot the pressure anisotropy follows the marginal stability condition for the mirror modes reasonably well , particularly at higher magnetization .
figures [ fig : beta20b2s]@xmath198-@xmath199 show @xmath200 and @xmath201 as a function of time .
while @xmath200 is about the same for the three magnetizations , @xmath201 is significantly smaller in the higher magnetization runs .
the relative magnitude of @xmath202 and @xmath203 is a measure of the orientation of the mirror s dominant wave vector , @xmath13 , with respect to @xmath28 : a large @xmath204 ratio implies that @xmath13 and @xmath28 are nearly perpendicular .
the results of figures [ fig : beta20b2s]@xmath198-@xmath199 imply that for larger magnetization , @xmath13 and @xmath28 are more perpendicular , which is consistent with linear calculations @xcite .     and
@xmath94 show the evolution of the ion magnetic moment @xmath116 for @xmath190 runs with @xmath114 ( beta20mag93a ; panel @xmath93 ) and @xmath149 and @xmath205 ( beta20mag670a and beta20mag2000 ; panel @xmath94 ) .
panels @xmath50 and @xmath101 compare the corresponding rate of change of @xmath116 with that of the mean magnetic field for runs beta20mag93a and beta20mag670a , respectively .
all plots show the evolution of the true average magnetic moment @xmath116 @xmath117 ( black line ) and an effective global magnetic moment , @xmath118 @xmath119 ( red line ) , where the subscript @xmath120 denotes an average over all particles .
the magnetic moment begins to decrease when the fluctuations reach @xmath206 ( see fig . [
fig : beta20b2s ] ) .
the difference between @xmath116 and @xmath118 is more significant than for @xmath160 ( cf fig .
[ fig : beta6b ] ) .
this is a consequence of particles bunching in mirrors , which leads to a correlation between @xmath122 and @xmath123 .
nonetheless , the ion magnetic moment changes at about the same rate as the background magnetic field after saturation at @xmath207 , indicating that pitch angle scattering regulates the nonlinear saturation on sufficiently long timescales.,width=340 ]    the orientation of the modes in the nonlinear stage can also be seen directly in figure [ fig : fieldsbeta20 ] , which shows the spatial distribution of @xmath208 at two different times for the runs with @xmath142 and @xmath143 .
figures [ fig : fieldsbeta20]@xmath93 and [ fig : fieldsbeta20]@xmath94 correspond to the time @xmath136 , while figures [ fig : fieldsbeta20]@xmath50 and [ fig : fieldsbeta20]@xmath101 correspond to @xmath137 . for each magnetization ,
the mirror fluctuations are dominated by two modes that are symmetric with respect to the magnetic field , and that are more oblique at higher magnetization .
note also that the wavelength of the dominant modes ( in units of the larmor radius of the particles , @xmath58 ) is larger for larger magnetization and that the modes tend to grow in wavelength as time goes on .
figure [ fig : mubeta20 ] shows the evolution of the ion magnetic moment @xmath116 and the rate of change of @xmath116 relative to that of the mean magnetic field for our @xmath190 calculations .
a comparison of figures [ fig : mubeta20 ] & [ fig : beta20b2s ] shows that the magnetic moment is reasonably conserved until the magnetic fluctuation amplitude is close to ( though somewhat less than ) its saturated value . as in our @xmath160 calculations
, this implies that there is a time lag of @xmath209 between the onset of the mirror instability and the onset of significant pitch - angle scattering ( and its associated decrease in @xmath116 ) .
note , however , that for @xmath129 the saturated amplitude of the fluctuations in the secular regime , @xmath210 , is larger than in the @xmath91 case , for which @xmath211 .
this implies that when _ only _ the mirror instability plays a significant role , a somewhat larger fluctuation amplitude is necessary for efficient pitch - angle scattering .
we will see below that @xmath210 continues even for @xmath212 ; the saturated amplitude of the mirror modes is thus fairly independent of the beta of the plasma as long as the ic instability is not significant ( @xmath213 6 ) .
figure [ fig : mubeta20 ] also shows that for @xmath190 , the true ion magnetic moment @xmath116 decreases somewhat more slowly than the effective global magnetic moment @xmath118 ( by @xmath214 ) .
this was not seen in our @xmath160 calculations ( figs .
[ fig : beta6b ] ) .
the modest difference between @xmath116 and @xmath118 indicates that the nonlinear saturation of the mirror instability involves correlations between @xmath122 and @xmath123 , i.e. , particles bunching in mirrors .
this effect is independent of magnetization , as can be seen from the solid and dotted lines in figure [ fig : mubeta20]@xmath94 ( corresponding to @xmath215 670 and 2000 , respectively ) .
run beta20mag670a at two different times during the nonlinear regime .
_ top : _ energy spectra as a function of the energies perpendicular and parallel to the magnetic field , @xmath150 ( black - dotted line ) and @xmath151 ( black - solid line ) .
_ bottom : _ energy spectra as a function of @xmath152 . in all cases , a bi - maxwellian energy spectrum ( red ) provides a good approximation to the numerical simulations.,width=340 ]    figure [ fig : spectrabeta20 ] shows that the particle energy spectrum stays fairly close to bi - maxwellian in our @xmath190 runs , regardless of the magnetization .
since these simulations are almost purely mirror dominated , this suggests that the mirror modes provide rather momentum - independent pitch - angle scattering to the particles . as a result ,
it is likely that ic modes produce the deviations from a bi - maxwellian distribution seen in figure [ fig : spectrabeta6 ] .
figure [ fig : spectrabeta20 ] also shows that the parallel velocity spectrum @xmath216 remains nearly maxwellian during most of the saturated regime ( with only a transient @xmath217 decrease for @xmath218 especially at the end of the exponential growth regime ) .
the lack of significant flattening of @xmath216 ( which contrasts with the prediction of * ? ? ?
* ) is analyzed in further detail in the appendix .
for @xmath212 we focus on simulations with @xmath143 and 1340 , respectively ( runs beta80mag670 and beta80mag1340 in table 1 ) .
the results are essentially the same as for the @xmath129 case so we do not discuss them in detail .
figure [ fig : beta80b2s ] shows that the mirror instability dominates the fluctuations in @xmath7 .
the evolution of the pressure anisotropy @xmath102 contains the initial overshoot followed by a rather flat behavior .
the amplitude of the overshoot is controlled by the magnetization , while the subsequent evolution is fairly well described by the threshold condition for mirror modes
.     runs with @xmath149 ( beta80mag670 ; left column ) and @xmath219 ( beta80mag1340 ; right column ) .
panels @xmath93 and @xmath94 show the evolution of @xmath95 ( @xmath96 ; solid ) and @xmath97 ( @xmath98 ; dotted ) , respectively ( @xmath99 and @xmath100 correspond to black , red , and green , respectively ) .
for the same runs , panels @xmath50 and @xmath101 show @xmath102 ( @xmath103 ; black - solid ) , compared with the linear mirror threshold in red ( @xmath104 ) for a pair plasma and growth rates @xmath107 and @xmath220 , respectively .
panels @xmath108 and @xmath109 show @xmath110 and @xmath111 .
the results shown here are very similar to the analogous @xmath129 results in figure [ fig : beta20b2s ] , consistent with the dominance of the mirror instability at high @xmath154.,width=340 ]    as in the @xmath129 case , the saturation @xmath221 is @xmath222 .
figure [ fig : mubeta80 ] shows the evolution of the ion magnetic moment for the @xmath223 simulations .
the evolutions of @xmath116 and @xmath224 are the same for the two magnetizations presented ( @xmath149 and @xmath219 ) , and show the same behavior seen in the @xmath129 case , which is discussed in detail in the previous subsection .     and
@xmath94 show the evolution of the ion magnetic moment @xmath116 for our @xmath223 runs with @xmath149 and @xmath219 ( runs beta80mag670 and beta80mag1340 , respectively ) .
panels @xmath50 and @xmath101 compare the rate of change of @xmath116 with that of the mean magnetic field . in all panels ,
we show the evolution of the true average magnetic moment @xmath116 @xmath117 ( black line ) and an effective global magnetic moment , @xmath118 @xmath119 ( red line ) , where the subscript @xmath120 denotes an average over all particles .
the magnetic moment begins to decrease when the fluctuations reach @xmath206 ( see fig . [
fig : beta80b2s ] ) .
all of the key results here are very similar to the @xmath129 results shown in figure [ fig : mubeta20].,width=340 ]      , @xmath225 , and @xmath226 ( run betai20magi240mass10 ) .
the ion results here are similar to the corresponding @xmath227 results in figures [ fig : beta20b2s ] and [ fig : mubeta20 ] , indicating that the pair plasma results capture most of the key ion physics .
panel @xmath93 shows the evolution of @xmath228 and @xmath229 , in solid and dotted lines , respectively .
the small contribution of @xmath195 implies the dominance of mirror modes .
panel @xmath94 shows the evolution of @xmath111 ( red ) and @xmath110 ( black ) with @xmath230 , implying that the mirror modes have @xmath231 .
panel @xmath50 shows the evolution of the ion pressure anisotropy @xmath232 ( black - solid ) . the linear mirror threshold for mirror modes with @xmath233 in a @xmath234 ,
@xmath235 plasma is also shown ( red line ) along with the analogous threshold for @xmath236 for comparison ( red dotted line ) .
panel @xmath101 shows the evolution of @xmath116 ( @xmath237 , black - solid ) and @xmath118 ( @xmath238 , red - solid ) ; see eq .
[ eq : mu ] and associated text for details.,width=340 ]    in this section we show that our use of @xmath83 in the previous sections does not have a significant effect on the evolution of the ion pressure anisotropy under the influence of the ic and mirror instabilities . in figure
[ fig : beta20b2smime4 ] we present the results for a simulation with @xmath235 , @xmath225 , and @xmath226 ( run betai20magi240mass10c ) , which is analogous to the previous @xmath83 simulations that use @xmath239 .
a comparison of figure [ fig : beta20b2smime4]@xmath93 with figures [ fig : beta20b2s]@xmath93 , [ fig : beta20b2s]@xmath94 , and [ fig : beta20b2s]@xmath50 shows that for both mass ratios the evolution of @xmath135 is dominated by the mirror instability , with the initial exponential regime being followed by secular growth .
figure [ fig : beta20b2smime4]@xmath94 shows the evolution of @xmath111 and @xmath110 . as in the @xmath83 case ( see figures [ fig : beta20b2s]@xmath198 , [ fig : beta20b2s]@xmath240 , and [ fig : beta20b2s]@xmath199 ) , @xmath230 in the secular stage , implying that mirror modes with @xmath231 are dominant .
the maximum amplitude of @xmath241 is also in good agreement with the @xmath83 simulations .
the evolution of the ion pressure anisotropy @xmath232 is shown in figure [ fig : beta20b2smime4]@xmath50 .
as in the @xmath83 runs ( see figures [ fig : beta20b2s]@xmath242 ) , @xmath243 follows the linear threshold of the mirror instability fairly closely .
this is seen in figure [ fig : beta20b2smime4]@xmath50 , which shows the linear mirror threshold expected for modes with @xmath233 , @xmath225 , and @xmath235 and 1836 ( red and dotted - red lines , respectively ) .
the small difference between our results and the @xmath235 and 1836 thresholds strongly suggests that the use of @xmath83 for most of our study captures reasonably well the key physics of the mirror and ic instabilities in a real electron - ion plasma .
finally , figure [ fig : beta20b2smime4]@xmath101 shows the evolution of @xmath116 ( black solid line ) and @xmath118 ( red solid line ) .
the results are quantitatively similar to the @xmath227 results in figures [ fig : mubeta20]@xmath93 and [ fig : mubeta20]@xmath94 .
this shows that the role played by pitch - angle scattering in the saturation of the mirror instability is well captured by our @xmath83 simulations .
we have used pic simulations to study the nonlinear , long - term evolution of ion velocity space instabilities in collisionless plasmas .
we have focused on instabilities driven by pressure anisotropy with @xmath244 : the ion cyclotron ( ic ) and mirror instabilities .
these instabilities are expected to arise when turbulence and/or mhd instabilities amplify @xmath7 , rendering @xmath244 due to conservation of the ion magnetic moment @xmath4 .
alternatively , in the solar wind @xmath245 can increase due to perpendicular heating and/or parallel cooling of the plasma .
in contrast to previous studies , we do not consider the initial value problem of the evolution of a given initial ion pressure anisotropy .
instead , we self - consistently induce the growth of an ion pressure anisotropy by continuously amplifying the mean magnetic field in our computational domain .
this setup allows us to study the long - term , saturated state of the plasma , rather than the ( much shorter ) initial , exponential stage .
the former is of most interest for the heliospheric and astrophysical applications of our work .
we have focused on the regime @xmath246 and plasma magnetizations @xmath43 in the range @xmath247 ( the magnetization is the ratio of the initial ion cyclotron frequency to the background shear frequency ) .
most of our simulations used @xmath83 , but we showed that our results for the ion scale physics are fairly independent of the mass ratio ( see  [ sec : compmime ] ) .
so long as @xmath45 , our results are also relatively independent of the precise value of the magnetization ( see , in particular , figs .
[ fig : beta20b2s ] & [ fig : mubeta20 ] ) .
thus we believe that our basic conclusions about the saturation of the mirror and ic instabilities in continuously driven systems can be applied to physical problems where @xmath21 is many orders of magnitude larger than @xmath37 .
our primary results are as follows :    * for the @xmath246 regime we have studied , the turbulent fluctuations in the saturated state are dominated by mirror modes .
this is despite the fact that ic and mirror instabilities are both present at early phases ( and have similar linear instability thresholds ) when @xmath248 a few .
a plausible explanation for this dominance of the mirror instability in the saturated state even at @xmath249 a few is that the ion velocity distribution departs from a bi - maxwellian spectrum in the presence of the ic instability .
this can in turn increase the threshold for the ic instability @xcite making it less important than the mirror instability in the nonlinear regime .
we indeed see more significant departures from a bi - maxwellian in the low @xmath30 regime where the ic instability is present ( figs .
[ fig : spectrabeta6 ] & [ fig : spectrabeta20 ] ) .
by contrast , the ion distribution function remains relatively bi - maxwellian in the high @xmath30 regime in which the mirror instability dominates over the ic instability at all times ( fig .
[ fig : spectrabeta20 ] ) .
the dominance of the mirror instability at @xmath248 a few is consistent with the nonlinear study of @xcite , which makes use of hybrid simulations of an expanding box .
* small - scale fluctuations driven by the mirror instability initially grow exponentially but saturate in a secular phase in which @xmath250 ( figures [ fig : beta6 ] & [ fig : beta20b2s ] ) .
the fluctuations in the saturated state have @xmath231 ( see , e.g. , figure [ fig : fieldsbeta6 ] ) consistent with linear theory expectations for mirror modes . *
the ion pressure anisotropy @xmath232 in the saturated state is reasonably well described by the marginal stability condition for the mirror instability ( e.g. , figs [ fig : beta6 ] & [ fig : beta20b2s ] ) .
note , however , that the marginal stability state of the mirror mode in our simulations is not the same as in real systems because of the comparable ion and electron masses and the smaller @xmath43 used in our calculations .
this is important to bear in mind when comparing our results to heliospheric measurements . *
the total thermal energy of the plasma increases in time at the theoretically predicted rate given the background shear in the plasma and the anisotropic stress associated with the pressure anisotropy ( fig .
[ fig : anistress ] ) .
this is consistent with the idea that the `` viscous '' heating rate in a collisionless plasma is set by how velocity space instabilities regulate the pressure anisotropy @xcite . * when the growth time of the mirror instability is long compared to the ion cyclotron period  as is the case in our simulations where the shear slowly amplifies the mean magnetic field
the ion magnetic moment is initially reasonably well conserved .
this is fundamentally because there is no mechanism to violate magnetic moment conservation during the linear phase of mirror growth .
so the modes continue to grow to larger and larger amplitudes , given the sustained free energy source created by the background velocity shear amplifying the magnetic field and continuously generating pressure anisotropy .
this regime lasts until @xmath3 and @xmath251 , and is consistent with the  secular regime " found by @xcite where @xmath116 is nearly conserved .
* saturation of the mirror modes starts to happen once the fluctuations attain @xmath3 . in this stage
the variation in magnetic field on the scale of the larmor radius is sufficiently large that the ion magnetic moment is no longer conserved .
this is consistent with theories of stochastic ion motion in large amplitude turbulent fluctuations ( e.g. , * ? ? ?
* ; * ? ? ?
in this phase the ion magnetic moment changes on nearly the same timescale as the mean magnetic field ( fig . [
fig : muderivs ] ) , which is consistent with pitch angle scattering maintaining marginal stability . in this stage
the magnetic fluctuations only experience a modest additional growth , reaching a maximum amplitude of @xmath252 . *
the need for finite amplitude fluctuations with @xmath3 in order to stop magnetic moment conservation implies the existence of a the temporal delay of @xmath185 between the onset of the ic and mirror instabilities and the onset of efficient pitch angle scattering ( @xmath37 is the rate at which the mean magnetic field is growing ) .
our results have a number of consequences for modeling the large scale dynamics of nearly collisionless astrophysical plasmas and for interpreting heliospheric measurements .
we briefly mention a few of these applications but defer a more detailed analysis to future work .
our results are consistent with a number of in situ heliospheric measurements .
for instance , our maximum saturated value of @xmath253 is in reasonable agreement with mirror modes observations in different environments like planetary magnetosheaths ( where large amplitude mirror modes are frequently found .
see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , the solar wind ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , and the heliosheath ( e.g. , * ? ? ?
* ; * ? ? ?
. corresponds to @xmath254 , where @xmath132 stands for volume average .
this means that our results still allow for _ local _ magnetic fluctuations of even larger amplitude .
it is also important to stress that if a plasma is rapidly driven to have a temperature anisotropy well in excess of the mirror threshold ( e.g. , at the earth s bow shock ; * ? ? ? * ) then the resulting amplitude of the mirror modes could significantly exceed what we find in our calculations .
in this case , a standard initial value calculation of the mirror evolution is likely to better capture the resulting dynamics than our model in which the anisotropy slowly increases in time . ]
our results also show that the dominant wavelength is initially close to a few ion larmor radii , with a subsequent migration to longer wavelengths through a coalescence process .
this wavelength growth is consistent with observations ( e.g. , * ? ? ?
* ; * ? ? ?
* ) and with theoretical models that predict the growth of mirror mode wavelengths due to magnetic diffusion @xcite .
the fact that our @xmath255 calculations find that the mirror instability dominates the nonlinear state of the fluctuations ( even though the ic and mirror instabilities have comparable linear thresholds and growth rates for a bi - maxwellian distribution ) is also consistent with heliospheric observations .
in particular , solar wind measurements at @xmath256 show that @xmath257 is limited by the linear threshold of the mirror instability ( instead of the ic instability .
see * ? ? ?
* ; * ? ? ?
the observations of @xcite also show enhanced magnetic fluctuations and magnetic compressibility that follow the linear mirror threshold fairly well .
it is important to emphasize , however , that the observations presented by @xcite and @xcite showing dominance of the mirror instability also include the regime @xmath258 , which is not considered in our study .
also , the dominance of the mirror instability could be affected by properties of the solar wind that are usually not included in calculations of instability thresholds .
these include the presence of @xmath259 of alpha particles ( see , e.g. , figure 7 of * ? ? ?
* ) , the drift of alpha particles with respect to protons @xcite , and the frequent presence of a proton beam @xcite .
we interpret the dominance of the mirror modes in our simulations as being due to the suppression of the ic waves , which is caused by the ion velocity distribution function departing from bi - maxwellian @xcite .
this interpretation also requires the ic modes to have finite amplitudes ( although smaller than the one of the mirror modes ) so that they can maintain , via pitch - angle scattering , an ion distribution function that is marginally stable to the growth of ic waves . to the best of our knowledge
, these ic modes have not been found by in - situ heliospheric measurements during mirror dominated events .
this is consistent with the fact that in our simulations the mirror and ic modes only have comparable amplitudes during the initial , exponential regime ( fig .
[ fig : beta6]@xmath93 and @xmath94 ) . in the saturated regime ( which is likely to be more representative of the observed state of the plasma in the heliospheric and astrophysical systems ) , we find that the mirror modes have significantly larger amplitudes
. moreover , this dominance increases with the magnetization @xmath43 , as can be seen by comparing figures [ fig : beta6]@xmath93 and @xmath94 . can be accounted for by the fact that the scattering rate of cyclotron - frequency fluctuations is expected to be @xmath187 . if we equate this to the shear parameter @xmath37 ( which is , roughly , the rate at which pitch - angle scattering takes place , as can be seen from the changing rate of @xmath116 in fig .
[ fig : muderivs ] ) , one obtains that @xmath260 . ]
thus , for realistic values of @xmath43 , we do not anticipate easily detectable ic waves by in - situ measurements in cases where the mirror modes dominate .
the departure from the bi - maxwellian ion distribution function has not been confirmed by in - situ heliospheric measurements of mirror - dominated events either .
we believe that this is due to the high sensitivity of resonant instabilities to gradients in the velocity distribution function .
thus , in cases where ic and mirror modes can have similar growth rates , only small deviations from the bi - maxwellian distribution can significantly reduce the importance of the ic modes and give way to the dominance of the mirror modes .
we , therefore , do not anticipate easily detectable departures from bi - maxwellian associated to this effect . nevertheless , future careful observations comparing the nature of the fluctuations and the ion distribution functions in the solar wind can help to unveil these effects .
a second application of our results is to low - collisionality accretion flows onto black holes . in these systems pressure anisotropies
can produce an anisotropic stress that can increase the transport of angular momentum and the heating of the plasma @xcite .
our results are consistent with @xcite s hypothesis that the mirror and/or ic instabilities would provide enhanced pitch angle scattering and regulate the evolution of the plasma pressure tensor .
axisymmetric kinetic simulations have also demonstrated this directly @xcite though the restriction to axisymmetry limits the quantitative applicability of these initial simulations to real disks .
one of the key conclusions of our study in this paper is that even modest plasma magnetization is sufficient to quantitatively capture the nonlinear saturation of the mirror instability .
for example , our @xmath261 results in figures [ fig : beta20b2s ] & [ fig : mubeta20 ] show that the evolution of the ion pressure anisotropy , ion magnetic moment , and turbulent fluctuations are quite similar for @xmath114 and @xmath149 .
the primary difference is that in the lower magnetization run the plasma has a larger `` overshoot '' of the saturation pressure anisotropy set by the linear mirror instability threshold .
this weak dependence on magnetization is encouraging because it is computationally infeasible to reach @xmath149 in three - dimensional kinetic simulations of astrophysical plasmas ( e.g. , accretion disks ) .
one of the interesting questions not fully addressed by our work is the saturation of velocity space instabilities in turbulent low - collisionality plasmas where the shear rate of the mean magnetic field itself fluctuates in time and the magnitude of the large - scale turbulent fluctuations satisfies @xmath262 .
our calculations in this paper have focused on the case where an imposed velocity shear amplifies a background magnetic field by order unity or more . if , on the other hand , the shear acts for less time such that the background field is amplified by significantly less , the plasma may never become mirror / ic unstable ( depending on @xmath30 and @xmath221 ) or the background shear may reverse sign before the mirror / ic fluctuations have reached sufficiently large amplitudes for pitch angle scattering to set in .
this regime will be explored in future work .
another key question not addressed by our study is isotropization of the electron pressure anisotropy . in future work
we will study how fast electrons interact with the large - scale mirrors generated by the ion pressure anisotropy , and whether there are separate electron - scale fluctuations that produce significant electron pitch - angle scattering ( e.g. , the whistler instability ) .
we thank tobi heinemann , ben chandran , phil isenberg and , in particular , matt kunz for fruitful discussions .
eq is supported in part by nasa htp grant nnx11aj37 g , nsf grant ast-1333682 , a simons investigator award from the simons foundation , the david and lucile packard foundation , and the thomas alison schneider chair in physics .
mr thanks support from the chilean comisin nacional de investigacin cientfica y tecnolgica ( conicyt ; proyecto fondecyt iniciacin n@xmath263 11121145 ) .
dv is supported in part by nasa htp grant nnx11aj37 g and nasa grant nnx12ab27 g .
99    anderson , b. j. , fuselier , s. a. , & murr , d. 1991 , geophysical research letter , 18 , 1955 balbus , s. a. , hawley , j. f. 1991 , , 376 , 214 balbus , s. a. 2001 , , 562 , 909 bale , s. d. , kasper , j. c. , howes , g. g. , quataert , e. , salem , c. , & sundkvist , d. 2009 , phys .
103 , 211101 baumgartel , k. , sauer , k. , & dubinin , e. 2003 , geophys .
, 30 , 1761 califano , f. , hellinger , p. , kuznetzov , e. a. , passot , t. , sulem , p. l. , & travnicek , p. m. 2008 , j. geophys . res . , 113 , a08219 chen , l. , lin , z. , & white , r. 2001 , phys .
of plasma , 8 4713 johnson , j. r. , & cheng , c. z. 2001 , geophys .
lett . , 28 , 4421 buneman , o. 1993 , `` computer space plasma physics '' , terra scientific , tokyo , 67 enriquez - rivera , o. , blanco - cano , x. , russell , c. t. , jian , l. k. , & luhmann , j. g. 2010 , aip conf . proc . , 1216 , 276 gary , s. p. 1992
, j. geophys .
res . , 97 , 8519 gary , s. p. & lee , m. a. 1994 , j. geophys .
res . , 99 , 11297 guo , j. , yang , z. , lu , q. , & wang , s. 2009 , plasma sci .
technol . , 11 , 274 hasegawa , a. 1969 , phys .
fluids , 12 , 2642 hasegawa , a. & tsurutani , b. t. 2011 , , 107 , 245005 hawley , j. f. , gammie , c. f. , & balbus , s. a. 1995 , , 440 , 742 hellinger , p. , travnicek , p. , kasper , j. c. , & lazarus , a. j. 2006 , geophys .
lett . , 33 , 9 hellinger , p. & travnicek , p. m. 2008 , j. geophys .
res . , 113 , a10109 hellinger , p. , kuznetzov , e. a. , passot , t. , sulem , p. l. , & travnicek , p. m. 2009 , geophys . res .
lett . , 36 , l06103 heuer , m. & marsch , e. 2007 , , 112 , a03102 horbury , t. s. & lucek , e. a. 2009 , , 114 , a05217 isenberg , p. a. 2012 , phys .
plasmas 19 , 032116 isenberg , p. a. , maruca , b. a. , & kasper , j. c. 2013 , , 773 , 164 joy , s. p. , kivelson , m. g. , walker , r. j. , khurana , k. k. , russell , c. t. , & paterson , w. r. 2006 , j. geophys .
res . , 111 , a12212 kunz , m. , schekochihin , a. a. , & stone , j. m. 2014 , phys . rev
112 , 205003 liu , y. , richardson , j. d. , belcher , j. w. , kasper , j. c. , & skoug , r. m. 2006 , j. geophys .
res . , 111 , a09108 lyutikov , m. 2007 , , 668 , l1 marsch , e. , muhlhauser , k. -h . ,
rosenbauer , h. , schewnn , r. , & neubauer , f. m. 1982 , j. geophys .
res . , 87 , 35 marsch , e. , muhlhauser , k. -h . ,
schewnn , r. , rosenbauer , h. , pilipp , w. & neubauer , f. m. 1982 , j. geophys .
res . , 87 , 52 maruca , b. a. , kasper , j. c. & bale , s. d. 2011 , phys .
107 , 201101 matteini , l. , landi , s. , hellinger , p. , pantellini , f. , maksimovic , m. , velli , m. , goldstein , b. , & marsch , e. 2007 , geophys .
lett . , 34 , l20105 matteini , l. , hellinger , p. , landi , s. , travnicek , p. m. , & velli , m. 2007 , space sci .
172 , 373 narayan , r. , mahadevan , r. , & quataert , e. 1998 , theory of black holes acretion disks pokhotelov , o. a. , sagdeev , r. z. , balikhin , m. a. , & treuman , r. a. 2004 , j. geophys .
, 109 , a09213 quataert , e. 2008 , , 673 , 758 remya , b. , reddy , r. v. , tsurutani , b. t. , lakhina , g. s. , & echer , e. 2013 , jgra , 118 , 785 riquelme , m. a. , quataert , e. , sharma , p. , & spitkovsky , a. 2012 , , nn , nn russell , c. t. , blanco - cano , x. , jian , l. k. , & luhmann , j. g. 2009 , geophys .
lett . , 36 , l05106 schekochihin , a. a. , cowley , s. c. , kulsrud , r. m. , hammett , g. w. , & sharma , p. 2005
, , 629 , 139 schekochihin , a. a. , cowley , s. c. , kulsrud , r. m. , rosin , m. s. , & heinemann , t. 2008 , phys .
lett . , 100 , 081301 schmid , d. , volwerk , m. , plaschke , f. , voros , z. , zhang , t.l . ,
baumjohann , w. , & narita , y. 2014 , ann .
geophys . , 32 , 651 schwartz , s. j. , burgess , d. , & moses , j. j. , 1996 , ann .
, 14 , 1134 sharma , p. , hammett , g. w. , quataert , e. , & stone , j. 2006 , , 637 , 952 shoji , m. , omura , y. , tsurutani , b. t. , verkhoglyadova , o. p. , & lembege , b. 2009 , , 114 , a10203 soucek , j. , lucek , e. , & dandouras , i. 2008 , j. geophys . res . , 113 , a04203 southwood , d.  j. , & kivelson , m.  g. 1993 , , 98 , 9181 spitkovsky , a. 2005 , aip conf .
proc , 801 , 345 , astro - ph/0603211 travnicek , p. , hellinger , p. , taylor , m. , escoubet , c. , dandouras , i. , & lucek , e. 2007 , geophys .
lett . , 34 , l15104 tu , c .- y . & marsch , e. 2002 , , 107 , a9 tsurutani , b. t. , smith , e. j. , anderson , r. r. , ogilvie , k. w. , scudder , j. d. , baker , d. n. , & bame , s. j. 1982 , , 87 , 6060 tsurutani , b. t. , lakhina , g. s. , verkhoglyadova , o. p. , echer , e. , guarnieri , f. l. , narita , y. , constantinescu , d. o. 2011 , , 116 , a02103 tsurutani , b. t. , echer , e. , verkhoglyadova , o. p. , lakhina , g. s. , guarnieri , f. l. 2011 , journal of atmospheric and solar - terrestrial physics , 73 , 1398 verscharen , d. , bourouaine , s. , chandran , b.  d.  g. , & maruca , b.  a.  2013 , , 773 , 8 volwerk , m. , zhang , t. l. , delva , m. , voros , z. , baumjohann , w. , & glassmeier , k.h .
2008 , geophys .
lett . , 35 ,
previous 1d studies of the mirror instability suggest that its growth and saturation are dominated by the  resonant " scattering of @xmath218 particles ( where @xmath264 is the particle velocity component parallel to the magnetic field ; * ? ? ?
* ; * ? ? ?
* ) , leading to a flattening of the distribution function near @xmath218 .
in contrast , our 2d runs discussed in the main text show that scattering by mirror modes is fairly momentum independent , with no significant flattening of @xmath265 . to compare our results more directly to the previous literature , in this appendix we compare the mirror saturation in 1d and 2d for different values of @xmath154 , paying especial attention to the behavior of @xmath218 particles .
we show that , compared with the 1d case , 2d runs significantly change the way pitch - angle scattering occurs , modifying the saturation mechanism .
to best compare with the previous literature , we use initial - value simulations in this appendix , where a linearly unstable pressure anisotropy is imposed as an initial condition .
figure [ fig:1d2d2.3 ] compares simulations with @xmath266 , @xmath267 , @xmath268 in 1d and 2d . in the 1d simulation , the initial magnetic field @xmath34 forms an angle of @xmath269 with the resolved dimension of the simulation . apart from the different number of dimensions
, these simulations also have the same numerical parameters .
panels [ fig:1d2d2.3]@xmath93 and [ fig:1d2d2.3]@xmath94 show @xmath116 and @xmath224 as a function of time for the 1d and 2d cases , respectively . in 1d
, @xmath116 decreases significantly slower than @xmath224 , which means that there is a significant correlation between @xmath163 and @xmath123 ( as expected for the mirror instability ) . in 2d , on the other hand , the difference between @xmath116 and @xmath224 is very small , implying that pitch - angle scattering occurs more efficiently in 2d .
this difference in the scattering efficiency can be explained considering the modes that grow in the 1d and 2d simulations .
indeed , in the 1d case , mirror modes with wave vectors pointing in only one direction can grow . in the 2d case ,
on the other hand , both ic and mirror modes appear ( with the mirror modes forming both positive and negative angles with respect to @xmath34 ) .
this can be seen from panels [ fig:1d2d2.3]@xmath50 and [ fig:1d2d2.3]@xmath101 , which show the evolution of the magnetic energy contained in fluctuations parallel ( @xmath167 ) and perpendicular ( @xmath270 ) to @xmath34 for the 1d and 2d cases . whereas the 1d case is dominated by @xmath167 ( characteristic of mirror modes ) , the 2d case is dominated by @xmath270 , showing the dominant presence of the ic modes .
panels [ fig:1d2d2.3]@xmath108 and [ fig:1d2d2.3]@xmath109 show snapshots of the parallel velocity distribution @xmath265 for the 1d and 2d runs at the saturated state @xmath271 ( black lines ) , and compare them with a bi - maxwellian distribution ( red line ) . in the 1d case
, we can see flattening of @xmath265 for @xmath218 .
this flattening is consistent with the results of @xcite , who argue that it is due to the resonant interaction of @xmath218 particles with the mirror modes .
the 2d run , however , does not have this flattening , implying that the saturation process in two dimensions is quite different .
this is indeed expected from the dominant role of the ic instability .        since at higher @xmath272 the ic instability
is expected to play a subdominant role , we also explored the case where initially @xmath273 .
figure [ fig:1d2d20 ] compares simulations with @xmath266 , @xmath274 , and @xmath275 in 1d and 2d , with both cases having the same numerical parameters .
panels [ fig:1d2d20]@xmath93 and [ fig:1d2d20]@xmath94 show the evolution of @xmath116 and @xmath224 for both cases . in the 1d case
these two quantities differ significantly , implying an important correlation between @xmath123 and @xmath163 .
this correlation almost disappears in the 2d case , which shows the presence of more efficient pitch - angle scattering in the 2d case .
since for @xmath273 the ic instability plays a subdominant role , this is likely a property of the mirror instability itself rather than the presence of the ic instability as was the case for the lower @xmath154 results shown in figure [ fig:1d2d2.3 ] . the higher rate of scattering in 2d also affects the maximum amplitude of the fluctuations , which are shown in panels [ fig:1d2d20]@xmath50 and [ fig:1d2d20]@xmath101 for the 1d and 2d cases .
although in both cases @xmath276 ( implying that mirror dominates ) , in 2d the amplitude of the modes is smaller , which means that field rearrangement plays a less important role in the approach to marginal stability .    finally , figure [ fig:1d2d20 ] also compares @xmath265 at the saturated state ( @xmath277 ) in the 1d and 2d cases .
there is very little flattening for @xmath218 in both cases .
this confirms our result from the main text that the pitch - angle scattering in the mirror dominated case is fairly independent of particle momentum .
it is interesting , however , that we do not see flattening even in 1d . since @xcite reported the appearance of flattening only near the threshold of the mirror instability ( and only for certain box sizes ) , we tested simulations with the same plasma parameters having anisotropies as low as @xmath278 and found similar results .