however , the role played by each of several acceleration and loss mechanisms is not yet established observationally [ millan and thorne , 2007 ] .
precipitation into the atmosphere is considered to be one of the major electron loss mechanisms , which can completely deplete the radiation belts of electrons during the main phase of some geomagnetic storms [ o'brien et al . , 2004 ;
selesnick , 2006 ] . using data from balloon - borne x - ray detectors , foat et al . , lorentzen et al . , and millan et al .   have reported precipitating relativistic electrons in the dusk sector .
these duskside precipitation events occur over a variety of magnetic activity levels [ lorentzen et al . , 2000 ;
kokorowski et al . , 2008 ] and a broad radial distribution ranging from l=38 [ millan et al . , 2013 ] .
the energy spectrum has been found to be well fit by an exponential distribution with an e - folding energy ranging from 0.5 to 3.6 mev [ millan et al . ,
2002 ] . based primarily on the fact that all these events were found at dusk ,
the cause has been suggested to be the gyroresonant scattering by electromagnetic ion cyclotron ( emic ) waves .
emic waves are observed throughout the inner magnetosphere but predominantly on the duskside and dayside [ anderson et al .
, 1992a , 1992b ; meredith , 2003 ; erlandson and ukhorskiy , 2001 ; fraser et al . , 1996 , 2006 ; usanova et al . , 2012 ] .
they are excited by anisotropic ring current ions injected into the inner magnetosphere [ e.g. , jordanova et al . , 2008 ] or by compressions of the magnetopause [ e.g. , anderson and hamilton , 1993 ] .
the linear growth rate of emic waves maximizes in high - density regions such as the duskside plasmapause or plasmaspheric drainage plume due to reduced resonant energies [ cornwall et al . , 1970 ;
horne and thorne , 1993 ] and wave guiding by steep density gradients near the plasmapause [ horne and thorne , 1993 ; jordanova et al .
field - aligned emic waves interact with relativistic electrons through the gyroresonance condition 1where e is electron gyrofrequency ,  is the relativistic factor , k|| and v|| are components of the wave propagation vector particle velocity along the direction of the ambient magnetic field .
emic waves are expected to effectively scatter relativistic electrons of geophysically interesting energies , preferably in duskside high - density regions , where the minimum energy limit is relatively low and the diffusion rate is close to the strong diffusion limit [ thorne and kennel , 1971 ; albert , 2003 ; summers , 2003 ] .
duskside simultaneous proton and relativistic electron precipitation has also been observed , supporting the theory of precipitation caused by emic wave scattering [ e.g. , bortnik et al . ,
quasi - linear theory [ kennel and petschek , 1966 ] has been the dominant treatment ; however , recent work shows that nonlinear effects can be very significant and even reverse the conclusions [ albert and bortnik , 2009 ] .
therefore , it is important to establish when and where each approach is applicable and to test these models with observations . existing theoretical discussions on emic waves as a possible loss mechanism mainly focus on the evolution of the trapped electron flux and its timescale . however , directly measured by balloon detectors is the energy spectrum of the bremsstrahlung x - rays produced by the precipitating electrons , and many satellites also measure the energy of the precipitating electrons [ millan et al . , 2007 ] .
it is very important to understand the energy dependence of the precipitating flux driven by emic waves , as well as the time evolution of this dependence in order to test the theory with observations .
this paper uses the quasi - linear formulation to evaluate the energy and time dependence of rep , and investigates how different parameters affect the diffusion coefficients and the energy spectrum .
when applied with input wave and particle data from satellites ( e.g. , van allen probes , goes ) , the results from our model can be directly compared with balloon ( e.g. , balloon array for rbsp ( radiation belt storm probes ) relativistic electron losses ) and low - altitude satellite ( e.g. , solar anomalous and magnetospheric particle explorer ( sampex ) and polar - orbiting operational environmental satellites ( poes ) ) measurements to investigate the role of emic waves in causing rep as well as the effectiveness of the adopted theoretical model .
we use quasi - linear diffusion theory to model the evolution of the distribution of electrons due to interactions with emic waves .
radial and energy diffusion are ignorable because the frequency of the emic waves is well above the drift frequency and well below the gyrofrequency of the resonant particles [ kennel , 1966 ] .
the bounce - averaged diffusion equation for pure pitch angle scattering can be written as [ davidson and walt , 1977 ; lyons , 1973 ; lyons and williams , 1984 ; shao et al . ,
2009 ] 
 2where t(0)=1.38020.3198(sin(0)+ sin1/2(0 ) ) is the normalized bounce time , f0 is the trapped electron phase space density , 0 is the equatorial pitch angle , d(0,e) is the bounce - averaged pitch angle diffusion coefficient , and atm is the timescale for losses to the atmosphere . assuming that losses occur only from within the loss cone and that the loss cone is emptied twice per bounce period , we take atm to be half of the bounce period inside the loss cone and infinite outside the loss cone [ lyons , 1973 ; davidson and walt , 1977 ; lyons and williams , 1984 ; shao et al . ,
the differential flux is related to the phase space density through the electron momentum p as j0=pf0 . to capture the feature of the isotropic flux distribution in the loss cone in the many strong diffusion cases in our parameter studies
, we set the boundary conditions to be . for weak diffusion cases , although does not have to be zero , the loss cone is essentially empty and a small number of electrons in the loss cone does not significantly affect the diffusion [ shprits et al . ,
, 2009 ] 
 3with an arbitrary scaling outside the loss cone and j0(t=0)=0 inside the loss cone . here
e is the particle energy in mev and 0l is the equatorial loss cone angle given by sin2(0l)=(4l3l ) at a particular l shell .
we calculate the bounce - averaged diffusion coefficient using the same method as summers   and summers et al .   for parallel - propagating emic waves in a multi - ion ( h , he , and o ) plasma .
the diffusion coefficient is corrected with a factor of 2 [ albert , 2007 ] .
the earth 's magnetic field is assumed dipolar , the wave amplitude is taken to be 1 nt , and the wave frequency spectrum is assumed to be a truncated gaussian , namely , 
 4where m is the center frequency , with lower frequency limit 1=m and upper frequency limit 2=m+.  is a measure of the bandwidth .
the emic waves are assumed to be confined to 15% in latitude and 10% of the electron drift orbit , propagating in a cold plasma of a storm time ion composition 70% h , 20% he , and 10% o [ meredith , 2003 ] .
the wave dispersion relation is 
 5where k is the wave number k1<k < k2 . the suffix j denotes the ion species ;
the values j=1,2 , and 3 refer to h , he , and o , respectively .
the minimum resonant energy can be approximated as [ summers , 2003 ] 
 6 figure  1 shows the bounce - averaged diffusion coefficient as a function of equatorial pitch angles we calculated for 0.15 mev electrons with four sets of equatorial background magnetic field strength b0 , cold plasma density n0 , and wave frequency . each diffusion coefficient spans from 0 to an upper resonant pitch angle limit 0u and maximizes to a value d(0,e)m at a pitch angle 0 m .
bounce - averaged pitch angle diffusion coefficients for emic waves interacting with electrons of energies 0.15 mev with an increment of 0.2 mev and four different sets of wave and plasma parameters .
we next investigate how the diffusion coefficient varies with energy and pitch angles for a full range of parameters b0 , n0 , and . the varying diffusion coefficients will determine how the energy spectrum evolves ( section  4 ) .
we choose b0 , n0 , and  to vary from 75 to 1150 nt ( corresponding to l38 ) , 1 to 1000 cm ( assumed to be constant over the entire interaction region ) , and 1.9 to 3.9 o+ ( covering roughly the entire helium band , where the spectral intensity of emic waves is enhanced [ thorne et al . , 2006 ; anderson et al . , 1992a ; hu et al . ,
2010 ] and pitch angle diffusion of geophysically interesting relativistic electrons is more likely to occur [ li et al . , 2007 ] ) , respectively . since emic waves are often observed near the geosynchronous orbit where n0100 cm [ halford et al . ,
2010 ] , we set the control b0 and n0 to be 100 nt and 100 cm , respectively .
we also choose the control wave frequency to be =2.43.4  o+ , estimated from wave observations [ meredith , 2003 ; ukhorskiy et al . , 2010 ] . in figure  2
( first to fourth rows ) ( the fifth row will be discussed in section  4 ) , we plot d(00,e),d(0,e)m , 0 m , and 0u as a function of e to characterize the diffusion coefficient when one parameter is varied while the others are held fixed . satisfying the gyroresonance condition equation  1 , the dispersion relation equation  5 , and the gaussian wave frequency distribution equation  4 ,
the plots show that for all b0 , n0 , 1 , and 2 , the pitch angles 0u and 0 m increase monotonically with energy , whereas d(00,e) and d(0,e)m increase from the minimum resonant energy up to a certain  peak  energy and then decrease . the value of the peak energy decreases with decreasing b0 or increasing n0 , 1 , or 2 .
the width of the peak significantly narrows with increasing 1 and 2 , meaning that at high frequencies , the diffusion coefficient sharply increases at the peak energy from surrounding energies . keeping energy fixed , as b0 decreases or n0 , 1 , or 2 increases , 0 m increases monotonically ; the diffusion coefficient increases at low energies yet decreases at high energies ( variations at > 5 mev are not shown in the figure ) ; d(0,e)m increases at low energies and decreases at high energies only with increasing 2 while increases at low energies and approaches the same values at high energies with decreasing b0 or increasing n0 or 1 ( variations at > 5 mev are not shown in the figure ) .
lastly , the minimum resonant energy emin decreases while 0u increases with decreasing b0 or increasing n0 or 2 , but emin and 0u are unaffected by 1 . ( first to fourth rows ) d(00,e) , d(0,e)m , 0 m , and 0u are shown , respectively .
( fifth row ) precipitating flux at t=1 s derived from initial distribution equation  9 .
( left column ) varying b0 , n0=100 cm , 1=2.4  o+ , 2=3.4  o+ ; ( middle column ) varying n0 , b0=100 nt ,
1=2.4  o+ , 2=3.4  o+ ; and ( right column ) varying  , b0=100 nt , n0=100 cm are shown .
applying the results for the diffusion coefficients from section  3 , we solve equation  2 for the trapped equatorial electron flux j0(0,e , t ) .
as we can see , the loss cone fills up quickly ( e.g. , the first 60 s in figure  3 ) before gradually being depleted as the total flux drops due to loss to the loss cone .
evolution of the pitch angle distribution of the trapped flux of 4 mev electrons in the simulation time 0 to 5 min .
b0 , n0 , 1 , and 2 are the same as figure  1(c )
. at b0=200 nt ( l6.7 ) , the loss cone is 2.3. at 0 min , the trapped flux distribution is the initial maxwellian flux distribution equation  3 . the omnidirectional flux at any latitude  on a dipole field line can be expressed in terms of equatorial flux and pitch angles as [ lyons and williams , 1984 ] , 
 7where . if we take  to be the latitude where the field line intersects the atmosphere boundary , then b/b0=(4l3l ) , =0lc and no particle bounces back to the equator .
the omnidirectional relativistic electron precipitation ( rep ) flux measured at the atmosphere boundary will then be 
 8jrep plotted versus e is then the energy spectrum of the precipitating electrons . to study how the jrep(e , t ) energy spectra are affected by the diffusion coefficients , we first apply an energy - independent initial distribution 
 9and examine the precipitation spectrum in the first second , i.e. , jrep(e , t=1s ) ( figure  2 , fifth row ) .
these plots strongly resemble those of d(00,e) ( figure  2 , first row ) , indicating that the precipitation at the outset is largely controlled by the diffusion coefficients at small pitch angles .
next , we again use the maxwellian initial distribution equation  3 and study how the spectrum is affected by time .
we found two types of spectra ( figure  4 ) : in most cases , the spectrum is singly peaked ( figure  4a ) , but occasionally , increasing wave frequency can induce another peak ( figure  4b ) . as shown in the figures , initially , both spectra resemble maxwellian distributions with a lower cutoff at the minimum resonant energy .
the precipitation quickly builds up from zero within the first second as the loss cone is being filled .
shortly after , the rep flux near the minimum resonant energies largely drops and the rep flux starts to gradually decrease as the loss cone is being depleted .
this rapid increase and slow decrease of the precipitation is consistent with many balloon observations [ millan et al . , 2007 ] .
however , in figure  4a , the peak of the rep flux keeps moving toward higher energies and the curves flatten out . in figure  4b , the peak stays at a roughly constant energy , followed by another peak forming at a later time at a higher energy and moving toward even higher energies .
the hardening of both spectra with time is due to the fact that 1 with increasing energy , diffusion becomes increasingly dominated by particles with larger pitch angles ( 0 m and 0u both increase with energy ) and they are scattered into the loss cone more slowly than the ones with smaller pitch angles ; 2 over time , with fewer particles left to interact with the waves , the precipitation is largely reduced at lower energies associated with shorter lifetimes ( larger diffusion coefficients ) and therefore the spectrum becomes harder and more isotropic ( flatter ) in energy .
furthermore , as we discussed in the previous section , the increase of the diffusion coefficient at the peak energy from surrounding energies is significantly enhanced with frequency . when the increase is sharp enough
, the precipitation of the particles with the highest diffusion coefficients can drastically reduce the flux and a trough can form in the middle of the spectrum . in figure  4b ,
no rep flux is produced above 4 mev since the diffusion coefficients in the loss cone are zero ( figure  1d ) .
the energy distribution of the precipitating electrons at time 1 s to 50 min and the time evolution of the energies of the rep peaks assuming maxwellian initial distribution equation  3 .
( a and b ) b0 , n0 , 1 , and 2 are the same as figures  1c and 1d , respectively .
note the x axis left end in figure  4a is higher than 100 kev for better visualization . in figure  5 , we show the variation of the spectra of rep at 30 min for the full range of parameters . in the first few seconds ( not shown ) , similar to the case of an initial distribution uniform in energy , with decreasing b0 or increasing n0 , 1 , or 2 , jrep increases at low energies but decreases at high energies .
but at 30 min , in the case of lower b0 or higher n0 , 1 , or 2 , the rep flux curves become fairly flat or doubly peaked , and sometimes intersect and fall below the curves of higher b0 or lower n0 , 1 , or 2
( e.g. , the precipitation with b0=125 nt is lower than that of 175 nt at 2.6 mev in figure  5a ; the precipitation with =3.4 - 3.9 o+ is lower than any other precipitation curve of 2=3.9  o+ at below 1 mev in figure  5c ) . when this happens , the role played by the changing parameter at these energy ranges is different than at early times .
our simulation also shows that increasing n0 ,  or decreasing b0 lowers the energies of the peaks in both single- or double - peak spectra .
energy spectra of precipitating electron flux at t=30 min with maxwellian initial distribution , plotted with the same parameters and color codes adopted in figure  2 .
( c ) the rep flux of 3.43.9 o+ is zero above 3.6 mev where diffusion coefficients are zero in the loss cone . in figure  6 , we integrate jrep in the energy ranges 0.10.3 , 0.31 , and 15 mev and vary two of the three parameters b0 , n0 , and  simultaneously . these spectra can be compared with satellites with similar energy channels ( e.g. , poes and sampex ) [ yando et al . ,
the curves of the same color ( for  the same fixed parameter ) get shorter and shorter , and can completely disappear when the energy of the channel is below the minimum resonant energies . with an increase in n0 or  or a decrease in b0 , the rep flux increases , except occasionally at late times when
the roles played by the parameters reverse ( explained in the previous paragraph ) , in which case strong initial scattering results in the significant decrease of the precipitation as the loss cone is depleted . precipitating electron flux in energy channels
( a ) 1=2.4  o+ , 2=3.4  o+ , b0=751150  nt , color - coded in n0 ranging from 25 to 1000  cm .
( b ) n0=100  cm , b0=751150  nt , color - coded in  ranging from 1.93.9  o+ .
( c )  b0=100  nt , n0=11000  cm , color - coded in  ranging from 1.9 to 3.9  o+ .
the goal of this paper is to investigate the shape of the energy spectrum of relativistic electron precipitation ( rep ) due to quasi - linear interactions with emic waves , and how it varies with time and changing parameters such as background magnetic field strength b0 , cold plasma density n0 , and wave frequency . the key results can be summarized as follows : 
 1 .
the rep energy spectrum is generally peaked , with a lower cutoff at the minimum resonant energy . over time
the second peak appears later at a higher energy and moves toward even higher energies.3 . in both single- and double - peak cases ,
the precipitating flux first rapidly increases as the loss cone is being filled , and then slowly decreases as the loss cone is being depleted.4 .
increasing n0 ,  , or decreasing b0 lowers the minimum resonant energy and the energies of the peak(s ) .
it causes the precipitation to increase at low energies and decrease at high energies . over time ,
when strong scattering slows down , the role played by the changing parameter is altered.5 .
we study how it is affected by changing parameters through varying two parameters at the same time .
the integrated flux monotonically increases with increasing n0 and  and decreases with increasing b0 , except at the energies where the role of the parameter reverses when strong scattering causes a large reduction of the particles interacting with the waves ( figure  6 ) .
the rep energy spectrum is generally peaked , with a lower cutoff at the minimum resonant energy . over time
, the peak moves toward higher energies and the spectrum flattens ( gets harder ) .
the second peak appears later at a higher energy and moves toward even higher energies .
the precipitating flux first rapidly increases as the loss cone is being filled , and then slowly decreases as the loss cone is being depleted .
4 . increasing n0 ,  , or decreasing b0 lowers the minimum resonant energy and the energies of the peak(s ) .
it causes the precipitation to increase at low energies and decrease at high energies . over time ,
when strong scattering slows down , the role played by the changing parameter is altered .
we study how it is affected by changing parameters through varying two parameters at the same time .
the integrated flux monotonically increases with increasing n0 and  and decreases with increasing b0 , except at the energies where the role of the parameter reverses when strong scattering causes a large reduction of the particles interacting with the waves ( figure  6 ) . to better explain the variation of the precipitation energy spectrum , we show how the spectrum is affected by three deterministic factors  the diffusion coefficient , the initial trapped flux , and time .
the diffusion coefficient is determined by the gyroresonance condition , the wave dispersion relation , and the wave frequency distribution and is a function of pitch angle .
we calculate the diffusion coefficient for a range of input parameters b0 , n0 , and  and characterize the pitch angle dependence of the diffusion coefficient with d(00,e) , the upper pitch angle limit 0u , and the maximum diffusion coefficient d(0,e)m along with its corresponding pitch angle 0 m .
we show that the diffusion coefficient maximizes at increasing pitch angles with increasing energy and the upper limit 0u also increases .
we also show that low b0 , high n0 and high  lower the minimum resonant energy , increase the diffusion coefficients of low energy particles and reduce the diffusion coefficients of high energy particles , and further increase the precipitation at low energies and decrease the precipitation at high energies at the beginning .
the number of particles scattered into the loss cone per unit of time increases with the number of initially trapped particles .
therefore , when we switch the initial trapped flux energy distribution from uniform to maxwellian , the spectrum becomes slanted with an enhancement at the lower energies , and the peak of the spectrum shifts to the left .
therefore , in the beginning , the precipitation is strongly affected by the diffusion coefficients at low pitch angles .
since the diffusion coefficients tend to extend out and maximize at higher pitch angles with increasing energy , the precipitating flux of high energies has a growing relative significance over time , and the precipitation energy spectrum gets harder .
in addition , if the scattering at a certain energy is initially strong , the particles at that energy are quickly lost and the precipitation will be largely reduced over time , and this also results in the rapid decrease of the precipitation at lower energies and the formation of the trough regions in the doubly peaked spectra .
it is worth noting that thermal heating [ anderson and fuselier , 1994 ; thorne et al . ,
2006 ] and the high plasma beta in the outer edge of the ring current during storm times [ lui et al . , 1987 ] may render the cold dispersion relation impractical .
several studies [ isenberg , 1984 ; chen et al . , 2011 ; silin et al . , 2011 ; chen et al . ,
2013 ] in which the hot dispersion relation is adopted suggest that at the vicinity of the ion cyclotron frequencies , the finite - beta effect may lead to the damping of the emic waves and the increase of the minimum resonant energies of relativistic electron scattering . based on their conclusions , in the situations where the warm / hot ions are abundant
, we should expect to see a higher minimum resonant energy cutoff in our rep energy spectrum if the wave frequencies are just below the helium gyrofrequency .
whether the doubly peaked spectra will still exist is uncertain because in our simulations , they usually only happen at high frequencies close to the helium gyrofrequency where the diffusion coefficient is significantly modified by the finite - beta effect .
these model energy spectra show what we should expect to observe given various wave and geomagnetic environmental conditions if the precipitation is caused by emic wave scattering and can be simulated by the adopted diffusion model .
further work will include event studies in which we will take the input data from satellites ( e.g. , van allen probes and goes ) and compare the simulated precipitating flux spectra with those detected by conjugate balloons ( e.g. , barrel ) and low altitude satellites ( e.g. , sampex and poes ) .
 [ 1]previous studies on electromagnetic ion cyclotron ( emic ) waves as a possible cause of relativistic electron precipitation ( rep ) mainly focus on the time evolution of the trapped electron flux 
 . however , directly measured by balloons and many satellites is the precipitating flux as well as its dependence on both time and energy . 
 therefore , to better understand whether pitch angle scattering by emic waves is an important radiation belt electron loss mechanism and whether quasi - linear theory is a sufficient theoretical treatment , we simulate the quasi - linear wave - particle interactions for a range of parameters and generate energy spectra , laying the foundation for modeling specific events that can be compared with balloon and spacecraft observations . we show that the rep energy spectrum has a peaked structure , with a lower cutoff at the minimum resonant energy 
 . the peak moves with time toward higher energies and the spectrum flattens . the precipitating flux , on the other hand , first rapidly increases and then gradually decreases . 
 we also show that increasing wave frequency can lead to the occurrence of a second peak . in both single- and double - peak cases 
 , increasing wave frequency , cold plasma density or decreasing background magnetic field strength lowers the energies of the peak(s ) and causes the precipitation to increase at low energies and decrease at high energies at the start of the precipitation . 
1. Introduction
2. Diffusion Equation
3. Diffusion Coefficients
4. Time Evolution of Trapped and Precipitating Electron Fluxes
5. Summary and Discussion