the radiation mechanisms of many high energy astrophysical objects are still an active issue since they often contain features which are not easily explained by the conventional synchrotron and inverse compton emissions .
recently , much attention has been paid to the radiation signatures from the turbulent electromagnetic fields ( e.g. , medvedev 2000 , fleishman 2006 , kelner , aharonian , & khangulyan 2013 , mao & wang 2013 , teraki & takahara 2013 ) .
the main scene of the emission regions of high energy astrophysical objects is collisionless shocks , and the turbulent electromagnetic fields would be generated in the shock region .
therefore , the electromagnetic turbulence should be taken into account when we consider the radiation .
however , for the major emission mechanisms of synchrotron radiation and inverse compton scattering , effects of small scale turbulence are not taken into account .
there is a room for novel emission signatures in the consideration of turbulence , which may be of relevance to observations . by reproducing the observed spectra
, we can extract physical parameters of astrophysical objects .
thus , researches of radiation signatures from the turbulent field would play a key role for the understanding of physical mechanisms of the high energy astrophysical objects .
radiation spectra from a small scale turbulent magnetic field have been well studied as the `` jitter radiation '' or `` diffusive synchrotron radiation '' ( medvedev et al .
2011 , fleishman & urtiev 2010 and references therein ) .
differences from the synchrotron radiation become significant when the typical spatial scale of an eddy @xmath20 is smaller than the photon formation length ( pfl ) of the synchrotron photons @xmath21 , where @xmath22 is the larmor radius and @xmath23 is the lorentz factor of the electron .
such a small scale magnetic field is thought to be generated by weibel instability around the shock front .
when the strength of this small scale turbulent magnetic field is dominant , the radiation spectra are determined by the turbulence and reveal various signatures of the turbulent field . for example , when @xmath24 the peak frequency in @xmath25 spectrum becomes @xmath26 , and the spectrum in the highest frequency region shows a power law @xmath27 when the turbulence exists up to the maximum wave number @xmath28 .
the power law index @xmath29 reflects that of the turbulent magnetic field @xmath30 .
the spectra show more complex signatures when @xmath31 ( medvedev 2011 , reville & kirk 2010 , teraki & takahara 2011 ) .
we note that the anisotropy of turbulence also affects the radiation spectra ( e.g. kelner et al .
2013 , reynolds & medvedev 2012 , medvedev 2006 ) .    the radiation from an electron which moves in non - uniform magnetic field in a laboratory
is well studied using the insertion device of synchrotron orbital radiation factory , where a series of magnets are line - upped to make the particle deflect periodically .
it is called `` wiggler '' or `` undulator '' ( jackson 1999 ) . for undulator , the strength of magnets @xmath32 and gaps between them @xmath0 are chosen to satisfy the condition that the observer is always in the beaming cone . on the other hand , for wiggler
, the observer is periodically in and off the beaming cone .
we estimate the critical distance @xmath33 which divides wiggler and undulator .
the deflection angle in one deflection is @xmath34 , where @xmath35 is the typical curvature radius of the orbit .
the radiation from a relativistic particle is concentrated into small cone with opening angle @xmath19 .
therefore , the critical condition dividing wiggler and undulator is @xmath36 , which is rewritten as @xmath37 .
thus , the device is called undulator when @xmath38 , while it is called wiggler when @xmath39 .
the radiation spectrum of undulator shows a sharp peak at @xmath40 , while wiggler shows a broad spectrum with peak frequency @xmath41 .
the relation between typical frequencies and deflection angle is a key point for understanding of the radiation spectra .
perturbative jitter radiation or perturbative dsr is recognized as extensions of the undulator radiation , since the spatial scale of turbulence @xmath0 is assumed to be much smaller than @xmath42 .
the fact that @xmath43 determines the radiation feature was also noted in astrophysics in early seventies by , e.g. , rees , ( 1971 ) , and getmantsev & tokarev ( 1972 ) .
rees assumed that the strong electromagnetic wave is emitted from the crab pulsar with frequency @xmath44 ( 30hz ) according to the `` oblique rotator '' model by ostriker & gunn ( 1969 ) , and argued that the radiation from an electron moving in the strong wave is not inverse compton scattering ( with frequency @xmath45 ) but synchrotron - like ( synchro - compton ) radiation , because the deflection angle @xmath43 is estimated as @xmath46 , where the factor @xmath47 comes from the ratio of cyclotron frequency to the wave frequency @xmath48 .
the magnetic field strength @xmath49 is estimated by equating the spindown luminosity and power of magnetic dipole radiation of the crab pulsar and by considering the distance from the pulsar .
he argued that the radiation signature of synchro - compton radiation resembles that of the synchrotron radiation .
the typical frequency is determined by the field strength as @xmath50 , and the low frequency spectral index is roughly @xmath18 .
getmantsev & tokarev ( 1972 ) studied the radiation spectra under general electromagnetic fields and stated that the radiation signature from a single charged particle is determined by the frequency or wavelength of the field for @xmath51 , while it is determined by the field strength for @xmath52 in the same way as synchrotron radiation .
note that they did not discussed explicit expressions of the radiation spectra from a single electron for @xmath53 , they presented a general expression of radiation spectrum from an ensemble of electrons with a power spectrum .
thirty years after these pioneering works , radiation spectra from a relativistic particle interacting with turbulent fields are gathering a renewed attention .
for example , the radiation mechanism from langmuir turbulence treated in the present paper has become an interesting topic .
the langmuir turbulence has been pointed out to be generated around the shock front of the relativistic shocks ( silva 2006 , dieckmann 2005 , bret , dieckmann , & deutch 2006 ) , so it should be as important as the radiation from turbulent magnetic field .
fleishman & toptygin ( 2007a , b ) have made a systematic treatment of diffusive radiation in langmuir turbulence ( see references cited in fleishman & toptyigin 2007b for other earlier relevant works ) .
their method is the most sophisticated treatment for langmuir turbulence , which is based on toptygin & fleishman ( 1987 ) .    for later discussions , we shortly review their treatment .
they treat the electron motion by a statistical approach and use a perturbative treatment for calculation of the radiation . the calculation formula for the radiation spectra
is based on the one written in landau & lifshitz ( 1971 ) . for @xmath54 , a rectilinear trajectory with constant velocity is assumed , but non - zero acceleration from the external field is taken into account .
the wavenumber of the langmuir waves is assumed to satisfy the condition @xmath55 , where @xmath4 is the plasma frequency .
they calculate the correlation between the acceleration and the langmuir waves .
the peak frequency is @xmath11 .
the spectrum shows an abrupt cutoff above the peak , and becomes @xmath56 in higher frequencies when the turbulence exists up to the maximum wavenumber @xmath57 for a power law turbulent spectrum @xmath58 .
the spectrum just below the peak is @xmath59 , and becomes @xmath60 in lower frequencies and @xmath61 in even lower frequency region .
spectrum is predicted in the lowest frequency region , which comes from the effect of wave dispersion in the plasma .
we do not discuss such effects in this paper . ]
@xmath61 spectrum comes from the effect of multiple scattering .
the angle between the velocity and observer direction becomes larger than the beaming cone after many deflections even when the deflection angle in one deflection @xmath43 is much smaller than @xmath62 , and the approximation of rectilinear trajectory is broken .
therefore , this treatment is beyond the perturbative treatment , and they call it `` non - perturbative treatment '' .
they treat the changing of direction of motion in many deflections by diffusion approximation . in consequence
, a spectral break emerges in the low frequency region with a suppression of low frequency photons .
the spectrum becomes @xmath61 from this effect and the index of @xmath63 comes from the diffusivity .
they claimed that the break frequency approaches the peak frequency as @xmath43 becomes large and the break and peak merge for @xmath53 .
they stated that even for @xmath52the spectrum in frequency region just below the peak of @xmath64 is @xmath61 .
this is inconsistent with the statement by getmantsev & tokarev ( 1972 ) .
this is one of the objectives of investigations in the present paper .
the effect of large angle deflection would come into play in forming the radiation spectra for langmuir turbulence as in the wiggler radiation for @xmath65 .
this point has been discussed for magnetic turbulence in kelner et al .
( 2013 ) , medevev et al . (
2011 ) , and teraki & takahara ( 2011 ) . in this paper
, we investigate general properties of the radiation spectra from relativistic electrons in a langmuir turbulence for various cases including this regime .    before proceeding to the formulation of the calculation employed in the present paper , we point out the differences between the magnetic field generated by weibel instability ( entropy modes ) and the langmuir waves .
the first is time variability .
the entropy mode does not oscillate , so the typical timescale is determined by the turnover time of an eddy .
it is longer than the crossing time of a relativistic electron if we assume that the background plasma is sub - relativistic .
therefore , we can treat the magnetic field as a static field when we calculate the radiation spectra for the zeroth - order approximation . on the contrary
, we should not treat the langmuir turbulence as a static field even for the zeroth - order approximation .
the crossing time can be comparable to or longer than the period of the langmuir waves , because the typical spatial scale is about inertial length @xmath66 , governed by the plasma frequency @xmath4 ( diekmann 2005 ) .
therefore , the time variability of the electric field can not be ignored .
this effect is well studied by fleishman & toptygin ( 2007a , b ) .
the second is the energy change of the radiating electrons .
for the case of turbulent magnetic field , the energy of electrons is conserved if we ignore the radiation back reaction .
however , the energy change can not be ignored for the langmuir turbulence , because the electric field can accelerate the electrons parallel to their velocity .
the lorentz factor of the electron can change in a short time by strong langmuir waves .
this effect may play a role for calculation of the radiation spectra .
the calculation of electron trajectory for the strong turbulence by analytical approach is hard to perform .
thus , we calculate the radiation spectra by numerical approach from first principle .
we calculate radiation spectra for a wide range of the field parameters .
we study about three factors , the scale length , time variability , and strength . by sweeping the parameter plane , we generally investigate the radiation spectra from a relativistic electron moving in langmuir turbulence .
in section 2 , we describe calculation method , and we show the results in section 3 . in section 4 , we give the physical interpretations of the discovered spectral features using radiation for a spatially uniform plasma oscillation . in section 5 , we make a summary and some discussions .
we use lienard - wiechert potential directly to calculate the radiation spectra .
it is the same approach as employed in teraki & takahara 2011 for magnetic turbulence .
this method is computationally expensive because the time step is restricted by observed frequency and the integration time has to be longer than the photon formation time ( pft ) .
pft is defined as the coherence time for forming a photon with a given frequency ( reville & kirk 2010 , akhiezer & shulga 1987 ) .
therefore , pft is a function of the frequency .
since pft is @xmath67 times larger than the inverse of the observed frequency and since the time step should be shorter than pft , the integration time is at least @xmath67 times larger than the inverse of the observed frequency . as a result
, large number of time steps are needed to calculate radiation spectra for highly relativistic case .
although less heavy approach is proposed by reville & kirk 2010 to overcome this problem by using appropriate approximations , we use the first principle method we show below , because it does not include any approximations in the calculation formula of the radiation spectra from prescribed trajectory . therefore , this method is suitable for studying the cases for which any approximations are difficult to adopt . by using moderately relativistic particles ,
we reduce the computational expenses .
we neglect radiation back reaction of electrons , because the strength of electric field is sufficiently weaker than @xmath68 ( cf .
jackson 1999 ) .
we use 3d isotropic langmuir turbulence in this paper .
if the langmuir waves are generated near the shock front , they may be highly anisotropic on the spot .
however , we assume the isotropic turbulence for two reasons .
first , the langmuir turbulence in large part of the emission region would be isotropic , since langmuir waves interact with background ions and form the kolmogorov type isotropic turbulence ( treumann & baumjohan 1997 ) .
second , the radiation spectra for 3d and 1d turbulence are not so different except for a particular case that the direction of radiating electrons and the wave vector are parallel . for this case ,
the radiation is linear acceleration emission .
however , emission from the perpendicular acceleration dominates as long as the direction differs from the wave vector direction even slightly ( fleishman & toptygin 2007a )
.    if the shock is ultra - relativistic , the downstream plasma may be relativistically hot . for simplicity
, we assume a mildly relativistic shock , so that the downstream plasma is sub - relativistic .
we ignore the thermal velocity of background plasma in the dispersion relation of langmuir waves @xmath69 .
thus , we use the propagating langmuir waves with the same frequency @xmath4 .    we generate 3d isotropic langmuir turbulence by using fourier transform description , which is slightly modified from the description for magnetic turbulence developed by giacalone & jokipii ( 1999 ) .
it is described by superpositions of @xmath70 fourier modes , each with random phase , and direction @xmath71 here , @xmath72 , @xmath73 , @xmath74 , and @xmath4 are the amplitude , phase , wave vector , frequency of the @xmath75th mode , respectively .
the amplitude @xmath72 of each mode is defined as @xmath76^{-1},\ ] ] where the variable @xmath1 represents the amplitude of turbulent field .
we use the following form for the power spectrum @xmath77 where @xmath78 is the correlation length of the field . here
, @xmath79 is chosen such that there is an equal spacing in logarithmic @xmath80-space , over the finite interval @xmath81 and @xmath82 , where @xmath83 , @xmath84 and @xmath85 is chosen to be @xmath86 or @xmath87 .
the spectrum has a peak at @xmath88 and the spectral index is for 3-dimensional isotropic langmuir turbulence .
then we define two parameters which characterize radiation spectra .
the first one is @xmath89 where @xmath90 and @xmath91 .
we call @xmath5 `` strength omega '' , and @xmath7 `` spatial omega '' .
the strength omega @xmath5 accounts for the effect of the field strength .
the meaning of strength omega can be understood by considering the curvature of orbit and beaming effects for the relativistic particles as follows .
for @xmath92 , the local curvature radius of the orbit suffering from perpendicular acceleration by the electric field is @xmath93 .
the radiation is concentrated in the beaming cone with an angle @xmath19 , so the searchlight sweeps the observer in the timescale of @xmath94 .
it is an analogy of the cyclotron frequency in the mechanism of synchrotron radiation .
therefore , @xmath95 represents the effect of the field strength on the radiation spectra . as for the spatial omega @xmath7 ,
since the electron moves nearly at the light speed , the rate of change in the force direction for the electron is @xmath96 .
the ratio @xmath97 parametrizes the field feature , which represents the ratio of the deflection angle @xmath43 to @xmath62 .
therefore , the spectra from the turbulent magnetic field changes for @xmath98 to @xmath99 drastically ( kelner et al .
2013 , medvedev et al .
2011 , and teraki & takahara 2011 ) .
note that @xmath14 is well known as the `` strength parameter '' .
it means the change of the lorentz factor of an electron which is accelerated along the electric field on the spatial scale @xmath100 .
we can understand it clearly by using the work by an electric field on an electron as @xmath101 where we assume the electric field is parallel to the velocity .
we should note that it is different from the one we used in teraki & takahara 2011 by a factor of @xmath102 , and it agrees with the definition of fleishman & toptygin ( 2007a , b ) .    the second one is the ratio of the frequency of the langmuir waves to @xmath7 , @xmath103 for @xmath104 , the force direction changes with frequency @xmath4 for all electrons . for @xmath105 ,
a change of the force direction is mainly from the spatial fluctuation . summarizing above ,
although there are three parameters of the turbulent field ( @xmath7 , @xmath5 , @xmath4 ) , we can reduce them to two parameters of ( @xmath106 ) when our interest is in the spectral signature .
we investigate the spectral features in the parameter plane of ( @xmath106 ) .
we inject monoenergetic electrons with lorentz factor @xmath107 into the generated turbulent field . the number of electrons used for each calculation is more than @xmath108 and
is written in the caption of each figure .
the initial velocity is randomly chosen to achieve an isotropic distribution .
next we solve the equation of motion @xmath109 by using the method which has second order accuracy for each time step .
the electrons get energy in the langmuir turbulence from the parallel ( to the velocity ) component .
if we pursue electrons much longer time than the pft of a given frequency , the radiation spectrum in high frequency region corresponds to an `` integrated spectrum '' .
this spectrum would be understood by a superposition of `` instantaneous spectra '' . here , the term of `` instantaneous spectrum '' means a spectrum which does not contain the effect of secular energy change of the radiating electrons .
if we do not obtain the instantaneous spectra , we can not discern either the intrinsic feature of instantaneous spectra or energy change of radiating electrons determines the spectral features .
therefore , in this paper we concentrate on `` instantaneous spectra '' . on the other hand ,
spectral broadening due to finite integration time of particle orbit is inevitable . to compromise the accuracy and the instantaneousness
, we choose the integration time as @xmath110 times pft of the peak or break frequency .
let us take an example of the spectrum for well known jitter radiation , for @xmath111 .
the break frequency of the jitter radiation is @xmath112 , and corresponding pft is @xmath113 .
we can understand the origin of this frequency by using `` the method of virtual quanta '' , ( cf .
jackson 1999 ) .
the spatial fluctuation of the electric field with a spatial scale @xmath114 . this fluctuation can be regarded as a photon which has the frequency with @xmath115 in the electron rest frame .
the electron scatters the photon , and transformed back to the observer frame , the frequency of the scattered photon is written as @xmath116 .
in other words , we can understand this process as an analogy of the inverse compton scattering for the photon with frequency @xmath7 .
we note that this analogy is based on the assumption that the observer is in almost the center of the beaming cone over the pft .
it is valid when @xmath117 .
the deflection angle in the time scale of @xmath118 is @xmath119 , and when we require that it is smaller than @xmath62 , we get the condition @xmath120 .
it can be transformed to @xmath121 .
an alternative explanation can be done by using `` the photon chasing effect '' ( rybiki & lightmani 1979 ) . since the radiation from an ultra - relativistic particle
is concentrated into the beaming cone , the particle chases the photon .
therefore , the observer sees the emitted radiation in a time span @xmath122 , so the frequency is @xmath123 . as a consequence ,
the typical frequency is @xmath112 . as we see in the next section , pft for the peak frequencies of the spectra is the one of the typical timescales of @xmath118 , @xmath124 , and @xmath125 .
we choose a suitable one for each case .
the integration time @xmath110 times the pft of the peak frequency is sufficiently long to resolve the spectral shape in the low frequency regions .
we calculate the radiation spectra using the acceleration , velocity , and position of electrons .
the energy @xmath126 emitted per unit solid angle @xmath127 ( around the direction @xmath128 ) and per unit frequency @xmath129 to the direction @xmath128 is computed as @xmath130 }     { ( 1 - \bm{\beta } \cdot \bm{n } ) ^2 } \exp\bigl\{{i\omega ( t^{\prime } - \frac{\bm{n } \cdot \bm{r}(t^{\prime})}{c})}\bigr\ } \bigr|^2,\ ] ] where @xmath131 is the electron trajectory , @xmath132 is retarded time ( jackson 1999 ) .
since we have not any specific constraints for @xmath14 and @xmath133 which are realized in the high energy astrophysical object , we sweep wide parameter range of @xmath14 and @xmath133 .
first , we show the radiation spectra for @xmath134 ( figure [ jit ] ) , i.e. , for the situation where typical spatial scale is shorter than the inertial length .
we set @xmath135 to @xmath136 , and fix @xmath137 ; specifically we set @xmath138 and @xmath139 , and change @xmath5 from @xmath140 to @xmath136 , and take @xmath141 .
the inequality @xmath105 can be transformed to @xmath142 , which means that the fluctuation scale is much smaller than the inertial length .
these fluctuations would be damped by landau damping and may not be realized in high energy astrophysical objects ( treumann & baumjohan 1997 ) .
however , we can not reject the possibility that @xmath143 is realized for some time span in relativistic plasmas , so that we study the spectra for @xmath144 . to explore the regime for @xmath144
clearly , we set an extreme value @xmath145 .
an example of the temporal variation of the lorentz factor of an electron is depicted in figure [ jit_ene ] .    for @xmath146 ,
the low frequency spectrum is as flat as @xmath60 , and there is a break at @xmath147 , and the spectrum declines with power law @xmath148 in the high frequency region . for @xmath149 , the break frequency becomes higher than that for @xmath146 , and the high frequency spectrum deviates from a power law . for @xmath150 , the spectrum in low frequency region becomes hard with an index @xmath151 , and the spectrum in higher frequencies reveals a cutoff feature . for @xmath152 , we see further different features . the spectrum in low frequency side of the peak becomes softer , with the spectral index of about @xmath18 , and we see a small deviation from a power law in the lowest frequency region .
the features for these spectra can be understood by using the analogy to the radiation theory from the stochastic magnetic field ( medvedev et al .
2011 , teraki & takahara 2011 ) .
since the wavelength of langmuir waves for @xmath105 is very short , the oscillation of the electric field can be neglected in the particle crossing time of the wavelength .
in fact , @xmath105 is also written by @xmath153 .
the crossing time @xmath154 corresponds to the pft of the typical frequency .    as explained above
, we can use the straight analogy to the radiation theory of the stochastic magnetic fields for the radiation spectra from langmuir turbulence for @xmath144 by substituting the electric field strength @xmath155 with magnetic field strength @xmath32 . for @xmath156 and @xmath144
, we can use the radiation theories of the dsr and jitter radiation theory ( fleishman 2006 , medvedev 2006 ) . on the contrary , for @xmath157 and @xmath144 jitter radiation of strong deflection regime can be applicable ( medvedev et al .
2011 , teraki & takahara 2011 ) .
we call @xmath156 and @xmath144 regime as `` jitter radiation '' regime because the jitter radiation is basically perturbative theory for @xmath156 .
we call the radiation for @xmath157 and @xmath144 regime as `` wiggler radiation in langmuir turbulence '' , or `` wrl '' in short .
although the wiggler radiation is not the radiation mechanism from the stochastic field but that for a fixed field , it has a common picture that the observer is in and off the beaming cone in the course of time .
first , we describe the signature of the radiation spectra of jitter radiation or diffusive synchrotron radiation .
the radiation signatures are determined by acceleration perpendicular to the motion , and the observer is always in the beaming cone , which is the same situation as the langmuir turbulence for @xmath98 and @xmath144 . for @xmath117 ,
the spectrum is written by a broken power law @xmath60 in the low frequency region and @xmath158 in the high frequency region .
the break frequency is @xmath159 . for @xmath160 , while the break frequency remains the same as @xmath161 , the multiple deflection effect comes into play in the spectral features near the break frequency .
the spectrum in the low frequency region becomes @xmath61 .
since the multiple deflection makes the angle between observer direction and velocity larger than the beaming cone angle @xmath62 , the observer sees the radiation over the timescale which is determined by the deflection condition .
fleishman supposed that the angle changes diffusively .
the break frequency of @xmath162 is calculated from the angle diffusion ( fleishman 2006 ) .
the spectral index of @xmath63 comes from the diffusivity , too .
next , we describe the signatures for strong deflection regime of @xmath163 . for magnetic turbulence ,
the spectral shape resembles synchrotron radiation in the middle frequency region and deviations from it would be seen in low frequency region and highest frequency region ( teraki & takahra 2011 ) .
the peak frequency of @xmath164 is from the sweeping of the beaming cone on the observer .
the picture is the same as the synchrotron radiation . for short wavelength regime @xmath165 of langmuir turbulence ,
the physical picture is the same as this , because the spatial fluctuation dominates the changing rate of the deflection angle .
a single deflection angle is @xmath166 , which is larger than @xmath62 for @xmath167 . as a result
, the beaming cone sweeps out of the observer within one deflection .
the intensity of the radiation off the beaming cone is weak .
therefore , the timescale which sweeps the observer @xmath168 corresponds to the pft of the typical frequency . considering the photon chasing effect
, we get the peak frequency @xmath169 .
the spectral break in the low frequency region would correspond to the deviation from local circular orbit , but the numerical error from finite integration time is also becoming large in the lower frequency region .
we do not discuss this point further here , since it is not the main point of contents in this paper .
the power law component in the highest frequency region comes from the smaller scale part of turbulence .
it is the same as the spectra of jitter regime .
we note here that the power law component in high frequency region can not be seen for @xmath150 and @xmath152 .
the reason may be as follows .
in contrast to the magnetic turbulence , the energy of radiating electrons changes for @xmath99 .
therefore , the high frequency region is determined by only the later part of integration time , because the peak frequency is @xmath170 for @xmath157 , and the electrons get energy and radiate higher energy photon in later time .
the power law component for the highest frequency region in our calculation is hidden by the component that small numbers of electrons with larger energy radiate by wiggler mechanism .
finally , we consider effects of the energy change in wrl regime , which is the one of the different points from the magnetic case . based on an example of the change of lorentz factor is depicted in figure [ jit_ene ] , we estimate the energy change in pft @xmath171 for peak frequency . the pft is @xmath124 .
the change of lorentz factor is estimated as @xmath172 therefore , the change of the lorentz factor in this timescale is smaller than @xmath17 .
thus , we do not need to consider the energy change in pft for peak frequency .
however , for sufficiently large @xmath14 , we have to take it into account for lower frequency spectrum .
we calculate the spectra for modest @xmath14 in this paper , so we omit this problem .
it will be studied in future works .
next we show the radiation spectra for @xmath156 and @xmath13 ( figure [ drl ] ) , i.e. , a situation where long wavelength turbulence with weak amplitudes dominates ; the spectra for @xmath173 and @xmath17 are also depicted for comparison .
an example of the temporal variation of the lorentz factor of an electron is depicted in figure [ drl_ene ] .
firstly , we consider the meaning of the parameters of @xmath156 and @xmath13 , which correspond to @xmath174 . in this regime ,
the changes of the direction of acceleration are mainly due to wave oscillation , rather than the spatial fluctuations because the crossing time is longer than the oscillation time .
moreover , @xmath175 is derived from @xmath176 .
therefore , we can regard the orbit as straight in the time scale of plasma oscillation @xmath177 .
the condition @xmath13 can be transformed to @xmath178 , which means that the inertial length is shorter than the wavelength .
therefore , this regime is likely to occur in the astrophysical objects .
we set @xmath138 , @xmath179 and change @xmath4 from @xmath140 to @xmath47 , so that @xmath180 and @xmath173 to @xmath47 , and take @xmath141 .
as was discussed in the previous subsection for @xmath173 , the spectrum shows jitter radiation signature .
the peak frequency is @xmath161 and @xmath181 in the low frequency region , and @xmath148 above the peak frequency reflecting the spectrum of the turbulent electric field @xmath182 .
as @xmath133 becomes larger , the peak shifts to higher frequency . for @xmath183 ,
the peak frequency is @xmath184 , which is identified with @xmath185 .
the spectral index of low frequency side is @xmath186 .
this feature coincides with the result of drl theory ( fleishman 2006 ) .
we call this regime `` drl regime '' .
the drl theory predicts the spectral shape for @xmath156 and @xmath13 as follows .
the peak frequency is @xmath185 with an abrupt cutoff in the higher frequency side and a power law component emerges in the highest frequency region for @xmath187 . in the low frequency side
, the spectrum becomes @xmath12 .
the peak frequency is determined by the timescale of plasma oscillation .
the shortest timescale of the electron motion is @xmath177 , and the observer located along the initial velocity direction can see this radiation , because the orbit is regarded as straight in this time scale as we showed above .
we consider doppler boosting , and we get the peak frequency @xmath188 in the observer frame . the origin of the power law component in the highest frequency region is the same as jitter radiation . the hard spectral index @xmath17 in the frequency region lower than the peak is from the photon chasing effect .
the electrons oscillate by the same oscillating frequency @xmath4 and move different angles to the observer in the observer frame .
the photon chasing effect shifts the photon frequency from each electron and makes the @xmath12 spectrum ( see fleishman 2007a , b ) .
it is regarded as the emission spectrum from an electron integrated over the solid angle , which can be understood in an analogy to the undulator theory ( jackson 1999 ) , although the force direction changes in this case not spatially but temporally .
when the particle mean velocity and the wavenumber are fixed , that makes no difference for the radiation spectra .
the peaky shape of the spectra is the most remarkable feature of the spectra for drl case .
a difference is that the wave number is a vector , while the frequency is a scalar . for dsr ,
the electron feels spatial fluctuation with wavenumber along the velocity @xmath189 , therefore the low frequency spectrum becomes flat . on the other hand , for drl
, all electrons feel the same frequency of @xmath4 for the langmuir turbulence ( c.f .
fleishman & toptygin 2007a , b and medvedev 2006 ) .
the remaining interesting regime of @xmath157 and @xmath13 is an open issue . as we explained in section [ int ] ,
the radiation spectrum for langmuir turbulence with @xmath15 predicted by fleishman & toptygin ( 2007b ) is @xmath61 below the peak of @xmath190 , which is inconsistent with the statements by getamantsev and tokarev ( 1972 ) .
here we investigate this regime and clarify the radiation spectra .
the inequations of @xmath157 and @xmath13 mean that the wavelength is longer than the inertial length and that the typical scale of pft is @xmath191 .
we set @xmath192 , @xmath193 to @xmath194 , so that @xmath195 to @xmath196 , and @xmath138 .
we set @xmath197 , and we change @xmath4 from @xmath136 to @xmath194 , and @xmath85 is chosen to be @xmath87 here .
we show the interesting results for @xmath157 and @xmath13 ( figure [ wig ] ) .
an example of energy change for @xmath198 and @xmath193 is depicted in figure [ wig_ene ] .
for clear discussion of the radiation spectra in this regime , we discuss the radiation spectra for @xmath199 and for @xmath200 separately . for @xmath200 ,
the peak frequency is @xmath201 .
the spectral index of low frequency region is @xmath17 , and cutoff feature can be seen above the peak .
this region is regarded as the drl regime from these signatures .
for @xmath200 , i.e. , @xmath202 , the particle is not deflected by large angle because the direction of the electric field changes in a time shorter than the time for which the beaming cone sweeps the observer .
we set @xmath203 , therefore , a power law component in high frequency region is not seen at all .    on the contrary ,
different features emerge for @xmath199 . the peak frequency becomes larger than @xmath11 .
moreover , the spectra below the peak frequency become softer , the index changes from @xmath17 to @xmath18 .
the energy change of electrons may cause the change of the peak frequency , but it can not explain the soft spectrum .
rather , it would be naturally understood that the peak frequency is @xmath164 and @xmath204 by using the analogy of the wiggler radiation .
we consider that we should use wrl theory not only for @xmath205 , but also for @xmath15 , because the deflection angle in one deflection is also larger than @xmath62 for this case .
this is in contrast to the drl theory , which predicts the same spectral signatures for the parameter range of @xmath15 as for @xmath206 .
thus , our numerical calculations reveal new features which have not been known previously for langmuir turbulence .
we ascribe that the parameter regime @xmath15 is in the wrl regime in @xmath10 plane .
this radiation feature seems to be consistent with the statements by getmantsev & tokarev ( 1972 ) . to clarify the spectral features for @xmath15 in more detail , and to confirm their statement and our consideration , we examine the radiation from a relativistic electron moving in pure plasma oscillation in the next section .
in this section , we investigate the emission of a relativistic electron suffering from pure plasma oscillation in order to discuss the interpretation of the features of the radiation spectra for @xmath15 .
to clarify the origin of the peak frequency @xmath207 , we set a simple configuration of the electric field , where electron motion is deterministic compared to stochastic character in turbulent fields .
we calculate the electron velocity analytically and the radiation spectra numerically . by comparing the motion and spectra
, we interpret the mechanism which determines the peak frequency .
lastly we consider the radiation spectra from the turbulent field by using these results .
we use a single langmuir wave which has infinitely large wavelength @xmath208 , in other words , @xmath209 .
therefore , it is a pure plasma oscillation .
we set @xmath210 with @xmath211 we characterize the field by using a single parameter of @xmath212 we inject an electron along the @xmath213-axis at @xmath214 with the initial lorentz factor @xmath215 and solve the equation of motion ( [ eom ] ) .
therefore , the orbit is determined by @xmath215 and @xmath216 . solving the equation of motion
, we get the momentum @xmath217 since the velocity is a periodic function , we can define the mean velocity by @xmath218 .
the mean velocity @xmath219 can not be represented elementarily in a general form .
then , we take the parameter @xmath220 and approximate the motion hereafter .
we expand the lorentz factor and the velocity , and get the mean velocity and the mean lorentz factor in the lowest order of @xmath221 .
@xmath222 for @xmath223 , @xmath224 , while for @xmath225 , @xmath226 .
we note that @xmath227 can be much smaller than @xmath215 for @xmath225 . using this approximated velocity ,
we calculate the maximum lorentz factor in the mean velocity frame , to clarify the fact that the radiation signatures depend on @xmath216 as @xmath228.\ ] ] for @xmath223 , @xmath229 the motion in this frame is non - relativistic , therefore , the radiation in this frame is dipole radiation . on the other hand , for @xmath225 ,
the maximum lorentz factor is @xmath230 therefore , the motion is relativistic even in this frame and the radiation spectrum consists of higher harmonics , because @xmath231 approaches @xmath17 .
it should be noted that for @xmath232 , the motion in the mean velocity frame is mildly relativistic with lorentz factor @xmath233 , and @xmath234 .
we can see that the transition from non - relativistic harmonic motion to relativistic motion occurs around @xmath235 a few from this fact .
next we show numerically calculated radiation spectra from the electron and their features are interpreted in terms of the properties of the orbit .
we fix @xmath236 , and change @xmath4 to change the parameter @xmath216 . the observer is on the @xmath213-direction .
we calculate radiation spectra using much longer integration time than the pft , because the electron moves perfectly periodically . as a consequence ,
the spectra show very sharp features , which makes it easier to understand the relation between spectral features and orbit .
first , we show the spectrum for @xmath237 ( @xmath238 , figure [ peaks](a ) ) . we see a sharp peak like a delta function at the frequency @xmath239 .
this is understood in terms of the motion of the electron in the mean velocity frame . for @xmath223 , in the mean velocity frame , electron motion is a non - relativistic simple harmonic motion with the frequency @xmath240 .
therefore , the radiation is the dipole radiation with the frequency of @xmath241 . since @xmath242 for @xmath223 , the radiation frequency in the observer frame is @xmath243
thus , we ascribe the frequency @xmath243 in the radiation spectra of perturbative regime ( @xmath202 ) to the doppler shifted dipole radiation .
next we show the spectrum for @xmath244 ( @xmath245 , figure [ peaks](b ) ) .
we can see the higher harmonics of @xmath246 .
it is from effects of retarded time , as is clearly seen in the mean velocity frame .
however , the effect is very weak , as the power of the second harmonics is about @xmath247 times smaller than the fundamental mode in the frequency resolution of this calculation .
the ratio of the power of the second harmonics to the fundamental mode is proportional to @xmath248 , so the second harmonics is much smaller than the fundamental mode in this case .
for @xmath249 , the spectral shape changes significantly .
first , the higher harmonics stand more strongly , because @xmath231 approaches @xmath17 .
many harmonics are as strong as the fundamental mode for @xmath250 ( figure [ peaks](c ) ) , and the envelope of the peaks of the harmonics shows an exponential cutoff .
we note that the spectrum in the frequency region higher than @xmath251 comes from numerical error .
second , the frequency of the fundamental mode deviates from @xmath243 , because the difference between @xmath227 and @xmath215 becomes larger .
the mean lorentz factor @xmath227 is @xmath252 for @xmath232 , thus the frequency of the fundamental mode in the original frame is @xmath253 .
the difference between @xmath254 and @xmath255 is small , but we can discern it in figure [ peaks](c ) .
next we discuss the peak frequency ( cutoff frequency ) for @xmath256 .
we show the spectra for @xmath257 , and @xmath196 in figure [ nc_peak ] .
the fundamental frequency is @xmath254 for @xmath232 , @xmath258 for @xmath259 , and @xmath260 for @xmath261 , but the cutoff frequency around a few hundreds does not change .
we see the cutoff frequency is not of the fundamental mode , but it is determined by the higher harmonics for @xmath256 .
the radiation spectra in the observer frame also can be derived by regarding it as a doppler boosted emission . however , since the mechanism which determines the peak frequency is the same as the wiggler radiation we can understand the peak frequency more easily by considering the pft in the observer frame .
the condition @xmath262 is equivalent to that pft of the typical frequency in wiggler mechanism , where @xmath124 is shorter than the oscillating time @xmath125 . on the other hand , the lorentz factor relevant for the peak frequency is not @xmath227 , but @xmath215 , because the beaming cone sweeps the observer in the phase around @xmath263 , where @xmath75 is a natural number .
we note that the change of the lorentz factor in the time scale of @xmath124 is @xmath17 at most , as seen in equation ( [ lfchange ] ) . as a result ,
the cutoff frequency is @xmath264 . in this way
, we get a clear understanding of the mechanism of the peak frequency shift around @xmath265 .    here
, we compare the results in this section with the radiation spectra obtained in the preceding section .
the case of langmuir turbulence with @xmath15 ( @xmath266 ) corresponds to the case of pure plasma oscillation with @xmath256 , since @xmath267 and @xmath209 for pure plasma oscillation .
moreover , the approximation we used in pure plasma oscillation of @xmath268 is also applicable for the case of langmuir turbulence , since @xmath269 and @xmath107 for the spectra shown in figure [ wig ] .
for pure plasma oscillation with @xmath256 , we have shown that the peak frequency is @xmath270 , and it consists of the higher harmonics of @xmath271 .
therefore , the peak frequency for the langmuir turbulence with @xmath272 in figure [ wig ] is interpreted as @xmath164 , and it naturally explains the fact that the peak frequency is larger than @xmath11 .
lastly we consider the spectral index for the langmuir turbulence with @xmath15 . as we see above
, the spectral index for @xmath15 is neither @xmath17 predicted by perturbative drl nor @xmath63 predicted by the angle diffusion effect .
we regard that the spectral index is around @xmath18 , because the radiation mechanism is identified as wiggler mechanism , and the angle integrated spectral index is @xmath18 in wiggler mechanism .
the energy change in pft for the frequency lower than peak can be larger than @xmath17 for @xmath15 .
this is a different point from the wiggler , and it can affect forming the spectral index . however , we choose modest value of @xmath273 for langmuir turbulence in the present paper , so that this effect does not stand out . summarizing above , we confirm that the peak frequency is @xmath274 and the spectral index in low frequency side of the peak @xmath275 for the langmuir turbulence with @xmath15 and @xmath276 .
the mechanism which determine these radiation features is wiggler mechanism , which is consistent with the statements by getamantsev and tokarev ( 1972 ) .
lastly we show an example of the extreme case of @xmath277 . for @xmath278 ,
the motion becomes strongly nonlinear and can not be treated analytically .
thus , we show numerically calculated electron orbit . we show the radiation spectra and the orbit for @xmath279 ( figure [ peaks](d ) , figure [ tr500 ] ) . as we expected , the peak frequency is @xmath264 , because the observer sees this electron mainly in the phase that the electric field and the electron velocity is nearly perpendicular .
we note that the spectral index of @xmath280 is the same as the wiggler radiation when the observers located in a particular direction , i.e. , the spectrum is not the angle integrated spectrum .
this spectral index is an collateral evidence for the spectral index of @xmath18 for the case of turbulence with @xmath15 for which angle integrated spectrum is calculated .
the lorentz factor changes from @xmath47 to @xmath281 , but @xmath215 is relevant to the observer situated at @xmath213-axis .
in other words , the beaming cone sweeps the observer when @xmath282 in the present geometry .
the lorentz factor relevant to an observer oriented in different angle is significantly different to each other .
moreover , in general the energy change in the pft of peak frequency becomes larger than @xmath283 , because the electric force in some phase of oscillation accelerates the electron linearly , and the curvature radius becomes larger .
we have to consider the linear acceleration emission in this case .
thus , we draw the line on @xmath284 and @xmath285 , and divide the @xmath15 region .
we call the radiation for this parameter range as `` non - linear trajectory '' radiation .
this part of the spectra from an electron moving in the 3d turbulent electric field is determined by the chaotic trajectory .
the generalization of the features of this regime is a future work .    summarizing this section ,
we have considered the motion and radiation in a single mode plasma oscillation .
we clarify that the cutoff frequency for @xmath256 ( @xmath286 ) is @xmath270 , which consists of higher harmonics of the fundamental frequency of @xmath271 .
it is from the effect that the beaming cone sweeps the observer , in the same way as the wiggler radiation . using this result
, we interpret that the peak frequency for 3d langmuir turbulence for @xmath15 ( @xmath266 ) is @xmath164 , where @xmath23 is determined by the acceleration . the shallower spectrum for @xmath15
can be explained by wrl mechanism .
lastly , we show numerically calculated radiation spectra for the extreme case of @xmath277 .
it shows wiggler like spectra in the middle and high frequency region , while the flattening can be seen in the low frequency region .
it is from the effect of elongated trajectory to the electric field direction .
the radiation signatures are summarized as a chart in the @xmath10 plane in figure [ a - b_plane ] .
we have calculated the radiation spectra from relativistic electrons moving in a langmuir turbulence by using first principle numerical calculation .
we characterize the radiation spectra by two parameters .
the one is @xmath287 , where @xmath3 is the strength omega , and @xmath288 is the spatial omega . the strength omega accounts for the effect of the field strength to the radiation spectra , and the spatial omega accounts for the effect of spatial fluctuation with a typical scale of @xmath0 .
the other is @xmath289 , where @xmath4 is the plasma frequency , which accounts for effects of the time variability of the waves .
we investigate the spectral signatures in the @xmath10 plane ( figure [ a - b_plane ] ) . for @xmath156 and @xmath144 ,
the spectral features are the same as those of jitter radiation or diffusive synchrotron radiation . for @xmath206 and @xmath290 ,
the theory of the diffusive radiation in langmuir turbulence is confirmed , where time variability plays a primary role . for @xmath15 , the spectra show novel features which are not predicted by drl theory . in this regime ,
the peak frequency is @xmath291 , which is higher than the predicted frequency @xmath11 from the drl theory .
the spectral index of the frequency region lower than the peak is @xmath275 .
these features are explained by the wiggler mechanism . to clarify the radiation features in this regime ,
we calculate the radiation spectra from an electron moving in an oscillating electric field , i.e. , for vanishing spatial omega .
we analytically calculate the motion of the electron , and numerically calculate the radiation spectra form this electron .
we show that for @xmath292 , the spectrum around the peak frequency consists of the higher harmonics of the fundamental mode , by considering the radiation in the mean velocity frame .
the electron motion becomes relativistic for @xmath256 even in this frame , so that strong higher harmonics photons are emitted because of the retarded time effect . as a result ,
the spectra in the observer frame consists of the higher harmonics of @xmath64 .
the peak frequency is characterized by @xmath207 , which is understood by the analogy of the wiggler radiation .
the feature that the radiation spectra from langmuir turbulence have a wide range of spectral indices can be important for high energy astrophysical objects , in particular gamma ray bursts .
the emission mechanism of grb is not settled for now .
the spectral indices of low frequency side of the band function are distributed as a gaussian with the central value of @xmath293 .
non - negligible number of grbs have spectral index harder than the theoretical limit for synchrotron radiation @xmath18 .
the photospheric emission model can overcome this difficulty , but it also has a difficulty that intrinsic photosheric emission is too hard with a low energy spectral index of @xmath294 and it is difficult to make it softer . on the other hand , the radiation mechanism from langmuir turbulence in this paper has some advantages . not only the spectral index is harder than the synchrotron radiation and it can reproduce very hard spectra of observed grb ( fleishman 2007b ) , but also it may explain a wide range of spectral indices . because the parameters of @xmath14 and @xmath133 are likely to have a value around 1 near the shock front ( silva 2006 , dieckmann 2005 ) , the radiation spectra change drastically around these parameters .
we are grateful to the referee for his / her constructive and helpful comments and for informing us of the interesting paper by getmantsev & tokarev ( 1972 ) .
we also thank k. toma and s. kimura for useful discussions .
this work is supported by the jsps research fellowships for young scientists ( y.t .
, 24593 ) .    99 akhiezer , a. i. , & shulga , n. f. 1987 , sov .
phys.usp . , 30 , 197 bret , a. , dieckmann , m. e. , & deutch , c. 2006 , phpl , 13 , 082109 dieckmann , m. e. 2005 , phys .
94 , 155001 fleishman , g. d. 2006 , apj , 638 , 348 fleishman , g. d. , & toptygin , i. n. 2007a , phys .
e. , 76 , 017401 fleishman , g. d. , & toptygin , i. n. 2007b , mnras , 381 , 1473 fleishman , g.d . , & urtiev , f. a. 2010 , mnras , 406 , 644 giacalone , j. , & jokipii , j. r. 1999 , apj , 520 , 204 getmantsev , g. g. , & tokarev , yu . v. 1972 , ap&ss , 18 , 135 jackson , j. d. 1999 , classical electrodynamics ( 3rd ed . ;
new york : wiley ) kelner , s. r. , aharonian , f. a. , & khangulyan , d. 2013 , apj , 774 , 61 landau , l. d. , & lifshitz , e. m. 1980 , the classical theory of fields ( 4th ed.;butterworth - heinemann ) mao , j. , & wang , j. 2013 , apj , 776 , 17 medvedev , m. v. 2000 , apj , 540 , 714 medvedev , m. v. 2006 , apj , 637 , 869 medvedev , m.v . ,
frederiksen , j. t. , haugbolle , t. , & nordlund , a. 2011 , apj , 737 , 55 ostriker , j. p. & gunn , j. e. 1969 , apj , 157 , 1395 rees , m. j. 1971 , iaus , 46 , 407 reville , b. , & kirk , j. g. 2010 , apj , 724 , 1283 reynolds , s. j. , & medvedev m. v. 2012 , phpl , 19 , 023106 rybicki , r. d. , & lightman , j. g. 1979 , radiative processes in astrophysics ( new york : wiley ) silva , l. o. aipc , 856 , 109 treumann , r. a. , baumjohann , w. 1997 , advanced space plasma physics ( london : imperial college press ) teraki , y. , & takahara , f. 2011 , apj , 735 , l44 teraki , y. , & takahara , f. 2013 , apj , 763 , 131 toptygin , i. n. & fleishman , g. d. 1987 , ap&ss , 132 , 213         and @xmath296 , and @xmath47 from top to down in low frequency range .
the number of electrons used for these calculations is @xmath297 .
the straight lines show a power law spectrum with index @xmath17 and @xmath298 .
we see the transition from jitter radiation to drl as @xmath133 increases .
, width=453 ]       with @xmath193 , @xmath300 , @xmath301 , and @xmath194 from top to down .
the number of electrons used in these calculations is @xmath194 .
the straight lines show the power law spectra with index @xmath18 and @xmath17 .
we see the transition from the drl to wrl as @xmath133 decreases .
, width=453 ]        , and @xmath196 ( @xmath303 , and @xmath304 ) from top to down .
a factor of @xmath305 is multiplied to the spectrum for @xmath232 , and @xmath306 is multiplied for @xmath261 .
the fundamental frequency is @xmath307 .
it is @xmath254 for @xmath232 , @xmath258 for @xmath259 and @xmath260 for @xmath261 .
the cutoff frequency is a few times of @xmath308 , which does not depend on @xmath309 explicitly .
, width=453 ]      , and vertical axis is @xmath311 . for the jitter regime with @xmath156 and @xmath144 ,
the radiation spectra are determined by the spacial fluctuations , because @xmath7 is the largest of the three .
the typical frequency for this case is @xmath161 . for @xmath206 ,
i.e. , @xmath312 , the radiation spectra are represented by drl theory , and typical frequency is @xmath313 .
the line @xmath314 divides the drl region and wrl region , and the spectral features for the wrl regime @xmath15 are newly clarified in this paper .
the typical frequency is @xmath164 and the spectral index of frequency region lower than the peak is @xmath275 , in the same way as synchrotron radiation . for @xmath315 and @xmath316 ,
the orbit of a radiating electron depicts non - linear trajectory , and its signature appears at the low frequency region of the spectrum .
, width=340 ]