Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: electron scattering is one of the most important processes in determining the emergent thermal spectrum from models of optically thick accretion discs around black holes and neutron stars . electron scattering opacity generally dominates absorption opacity in the atmospheres of the innermost regions of such discs @xcite . in the case where coherent ( thomson ) scattering is a good approximation , the resulting local thermal spectrum of some annulus in the disc is generally harder than a blackbody with the same effective temperature , due to incomplete thermalization at the scattering photosphere . however , incoherent ( compton ) scattering in the disc surface layers by thermal electrons can reduce this spectral hardening by increasing the energy exchange between the photons and the plasma @xcite . in addition to the disc atmosphere itself , many models invoke a powerful corona above the disc consisting of high temperature or nonthermal electrons that compton upscatter disc photons to produce the energetically significant hard x - rays that exist in certain classes of active galactic nuclei and in certain black hole x - ray binary accretion states ( e.g. @xcite ) . these hard x - rays in turn interact with the relatively cool disc atmosphere to produce reflection spectra that are widely observed in many black hole sources ( e.g. @xcite ) .    in this work we explore turbulent comptonization , which is the effect of comptonization by bulk , turbulent electron motions on photon spectra @xcite . in sources with radiation pressure dominated accretion flows , bulk velocities may exceed thermal electron velocities , a phenomenon first pointed out in @xcite . the turbulent speeds @xmath0 on the outer scale of an mhd turbulent cascade will be of order the alfvn speed , and the ratio of this to the root mean square electron thermal velocity is therefore @xmath1 here @xmath2 , @xmath3 , and @xmath4 are the gas , radiation , and magnetic pressures , respectively , and @xmath5 is the ratio of the proton to electron mass . stratified shearing box simulations of magnetorotational turbulence generally have disc atmospheres that are supported by magnetic fields rather than thermal pressure @xcite . hence the first factor generally exceeds unity in an otherwise radiation pressure dominated disc . bulk speeds on the outer scale of the turbulence will therefore exceed the electron thermal speeds whenever the radiation pressure to gas pressure ratio exceeds the ratio of the proton to electron mass ratio , and even smaller depending on how magnetically supported is the disc atmosphere . in this regime , bulk comptonization by the turbulence may dominate thermal comptonization in determining the shape of the spectrum emitted by a local patch of the disc . this regime is commonly reached for near eddington accretion on black holes of all mass scales . indeed , the inner disc solution of the standard geometrically thin model of @xcite gives vertically averaged radiation to gas pressure ratios of approximately @xmath6 where @xmath7 is the ratio of turbulent stress to thermal pressure , @xmath8 is the black hole mass , @xmath9 is the radiative efficiency of the disc as a whole , @xmath10 is the luminosity in eddington units , and @xmath11 is the radius in the disc in units of the gravitational radius @xmath12 . hence the radiation to gas pressure ratio of the innermost disc will generally exceed the proton to electron mass ratio for near - eddington accretion even for stellar mass black holes , and certainly for supermassive black holes . on the other hand , energy exchange between the photons and the plasma is generally dominated by true absorption opacity in standard disc models for the most supermassive black holes @xcite , which may reduce bulk comptonization by turbulence in these sources . we are also motivated by observational evidence in many sources for a cooler , more optically thick medium that also compton upscatters disc photons in high luminosity sources , but not to the high energies usually associated with the more traditional external corona . the soft x - ray excess of the type 1 seyfert ngc  5548 was fit by @xcite with a warm ( @xmath13  kev ) optically thick ( @xmath14 ) comptonizing medium . similar parameters were obtained for the soft x - ray excess of the type 1 seyfert mrk 509 by @xcite : @xmath15  kev and scattering optical depth @xmath16 to @xmath17 . @xcite fit the soft x - ray excess of the extreme narrow line seyfert 1 re  j0134 + 396 assuming extra comptonization within the disc atmosphere itself with @xmath18  kev and @xmath19 . such a spectral decomposition is consistent with the lack of short time scale variability in the soft x - ray excess compared to the hard x - rays in this source @xcite . @xcite fit the soft x - ray excesses of eleven bright quasars with comptonization in a medium with @xmath20  kev and @xmath21 . @xcite fit soft intermediate state spectra of the black hole x - ray binary gx 339 - 4 with a comptonizing medium with @xmath22  kev and @xmath23 . thermal comptonization fits to the soft x - ray spectra of a number of ultraluminous x - ray sources have @xmath24  kev and @xmath25 @xcite . such fits have been nicely confirmed with _ nustar _ and _ xmm - newton _ spectra of the ultraluminous x - ray sources ngc 1313 x-1 and x-2 , which require warm @xmath26  kev , optically thick @xmath27 comptonizing media @xcite . @xcite and @xcite suggest that one way of forming these optically thick cool coronae is by a radiation pressure driven outflow . the alternative we consider here is bulk comptonization by turbulence in the disc atmosphere itself . comptonization by bulk motions in the accretion flow has also been considered by others . @xcite considered bulk comptonization in converging flows and shocks . starting from this seminal work , bulk comptonization by radial flows has been calculated in detail by numerous authors @xcite . @xcite included bulk comptonization in their monte carlo calculations of photon spectra from radiation mhd simulations of super - eddington accretion flows , and found that it produced significant spectral hardening which resembled spectra of ultra - luminous x - ray sources . here we focus on smaller scale bulk comptonization by turbulence within the disc atmosphere itself . turbulent comptonization has also been invoked in other areas of astrophysics . @xcite and @xcite used then current limits on cosmic microwave background temperature anisotropies to constrain possible turbulent energy on cosmological scales prior to recombination . @xcite considered comptonization by alfvnic turbulence in a relativistic outflow as a model for the spectrum of gamma - ray bursts .    the purpose of this paper is to determine how photon spectra produced by turbulent comptonization depend on properties of the turbulence itself , and how to resolve and interpret this effect in radiation mhd simulations . the structure of this paper is as follows . in section [ sec_general ] we show that the macroscopic physical origins of turbulent comptonization energy exchange are work due to radiation pressure and viscous dissipation due to the radiation viscous stress tensor , and we discuss why this requires us to treat divergenceless turbulence separately from turbulence with non - zero divergence . in section [ sec_energy_exchange ] we discuss the consequences of this for correctly implementing radiative transport in simulations , and derive the appropriate radiation energy equation in both lab and fluid frame variables . in section [ sec_incompressible ] we address the conjecture of @xcite that turbulent comptonization can be treated as thermal comptonization by solving the kompaneets equation with an equivalent  wave " temperature critically dependent on the photon mean free path . we show this is true only for divergenceless turbulence , derive the exact wave temperature with an analytic solution of the radiative transfer equation , and use this result to discuss how the wave temperature depends on the power spectrum of the turbulence . to provide physical insight , we also perform an intuitive , heuristic calculation of the wave temperature which well approximates the analytic solution . in section [ sec_compressible ] we consider bulk comptonization by turbulence with non - zero divergence . we show that comptonization by turbulence whose wavelengths are short relative to the photon mean free path can be treated as thermal comptonization with an equivalent temperature given by the full turbulent power . in the limit of extremely optically thick turbulence , we show how the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids . in section [ sec_discussion ] we discuss our results , and we summarize our findings in section [ sec_conclusion ] . in order to determine how photon spectra produced by turbulent comptonization depend on properties of the turbulence itself , it is useful to understand the distinct macroscopic physical origins of net energy exchange due to this effect . we show that these are simply the work done by radiation pressure and radiation viscous dissipation , and discuss the major consequences of this . we limit our consideration in this paper to non - relativistic velocities .    before proceeding , we define terms and quantities that will be used repeatedly . by bulk comptonization we mean comptonization by bulk motions in general , and by turbulent comptonization we mean bulk comptonization specifically due to turbulence . we denote the characteristic photon mean free path @xmath28 , where @xmath29 is the electron density and @xmath30 is the thompson cross section . we denote the typical length scale for bulk velocity variations @xmath31 , such as the wavelength if there is a well - defined spatial period . unless otherwise stated , by the terms optically thin and thick we mean @xmath32 and @xmath33 , where @xmath34 and @xmath35 are the minimum and maximum length scales in the turbulent cascade , respectively , not referring to the optical depth that a photon would need to travel to escape the medium . net energy exchange due to bulk comptonization is simply the net energy exchange between gas mechanical energy and radiation . inside the photosphere , the mechanical energy loss rate per unit volume to radiation is @xmath36 where @xmath37 is the lab frame radiation pressure tensor . this can also be written as @xmath38 where @xmath39 is the trace of the radiation pressure tensor , @xmath40 is the radiation viscous shear stress tensor , and @xmath41 is the velocity shear tensor . we see that the energy exchange is separated into two pieces , one due to only a diverging velocity field and another due to a shearing velocity field in the presence of a radiation viscous shear stress tensor . the first piece has contributions from two effects , ordinary work done by radiation pressure , and radiation viscous dissipation . the former effect is first order in velocity since it is due to the contribution to @xmath42 that is zeroth order in velocity . energy exchange due to viscous effects , on the other hand , is second order in the velocity , as @xmath43 and the relevant contribution to @xmath42 must themselves be at least first order in velocity since they are a consequence of the velocity field . @xcite conjectured that turbulent comptonization can be treated as thermal comptonization by solving the kompaneets equation with an equivalent  wave " temperature critically dependent on the photon mean free path . this is physically intuitive for divergenceless turbulence since in this case energy exchange is entirely due to radiation viscous dissipation and is therefore second order in velocity . in section [ sec_incompressible ] we prove this conjecture for an arbitrary , divergenceless velocity field of uniform density , derive the exact expression for the wave temperature , and show how it is related to the standard expression for radiation viscosity in the optically thick limit . the connection between bulk comptonization and the radiation viscous shear stress tensor is useful because it correctly suggests that the wave temperature can be written in terms of @xmath43 . since pressure work , on the other hand , is an effect first order in velocity , and since comptonization by a velocity field with non - zero divergence is a combination of pressure work and radiation viscous dissipation , it is not surprising that in this case bulk comptonization can not be treated as thermal comptonization , as we show in section [ sec_compressible ] . but in the optically thin limit , i.e. when the mean free path is significantly larger than the largest length scale in the turbulence , energy exchange that is first order in velocity vanishes since photons are equally likely to downscatter as they are to upscatter . bulk comptonization by a velocity field with non - zero divergence is then solely due to radiation viscous dissipation , and in section [ sec_compressible ] we show that it may be treated as thermal comptonization by solving the kompaneets equation with a wave temperature equivalent to the full bulk power .    in the optically thick case , i.e. when the mean free path is significantly smaller than the smallest scale in the turbulence , the lowest order energy exchange is the work done by radiation pressure to compress the gas , since it is first order in velocity and since radiation viscous effects are suppressed by a factor @xmath44 ( section [ sec_incompressible ] ) . @xcite assumed that effects first order in velocity always vanish on average for turbulent eddies , but in the optically thick limit photons trapped in a converging ( diverging ) region undergo systematic upscattering ( downscattering ) . in the extremely optically thick limit in which the photon and gas fluids are strongly coupled , velocity convergence corresponds to compression in which gas mechanical energy is transferred locally to the photons . in section [ sec_compressible ] we show that in this process a locally thermal photon distribution remains thermal and only changes temperature , completely analogous to the evolution of the cosmic microwave background radiation under the expansion of the universe . unlike energy exchange due to viscous dissipation , this process is reversible . the effect of this process on the emergent spectrum of the disc will depend primarily on how effectively photons are able to escape from such regions to the observer . self - consistent radiation mhd simulations of turbulent , radiation pressure dominated accretion flows now exist , both in local vertically stratified shearing box geometries @xcite and in global simulations @xcite . although these simulations use frequency integrated equations , the emergent radiation spectrum can be computed , including the effects of bulk comptonization , using post - processing monte carlo simulations . indeed , this has already been done by @xcite . however , in order for such calculations to be self - consistent , the frequency integrated radiation mhd equations used in the simulations themselves must include energy exchange due to bulk comptonization . we now discuss the consequences of the macroscopic physical origins of such energy exchange detailed in section [ sec_general ] for ensuring this effect is captured in simulations . we then proceed to derive the appropriate frequency integrated source terms due to compton scattering for the gas and radiation energy equations in both lab frame and fluid frame variables . using these results , we discuss the extent to which bulk comptonization is captured by existing radiation mhd simulation codes . the decomposition of bulk comptonization energy exchange into pressure work and radiation viscous dissipation shows that radiation mhd schemes that neglect contributions to the viscous stress tensor that are first order in velocity can not capture bulk comptonization energy exchange due to a shearing velocity field or any optically thin velocity field with non - zero divergence . as these effects are second order in velocity , we also note that a necessary , but not sufficient , condition for capturing these effects is inclusion of energy terms second order in velocity . without such terms , turbulence in this form , instead of exchanging energy with photons , will eventually cascade down to the gridscale ( or viscous or resistive scale if the code has explicit viscosity or resistivity ) , and increase the internal energy of the gas . gas internal energy may then be exchanged with photons through thermal comptonization . the omission of viscous dissipation by radiation therefore does not prevent the eventual transfer of turbulent energy to radiation , but it may have other physical effects that can in turn affect radiation spectra .    to derive the appropriate frequency integrated source terms due to comptonization , we start with the zeroth moment of the radiative transfer equation , correct to order @xmath45 and @xmath46 @xcite , @xmath47+\frac{v_i}{c } n^i\biggr\}\cr + \left({18\over5}+{17\over5}\epsilon{\partial\over\partial\epsilon } + { 11\over20}\epsilon^2{\partial^2\over\partial\epsilon^2}\right ) \left(n^{ij}{v_iv_j\over c^2}-{v^2\over3c^2}n\right ) . \label{eqpsaltis}\end{aligned}\ ] ] here @xmath48 is the photon energy , and the various angle - averaged moments are defined in terms of the energy and direction ( @xmath49 ) dependent photon occupation number @xmath50 by @xmath51 in principle , the energy equation is obtained by writing eq . ( [ eqpsaltis ] ) in terms of moments of the specific intensity and then integrating over all frequencies . unfortunately , we can not integrate over terms multiplied by @xmath52 without prior knowledge of the spectrum . for the purpose of simulations , then , we make two approximations . first , we observe that the fractional energy change per scattering off of non - relativistic electrons is small , so that only regions inside the photosphere contribute to comptonization . since the stimulated scattering terms are already order @xmath52 and in these regions departures from isotropy are order @xmath53 , we make the following approximation for these terms : @xmath54 second , we assume that the spectrum can be approximated by a bose - einstein distribution with temperature @xmath55 . with these approximations , eq . ( [ eqpsaltis ] ) yields @xmath56 this is the correct energy equation , but in order for it to capture second order energy exchange , the substituted value of @xmath57 must be calculated with a moment closure scheme that does not neglect contributions to the radiation viscous stress tensor that are first order in velocity . for example , it is not adequate to calculate @xmath57 by substituting @xmath58 into the first moment equation . here . ] this is equivalent to flux limited diffusion in the diffusion regime , such as that implemented in @xcite . these do , however , capture energy exchange due to pressure work in the optically thick regime . to show this , we substitute into eq . ( [ integrated1 ] ) the standard closure relation , @xmath59 which gives @xmath60 we see that energy exchange second order in velocity is not present . furthermore , we see that energy exchange due to a converging velocity field is indeed the work done by radiation pressure to compress the gas , @xmath61 . the m1 closure scheme @xcite , implemented in , e.g. , @xcite , also captures first order energy exchange but not second order energy exchange . this scheme assumes that there exists a frame in which @xmath62 . the lab frame radiation pressure tensor can then be expressed in terms of the energy density and flux @xcite : @xmath63 where @xmath64 and @xmath65 to show that second order energy exchange is not captured , we consider the case of a non - zero radiation viscous shear stress tensor due to a non - relativistic velocity field in an otherwise homogeneous medium . the lowest order contribution to the flux must be first order in velocity . in this scheme , then , the radiation viscous shear stress tensor is zero to first order in velocity , and hence second order energy exchange , which requires a contribution first order in velocity ( eq . [ energy_exchange ] ) , is not captured . another way to understand why both flux - limited diffusion and the m1 closure scheme fail to capture second order energy exchange is to observe that they both bridge generic optically thick conditions with optically thin conditions , while optically thin turbulence does not fall into either category . in optically thin turbulence , the turbulence length scales are optically thin ( @xmath66 ) , but conditions are otherwise optically thick ( @xmath67 ) , since we are far enough inside the photosphere that photons must scatter many times before escaping . it seems that only a more sophisticated approach , such as explicitly solving the transfer equation as done by @xcite , can capture this effect . we note that @xcite add an artificial viscosity to the m1 closure scheme in order to address a numerical problem associated with artificial shocks in their simulations . they assume a kinematic radiation viscosity given by @xmath68 which , as they acknowledge , underestimates the actual viscosity in the optically thick limit by a factor of @xmath69 ( eq . [ scattering_nu ] ) . since we show ( section [ sec_incompressible ] ) that the viscosity is smallest in the optically thick limit ( eq . [ scattering_nu_k ] ) , this artificial viscosity always underestimates the actual radiation shear stress tensor , and therefore only partially captures secord order energy exchange .    for completeness , we also write eq . ( [ integrated1 ] ) in terms of fluid frame radiation variables , indicated by subscript zero : @xmath70 since @xmath71 eq . ( [ integrated2 ] ) simplifies to @xmath72 @xcite ( hereafter s04 ) conjectured that comptonization by turbulence can be treated as thermal comptonization by solving the kompaneets equation with an equivalent  wave " temperature given by @xmath73 where @xmath74 is the temperature corresponding to the total electron energy density @xmath75 at wavenumber @xmath76 . equation ( [ twave1 ] ) generalizes s04 eq . ( 8) , @xmath77 which gives @xmath78 for a kolmogorov spectrum , @xmath79 , with maximum wavelength @xmath80 . equation ( [ twave1 ] ) is a weighting scheme of the form @xmath81 with the weighting function @xmath82 given by @xmath83 this scheme roughly gives full weight to wavelengths less than @xmath84 and zero weight to wavelengths greater than @xmath84 . for a periodic velocity field , since the modes are discrete , eq . ( [ eq_scheme_form ] ) is more clearly written @xmath85 where @xmath86 is the temperature of mode @xmath87 . for the remainder of this section we use eq . ( [ twave2 ] ) since it is more useful for applications to radiation mhd simulations , but note that all results also hold for eq . ( [ eq_scheme_form ] ) , as this is just the continuum limit . we also define @xmath88 = @xmath89 , the optical depth divided by @xmath90 across a mode with wavenumber @xmath76 . ( [ twave3 ] ) , for example , can then be written , @xmath91    the structure of this section is as follows . in section [ sec_exact ] we build on s04 by proving that comptonization by an arbitrary velocity field of uniform density can indeed be treated as thermal comptonization , provided that the field is divergenceless and the photon escape probability is sufficiently small . the exact equivalent wave temperature is of the form eq . ( [ twave2 ] ) , with the weighting function @xmath82 given by @xmath92 where @xmath93 . \label{eq : qdefinition}\ ] ] the limiting cases of @xmath82 are @xmath94 we plot @xmath95 and @xmath82 in fig . [ mcdata1 ] . since log scaling is used , the curve for @xmath95 disappears for @xmath96 . we see that @xmath95 does roughly approximate @xmath82 . in section [ sec_heur ] , to provide physical insight into the exact solution , we also find a heuristic , approximate expression for @xmath78 , also plotted in fig . [ mcdata1 ] , with a simple physical model motivated by ideas put forth in s04 . finally , in section [ sec_spec ] we use the correct expression for @xmath78 to discuss how the wave temperature depends on the power spectrum of the turbulence and determine which turbulent wavelengths contribute most to comptonization .      to derive eqs . ( [ twave2 ] ) and ( [ eq : analyticsol ] ) , we first find an expression for the wave temperature in terms of the second moment of the occupation number , since it is the latter which gives rise to the radiation shear stress tensor . we start with the zeroth moment of the radiative transfer equation , eq . ( [ eqpsaltis ] ) , correct to order @xmath45 and @xmath46 . note that terms of order @xmath46 are second order in velocity since for comptonized photons @xmath97 . consider a solution for the angle dependent occupation number @xmath98 constructed with the expansion @xmath99 where @xmath100 is formally @xmath87th order in velocity . we make the ansatz that @xmath101 is independent of position and angle , so that to second order in velocity eq . ( [ eq_expansion_1 ] ) is @xmath102 to second order in velocity the zeroth moment equation then simplifies to @xmath103 we also invoke the first moment of the transfer equation , @xcite eq . we multiply this by @xmath104 , so that in steady state , to second order it becomes @xmath105 substituting eq . ( [ eq_first_moment ] ) into eq . ( [ eq_zeroth_simplified ] ) gives @xmath106 we average this over space to obtain @xmath107 this is the kompaneets equation with a bulk comptonization contribution to the temperature given by @xmath108    to evaluate eq . ( [ eq_twave_n ] ) , we need to derive @xmath109 to first order in velocity . we start with the radiative transfer equation , eq . ( a1 ) of @xcite , correct to order @xmath45 and @xmath46 , and substitute in eq . ( [ eq_expansion_2 ] ) . in steady state , to first order in velocity , the transfer equation then simplifies to @xmath110 if the density is constant and the velocity field is divergenceless with sinusoidal mode decomposition @xmath111 then the solution is given by ( appendix [ transverse_solution ] ) @xmath112 where @xmath82 is given by eq . ( [ eq : analyticsol ] ) .    finally , we evaluate eq . ( [ eq_twave_n ] ) by plugging in eq . ( [ eq_nij ] ) to get @xmath113 since the modes are sinusoidal , eq . ( [ eq_twave4 ] ) simplifies to @xmath114 which is equivalent to eq . ( [ twave2 ] ) . if , however , one solves the kompaneets equation with an escape probability term , @xmath115 , we must also add such a term to the radiative transfer equation , which may affect the occupation number second moment . for an optically thick system , for example , we must require that this term be small compared to the term that sets the diffusion time scale . this term comes from the term @xmath116 , which is approximated by @xmath117 the diffusion timescale is set by the second derivative term , @xmath118 which gives the condition @xmath119    we now discuss the extent to which , in this case , radiation viscous dissipation can be described by a coefficient of kinematic viscosity @xmath120 , which is usually defined by @xmath121 where @xmath122 is the fluid mass density . first we find the radiation shear stress tensor corresponding to eq . ( [ eq_nij ] ) by frequency integrating and subtracting off the scalar radiation pressure , which gives @xmath123 incidentally , we note that eq . ( [ eq_twave_n ] ) can then alternatively be written in terms of the radiation viscous shear stress tensor , @xmath124 according to eqs . ( [ eq_kinematic ] ) and ( [ eq_vis_p ] ) , although we can not in general define a kinematic viscosity coefficient , for any single velocity mode the kinematic radiation viscosity is given by @xmath125 so that the radiation viscous shear stress tensor can be written @xmath126 however , in the optically thick limit ( i.e. @xmath127 ) , the kinematic viscosity is independent of @xmath128 , @xmath129 so that in this limit the kinematic viscosity is well - defined for an arbitrary ( divergenceless ) velocity field . we also note that in this limit eq . ( [ eq_twave4 ] ) agrees with the  heating temperature " in @xcite , @xmath130 they also mention that this can be understood in terms of the kinematic radiation viscosity given by eq . ( [ scattering_nu ] ) . from the more commonly cited value , @xmath131 , in , e.g. , @xcite , @xcite , and @xcite . the reason for the discrepancy is that the more commonly cited value , first derived by @xcite , assumes pure absorption , while eq . ( [ scattering_nu ] ) is correct for pure scattering ( @xcite ; @xcite ) . ]      to provide physical insight into our analytic solution given by eqs . ( [ twave2 ] ) and ( [ eq : analyticsol ] ) , we now find a heuristic , approximate expression for @xmath132 with a simple physical model motivated by ideas put forth in s04 . to do this , we first consider the wave temperature of a single mode with wavenumber @xmath133 , which we denote @xmath134 . s04 suggested that photons can only sample turbulent velocities on scales @xmath135 , since longer turbulent wavelengths will advect photons back and forth with the flow without allowing the photons to  sample " their velocities . a rough interpretation of this reasoning leads to @xmath136 if one then assumes that the wave temperature of an arbitrary field is the sum of the wave temperatures of its modes , then eq . ( [ twave1 ] ) follows . but a more subtle interpretation of this reasoning suggests that @xmath134 is determined by the second moment of the distribution of velocity differences between subsequent scatterings , @xmath137 . in the long wavelength limit , velocity differences between subsequent scatterings are negligible and so the wave power does not contribute to comptonization . in the short wavelength limit , velocity differences allow photons to sample the full power of the wave . this model of comptonization suggests that we define @xmath138 where @xmath139    before proceeding , we note that defining @xmath137 is potentially tricky because the distribution of @xmath140 for a photon is dependent on its current location , in effect introducing correlations into subsequent @xmath140 s . in other words , subsequent @xmath140 s are not independent . but if the escape probability is low enough , a condition we quantify below , then the set of @xmath140 s encountered by a photon before it escapes is indistinguishable from a set of @xmath140 s independently drawn from the position - averaged @xmath140 distribution . the order of @xmath140 s is different in the two cases , but the total photon energy change does not depend on the order because the fractional photon energy change per scattering is small for @xmath141 . with these potential problems accounted for , we proceed to calculate @xmath142 . first we find the proportionality constant between @xmath142 and @xmath137 by evaluating both sides in the short wavelength limit . in this limit , the full wave power must contribute , so @xmath143 . to evaluate @xmath137 in this limit , let @xmath144 be any normalized distribution of velocities . then , @xmath145 therefore , @xmath146    we now calculate @xmath147 for the position - averaged @xmath140 for a single , divergenceless ( i.e. transverse ) mode with wavelength @xmath148 , @xmath149 let @xmath150 be the probability density that a photon travels a displacement @xmath151 between scatterings . then , at a given position @xmath152 , @xmath153 averaging over all positions in a volume @xmath154 , this is @xmath155 for a single mode , eq . ( [ vel_field ] ) gives @xmath156 and @xmath157 the probability that a photon with mean free path @xmath158 travels a distance between @xmath159 and @xmath160 is @xmath161 . let @xmath162 , where @xmath163 is the angle between the photon propagation direction and the @xmath164-axis , so that @xmath165 . then , expressing @xmath166 in spherical polar coordinates and invoking axisymmetry about the @xmath164-axis , equation ( [ eq : dv2deltar ] ) becomes @xmath167 ^ 2 e^{-s / l}.\end{aligned}\ ] ] this is easily evaluated by performing the integral over @xmath164 first , giving @xmath168 . \\ & = 2 \sigma_{{\bf v}}^2\left[1 - \tau_k \tan^{-1}\left(\frac{1}{\tau_k}\right)\right].\end{aligned}\ ] ] eq . ( [ fheur ] ) then gives @xmath169 the limiting cases are @xmath170 by comparison , the exact solution for a single mode is determined by eq . ( [ eq : analyticsol ] ) . our heuristic result is plotted in fig . [ mcdata1 ] , and is remarkably close to the exact solution . in particular , it is a much better approximation than the rough weighting function , eq . ( [ twave3 ] ) , also shown in fig . [ mcdata1 ] . our model based on the second moment of the velocity difference distribution therefore captures the essential physics of comptonization by a single mode .    before proceeding to define @xmath171 for an arbitrary velocity field , we quantify the condition that the escape probability be low enough , presupposed by our derivation of @xmath147 . since the distribution of velocity differences , as a function of position , repeats every quarter wavelength , our results should be valid provided photons travel a distance @xmath172 in the @xmath164 direction that is greater than @xmath173 before escaping . for an ( optically thick ) random walk , @xmath174 , where @xmath175 is the average number of scatterings . since @xmath176 , where @xmath177 is the escape probability per unit time during the average time @xmath178 between subsequent scattering events , @xmath179 for optically thick modes . up to a factor of unity , this agrees with eq . ( [ eq_escape_constraint1 ] ) . we now use our model to compute @xmath171 for an arbitrary velocity field in terms of @xmath180 . proceeding analogously to the single mode case , we define @xmath181 where @xmath182 is the temperature corresponding to the full electron energy density of the velocity field . to simplify this , we compute @xmath183 ^ 2 p_{\delta { \bf r}}(\delta { \bf r } ) d^3 \delta { \bf r } d^3 { \bf r } \notag\\ & = \frac{1}{v}\int \sum_i \left({\bf v}_i({\bf r}+\delta { \bf r})-{\bf v}_i ( { \bf r})\right)^2 p_{\delta { \bf r}}(\delta { \bf r } ) d^3 \delta { \bf r } d^3 { \bf r } \notag\\ % & = \frac{1}{v}\int \sum_i \left ( \delta { \bf v}_i(\delta { \bf r},{\bf r } ) % \right)^2 p_{\delta { \bf r}}(\delta { \bf r } ) d^3 \delta { \bf r } % d^3 { \bf r}\notag \\ & = \sum_i \langle ( \delta { \bf v}_i)^2 \rangle . \label{deltavcalc}\end{aligned}\ ] ] to get from line 2 to line 3 we made use of the orthogonality for distinct sinusoidal modes . that is , for two distinct modes , @xmath184 and @xmath185 , @xmath186 , and any displacement @xmath187 , @xmath188 then , since @xmath189 , eqs . ( [ fheur ] ) , ( [ theur ] ) , and ( [ deltavcalc ] ) give @xmath190 or , alternatively , @xmath191 eq . ( [ theur2 ] ) is the same as the exact solution for the wave temperature , eq . ( [ twave2 ] ) , except that here the heuristic weighting function is used . note that in our heuristic derivation it is the orthogonality of distinct modes that allows us to express the wave temperature of an arbitrary velocity field as a sum over the wave temperatures of its modes . unsurprisingly , orthogonality was used analogously to prove eq . ( [ twave2 ] ) . therefore , our model based on the second moment of the velocity difference distribution captures the essential physics of comptonization by an arbitrary , divergenceless velocity field . we now analyse the dependence of @xmath132 on the power spectrum of the turbulence . for the remainder of this section , we write @xmath76 in units of @xmath192 for clarity ( i.e. @xmath193 and @xmath194 denote optically thick and thin scales , respectively ) . if the turbulence is completely optically thin on all scales ( @xmath195 ) , then @xmath196 , independent of the energy spectral index , @xmath197 . however , if some scales in the turbulent cascade are optically thick , then @xmath197 will affect @xmath132 .    for the case where all scales in the turbulence are optically thick ( @xmath198 ) , eq . ( [ eq : analyticsol ] ) implies that @xmath199 integrating this over an energy spectrum @xmath200 then gives @xmath201 for a broad power spectrum , i.e. @xmath202 , this simplifies to @xmath203 this is illustrated in fig . [ ueff_p ] . note that @xmath204 drops significantly for @xmath205 , because then the energy bearing modes are on the largest scales . these are the most optically thick and therefore the most downweighted in their contribution to bulk comptonization . we next analyse whether @xmath204 , for a given spectral index @xmath197 and range of modes , @xmath206 , is dominated by small or large scales . in other words , we examine which turbulent modes in a given spectrum contribute most to bulk comptonization . the relative contribution of a scale with wavenumber @xmath76 is @xmath207 where , from eq . ( [ eq : analyticsol ] ) , @xmath208 for optically thick ( @xmath209 ) scales , and @xmath210 for optically thin ( @xmath211 ) scales . now consider an underlying power spectrum with some @xmath212 and @xmath213 . we see that for @xmath214 the exponent in eq . ( [ contribution ] ) is always positive , and so small scales contribute most to bulk comptonization , regardless of @xmath212 and @xmath213 . this is physically intuitive ; for @xmath214 , the turbulent power is concentrated on small scales . since the weighting factor @xmath82 also favors small scales , they of course contribute most . for @xmath215 , the exponent is always negative , and so large scales always contribute most . in this case , the turbulent power is so concentrated on large scales that they contribute more even though @xmath82 favors small scales .    for @xmath216 , we first consider the part of the spectrum with @xmath194 ( if it exists ) . since @xmath210 , small scales contribute more than large scales for these modes . now consider the part of the spectrum with @xmath193 ( if it exists ) . here , large scales contribute more than small scales . therefore , looking at the entire power spectrum , it is intermediate scales that contribute most , assuming it is broad enough to include regions of both small and large @xmath76 . if it is not sufficiently broad , then whether small or large scales contribute most depends on @xmath212 and @xmath213 relative to @xmath217 ( the optically thin to thick transition wavenumber ) . these results are depicted in fig . [ ueff_dist ] . the curve in this figure shows the values of @xmath197 and @xmath76 such that the derivative of @xmath218 is zero , using the full analytic expression for @xmath82 from eq . ( [ eq : analyticsol ] ) . to connect this figure to our discussion , draw a horizontal line from @xmath212 to @xmath213 at a given value of @xmath197 . if the line lies in the lower ( upper ) region , then small ( large ) scales contribute most . if the line straddles the two regions , then for the part of the line that lies in the lower region small scales contribute most , and for the part that lies in the upper region large scales contribute the most . in this case , then , for the entire spectrum it is the scales that straddle the curve which contribute most . note that for @xmath214 and @xmath215 a spectrum can never straddle the curve , whereas for a kolmogorov spectrum ( @xmath219 ) , e.g. , it can , if @xmath220 . @xcite conjectured that comptonization by turbulence can be treated as thermal comptonization by solving the kompaneets equation with an equivalent  wave " temperature . in section [ sec_incompressible ] we proved this under certain conditions , one of which is that the turbulence be divergenceless . in this section we investigate comptonization by velocity fields with non - zero divergence , an effect that usually can not be treated as thermal comptonization , with the aim of understanding how it impacts radiation spectra in generic , turbulent regions of stratified accretion disc atmospheres . while we do not have an exact solution for comptonization by an arbitrary velocity field with non - zero divergence , we have solutions for two limiting cases . the trivial case is the optically thin case , i.e. when the mean free path is significantly larger than the largest length scale in the turbulence , @xmath221 . electron velocities  sampled " by photons are uncorrelated and so compton scattering should depend on only the total spatial average distribution of electron velocities . this is , therefore , the one case where comptonization by a velocity field with non - zero divergence can be treated as thermal comptonization , by solving the kompaneets equation with @xmath132 given by @xmath222 note that although we formally write @xmath132 as a function of @xmath84 , in this limit @xmath132 is independent of @xmath84 . energy exchange first order in velocity vanishes since photons are equally likely to downscatter as they are to upscatter . bulk comptonization is then solely due to radiation viscous dissipation ( see section [ sec_general ] ) . as viscous effects are second order in velocity , it is unsurprising that they can be characterized by a temperature . we also note that in this limit the wave temperature is the same as that for a divergenceless velocity field , eq . ( [ eq_limits ] ) . optically thin bulk comptonization is therefore a single phenomenon that depends on only the mean square speed of the velocity field . to arrive at this result with a more formal approach , we start with the zeroth moment of the radiative transfer equation , eq . ( [ eqpsaltis ] ) . in the limit of optically thin turbulence the radiation variables must be homogeneous and isotropic , so that @xmath223 , and @xmath224 . then , averaging eq . ( [ eqpsaltis ] ) over the largest scale @xmath35 gives @xmath225\biggr\}. \label{eqkomp}\end{aligned}\ ] ] this is the kompaneets equation , with the contribution from the velocity field to the comptonization temperature given by eq . ( [ compressible_thin ] ) .    in the optically thick case , i.e. when the photon mean free path is significantly smaller than the smallest scale in the turbulence , the lowest order energy exchange is the work done by radiation pressure to compress the gas , since it is first order in velocity and since radiation viscous effects are suppressed by a factor @xmath44 ( section [ sec_incompressible ] ) . we focus on the extremely optically thick case , which we define as the limit in which photon diffusion is negligible relative to photon advection , so that the photon and gas fluids are strongly coupled . if we define @xmath226 , the photon number density at energy @xmath48 , and @xmath227 , the total photon number density , then the advection and diffusion fluxes are given by @xmath228 and @xmath229 respectively . the extremely optically thick limit is then given by @xmath230 in this case , velocity convergence corresponds to compression in which gas mechanical energy is transferred locally to the photons . we expect that photons with wavelength @xmath231 are effectively compressed at a rate given by the velocity difference across @xmath231 , @xmath232 so that , for example , a locally thermal photon distribution remains thermal and only changes temperature , completely analogous to the evolution of the cosmic microwave background radiation under the expansion of the universe . this is equivalent to a fractional energy change per scattering given by @xmath233 since @xmath234 we now confirm that eq . ( [ upscattering_fraction ] ) correctly describes extremely optically thick compton scattering in a converging or diverging flow , both by providing a heuristic argument and by deriving it from the radiative transfer equation .    before proceeding , we note that the study of photon upscattering by a converging velocity field can be traced back to @xcite and @xcite , who , along with later authors , made detailed spectral calculations for specific velocity fields in shocks and spherically accreting systems . in fact , eq . ( [ upscattering_fraction ] ) can equivalently be stated as the upscattering timescale @xmath235 given in @xcite . in this section , by contrast , we have been investigating how this effect manifests itself locally in a generic , turbulent region of a stratified disc atmosphere , with the goal of resolving and interpreting it in spectral calculations of radiation mhd simulations . we have been focusing on the extremely optically thick case , in which the photon and gas fluids are strongly coupled , because the physics is both relevant and intuitive . in the moderately optically thick case , on the other hand , i.e. @xmath236 such as photon upscattering in a radiation pressure dominated shock @xcite , diffusion competes with advection so that photon distributions at neighboring fluid elements mix . photon upscattering in such a converging flow may not be viewed as simply the compression of a photon fluid strongly coupled to the gas , and photon upscattering in which a photon thermal distribution is not preserved can occur .    to heuristically derive eq . ( [ upscattering_fraction ] ) , consider a disturbance converging in the @xmath237 direction given by @xmath238 for a 3d random walk , the average distance between scatterings traveled by a photon in the direction of convergence is @xmath239 . since the fractional energy change per scattering for low energy photons is @xmath53 , at @xmath240 this gives @xmath241 in agreement with eq . ( [ upscattering_fraction ] ) . we now derive eq . ( [ upscattering_fraction ] ) with the radiative transfer equation , eq . ( [ eqpsaltis ] ) . if we ( 1 ) omit stimulated scattering terms to facilitate comparison with simulations , ( 2 ) substitute in the standard closure relations for the first moment in the optically thick limit , @xmath242 and @xmath243 and ( 3 ) substitute in the photon number density @xmath244 , then the radiative transfer equation to second order in velocity and first order in @xmath245 becomes @xmath246 neglecting stimulated scattering , this is @xcite eq . ( 18 ) , cast in the physically revealing form of a fokker - plank equation . the terms inside the divergence operator correspond to spatial drift ( i.e. photon advection ) and spatial diffusion , respectively . the next term corresponds to energy drift due to photon upscattering ( downscattering ) in the presence of a converging ( diverging ) velocity field . the remaining terms correspond to energy drift and diffusion due to thermal comptonization . note that even though @xcite start with a zeroth moment equation correct only to first order in @xmath53 , their resulting equation is the same because with the standard closure relation ( eq . [ closure1 ] ) the second order terms cancel . the fractional energy change per scattering given by eq . ( [ upscattering_fraction ] ) follows from the bulk upscattering term . in the extremely optically thick limit the spatial diffusion term is negligible , and so photons are advected with the velocity field and upscatter according to eq . ( [ upscattering_fraction ] ) . the effect of this process on the emergent spectrum of the disc will depend primarily on how effectively photons are able to escape from converging ( or diverging ) regions to the observer . we have investigated how photon spectra produced by turbulent comptonization depend on the properties of the turbulence itself , and how to resolve and interpret this effect in radiation mhd simulations . our principle results are as follows . first we address the conjecture of @xcite that turbulent comptonization can be treated as thermal componization by solving the kompaneets equation with an equivalent  wave " temperature . for comptonization by divergenceless turbulence , which is due to the radiation viscous shear stress tensor and is therefore an effect second order in velocity , this conjecture holds provided that the density is uniform and the escape probability is sufficiently small ( eq . [ eq_escape_constraint1 ] ) . the correct wave temperature can be expressed in terms of the radiation viscous shear stress tensor ( eqs . [ eq_twave_n ] and [ eq_twave_p ] ) and is in turn given by eq . ( [ twave2 ] ) . this is a weighting over the velocity power at each scale in the turbulence , with the weighting function given by eq . ( [ eq : analyticsol ] ) . optically thin scales have unit weight factors , implying that the power at these scales contributes fully to the wave temperature . optically thick scales , on the other hand , are downweighted by a factor proportional to one over the square of the optical depth across those scales . this downweighting arises because bulk comptonization is sensitive to the difference in bulk fluid velocity between subsequent scatterings , as demonstrated by our heuristic argument in section [ sec_incompressible ] . these velocity differences are much smaller than the velocity amplitude on some scale if the photon mean free path is much smaller than that scale . we note that our expression for the wave temperature of divergenceless turbulence assumes a uniform density . furthermore , it is non - local in that it is expressed in terms of the mode decomposition of the entire velocity field . although a generic velocity field has density variations , if the largest turbulent length scale @xmath35 is smaller than the scale on which the density varies then we can define the wave temperature for regions of order @xmath35 . if the density also varies on smaller scales , then we can define the wave temperature for smaller regions by omitting contributions from larger modes . assuming these modes would be significantly downweighted anyway , this procedure would likely capture nearly all comptonized power in the turbulence .    for turbulence with non - zero divergence , bulk comptonization energy exchange is due to both work done by radiation pressure and radiation viscous dissipation . since pressure work is an effect first order in velocity , it can not be treated with an equivalent temperature . but in the optically thin limit , i.e. when the mean free path is significantly larger than the largest length scale in the turbulence , energy exchange first order in velocity vanishes since photons are equally likely to downscatter as they are to upscatter . bulk comptonization by a velocity field with non - zero divergence is then solely due to radiation viscous dissipation , and it may be treated as thermal comptonization . the wave temperature is given by the full power in the turbulence , eq . ( [ compressible_thin ] ) , directly analogous to the optically thin limit of divergenceless turbulence .    in the optically thick case , i.e. when the mean free path is significantly smaller than the smallest scale in the turbulence , the lowest order energy exchange is the work done by radiation pressure to compress the gas , since it is first order in velocity and since radiation viscous effects are suppressed by a factor @xmath44 . the extremely optically thick limit , defined as the limit in which spatial photon diffusion is negligible relative to photon advection , is given by eq . ( [ eq_extremely_optically_thick ] ) . in this limit , the photon and gas fluids are strongly coupled , so that velocity convergence corresponds to compression in which gas mechanical energy is transferred locally to the photons . photons upscatter with fractional energy change given by eq . ( [ upscattering_fraction ] ) . in this process a locally thermal photon distribution remains thermal and only changes temperature , analogous to the evolution of the cosmic microwave background radiation under the expansion of the universe . unlike energy exchange due to viscous dissipation , this process is reversible . the effect of this process on the emergent spectrum of the disc will depend primarily on how effectively photons are able to escape from such regions to the observer .    for moderate optical depths , such as photon upscattering in a radiation pressure dominated shock @xcite , diffusion competes with advection so that photon distributions at neighboring fluid elements mix . the effect of photon upscattering in a converging region on local spectra is more complex and depends on the details of the flow . it is important to note that optically thick comptonization by a velocity field with non - zero divergence may be very sensitive to the time dependence of the velocity field . this is important because usually post - processing monte - carlo simulations invoke time - independent atmospheres with the assumption that they approximately capture the effects of interest . for example , consider comptonization in a converging region . a time - independent velocity field results in the accumulation of photons and subsequent upscattering to high energies at the point of zero velocity . if the region is near the photosphere , the emergent spectrum will be strongly upscattered . but in a time - dependent velocity field photons have a limited time to upscatter before converging regions become diverging regions , and so the effect on the emergent spectrum is significantly different . real accretion flows are spatially stratified . simulations of magnetorotational turbulence generally indicate that this turbulence dominates the fluid velocities in the midplane regions of the accretion flow . this turbulence is largely incompressible ( divergenceless ) , although it generally excites compressible spiral acoustic waves @xcite . sufficiently far from the midplane , magnetic forces always dominate thermal pressure forces , and support the flow vertically against the tidal gravitational field of the compact object . such regions are dominated by parker instability dynamics , and exhibit considerable compressive behavior ( i.e. the flow has non - zero divergence ) with large density fluctuations @xcite . in order for radiation mhd simulations to properly account for energy exchange due to turbulent comptonization so that post - processing monte carlo simulations of photon spectra are self - consistent , they must include energy terms second order in velocity and use a moment closure scheme that correctly captures contributions to the radiation stress tensor that are first order in velocity . the appropriate energy equation source terms in lab and fluid frame variables are given by eqs . ( [ integrated1 ] ) and ( [ eq_intergrated2 ] ) , respectively . flux - limited diffusion and the m1 closure scheme are insufficient because they neglect the lowest order contribution to the radiation stress tensor . in this case , instead of exchanging energy with photons , turbulence will eventually cascade down to the gridscale ( or viscous or resistive scale if the code has explicit viscosity or resistivity ) and increase the internal energy of the gas . gas internal energy may then be exchanged with photons through thermal comptonization . the omission of viscous dissipation by radiation therefore does not prevent the eventual transfer of turbulent energy to radiation , but it may have other physical effects which can in turn affect radiation spectra . radiation pressure dominated accretion flows may have bulk velocities that exceed thermal electron velocities . turbulent comptonization , which is the effect of comptonization by bulk , turbulent electron motions on photon spectra , should be significant in such sources .    for turbulence with zero divergence , bulk comptonization is equivalent to a radiation viscous shear stress tensor acting on the fluid . the resulting comptonized spectrum can be calculated by solving the kompaneets equation with an equivalent  wave " temperature , which is a weighted sum over the power present at each scale in the turbulent cascade ( eq . [ twave2 ] ) . bulk comptonization by turbulence with non - zero divergence is due to both pressure work and radiation viscous dissipation . in the limit of optically thin turbulence , pressure work has negligible effect on photon spectra and radiation viscous dissipation alone can be treated as thermal comptonization with a temperature equivalent to the full turbulent power ( eq . [ compressible_thin ] ) . in the limit of extremely optically thick turbulence ( eq . [ eq_extremely_optically_thick ] ) , radiation viscous dissipation is suppressed , and the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids . photons upscatter with fractional energy change given by eq . ( [ upscattering_fraction ] ) . in this process a locally thermal photon distribution remains thermal and only changes temperature . the effect of this process on the emergent spectrum of the disc will depend primarily on how effectively photons are able to escape from such regions to the observer . modeling turbulent comptonization will ultimately require detailed analysis of radiation mhd simulations and post - processing monte - carlo simulations . by exploring how photon spectra produced by turbulent comptonization depend on the properties of the turbulence itself , we have laid the groundwork necessary to make sure this effect is both captured and correctly interpreted in these simulations . we thank aristotle socrates for useful discussions . we also thank donald pakey for forwarding errata to the paper by @xcite , and li xue and feng yuan for hosting ob at xiamen university and shanghai astronomical observatory where part of this work was carried out . we gratefully acknowledge support by nasa astrophysics theory program grant nnx13ag61 g , and the international space science institute ( issi ) in bern . bachetti m. et al . , 2013 , apj , 778 , 163 blaes o. m. , hirose s. , krolik j. h. , 2007 , apj , 664 , 1057 blaes o. m. , krolik j. h. , hirose s. , shabaltas n. , 2011 , apj , 733 , 110 blandford r. d. , payne d. g. , 1981a , mnras , 194 , 1033 blandford r. d. , payne d. g. , 1981b , mnras , 194 , 1041 chan k.l . , jones b.j.t . , 1975 , apj , 200 , 454 colpi m. , 1988 , apj , 326 , 223 collinson j. s. , ward m. j. , done c. , landt h. , elvis m. , & mcdowell j. c. , 2015 , mnras , in press [ arxiv:1502.05850 ] davis s. w. , blaes o. m. , hirose s. , krolik j. h. , 2009 , apj , 703 , 569 davis s. w. , blaes o. m. , hubeny i. , turner n. j. , 2005 , apj , 621 , 372 davis s. w. , done c. , blaes o. m. , 2006 , apj , 647 , 525 done c. , davis s. w. , jin c. , blaes o. m. , ward m. , 2012 , mnras , 420 , 1848 gladstone j. c. , roberts t. p. , done c. , 2009 , mnras , 397 , 1836 guan x. , gammie c. f. , 2011 , apj , 728 , 130 haardt f. , maraschi l. , 1991 , apj , 380 , 51 heinemann t. , papaloizou j. c. b. , 2009a , mnras , 397 , 52 heinemann t. , papaloizou j. c. b. , 2009b , mnras , 397 , 64 hirose s. , krolik j. h. , stone j. m. , 2006 , apj , 640 , 901 hirose s. , krolik j. h. , blaes o. m. , 2009 , apj , 691 , 16 hubeny i. , blaes o. m. , krolik j. h. , agol e. , 2001 , apj , 559 , 680 jiang y .- f . , stone j. m. , davis s. w. , 2013 , apj , 778 , 65 jiang y .- f . , stone j. m. , davis s. w. , 2014a , apj , 784 , 169 jiang y .- f . , stone j. m. , davis s. w. , 2014b , apj , 796 , 106 kawashima t. , ohsuga k. , mineshige s. , yoshida t. , heinzeller d. , matsumoto r. , 2012 , apj , 752 , 18 kolehmainen m. , done c. , daz trigo m. , 2011 , mnras , 416 , 311 laor a. , netzer h. , 1989 , mnras , 238 , 897 lightman a. p. , white t. r. , 1988 , apj , 335 , 57 levermore c. d. , 1984 , j. quant . transfer , 31 , 149 magdziarz p. , blaes o. m. , zdziarski a. a. , johnson w. n. , & smith d. a. , 1998 , mnras , 301 , 179 masaki i. , 1971 , pub . japan , 23 , 425 mckinney j. c. , tchekhovskoy a. , sadowski a. , narayan r. , 2014 , mnras , 441 , 3177 mehdipour m. et al . , 2011 , a&a , 534 , a39 mihalas d. , mihalas b. , 1984 , foundations of radiation hydrodynamics ( new york : oxford univ . press ) middleton m. , done c. , ward m. , gierliski m. , schurch n. , 2009 , mnras , 394 , 250 miller k. a. , stone j. m. , 2000 , apj , 534 , 398 niedwiecki a. , zdziarski a. a. , 2006 , mnras , 365 , 606 ohsuga k. , mineshige s. s. , 2011 , apj , 736 , 2 payne d. g. , blandford r. d. , 1981 , mnras , 196 , 781 pozdnyakov l. a. , sobol i. m. , & sunyaev r. a. , 1983 , astrophys . space phys . , 2 , 189 psaltis d. , 2001 , apj , 555 , 786 psaltis d. , & lamb f. k. , 1997 , apj , 488 , 881 ross r. r. , fabian a. c. , 1993 , mnras , 261 , 74 sadowski a. , narayan r. , tchekhovskoy a. , abarca , d. , zhu y. , mckinney , j. , 2015 , mnras , 447 , 49 sadowski a. , narayan r. , tchekhovskoy a. , zhu y. , 2015 , mnras , 429 , 3533 shakura n. i. , sunyaev r. a. , 1973 , a&a , 24 337 shimura t. , takahara f. , 1995 , apj , 445 , 780 socrates a. , 2010 , apj , 719 , 784 socrates a. , davis s. w. , blaes o. m. , 2004 , apj , 601 , 405 stobbart a .- m . , roberts t. p. , wilms j. , 2006 , mnras , 368 , 397 straumann n. , 1976 , helv . acta , 49 , 269 svensson r. , zdziarski a. a. , 1994 , apj , 436 , 599 takeuchi s. , ohsuga k. , mineshige s. , 2013 , pasj , 65 , 88 thomas , l.h . , 1930 , q.j . math . , 1 , 329 thompson c. , 1994 , mnras , 270 , 480 titarchuk l. , mastichiadis a. , kylafis n. , 1997 , apj , 487 , 834 weinberg s. , 1971 , apj , 168 , 175 weinberg s. , 1972 , gravitation and cosmology ( new york : wiley ) yoshida t. , isobe n. , mineshige s. , kubota a. , mizuno t. , saitou k. , 2013 , pasj , 65 , 48 zeldovich ya . b. , illarionov a. f. , sunyaev r. a. , 1972 , sov . jetp , 35 , 643 beginning with eq . ( [ eq_transfer_first_order ] ) , we prove that the steady state occupation number second moment for a divergenceless velocity field of uniform density to first order in velocity is given by eq . ( [ eq_nij ] ) . first we find the solution for a single mode given by @xmath247 for this mode , the transfer equation is @xmath248 where we assume @xmath249 and @xmath250 , which we can check later . then , the transfer equation is @xmath251 first we address the @xmath164-dependence . because this equation is linear , it must be that @xmath252 is a superposition of a sine and a cosine , @xmath253 this gives two coupled integral equations for @xmath254 and @xmath255 , @xmath256 and @xmath257 it seems that both @xmath254 and @xmath255 are proportional to one power of @xmath258 . writing @xmath259 and @xmath260 , and @xmath261 , we then obtain @xmath262 and @xmath263 observing that letting @xmath264 and @xmath265 be odd and even , respectively , is a consistent solution , the @xmath266 integral of @xmath265 vanishes , and the two equations can be combined to give a single equation for @xmath264 , @xmath267 we then solve this equation with a series expansion of odd powers of @xmath268 : @xmath269 gives @xmath270 which is just the expansion of @xmath271 , so that @xmath272 substituting this back into the integral equation gives @xmath273 , which completes the solution . so far then , we have @xmath274\sin \left({2\pi z\over\lambda}\right)\biggr\},\end{aligned}\ ] ] where @xmath275 . \label{eq : qexact}\ ] ] note that @xmath276 in the optically thin limit , and @xmath277 in the optically thick limit . this solution is consistent with our assumptions and solves eq . ( [ eq_transfer_first_order ] ) to first order in velocity . the second moment is @xmath278 where @xmath279 is given by eq . ( [ eq : analyticsol ] ) . ( [ eq_transfer_first_order ] ) is linear , the solution for an arbitrary , divergenceless velocity field of uniform density is then given by eq . ( [ eq_nij ] ) .
Abstract:

Answer: radiation pressure dominated accretion discs around compact objects may have turbulent velocities that greatly exceed the electron thermal velocities within the disc . bulk comptonization by the turbulence may therefore dominate over thermal comptonization in determining the emergent spectrum . bulk comptonization by divergenceless turbulence is due to radiation viscous dissipation only . it can be treated as thermal comptonization by solving the kompaneets equation with an equivalent  wave " temperature , which is a weighted sum over the power present at each scale in the turbulent cascade . bulk comptonization by turbulence with non - zero divergence is due to both pressure work and radiation viscous dissipation . pressure work has negligible effect on photon spectra in the limit of optically thin turbulence , and in this limit radiation viscous dissipation alone can be treated as thermal comptonization with a temperature equivalent to the full turbulent power . in the limit of extremely optically thick turbulence , radiation viscous dissipation is suppressed , and the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids . we discuss the consequences of these effects for self - consistently resolving and interpreting turbulent comptonization in spectral calculations in radiation mhd simulations of high luminosity accretion flows . [ firstpage ]    accretion , accretion discs  radiation mechanisms : non - thermal  turbulence  galaxies : active  x - rays : binaries .