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: over the last few years , several methods have been developed to cool molecules to temperatures below 1 degree kelvin . these include buffer - gas cooling in cryogenic helium @xcite , stark deceleration with switched electric fields @xcite , velocity filtering @xcite , optical deceleration with laser fields @xcite , and crossed molecular beam scattering @xcite . the availability of cold molecules opens up a new field of low - energy collision studies , in a novel regime where collisions are dominated by long - range forces and resonances . stark decelerators can be used to _ control _ beam velocity as well as to reduce it @xcite . pulsed molecular beams can be generated with much smaller velocity spreads than is possible with conventional supersonic sources . this offers the opportunity to study scattering resonances at much higher resolution than has been possible in the past . _ @xcite have carried out a proof - of - concept experiment in which a velocity - controlled beam of oh radicals collided with xe atoms in a jet . in this case the energy resolution of the experiment was limited to 13 @xmath1 by the velocity spread of the xe atoms . however , experiments are under way to collide two velocity - controlled beams , which will provide much higher resolution . an alternative approach is to collide a beam with a sample of trapped atoms or molecules that are already nearly at rest . _ @xcite have recently measured collision cross sections for he atoms and h@xmath3 molecules in pulsed beams colliding with magnetically trapped oh radicals , and achieved a resolution of 9 @xmath1 . experiments to investigate the collisions of a velocity - controlled beam of nd@xmath0 with trapped ultracold rb atoms are also under way , and again should be able to provide much better velocity resolution . scattering resonances in molecular collisions have been studied for many years @xcite and can be very important in chemical reactions @xcite . in simple systems where only a few partial waves contribute , they may produce sharp structures in cross sections as a function of collision energy @xcite . however , in more complicated systems with dense energy level patterns , the resonant structures may get lost in the background . it is therefore important to explore whether well - defined resonant structures are expected for collisions of molecules that can be velocity - controlled ( such as nd@xmath0 and oh ) and atoms that can be laser - cooled ( such as alkali metal atoms ) .    in this paper we study collisions of rb atoms with nd@xmath0 and nh@xmath0 for collision energies between 0 and 100 @xmath1 . we carry out quantum - mechanical calculations of integral elastic and state - to - state inelastic cross sections using the coupled states approximation and observe numerous scattering resonances . to understand the nature of the scattering processes , we study the resonance structure of the scattering cross sections for individual partial waves . the resonances can be explained in terms of the bound states of the rb  nh@xmath0 complex . we calculate the pattern of bound states near the lowest dissociation limits of the complex as a function of the end - over - end angular momentum and identify the bound states responsible for the resonances . there is a further reason for being interested in rb  nh@xmath0 collisions . methods such as stark deceleration can slow molecules to velocities of a few metres per second , and the resulting molecules can then be trapped at temperatures of 10 to 100 mk @xcite . it has not yet proved possible to cool such molecules further , towards the temperatures and phase space densities at which they might undergo condensation to form quantum gases . there is great interest in methods that might be used to achieve this , and one of the most promising is _ sympathetic cooling _ @xcite , in which molecules are cooled by thermal contact with a laser - cooled gas of atoms such as rb . however , sympathetic cooling can work only if the collisions are predominantly elastic rather than inelastic : if the molecules are in excited states and undergo inelastic ( deexcitation ) collisions , the kinetic energy released is usually enough to eject both collision partners from the trap . we therefore also carry out calculations of elastic and inelastic cross sections at the temperatures relevant to sympathetic cooling . in previous work we obtained the potential energy surface for nh@xmath0 interacting with rb from highly correlated electronic structure calculations @xcite . spin - restricted coupled - cluster calculations including single , double and perturbative triple excitations [ rccsd(t ) ] were carried out on a grid of points in the intermolecular distance @xmath4 and intermolecular angles @xmath5 and @xmath6 . the angle @xmath6 corresponds to rotation of nh@xmath0 about its @xmath7 axis . the angles @xmath8 were chosen to be the points for 20-point gauss - lobatto quadrature , which include @xmath9 and @xmath10 ( with @xmath9 corresponding to approach toward the lone pair of nh@xmath0 ) . calculations were carried out on a grid of @xmath4 values from 3.5 to 25 @xmath11 and for two azimuthal angles @xmath12 and @xmath13 .    for collision calculations , we need to evaluate matrix elements of the interaction potential between angular momentum basis functions . this is most easily achieved by expanding the potential in renormalized spherical harmonics @xmath14 . after taking account of the @xmath15 symmetry of nh@xmath0 , the expansion takes the form @xmath16 . \label{pot_exp1}\end{aligned}\ ] ] the dominant terms in this expansion are those with @xmath17 and 3 , and the corresponding coefficients may be written @xmath18 where @xmath19 is an associated legendre polynomial and @xmath20 \nonumber \\ v_3(r,\theta ) & = & \textstyle{\frac{1}{2}}[v(r,\theta,0)- v(r,\theta,60^\circ ) ] . \label{v0v3}\end{aligned}\ ] ] @xmath21 can be viewed as the interaction potential averaged over @xmath6 , while @xmath22 describes the leading anisotropy of the potential with respect to rotation about the @xmath7 axis of nh@xmath0 . to obtain the expansion coefficients in eq.([pot_exp2 ] ) at a given value of @xmath4 for rb  nh@xmath0 , the potential functions @xmath23 at each grid point @xmath8 are first evaluated by reproducing kernel hilbert space ( rkhs ) interpolation @xcite and the angular integrations are then carried out by gauss - lobatto quadrature .    for calculations on rb  nd@xmath0 we need to reexpand the potential in a coordinate system based on the center of mass of nd@xmath0 instead of nh@xmath0 . this is done by first generating a set of points for rb  nd@xmath0 at the same grid of distances and angles as for rb  nh@xmath0 . the center of mass of nd@xmath0 is shifted by a distance @xmath24 towards the d atoms from that of nh@xmath0 . each point is obtained by as @xmath25 where @xmath26 and @xmath27 . the potential coefficients for rb  nd@xmath0 are then obtained by quadrature exactly as for rb  nh@xmath0 , but with the transformed set of points .    the @xmath21 and @xmath22 interaction potentials for rb  nh@xmath0 are shown in fig . [ potential ] . the potential is strongly anisotropic ; there is a very deep minimum , nearly 1900 @xmath1 deep , at the n side of nh@xmath0 . this is due to the interaction of the lone pair of nh@xmath0 with the singly occupied @xmath28 orbital of the rb atom . there is also a shallow secondary minimum , only about 100 @xmath1 deep , on the h side of nh@xmath0 , arising from dispersion interactions . both minima are on the @xmath7 axis of nh@xmath0 . a detailed discussion of the alkali - atom  nh@xmath0 interaction can be found in the ref . @xcite . the basis set which we used in our previous work ( augmented triple - zeta for n and h and a basis set of triple - zeta quality for rb ) gave depths for the minima at linear geometries of 1862 and 110 @xmath1 ( for n - side and h - side configurations , respectively ) . to improve this in the present work , we performed additional calculations for linear geometries near the minima with a significantly larger basis set . with the basis set extended to augmented quadruple - zeta for n and h , and with additional @xmath29 functions on the rb atom ( with exponents 0.45,0.167 ) , and an additional midbond @xmath30 function ( with exponent 0.4 ) , we obtained a global minimum well depth of 1881 @xmath1 with negligible ( less than 0.1% ) change in the equilibrium distance . all the potential points were then scaled by the ratio of the well depths calculated with quadruple - zeta and triple - zeta basis sets . the depth of the secondary minimum on the surface obtained by such scaling is 111 @xmath1 , which agrees well with the value obtained from _ ab initio _ calculations with the quadruple - zeta basis set . in the scattering calculations described below we have used the scaled triple - zeta potential .    for very large rb  nh@xmath0 separations the expanded potential is dominated by the @xmath31 and @xmath32 terms , which vanish as @xmath33 . since the inherent error of the supermolecular approach is quite significant at large @xmath4 , we represent the long - range behavior of @xmath31 and @xmath32 using @xmath34 and @xmath35 coefficients obtained from perturbation theory @xcite . at short range we switch over to the full expansion using a switching function @xmath36 \}^{-1}$ ] with @xmath37 and @xmath38 .     from ccsd(t ) calculations : @xmath39 component ( upper panel ) and @xmath40 component ( lower panel ) . contours are labeled in @xmath1 . to aid visualization , @xmath22 is plotted only in the energetically accessible region defined by @xmath41 . , title="fig : " ]   from ccsd(t ) calculations : @xmath39 component ( upper panel ) and @xmath40 component ( lower panel ) . contours are labeled in @xmath1 . to aid visualization , @xmath22 is plotted only in the energetically accessible region defined by @xmath41 . , title="fig : " ]      the hamiltonian that describes collisions between an atom and a molecule may be written @xmath42 + h_{\rm at } + h_{\rm mol } + v_{\rm inter},\ ] ] where @xmath43 is the reduced mass for the colliding pair , @xmath44 is the operator for rotation of the collision partners about one another , and @xmath45 and @xmath46 are the hamiltonians for the isolated atom and molecule respectively . in the special case where either the atom or the molecule is in a closed - shell singlet state , the main part of the interaction potential @xmath47 is a single potential energy surface @xmath48 that couples the molecular rotations strongly to the intermolecular distance . the molecular rotations are in turn coupled to the nuclear spins in the molecule by hyperfine terms . however , for an atom in an s state the only terms that couple the molecular degrees of freedom to the _ atomic _ electron and nuclear spins are ( i ) very weak magnetic dipole - dipole interactions between atomic and molecular spins and ( ii ) any dependence of the atomic hyperfine coupling on the intermolecular distance @xmath4 . if these very small terms are neglected , the collision problem may be treated as the scattering of the molecule from an _ unstructured _ atom . the hamiltonian of a rigid symmetric - top molecule may be written @xmath49 where @xmath50 is the operator representing the angular momentum of the molecule and @xmath51 denotes its projection onto the molecule - fixed symmetry axis @xmath52 , while @xmath53 and @xmath54 are the rotational constants for rotation about the symmetry axis and an axis perpendicular to it . the eigenfunctions of this hamiltonian are labeled by three quantum numbers : @xmath55 , the total angular momentum of the rotor ; @xmath56 , the projection of @xmath55 on the laboratory - frame @xmath57-axis and @xmath58 , the projection of @xmath55 on the body - fixed molecular symmetry axis @xmath52 , @xmath59    the n - h stretching vibrations are at sufficiently high frequency that they can be neglected in the present work , but nh@xmath0 also has a bending ( umbrella ) vibration with two minima that are equivalent in the free molecule . tunnelling between these two minima produces a low - frequency splitting that corresponds to inversion of the molecule through a planar geometry . to describe this we introduce an additional degree of freedom , the inversion coordinate @xmath29 . the vibration - inversion functions @xmath60 and @xmath61 ( corresponding to the lower and upper component of the inversion doublet ) may be written @xmath62,\ ] ] where @xmath63 is a function that is sharply peaked around @xmath64 , so that @xmath65 is a bending function that is sharply peaked around one of the equilibrium values @xmath66 ( corresponding to one of the individual `` umbrella '' states of nh@xmath0 ) . nh@xmath0 occurs in two different nuclear spin configurations , corresponding to @xmath67 and @xmath68 symmetries . in the first case ( referred to as ortho - nh@xmath0 ) the molecule can occupy rotational levels with @xmath69 , while in the latter ( referred to as para - nh@xmath0 ) only @xmath70 are allowed . for ortho - nh@xmath0 in @xmath71 states , only one component of each inversion doublet is allowed by symmetry . nd@xmath0 is slightly different and occurs in three nuclear spin symmetries , @xmath67 , @xmath72 and @xmath68 . molecules in nuclear spin states of @xmath67 and @xmath72 symmetries can occupy rotational levels with @xmath71 , 3 , 6 , etc . , referred to as ortho - nd@xmath0 , while those in the @xmath68 state can occupy levels with @xmath73 , 2 , 4 , 5 , etc . , referred to as para - nd@xmath0 . ortho- and para - species do not readily interconvert in collisions . by contrast with nh@xmath0 , however , nd@xmath0 molecules in @xmath71 states can exist in both upper and lower components of the inversion doublet . the rotational functions adapted for permutation - inversion symmetry may be written @xcite @xmath74 | \pm \rangle\ ] ] for para - nh@xmath0 and @xmath75 | \pm \rangle\ ] ] for para - nd@xmath0 . if @xmath76 and @xmath77 do not overlap significantly , evaluating of the matrix elements of the potential between the symmetry - adapted basis functions gives expressions that are isomorphic to those for the case of rigid symmetric top with parity - adapted monomer basis functions . it is therefore sufficient to carry out the scattering calculation exactly as for a rigid symmetric top , but with the monomer energies for the even and odd combinations of @xmath78 and @xmath79 shifted up and down from the rigid - rotor values by half the inversion splitting .    the approach of restricting the basis set of inversion functions to just the lowest pair of states was introduced by green @xcite and davis and boggs @xcite , but was later verified to be a good approximation by van der avoird and coworkers @xcite , who carried out calculations on ar  nh@xmath0 collisions to compare with theory that took the @xmath29 coordinate into account explicitly . the rotational constants and inversion parameters used for nh@xmath0 and nd@xmath0 in the present paper are given in table [ molpar ] . since inelastic transitions between inversion states are of key interest for the present work , we allow the inversion splitting to be a function of rotational state @xcite , @xmath80   - \nu_b k^2\ ] ] the energy levels for @xmath81 to 3 for both ortho- and para - nh@xmath0 and nd@xmath0 are given in table [ enlevels ] . in each case the zero of energy is the energy of the lowest allowed tunneling component for @xmath82 , @xmath71 .    in the following sections we identify tunneling components by @xmath83 ( upper ) or @xmath84 ( lower ) instead of symmetry labels , because the @xmath85 designation indicates the energetics more directly . .rotational constants @xcite and inversion splitting parameters ( with centrifugal distortions ) @xcite used for nh@xmath0 and nd@xmath0 . units are @xmath1 . [ cols="<,<,<",options="header " , ] the nd@xmath0 and nh@xmath0 molecules can be slowed by stark deceleration in their low - field - seeking @xmath81 , @xmath73 states , which correlate at low field with the upper level of the inversion doublet . since nd@xmath0 has a much smaller tunneling splitting than nh@xmath0 , its stark effect is more nearly linear and it is easier to decelerate . we have calculated the energy dependence of the state - to - state integral cross sections for molecules initially in the upper inversion state of nd@xmath0 and nh@xmath0 for collision energies between 0 and 100 @xmath1 . the results are shown in fig . [ totalcs ] . ( upper panel ) and rb  nd@xmath0 scattering ( lower panel ) from the upper component of the inversion doublet for @xmath81 , @xmath73 state . the smooth dashed lines show the result of the semiclassical background formula for elastic cross sections and of langevin capture theory for the total inelastic cross sections . , title="fig : " ]   ( upper panel ) and rb  nd@xmath0 scattering ( lower panel ) from the upper component of the inversion doublet for @xmath81 , @xmath73 state . the smooth dashed lines show the result of the semiclassical background formula for elastic cross sections and of langevin capture theory for the total inelastic cross sections . , title="fig : " ]    for the whole energy range considered , the elastic cross sections @xmath86 are large compared to the total inelastic cross sections for both nd@xmath0 and nh@xmath0 . fig.[totalcs ] includes the elastic cross section obtained from the semiclassical background formula @xcite for a pure @xmath33 potential ( dashed red line ) . for both rb  nh@xmath0 and rb  nd@xmath0 the elastic cross sections follow the background formula fairly closely , with slow glory oscillations superimposed on the background . the total cross sections ( elastic + inelastic ) are actually in even better agreement with the semiclassical model than the elastic cross sections . there is also structure due to scattering resonances at low collision energies , but it is not very strong with respect to the background for the elastic cross sections . at higher collision energies the resonances are lost in the background . the overall magnitude of the inelastic cross sections is at first sight less than expected . in a strongly coupled system where every collision that crosses the centrifugal barrier leads to inelasticity , the total inelastic cross section is given by the langevin capture formula @xcite , @xmath87 the langevin result is shown as a dashed green line in fig.[totalcs ] and it may be seen that the actual @xmath88 inelastic cross section is substantially below it . this is surprising in view of the very large anisotropy which directly couples the @xmath89 and @xmath90 states . further insight is given by the partial - wave contributions to the @xmath88 cross sections , which are shown for collision energies of 0.6 @xmath1 and 10 @xmath1 in fig . [ opacity ] . in this case the capture cross section corresponds to an elastic s - matrix element of zero for every @xmath91 below a cutoff due to the centrifugal potential . the partial inelastic cross section is far less than the capture value across most of the range of @xmath91 . this contrasts with the behavior observed in other strongly coupled systems , such as high - energy hf  hf collisions @xcite , rotationally inelastic collisions in ar  n@xmath3 at 300  k @xcite , and barrierless reactions in alkali metal atom + diatom collisions at energies above 1  mk @xcite , where the partial cross sections _ are _ close to the capture value for all @xmath91 . cross sections for rb  nd@xmath0 at collision energies of 0.6 @xmath1 ( upper panel ) and 10 @xmath1 ( lower panel ) . the straight dashed lines show the result corresponding to langevin capture theory , with a cutoff determined by the heights of the centrifugal barriers calculated for the full potential.,title="fig : " ]   cross sections for rb  nd@xmath0 at collision energies of 0.6 @xmath1 ( upper panel ) and 10 @xmath1 ( lower panel ) . the straight dashed lines show the result corresponding to langevin capture theory , with a cutoff determined by the heights of the centrifugal barriers calculated for the full potential.,title="fig : " ]    a possible explanation for the low inelasticity for rb  nd@xmath0 is that , despite the large anisotropy , the collisions are approximately adiabatic in the rotation - inversion coordinates . fig.[adiabats ] shows  adiabatic bender " curves for rb  nd@xmath0 , which are eigenvalues of the rotation - inversion hamiltonian for para - nd@xmath0 at fixed @xmath4 . it may be seen that the curves correlating with the @xmath90 and @xmath89 states stay far apart for all values of @xmath4 , with no avoided crossings between them where strong nonadiabatic transitions would be expected .    . the line types indicate the level that each curve correlates with at long range and not the character of the adiabatic state at short range . the cross indicates the position of the absolute minimum of the potential surface and the horizonal line indicates the energy of the lowest bound state , @xmath92 @xmath1 . ] the resonances are much stronger in the inelastic cross sections than the elastic cross sections . they are particularly strong for collision energies up to about 20 @xmath1 , as can be seen in the expanded view of the energy dependence of @xmath93 in fig . [ totalcszoom ] . for rb  nd@xmath0 at low energies the resonances can enhance inelastic cross sections by up to a factor of 2 . as the collision energy increases the resonances again become weaker and wash out . we see significantly more resonances for rb  nd@xmath0 than for rb  nh@xmath0 . a detailed description of the origin of the resonances for low kinetic energies will be given in section iii.b . cross sections for rb  nh@xmath0 and rb  nd@xmath0 at low collision energies.,title="fig : " ]   cross sections for rb  nh@xmath0 and rb  nd@xmath0 at low collision energies.,title="fig : " ]    the coupled states approximation is not expected to be accurate at very low energies , because it approximates the centrifugal terms in the potential and neglects couplings between different helicities . unfortunately , full close - coupling calculations for higher energies are prohibitively expensive for partial waves corresponding to large values of @xmath94 . for example , the basis set used above ( @xmath95 and @xmath96 ) , which gives 210 channels for para - nd@xmath0 in the cs approximation , produces 2745 channels for each parity in cc calculations for @xmath97 . however , it is possible to make the comparison for a restricted basis set , and the results are shown in fig.[ccvscs ] for a basis set with @xmath98 , @xmath99 . the terms that are neglected in the cs approximation cause small shifts in resonance positions and intensities , so that it is not possible to make a valid comparison at any individual collision energy . however , the general extent of the resonance structure is similar in the two calculations and the inelastic cross sections typically agree within 5% above 30 @xmath1 . we can therefore have confidence in the general features of the cs calculations for the full basis set . ( upper panel ) and inelastic @xmath100 ( lower panel ) cross sections for rb  nd@xmath0 calculated with cs and cc methods for a restricted basis set with @xmath98 , @xmath101.,title="fig : " ]   ( upper panel ) and inelastic @xmath100 ( lower panel ) cross sections for rb  nd@xmath0 calculated with cs and cc methods for a restricted basis set with @xmath98 , @xmath101.,title="fig : " ]    the state - to - state cross sections from the @xmath90 state to other @xmath102 states are shown as a function of collision energy in fig .  [ inelcs ] . as the transition to each additional rotational state becomes energetically allowed , the corresponding excitation cross section rises sharply from zero . this behavior could be seen readily in experiments that use state - selective detection for the scattered nd@xmath0 . however , since the cross sections for rotational excitation are very small compared to those for @xmath100 relaxation at these energies , we do not see any significant changes in the total inelastic cross section as new channels open up .     state of rb  nh@xmath0 and rb  nd@xmath0 as a function of collision energy , from coupled - states calculations.,title="fig : " ]   state of rb  nh@xmath0 and rb  nd@xmath0 as a function of collision energy , from coupled - states calculations.,title="fig : " ]    fig .  [ wigner_open ] shows the threshold behavior in more detail , for the contributions to the cross section @xmath103 obtained from cc calculations for different partial waves @xmath94 with the full basis set . cc calculations are feasible in this case because only low partial waves contribute just above threshold . at very low energy , each partial - wave contribution follows the wigner threshold law @xcite and is proportional to @xmath104 , where @xmath105 and @xmath106 are the collision energy and the threshold energy respectively . @xmath91 takes the lowest allowed value for the threshold level : for @xmath107 , @xmath108 for @xmath109 , @xmath110 for @xmath111 and 3 and @xmath112 for @xmath113 and 4 . for any final @xmath55 there is always one total angular momentum ( @xmath114 ) and parity combination that allows @xmath108 in the outgoing channel , and this dominates the integral cross section for excess energies below @xmath115 @xmath1 . above this , however , additional partial waves contribute and simple @xmath116 behavior is not expected for the integral cross section . inelastic excitation cross section , from close - coupling calculations , as a function of the excess energy above the threshold for excitation . ]    for rb  nh@xmath0 , the @xmath117 transition is much stronger than any other inelastic transition for the whole range of collision energies considered here . for rb  nd@xmath0 , @xmath103 becomes comparable to @xmath118 for kinetic energies around 100 @xmath1 . the resonances in individual state - to - state cross sections are in some cases quite strong compared to their backgrounds , especially closely above the thresholds . however , since their magnitude is much smaller than @xmath93 , these resonances are not visible in the total inelastic cross section . there are clear propensity rules that apply to the cross sections , as observed previously for he  nh@xmath0 @xcite and ar  nh@xmath0 @xcite collisions . the general effect has been discussed by alexander @xcite . for @xmath58-conserving collisions , inversion - changing transitions are preferred when @xmath119 is odd and inversion - conserving transitions are preferred when @xmath119 is even . thus @xmath120 is larger than @xmath121 , since the @xmath90 and @xmath122 channels are directly coupled by @xmath123 potential term . conversely , for collisions with @xmath124 , which are driven by the @xmath125 potential term , inversion - conserving transitions are preferred when @xmath119 is odd , whereas inversion - changing transitions are preferred when @xmath119 is even . thus @xmath126 is much larger than @xmath127 . it is interesting that these propensity rules survive in a system as strongly coupled as rb  nh@xmath0 . there are in principle two types of resonance that might be seen in low - energy collisions . shape resonances correspond to quasibound states that are confined behind a centrifugal barrier , whereas feshbach resonances correspond to quasibound states that reside principally in channels corresponding to internally excited states . the height of the centrifugal barrier for each @xmath91 is given approximately by @xcite @xmath128^{3/2 } ( 54c_6)^{-1/2}. \label{eq : vl}\ ] ] this function is compared with the actual barrier heights obtained from the adiabatic bender curves for rb  nd@xmath0 in fig .  [ barriers ] ; it may be seen that the long - range formula ( [ eq : vl ] ) is accurate for low @xmath91 but overestimates the barrier height by about 10% by @xmath129 . however , since rb  nd@xmath0 has a strongly attractive long - range potential ( @xmath130 ) and large reduced mass , the heights of the centrifugal barriers are much smaller than in lighter systems . as will be seen below , most of the resonances observed in the present work are feshbach resonances . predicted by long - range formula ( [ eq : vl ] ) and obtained from the full potential . ] ( top panel ) and @xmath131 ( center panel ) to the inelastic @xmath132 cross section for rb  nd@xmath0 . the resonances are labeled by roman numerals corresponding to the series of bound states shown in fig . the bottom panel shows eigenphase sums for @xmath131 . the dashed vertical lines in the bottom panel indicate the positions of centrifugal barriers.,title="fig : " ]   ( top panel ) and @xmath131 ( center panel ) to the inelastic @xmath132 cross section for rb  nd@xmath0 . the resonances are labeled by roman numerals corresponding to the series of bound states shown in fig . the bottom panel shows eigenphase sums for @xmath131 . the dashed vertical lines in the bottom panel indicate the positions of centrifugal barriers.,title="fig : " ]   ( top panel ) and @xmath131 ( center panel ) to the inelastic @xmath132 cross section for rb  nd@xmath0 . the resonances are labeled by roman numerals corresponding to the series of bound states shown in fig . the bottom panel shows eigenphase sums for @xmath131 . the dashed vertical lines in the bottom panel indicate the positions of centrifugal barriers.,title="fig : " ]     ( full circles ) for helicity @xmath133 as a function of @xmath134 . open circles indicate the positions of the strongest resonances in fig . the arrow points out an example of an avoided crossing between bound states . the horizontal solid lines at 8.257 and 8.310 @xmath1 show the @xmath89 and @xmath90 thresholds , respectively . the dashed line shows the position of centrifugal barrier . ]    the origin of the resonances can be understood if we study the contribution to the inelastic cross sections from individual terms in the partial - wave expansion . the top two panels of fig . [ jexp ] show the @xmath100 partial cross sections for @xmath135 and @xmath131 for rb  nd@xmath0 collision energies up to 2 @xmath1 , where there are very strong resonances . the bottom panel shows the s - matrix eigenphase sums @xcite for @xmath131 , which show a sharp rise through @xmath136 as the energy passes through a resonance .    in general each partial inelastic cross section shows a non - resonant peak near the corresponding barrier maximum and then dies off at higher energies . the non - resonant peak is quite sharp for low @xmath91 but broadens as @xmath91 increases . in particular , the peaks below @xmath137 @xmath1 for @xmath138 are non - resonant . superimposed on the non - resonant background are peaks and troughs due to resonances . below the barrier maximum , both shape and feshbach resonances may occur ; however peaks due to feshbach resonances are suppressed when they are below barriers and it is usually shape resonances that give large peaks in this region .    above the barrier maximum for each @xmath91 , langevin capture theory would predict a partial inelastic cross section @xmath139 except at resonances and occasionally near the barrier maximum , the partial cross sections are generally much smaller than this . it is thus evident that low - energy inelastic scattering in rb  nd@xmath0 is _ dominated _ by resonant effects . we also carried out bound - state calculations close to the @xmath2 thresholds to help identify the states responsible for the resonances . the bound - state calculations used the coupled - states approximation with the same basis set as for the collision calculations . the helicity was restricted to @xmath133 since this is the only value that contributes to the @xmath140 cross section . the resulting bound states are shown in fig . [ bnd ] . since the rotational energy of the bound and quasibound states of the rb  nd@xmath0 complex is approximately proportional to @xmath134 , it is convenient to plot the energies as a function of this quantity . this allows rotational progressions to be identified as nearly straight lines , and rotational constants can be extracted from the slopes of the lines . however , if there are two or more bound states for the same @xmath91 with energies very close together , they can mix and avoid one another ; an example of this is shown by the dotted lines in fig . [ bnd ] .    most of the peaks in fig . [ jexp ] can be identified with series of bound states . for example , the cross sections for @xmath141 shows three strong peaks . the first of these , at about 0.5 @xmath1 , is a nonresonant peak associated with the barrier maximum . at higher energies are two resonant peaks , labeled as i , and ii in fig . they are associated with characteristic features in the eigenphase sum and can be identified as feshbach resonances associated with the series of bound states along the dashed lines i and ii in fig . [ bnd ] ; series ii also explains the strong peak in the cross section for @xmath142 near 0.6 @xmath1 . lines i and ii both have a large slope , corresponding to a rotational constant @xmath143 @xmath1 and an effective distance between rb and nd@xmath0 around @xmath144 . it is also not difficult to assign the resonances labeled by iv , v , vi and vii to the appropriate series of bound states ; all these are feshbach resonances that arise from bound states with small effective intermolecular distances @xmath145 . for all these bound states , inspection of the wavefunction near the van der waals minimum confirms that they are dominated by rotationally excited basis functions with @xmath146 . the levels along the line iii in fig . [ bnd ] , which give strong resonances for @xmath147 and 9 , are associated with quasibound states of significantly smaller rotational constant , @xmath148 @xmath1 , and an effective intermolecular distance around 20 @xmath11 . the dominant contributions to the wavefunction in this case are the two @xmath2 states ; the @xmath83 and @xmath84 states contribute almost equally and with opposite signs , corresponding to a single umbrella state of non - inverting nd@xmath0 . this suggests that these are shape resonances , and indeed for @xmath147 and 9 they lie below the energy of the centrifugal barrier , shown as a dashed curve in fig .  [ bnd ] . for @xmath149 the extrapolated energy of line iii is above the barrier maximum and no well - defined shape resonances occur . it is important to appreciate that the absolute energies of the resonances and bound states close to threshold are strongly dependent on the potential energy surface . thus , the studies reported in this section should be treated as a paradigm for understanding the resonant behavior of cross sections , rather than as predictions of the actual energies at which resonances will appear for rb  nd@xmath0 . as described in the introduction , there is great interest in the possibility of _ sympathetic cooling _ , in which molecules are cooled by thermal contact with a gas of laser - cooled atoms such as rb . the simplest sympathetic cooling experiment would overlap an electrostatic trap containing a sample of cold molecules such as nd@xmath0 with a gas of atoms in a magnetic or magneto - optical trap . it is therefore of considerable interest to explore collision cross sections for promising candidate atoms and molecules in the temperature regime between 1 @xmath43k and 100 mk . an electrostatic trap works only for molecules in low - field - seeking states , and for nh@xmath0 and nd@xmath0 the low - field - seeking states correlate at low field with the upper state of the tunneling doublet . it is thus very important to know whether collisions with ultracold rb will cause relaxation to the lower tunneling state , which is a high - field - seeking state and can not be trapped electrostatically . more specifically , we need to know whether the cross section for the relaxation from @xmath90 to @xmath89 state is sufficiently small compared to the elastic cross section .    to explore this we have performed close - coupling calculations for very low energies using the same rotational basis set as for the cs calculations described above . in fig . [ ultracold ] we show the elastic and total inelastic integral cross sections calculated for partial waves @xmath113 to 8 , which gives reasonable convergence for collision energies up to 100 mk . the contribution to the total cross sections from @xmath150-wave scattering ( @xmath111 partial wave ) is dominant in the @xmath43k regime . the wigner threshold behavior for @xmath150-wave scattering is followed for kinetic energies up to about 100@xmath43k . elastic cross sections are larger than inelastic cross sections for collision energies above 100 @xmath43k , but never by much more than a factor of 10 . for sympathetic cooling , the commonly stated rule of thumb is that the ratio of elastic to inelastic cross sections should be at least 100 . thus it appears that the inelastic cross sections are too large to allow nd@xmath0 molecules in low - field - seeking states to be cooled to sub - mk temperatures by collision with ultracold rb atoms .     and inelastic @xmath100 cross sections and the contribution from @xmath150-wave scattering ( @xmath111 ) , from close - coupling calculations . ] there remains the possibility of using sympathetic cooling for molecules in high - field - seeking states , which can be confined using an alternating current trap @xcite . confinement of rb atoms in such a trap has also been demonstrated @xcite , but the ac frequencies required are quite different in the two cases . it is therefore important to know whether sympathetic cooling of nd@xmath0 molecules in high - field - seeking states might be feasible using magnetically trapped rb atoms , which are themselves in low - field - seeking states that are not the ground state in the applied magnetic field . in previous work on rb  oh collisions @xcite , we found large cross sections for low - energy collisions that changed the hyperfine state of rb . we initially anticipated @xcite that this would occur for rb  nd@xmath0 collisions as well . however , the more thorough analysis of the collision hamiltonian in section ii.b above suggests that for molecules in closed - shell singlet states the atomic spins are likely to be almost unaffected by collisions . we therefore now consider that sympathetic cooling of high - field - seeking states of nd@xmath0 ( or nh@xmath0 ) by magnetically trapped rb ( or another laser - cooled atomic gas ) has a good prospect of success . collision calculations at ultralow kinetic energies are very sensitive to details of the potential used in the calculations . quantitative theoretical predictions of parameters related to scattering cross sections , such as the highest bound states or scattering lengths , are possible only for the simplest , lightest systems such as metastable helium @xcite or hydrogen - lithium mixtures @xcite . the calculated rb  nh@xmath0 interaction potential we use in this paper suffers from many uncertainties resulting from incompleteness of the electronic basis set , the approximate treatment of electronic correlation , relativistic effects , etc . , which limit its accuracy to at best a few percent .    to explore this sensitivity , we introduce an additional scaling factor @xmath151 into the definition of the interaction potential . we varied the scaling factor between 0.90 and 1.05 and explored the variation of the @xmath150-wave elastic @xmath86 and inelastic @xmath93 cross sections ( using close - coupling calculations ) for a collision energy of 10 @xmath43k . the results for a representative slice of the range of @xmath151 studied are shown in fig . [ scalingparam ] . the elastic and inelastic cross sections vary dramatically as bound states cross threshold and appear as scattering resonances . in the range of @xmath151 explored , the elastic and inelastic cross sections pass through more than ten resonances . there are some values of @xmath151 , such as near @xmath152 , where the inelastic cross section @xmath93 is significantly suppressed by the presence of a strong resonance , but for most other values of the scaling parameter the elastic - to - inelastic ratio ranges lies between 10 and 0.1 . although it is possible that the real potential might produce inelastic cross sections low enough to allow sympathetic cooling for low - field - seeking states , it is quite unlikely .     and inelastic @xmath100 cross sections for rb  nd@xmath0 at a collision energy of 10 @xmath43k to a scaling factor @xmath151 applied to the interaction potential . ] we have studied collisions between rb atoms and nh@xmath0/nd@xmath0 molecules , motivated by recent progress @xcite in collision experiments involving velocity - controlled and stark - decelerated beams of molecules and by interest in the possibility of sympathetic cooling . we focused principally on nh@xmath0 or nd@xmath0 molecules initially in the upper component of the inversion doublet for @xmath81 , @xmath73 . this level correlates with the low - field - seeking state in an electric field , which can be slowed by stark deceleration and trapped in an electrostatic trap . we considered collision energies between 0 and 100 @xmath1 . using the coupled - states approximation , we calculated the elastic cross sections and state - to - state inelastic cross sections from the low - field - seeking @xmath2 state to other rotation - inversion levels . the inelastic cross sections are smaller than expected for such a strongly coupled system , but are still only about a factor of 10 smaller than the elastic cross section over most of the energy range considered . both the elastic and inelastic cross sections show dense structure due to scattering resonances . the resonances in the elastic cross sections are diffuse and rather weak compared to the background . the total inelastic cross sections have much stronger resonances compared to their background , especially at collision energies below about 20 @xmath1 and one can consider the inelastic scattering as _ mostly _ resonant in nature . these resonances are washed out for larger collision energies . for energies below 90 @xmath1 the @xmath100 inelastic cross section makes the largest contribution to the total inelastic cross section . transitions to @xmath153 and 3 do not change the total inelastic cross section significantly as the new channels become energetically accessible at higher collision energies . we have considered the origin of the scattering resonances , using calculations of the bound states of rb  nd@xmath0 near threshold . since the long - range attraction between rb and nd@xmath0 is very strong and the reduced mass is fairly large , the resulting centrifugal barriers are much smaller than in lighter systems . most of the resonances appearing in the cross sections arise from quasibound states in which rb interacts with rotationally excited nh@xmath0 , and can be classified as feshbach resonances . finally we studied the ultracold limit of the elastic and inelastic cross sections , in order to test whether low - field - seeking molecules could be cooled to microkelvin temperatures by contact with an ultracold gas of rb atoms . hyperfine effects were neglected . we found that inelastic collisions are strong , and the elastic / inelastic ratio is unlikely to be sufficient to achieve sympathetic cooling of nd@xmath0 molecules in the @xmath90 state . however , there is a good prospect that sympathetic cooling will be possible for high - field - seeking states of nd@xmath0 or nh@xmath0 , even if the atoms used as a coolant are in low - field - seeking states . the authors are grateful to heather lewandowski for sparking their interest in rb  nh@xmath0 collisions and to ruth lesueur for discussions on the effect of the ammonia inversion on the collisions . the authors are grateful to epsrc for funding of the collaborative project copomol under the esf eurocores programme euroquam .
Abstract:

Answer: we carry out quantum inelastic scattering calculations of collisions of rb atoms with inverting nh@xmath0 and nd@xmath0 molecules in the energy range between 0 and 100 @xmath1 , which are important for experiments using velocity - controlled molecular beams to probe scattering resonances . we focus on molecules initially in the upper level of the ammonia inversion doublet for @xmath2 , which is low - field - seeking and can be controlled in a stark decelerator . we calculate the integral elastic and state - to - state inelastic cross sections in the coupled states approximation . we demonstrate the presence of both shape and feshbach resonances in the elastic and inelastic cross sections at low collision energies and discuss their origin in terms of the bound states of rb nd@xmath0 complex . we also consider elastic and inelastic cross sections in the ultracold regime , using close - coupling calculations , in order to assess the viability of sympathetic cooling of nd@xmath0 by rb . the inelastic cross section for relaxation to the lower level of the inversion doublet is smaller than expected for such a strongly coupled system , but is still likely to be too large to allow sympathetic cooling for nd@xmath0 in low - field - seeking states . however , there is a good prospect that sympathetic cooling will be possible for molecules in high - field - seeking states , even when the collision partner is a magnetically trapped atom in a low - field - seeking state .


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: compton scattering in quasi - spherical accretion flows is thought to be responsible for the x - ray spectra of many accreting compact objects . for example , the hard x - ray spectra of weakly - magnetic accreting neutron stars can be accounted for if a geometrically thick , hot medium surrounds the neutron star and inner accretion disk ( see , e.g. , psaltis , lamb , & miller 1995 ) . the spectral signature of isolated neutron stars accreting from the interstellar medium , which is required for their positive identification , is determined mainly by the efficiency of compton upscattering of thermal photons emitted from the stellar surfaces ( see , e.g. , zampieri et al . moreover , compton scattering plays a major role in the formation of the high - energy spectra of advection - dominated accretion flows , which may be the relevant mode of accretion onto galactic and supermassive black holes ( see , e.g. , esin et al . 1997 ) . compton scattering in static and moving media , even in plane - parallel or spherical symmetry , couples strongly all three independent phase - space coordinates of the photons , namely their energy , direction of propagation , and spatial coordinate . for this reason , monte carlo methods have been used extensively in calculating comptonization spectra , because of their easy implementation and natural treatment of the diverse length scales in an accretion flow ( see , e.g. , laurent & titarchuk 1999 for a recent study ; methods and results based on monte carlo treatments have been reviewed by pozdnyakov , sobol , & sunyaev 1983 and will not be discussed here ) . alternatively , compton scattering can be studied by solving the energy - dependent radiative transfer equation or its moments . implementing such a method is more difficult and often relies on a number of approximations for describing the interaction of radiation with matter ( see , however , poutanen & svensson 1996 and hsu & blaes 1998 for complete solutions in static media ) . however , at the same time it is substantially faster and offers significant advantages when the radiative transfer problem is coupled to a solution of the fluid dynamics , as well as in treating steep spectra and processes that do not conserve photons .    over the last twenty five years several authors have derived the photon kinetic or radiative transfer equation and their moments for compton scattering in static and moving media ( see , e.g. , kompaneets 1957 ; babuel - peyrissac & rouvillois 1969 ; pomraning 1973 ; chan & jones 1975 ; payne 1980 ; blandford & payne 1981a ; thorne 1981 ; fukue , kato , & matsumoto 1985 ; madej 1989 , 1991 ; titarchuk 1994 ) . the moment equations derived by these various sets of authors have been widely used in studies of comptonization by static media ( see , e.g. , katz 1976 ; shapiro , lightman , & eardly 1976 ; sunyaev & titarchuk 1980 ) or by strong shocks and accretion flows onto compact objects ( see , e.g. , blandford & payne 1981b ; payne & blandford 1981 ; lyubarskij & sunyaev 1982 ; colpi 1988 ; riffert 1988 ; mastichiadis & kylafis 1992 ; titarchuk & lyubarskij 1995 ; turolla et al . 1996 ; titarchuk , mastichiadis , & kylafis 1997 ) , either under the diffusion approximation or with closure relations that were specified a priori . most of the above analyses of compton scattering in moving media were performed in the system frame and to first order in the electron bulk velocity . on the other hand , schmid - burgk ( 1978 ) , zane et al . ( 1996 ) , and titarchuk & zannias ( 1998 ) have solved the general - relativistic radiative transfer equation for compton scattering without the need for specifying a priori closure relations .    in paper i ( psaltis & lamb 1997 ) of this series , the time - dependent photon kinetic and radiative transfer equations were derived , as well as their zeroth and first moments that describe absorption , emission , and spontaneous and induced compton scattering in static and moving media , correcting various errors in the literature . the system - frame equations that were derived are valid to first order in @xmath2 and @xmath3 , and to _ second _ order in @xmath1 , where @xmath4 and @xmath5 are the electron mass and temperature , @xmath6 is the photon energy , and @xmath1 is the bulk velocity of the electrons in units of the speed of light ; the fluid - frame equations that were derived are valid for arbitrary values of the bulk velocity @xmath1 . using these equations it was argued in paper i that the effects of comptonization by the bulk electron velocity that are described by the terms that are second - order in @xmath1 are qualitatively different than and can become at least as important as the effects described by the terms that are first - order in @xmath1 , even when @xmath1 is small . as a result , these second - order terms should generally be retained ( see also yin & miller 1995 ) .    in this paper , i describe a numerical algorithm for the solution of time - independent radiative transfer problems in systems with spherical symmetry . even though the natural reference frame for solving transport problems is the one comoving with the flow , i work in the system frame , which can be any inertial frame in which the central object is at rest . this is important for the fast algorithm described below , because the moment equations in the system frame are always parabolic ( as shown in  3 ) in contrast to the fluid - frame equations ( turolla , zampieri , & nobili 1995 ; smit , cernohorsky , & dullemond 1997 ; see also krner & janka 1992 ) . moreover , in problems of mass accretion onto a compact star , the boundary conditions are often more easily implemented in the system frame . i also neglect any general relativistic effects , which can be shown to introduce only small quantitative corrections to the emerging spectra ( papathanassiou & psaltis 2000 ) , and truncate the transfer equation keeping only terms up to second order in @xmath1 . the algorithm used is a generalization of the method of variable eddington factors ( mihalas 1980 ; see also mihalas 1978 , p.  157158 ) . it is therefore accurate for systems of arbitrary optical depth and does not require a priori specification of closure relations . it is based on the iterative solution of both the transfer equation and the systems of its first two moments and hence the validity and accuracy of the solution can be verified explicitly at the end of each run . the method also requires only a small number of iterations in order to converge to the solution and is therefore ideal for future extensions to time - dependent and multi - dimensional transport problems .    in  2 , i summarize the results from paper i and introduce the notation . in  3 , i discuss the mathematical character of the system of equations and in  4 , i present the numerical algorithm . in  5 , i describe the results of numerical solutions of the transfer equations for a variety of situation . finally , in  6 , i briefly summarize the implications of the above for comptonization spectra in static and moving media . _ units_throughout this paper i set @xmath7 , where @xmath8 is planck s constant , @xmath9 is boltzmann s constant , and @xmath10 is the speed of light . i also normalize all the spatial coordinates to the inner radius of the flow , the electron temperature @xmath5 and photon energy @xmath6 to the electron rest mass @xmath4 , and the electron density @xmath11 to its value at the inner radius of the flow . finally , i normalize the emission and absorption coefficients @xmath12 and @xmath13 to the inverse of the electron scattering mean - free path in the thomson limit . _ the electron gas_let @xmath14 be the system - frame three - velocity of a given electron . i define the local bulk velocity @xmath15 of the electrons measured in the system frame as @xmath16 where the sharp brackets indicate the average over the local electron velocity distribution . i also assume that the electron velocity distribution in the fluid frame is a relativistic maxwellian and therefore @xmath17 where i have neglected terms of order @xmath18 and higher . _ the radiation field_i describe the radiation field at any position in the system frame with coordinate vector @xmath19 using the monochromatic specific intensity @xmath20 , where @xmath21 is the direction of propagation and @xmath6 is the photon energy . here i consider only unpolarized radiation and hence i have suppressed the dependence of the specific intensity on polarization . because of the assumed spherical symmetry of the problem , the radiation field depends only on the radial distance from the center of symmetry ( i.e. , the @xmath22component of the vector @xmath19 ) , on the angle @xmath23 between the radial direction and the direction of propagation ( i.e. , @xmath24 ) , and on the photon energy @xmath6 . i define the first five moments of the monochromatic specific intensity as @xmath25 the non - zero components of the first five moments of the specific intensity for systems with spherical symmetry are given in appendix  a.    _ the transfer equation and its moments_following mihalas ( 1980 ; see also mihalas 1978 ) , i solve the radiative transfer equation in a plane that contains the center of symmetry of the system , between the inner and outer radii @xmath26 and @xmath27 of the flow , using the coordinates @xmath28 and @xmath29 . the radiative transfer equation along a ray of constant impact parameter @xmath30 is @xmath31 where the signs ` @xmath32 ' and @xmath33 correspond to the equations for the outgoing and incoming rays respectively and i have neglected the effects of induced compton scattering . in equation  ( [ rteq ] ) i have used the optical depth along the ray defined as @xmath34 dz$ ] , where @xmath35 at the outer boundary and @xmath36 is the angle - integrated thomson scattering cross section . @xmath37 is the generalized source function for absorption , emission , and compton scattering in systems with spherical symmetry and is given in appendix  b.    the equations for the moments of the specific intensity can be derived either directly from equation ( [ rteq ] ) or from equations ( 34 ) and ( 40 ) of paper i. written in compact form and using the dimensionless quantities defined above , the moment equations are term reported in paper i ( psaltis & lamb 1997 ) contained a typo that has been corrected here . ] @xmath38 where @xmath39 , @xmath40 , @xmath41 , and @xmath42 . when the electron density is uniform , then the dimensionless quantity @xmath43 is equal to the electron scattering optical depth in the radial direction measured from @xmath27 . the coefficients in equations ( [ 0th3 ] ) and ( [ 1st ] ) are given in appendix  c.    _ boundary conditions_the boundary conditions for both the transfer equation and the system of its moments depend on the specific problem under study . for the model problems discussed in  5 , i shall assume that the flow is not illuminated from the outside , i.e. , @xmath44 at the outer boundary , and specify the radiation flux at the inner boundary as @xmath45    for solving the moments of the radiative transfer equation i shall specify , for all photon energies @xmath6 , the first moment of the monochromatic specific intensity at the inner radius , i.e. , @xmath46 as well as the ratio of the first to the zeroth moment of the monochromatic specific intensity at the outer radius @xmath47 the quantity @xmath48 is not known a priori but is calculated in a self - consistent way in the iterative method described below . i shall also set at all radii @xmath49 and @xmath50 which reflect the fact that the photon occupation number and the radiation flux are finite quantities . in the method of variable eddington factors , the system of partial differential equations ( [ 0th3 ] ) and ( [ 1st ] ) is solved iteratively with the radiative transfer equation  ( [ rteq ] ) , given a set of eddington factors computed from the previous iteration . the boundary conditions and the method used for the solution of the system depends on the character of the partial differential operator . for example , the system of moment equations in the frame _ comoving _ with the flow ( or the local lorentz frame when gravitational effects are taken into account ) may be hyperbolic or elliptic , depending on the velocity field ( turolla et al . 1995 ) . in this section i follow the procedure outlined by ames ( 1992 , pp . 812 ) to show that the _ system - frame _ , energy dependent moment equations always form a parabolic system of partial differential equations , thus simplifying the implementation of the numerical algorithm . i define the quantities @xmath51 and @xmath52 , which together with the partial differential equations ( [ 0th3 ] ) and ( [ 1st ] ) , and the equations for the differentials @xmath53 , @xmath54 , @xmath55 , and @xmath56 , form an @xmath57 system of algebraic equations for the eight derivatives @xmath58 \cdot \left [ \begin{array}{c } \partial_\e j\\ \partial_{\tau_r } j\\ \partial_\e h^r\\ \partial_{\tau_r } h^r\\ \partial_\e l\\ \partial_{\tau_r } l\\ \partial_\e m^r\\ \partial_{\tau_r } m^r\\ \end{array } \right ] & = & \qquad\qquad\qquad\qquad\nonumber\end{aligned}\ ] ] @xmath59\;.\ ] ] equating the determinant of the coefficient matrix to zero , i obtain the characteristic equation for the system of equations , @xmath60 the coefficients @xmath61 and @xmath62 are proportional to the electron temperature and the square of the bulk velocity ( see appendix  c ) and hence for the problems under consideration they are never zero . therefore , for a finite electron energy @xmath6 , the system of partial differential equations ( [ 0th3 ] ) and ( [ 1st ] ) has only one multiple , real characteristic direction , i.e. , @xmath63 , and hence is always parabolic . the difference with respect to the moment equations in the comoving frame arises from the presence of second order energy derivatives in the differential operator . equation ( [ rteq ] ) is an integro - differential equation for the specific intensity of the radiation field and therefore any numerical solution of this equation is not trivial to obtain . on the other hand , the system of equations ( [ 0th3 ] ) and ( [ 1st ] ) depends on two variable eddington factors and an unknown outer boundary condition , which can be computed only when the solution of equation ( [ rteq ] ) for the specific intensity is obtained . here , i employ a numerical algorithm that is a generalization of the method of variable eddington factors ( mihalas 1980 ; see also mihalas 1978 , pp . 157158 ) in order to solve iteratively equation ( [ rteq ] ) for the specific intensity of the radiation field and the system of partial differential equations ( [ 0th3 ] ) and ( [ 1st ] ) for its zeroth and first moments . this algorithm has been proven to be very efficient and gives a solution to the complete transfer equation that is limited only by numerical accuracy , independent of the optical depth of the medium .    in this method , i first solve the moments of the transfer equation using an initial guess for the variable eddington factors and for the outer boundary condition . i then compute the generalized source function using the calculated moments of the specific intensity , solve the radiative transfer equation , and use this solution to update the eddington factors . i repeat the whole procedure until convergence is achieved . the moments of the transfer equation depend only on radius and photon energy , when the variable eddington factors are known . therefore , i discretize the domain of solution and the differential operators of equations ( [ 0th3 ] ) and ( [ 1st ] ) on a two  dimensional mesh of @xmath64 grid points over the variable @xmath43 and @xmath65 grid points over the photon energy @xmath6 .    the mesh of discrete grid points over the variable @xmath43 must resolve the rapid change of the variable eddington factors with optical depth near both boundaries of the domain of solution . in order to achieve this , i define the mesh points to be equidistant in the quantity @xmath66 where @xmath67 is the total optical depth ( and not the quantity @xmath43 ) in the radial direction measured from @xmath27 , @xmath68 is the maximum radial optical depth in the medium , and @xmath69 is a parameter that allows a combination of a logarithmic grid near the boundaries with a linear grid in the interior . because the spectra that emerge from comptonizing media are usually power laws , i use a logarithmic mesh of discrete grid points over photon energy @xmath6 . the system of moment equations involve differentiation over two density - like quantities , the radiation energy density @xmath70 and the first eddington factor @xmath71 , and two flux - like quantities , the radiation energy flux @xmath72 and the second eddington factor @xmath73 . i discretize density - like quantities in a shell - centered fashion and flux - like quantities on the grid points . let @xmath74 be the value of a physical quantity at the @xmath75th grid point in the variable @xmath43 and the @xmath76th grid point in photon energy . away from the boundaries , i differentiate density - like quantities over the variable @xmath43 on the grid points as @xmath77 and flux - like quantities as @xmath78 in an explicit or fully implicit scheme .    in both equations i approximate the first derivative of a physical quantity @xmath79 with respect to photon energy by @xmath80 and the second derivative with respect to photon energy by @xmath81 where the index ` @xmath75 ' in the quantity @xmath74 may also be fractional , corresponding to a shell center , as needed by an implicit differencing scheme over the variable @xmath43 . the energy differencing scheme  ( [ 1stder])([2ndder ] ) can not be used directly in differencing the radiation energy density and flux in equations  ( [ 0th3 ] ) and ( [ 1st ] ) because it does not guarantee particle conservation when only scattering is taken into account ( chang & cooper 1970 ) . this is important when the energy spectrum has reached quasi - equlibrium and small numerical errors lead to significant deviation from particle conservation ( see eq .  [ 10 ] in chang & cooper 1970 ) . a better differencing scheme can be deduced by examining the properties of the differential operator at that limit . the solution of the moment equations reaches quasi - equlibrium when the photon mean free path is very small compared to any characteristic length scale in the system . in that limit , @xmath82 , and when only processes that conserve photon number are operating , the moment equation of zeroth order becomes @xmath83\ ; ,    \label{cons}\ ] ] where @xmath11 is the dimensional electron density . the photon number flux is simply equal to @xmath84 and integrating equation  ( [ cons ] ) over photon energy gives @xmath85 i.e. , photon conservation . following chang & cooper ( 1970 ) , i define a generalized photon current @xmath86 as @xmath87 and write the zeroth and first moment equations as @xmath88 where the coefficient are given in appendix  d. differencing equations  ( [ 0th4])([1st4 ] ) according to the scheme  ( [ back])([2ndder ] ) results in strict photon number conservation in the limit of small photon mean free paths .    a final point with the energy differencing is related to the requirement that the photon energy density is everywhere positive . this introduces a restriction on the spacing of the energy grid , depending on the differencing scheme . for the scheme used above it can be shown that the requirement for a positive photon energy density results in the restriction ( cf . .  [ [ cons ] ] and chang & cooper 1970 , eq .  [ 14 ] ) @xmath89 which can be easily met . i solve the system of equations  ( [ 0th4 ] ) and ( [ 1st4 ] ) for the zeroth and first moments of the specific intensity of the radiation field , @xmath70 and @xmath72 , in @xmath64 grid points in optical depth and @xmath90 grid points in photon energy . therefore the system of equations has @xmath91 unknowns . in all the interior grid points ( @xmath92 ; @xmath93 ) i solve both equations  ( [ 0th4 ] ) and ( [ 1st4 ] ) using the differencing scheme of equations  ( [ back])([2ndder ] ) . for the inner spatial boundary , i.e. , for @xmath94 and for all @xmath76 , i solve equations  ( [ inboundmom ] ) and ( [ 1st4 ] ) , whereas for the outer spatial boundary , i.e. , for @xmath95 and for all @xmath76 , i solve equations  ( [ outboundmom ] ) and ( [ 0th4 ] ) . finally , at both energy boundaries , i.e. , for all @xmath75 and for @xmath96 , i set both the zeroth and first moments of the specific intensity equal to zero , according to boundary conditions ( [ ein ] ) and ( [ eout ] ) . the resulting system of equations is linear in both the zeroth and first moments of the specific intensity . defining the quantity @xmath97 as @xmath98 the system of equations takes the form shown in figure  1 . the matrix of coefficients is a sparse matrix and can be solved efficiently with the biconjugate gradient method ( see , e.g. , press et al . 1992 , pg . 209 ) with a preconditioner matrix . i choose the preconditioner matrix to be the coefficient matrix , with the elements at the outer two diagonal bands set equal to zero . the preconditioner matrix can then be inverted efficiently with the lu - decomposition method ( see , e.g. , press et al . 1992 , pg . 202 ) . in this method , in order to avoid round - off errors , i scale all the rows in the coefficient matrix so that the highest value of the elements in each row is unity . depending on the complexity of the particular problem , the system of difference equations can be solved to a fractional accuracy of @xmath99 within less than fifty iterations , which can be performed in a few seconds on a current workstation .      in the iterative method of the variable eddington factors , the transfer equations  ( [ rteq ] ) is solved at every grid point in the @xmath100 plane for the boundary conditions  ( [ inbound ] ) and ( [ outbound ] ) . because the generalized source function is calculated using the results from the previous iterations , the solution of equation ( [ rteq ] ) is just the formal solution that can be obtained analytically , i.e. , @xmath101 where @xmath102 and @xmath103 or @xmath104 when @xmath105 or @xmath106 respectively . in calculating the integrals in equation  ( [ formal ] ) , i use a piece - wise linear interpolation of the source function through each shell and evaluate the energy derivatives according to equations  ( [ 1stder])([2ndder ] ) . the rapid convergence of the method of variable eddington factors is preserved in the current implementation for the problem of compton scattering . this is shown in figure  2 , where the maximum fractional error in the radiation energy density is plotted for each iteration , for two uniform media with different scattering optical depths . for the model problems shown , the outer radius of the scattering medium is set to @xmath107 , its electron temperature to 10  kev , and the energy - dependence of the radiation flux at the inner boundary to be that of a blackbody with a temperature of 0.5  kev . the converge is typically faster for problems with higher optical depths , because i have used @xmath108 as an initial value for the first variable eddington factor , which is closer to the real value at the limit of high optical depth . the numerical method described in this section also allows for a consistency check between the solution of the transfer equation and the solution of the moment equations that can be performed for each individual problem . for this test , equation  ( [ rteq ] ) can be solved iteratively without using the solution of the system of equations  ( [ 0th3])([1st ] ) , by starting with an assumed generalized source function [ e.g. , @xmath109 , solving the transfer equation for the specific intensity , updating the generalized source function , and repeating until convergence is achieved . the moments of the specific intensity can then be obtained in this way and compared directly with the moments obtained by solving equations  ( [ 0th3 ] ) and ( [ 1st ] ) . the desired accuracy of the calculations , which is affected mainly by the choice of the discrete mesh , depends on the particular problem . because the main goal of this study is the calculation of the spectra of accreting compact objects , the x - ray colors of which are observed to change by a few percent when the mass accretion rate changes by a factor of @xmath110 , i require for the solutions a fractional accuracy of @xmath111 at each grid point . in this section i solve the radiative transfer equation for a number of model problems with spherical symmetry that are related to various astrophysical systems and in particular to media around compact objects . in  5.1 , i discuss comptonization of soft x - ray photons by a static , uniform medium of hot electrons . in  5.2 , i solve the radiative transfer equation in a cold , divergence - free inflow ; although such a configuration requires fine tuning of the physical conditions around an accreting object and may not occur in nature , it is used here to provide physical understanding of the effect of comptonization by the bulk electron velocity . finally , in  5.3 , i discuss more realistic problems of spherical accretion onto and outflows from compact stars . for the first model problem , i consider a static , spherically symmetric , purely scattering medium and set the parameters to values that are typical for weakly magnetic , accreting neutron stars ( see psaltis et al . i set the outer radius of the scattering medium to @xmath107 , its electron temperature to 10  kev , its electron scattering optical depth to 4 , and the energy - dependence of the radiation flux at the inner boundary to be that of a blackbody with a temperature of 0.5  kev ; the normalization of the input flux is arbitrary because i have assumed that there are no sources of photons in the medium . this is a frequently solved problem that has been analyzed in detail by , e.g. , katz ( 1976 ) , shapiro et al . ( 1976 ) , and sunyaev & titarchuk ( 1980 ) . all these authors showed that the emerging radiation spectrum at energies much larger than the energy of the injected photons is largely independent of the details of the input spectrum and can be described by a power - law with an exponential cut - off at @xmath112 . here i solve this simple problem to study the fact that the variable eddington factors are not independent of energy even when the photon mean - free path is independent of energy and the medium is static ( see also the discussion in paper  i ) . figure  3 shows the zeroth and first moments of the specific intensity of the radiation field , @xmath70 and @xmath72 , as a function of photon energy and radius . both quantities decrease overall with increasing radius because of the dilution of the radiation field . the energy dependence of the radiation energy density changes within one photon mean - free - path from the inner boundary but remains unchanged at larger radii . on the other hand , the radiation flux evolves throughout the scattering medium . this is a consequence of the fact that the system of moment equations has been solved by imposing an inner boundary condition on @xmath113 and an outer boundary condition on the ratio @xmath114 and not on @xmath70 itself . figure  4 shows the first eddington factor @xmath115 at two photon energies as well as the ratio of the energy - integrated second to zeroth moments of the specific intensity of the radiation field . the eddington factor depends on photon energy even though the scattering cross section and hence the photon mean - free - path are mostly energy independent at these low photon energies . this is because photons were injected into the medium with energies @xmath116 and on average gain energy at each scattering . as a result , photons emerging with low energy have experienced on average a smaller number of scatterings than photons with higher energy and therefore their distribution is less isotropic . figure  5 compares the emerging radiation spectrum obtained using the method of variable eddington factors with the spectra obtained using the same boundary conditions but two energy - independent prescriptions of the first eddington factor that are commonly used ; the second eddington factor does not enter the calculation when the bulk velocity of the electrons is zero . note that neither prescription is the correct diffusion approximation for compton scattering because they both neglect the energy dependence of @xmath71 on photon energy shown in figure  4 . however , the discrepancy between the self - consistent solution and the ones obtained using prescribed eddington factors is small for the prescription that depends on optical depth and can be as large as @xmath117% at high photon energies for the prescription that is independent of optical depth . this discrepancy increases when the ratio of the photon mean - free path to the smallest characteristic length scale in the system increases .      in a moving medium , the photons gain energy by scattering off the fast moving electrons . as measured by an observer _ comoving with the medium _ , the photon energy density increases mostly due to the convergence of the flow , which is described by terms that are proportional to @xmath118 ( blandford & payne 1981 ) . on the other hand , the flux that is measured by an observer _ at infinity _ increases because of the systematic energy exchange between photons and electrons , which is described mostly by terms that are second - order in @xmath1 ( psaltis & lamb 1997 ) .    in order to investigate the effects of comptonization by the bulk electron flow and distinguish them from the effects of the convergence of the flow , i study in this section the photon - electron interaction and the formation of spectra in fast , divergence - free flows . for the model problems shown below , i set the outer radius of the scattering medium to @xmath119 , the electron temperature to zero , and the electron velocity to the divergence - free profile @xmath120 . finally , i set the electron density to a constant value , which arises from mass conservation in the flow , and characterize each flow by the total electron scattering optical depth in the radial direction . figure  6 shows the radiation spectrum at the inner boundary , the emerging radiation spectrum when only terms that are first - order in the electron bulk velocity are taken into account , as well as the emerging radiation spectrum that is correct to second order in the electron bulk velocity , for a flow with a scattering optical depth of 2 . the terms that are first - order in the electron bulk velocity describe mainly the effects of non - relativistic doppler shift and hence produce a displacement in @xmath121 of the input spectrum towards higher photon energies ( see also psaltis & lamb 2000 ) . on the other hand , the terms that are second - order in the electron bulk velocity result in a systematic photon upscattering and produce a power - law tail at photon energies higher than the average energy of the injected photons . note here that the approximation employed here of keeping , in the scattering kernel , only terms that are first order in @xmath6 limits the validity of the calculation to energies @xmath122  kev .    when the total optical depth of the scattering medium , and hence the average number of scatterings that each photon experiences , are increased , then on average photons gain more energy by scattering off moving electrons . as a result , the emerging radiation spectrum is displaced towards higher photon energies , and the power - law tail becomes flatter ( see fig . when the bulk electron velocity in the flow increases , the average energy gain per scattering increases as well , and hence the power - law tail also becomes flatter ( fig . 7b ) . as a result , scattering of photons by fast moving electrons produces power - law tails at high photon energies , even in divergence - free flows , in a way that is similar to the well understood process of thermal comptonization . i finally perform simple model calculations of accretion onto and outflows from a central object by assuming a free - fall density profile and setting everywhere in the flow the electron temperature to zero and the inward radial velocity equal to the free - fall velocity from infinity onto a compact object with a radius of 10  km and a gravitational mass of @xmath123 . i also set the outer radius of the medium to @xmath119 ; at these large radii the photon mean - free path is large compared to radius and the electron bulk velocity is small compared to the speed of light that the flow does not effect significantly the photon distribution . finally , instead of specifying the mass accretion rate onto the compact object , i specify the total electron - scattering optical depth of the flow . figure  8a shows the radiation spectrum at the inner boundary , the emerging radiation spectrum when only terms that are first - order in the electron bulk velocity are taken into account , as well as the emerging radiation spectrum that is correct to second order in the electron bulk velocity . as in the model calculation discussed in  5.2 , the terms that are first - order in the electron bulk velocity describe mostly the non - relativistic doppler shifts and produce a displacement in @xmath121 of the input radiation spectrum towards higher photon energies . on the other hand , the terms that are second - order in the electron bulk velocity systematically upscatter photons and produce a power - law tail at high photon energies . when the optical depth of the flow is increased , the systematic upscattering of photons becomes more efficient and the high - energy tail becomes flatter . figure  9 shows the same information as figure  8 but for the case of an outflow . for this model problems i have used the same parameters as in the inflow calculation but changed the sign of the velocity at all radii . in an outflow , the terms that are first - order in the electron bulk velocity produce a displacement in @xmath121 of the input radiation spectrum towards _ lower _ photon energies . on the other hand , the effect of the terms that are second - order in the electron bulk velocity is the same as in the case of inflow , i.e. , they systematically upscatter the photons and produce a tail at high photon energies . for the optical depths and bulk velocities considered here , the effect of the second - order terms nearly cancels the effect of the first - order terms and the major difference between the input and emerging spectra is the power - law tail at high photon energies . as in the case of the inflow , when the optical depth of the flow increases , the power - law tail becomes flatter , but the overall effects is significantly reduced . in the previous sections i described a method for solving numerically the radiative transfer equation that describes compton scattering in spherically symmetric , static and moving media , and applied it to a variety of cases related to flows around compact objects . i showed that the first eddington factor , which describes the degree of isotropy of the radiation field at a given point in space , may depend on photon energy , even though the photon mean - free - path is energy independent , if the photons are preferably upscattered or downscattered in energy at each scattering . i also demonstrated that the terms that are first - order in the electron bulk velocity alter the input spectrum in a way that is qualitatively and quantitatively different than the effect of the terms that are of second - order . the systematic upscattering of photons caused by bulk comptonization is described by the terms that are second - order in @xmath1 . this becomes apparent by the fact that , at each scattering with a moving electron , the average fractional energy change of a photon in the system frame is always positive , independent of the photon energy , and equal to @xmath124 ( see , e.g. , rybicki & lightman 1979 , ch . 7 ; see also the discussion in psaltis & lamb 1997 ) . if the distribution of photon residence times in the comptonizing medium drops exponentially to zero at late times , then a power - law tail is produced at high photon energies ( for total scattering optical depths @xmath125 ) , in analogy to the generation of power - law tails by thermal comptonization . the velocity at which the terms that are second - order in @xmath1 become important depends on the ratio of the photon mean - free - path @xmath126 to the shortest characteristic length scale @xmath127 in the problem . in general , when @xmath128 , i.e. , in regions close to and within the photon - trapping radius , the terms that are first- and second - order in @xmath1 become comparable in size ( yin & miller 1995 ; psaltis & lamb 1997a ) and hence the latter should not be neglected . the generalization of the method to time - independent multi - dimensional configurations is conceptually easy . however , a problem with axial symmetry requires the solution of the radiative transfer equation in five dimensions in phase space and the introduction of a large number of independent variable eddington factors . moreover , the presence of high flow velocities in multi - dimensions introduces sharp gradients in the variable eddington factors , which can not be guessed a priori , as required for the first iteration ( dykema , klein , & castor 1996 ) . even with these problems , however , the method of variable eddington factors will still offer a significant advantage , when compared to any other iterative method , which arises from the fact that the eddington factors are bound to lie over a very narrow range of values . the fast convergence of the numerical algorithm presented here makes it also ideal for future extensions to transport problems coupled to time - dependent hydrodynamics . for the case of accretion onto and outflows from compact objects , the characteristic photon diffusion timescale from a flow of dimension @xmath129 is of order @xmath130 which is significantly smaller than the corresponding @xmath131  ms dynamical timescale . as a result , for most cases the structure of the radiation field can be obtained by solving the time - independent radiative transfer problem for each snapshot of the fluid properties . given that the variable eddington factors will be only marginally different between successive timesteps , the convergence of the algorithm will be even faster than what is shown in figure  2 . i am grateful to fred lamb for his support and numerous useful discussions , to dimitri mihalas for his help on solving numerically radiative transport equations , and to d.  swesty for bringing to my attention the biconjugate gradient method for solving linear systems of equations . it is also a pleasure to thank luca zampieri and feryal zel for many helpful discussions and comments and an anonymous referee for a very careful review of the manuscript . i acknowledge the support of a postdoctoral fellowship of the smithsonian institution . this research was also supported in part by nsf grant ast  93 - 15133 and nasa grant nag  5 - 2925 .    * appendices * in defining the moments of the specific intensity in the system frame i shall use the quantities @xmath132\ ] ] and @xmath133\;.\ ] ] the zeroth moment of the specific intensity is then @xmath134 where @xmath22 depends on @xmath135 and @xmath30 through the relation @xmath136 . because of the assumed spherical symmetry , the only non - zero component of the first moment of the specific intensity is the @xmath22component : @xmath137 similarly , the non - zero components of the second , third , and fourth moments of the specific intensity are @xmath138 finally , i shall also use the first two eddington factors defined by @xmath139 using equation ( a8)(a12 ) of paper i , i derive the generalized source function for a spherically symmetric system : @xmath140 where @xmath141 j\\ s_2^\pm & = & \frac{3}{8 } \left\ {      \left(-1 + \e\partial_\e\right )     \left [ 3 - \cos^2\theta+\left(3\cos^2\theta-1\right)f^{rr }           \right]\right . j \nonumber\\ & & \quad\quad+\left.\left(1-\e\partial_\e\right )      \left[5 - 3\cos^2\theta+\left(5\cos^2\theta-3\right )           g^{rr}\right]h^r\cos\theta\right\}\\ s_3^\pm & = & \frac{3}{4}\left\ {     \left[3\cos^2\theta-1-      3\left(3\cos^2\theta-1\right)f^{rr}\right ]     j\right.\nonumber\\   & & \quad\quad+ 2\left[1 - 3g^{rr}+     \left(5 g^{rr}-3\right)\cos^2\theta \right]h^r     \cos\theta\nonumber\\   & & \quad\quad+   \frac{1}{2}\left(-2\e \partial_\e     + \e^2\partial^2_\e\right )     \left[3-\cos^2\theta+\left(3\cos^2\theta-1\right)f^{rr }     \right]j\nonumber\\ & & \quad\quad+ \frac{1}{2}\left(2\e \partial_\e     -\e^2\partial^2_\e\right )     \left[5 - 3\cos^2\theta+\left(5\cos^2\theta-3\right)g^{rr }     \right]h^r\cos\theta\\ s_4^\pm & = &   v\cos\theta i + \frac{3}{8}v\left\ {     4\cos\theta\left[2-\cos^2\theta        + \left(3\cos^2\theta-2\right)f^{rr}\right]j      \right.\nonumber\\   & & \quad\quad    + \left[-1 - 5\cos^2\theta+\left(3\cos^2\theta-         1\right)g^{rr}\right]h^r\nonumber\\   & & \quad\quad    -\cos\theta\partial_\e \left[3-\cos^2\theta+      \left(3\cos^2\theta-1\right)f^{rr}\right]j\nonumber\\   & & \quad\quad\left . + \partial_\e    \left[\frac{1}{2}\left(3-\cos^2\theta\right)+     \left(3\cos^2\theta-1\right)g^{rr}\right]h^{r}\right\}\\ s_5^\pm & = & \frac{3}{4}v^2 \left\ { \left [   -5\cos^4\theta+\frac{21}{2}\cos^2\theta-1    -2\cos^2\theta\left(2-\cos^2\theta\right )   \e\partial_\e \right.\right.\nonumber\\ & & \quad\quad\quad\quad   \left.+\frac{1}{4}\cos^2\theta\left(3-\cos^2\theta\right )    \e^2 \partial_\e^2\right]j\nonumber\\ & & \quad\quad+\left[4 - 10\cos^2\theta+     \left(3+\cos^2\theta\right)\e \partial_\e    -\frac{1}{2}\left(3-\cos^2\theta\right )    \e^2\partial_\e^2\right]\cos\theta h\nonumber\\ & & \quad\quad + \left [    \frac{1}{2}\left(30\cos^4\theta-37\cos^2\theta+7\right )    + 2\left(1 - 3\cos^4\theta\right)\e\partial_\e   \right.\nonumber\\ & & \quad\quad\quad\quad \left . + \frac{1}{4}\left(3\cos^4\theta-3\cos^2\theta+4\right )    \e^2\partial_\e^2\right]f^{rr}k^{rr}\nonumber\\ & & \quad\quad + \left [   \left(6\cos^2\theta-4\right)+\left(3\cos^2\theta-3\right )   \e\partial_\e   -\frac{1}{2}\left(3\cos^2\theta-1\right)\e^2   \partial_\e^2\right]\cos\theta g^{rr } h^r   \nonumber\\ & & \quad\quad \left . + \left(2\cos^2\theta-1\right)\left ( 1 + 2\e    \partial_\e + \frac{1}{2}\e^2\partial_\e^2    \right)r^{rrrr}\right\}\end{aligned}\ ] ] where @xmath142 and the variable eddington factors and the moments of the specific intensity are defined as in appendix a. using equations ( 34 ) and ( 40 ) of paper i , i derive @xmath143     \nonumber\\ & & \quad\quad-\left[1+\frac{1}{2}v^2      + \left(f^{rr}+\frac{1}{2}\right)v^2 \e\partial_\e      + \frac{1}{2}f^{rr}v^2 \e^2\partial^2_\e\right]\chi\\ a_2 & = & n_e \left\{\e-2t_e-\frac{2}{3}v^2     + \left[\frac{1}{10}\left(f^{rr}-\frac{1}{3}\right )        + \frac{22}{20}\e\partial_\e f^{rr}\right]v^2\right\}\\ a_3 & = & n_e \left[t_e+\frac{1}{3}v^2 +     \frac{11}{20}\left(f^{rr}-\frac{1}{3}\right)v^2\right]\\ a_4 & = & -\frac{2}{\tau_0-\tau_{\rm r}}+v\chi + v \e\partial_\e       \chi\\   a_5 & = & n_e v\\ c_1 & = & -\chi s_e + \frac{1}{6}v^2 s_e      \left(\e\partial_\e - \e^2 \partial^2_\e\right)\chi     + \frac{1}{6}v^2 \e \left(\partial_\e s_e\right )     \left(1 - 2\e\partial_\e\right )     \chi\nonumber\\ & & - \frac{1}{6 } v^2 \e^2 \left(\partial^2_\e s_e\right ) \chi\end{aligned}\ ] ] and @xmath144\nonumber\\ & & \quad\quad      -\left[1+\frac{1}{2}v^2+v^2\left(g^{rr }     + \frac{1}{2}\right ) \e\partial_\e     + \frac{1}{2}v^2g^{rr}\e^2\partial^2_\e\right]\chi\\ b_4 & = & n_e \left[-\frac{2}{5}\e+\frac{4}{5}t_e     + \frac{9}{10}v^2     -\frac{3}{10}v^2\left(1+\frac{2}{3}\e\partial_\e\right )       g^{rr}\right]\\ b_5 & = & n_e\left[-\frac{2}{5}t_e -\frac{1}{10}v^2     \left(g^{rr}+3\right)\right]\\ c_2 & = & -\frac{2}{3}v\chi s_e+     \frac{1}{3}v \chi \e\partial_\e s_e + \frac{1}{3}v s_e \e\partial_\e \chi\;,\end{aligned}\ ] ] where @xmath145 and i have suppressed the dependence of the various quantities on spatial position and photon energy . using equations ( 34 ) and ( 40 ) of paper i i derive @xmath146\chi\\ a_2 ' & = & n_e \left[\frac{1}{10}\left(f^{rr}-\frac{1}{3}\right )        + \frac{22}{20}\e\partial_\e f^{rr}\right]v^2\\ a_3 ' & = & n_e      \frac{11}{20}\left(f^{rr}-\frac{1}{3}\right)v^2\\ a_4 ' & = & -\frac{2}{\tau_0-\tau_{\rm r}}+v\chi + v \e\partial_\e       \chi\\   c_1 ' & = & -\chi s_e + \frac{1}{6}v^2 s_e      \left(\e\partial_\e - \e^2 \partial^2_\e\right)\chi     + \frac{1}{6}v^2 \e \left(\partial_\e s_e\right )     \left(1 - 2\e\partial_\e\right )     \chi\nonumber\\ & & - \frac{1}{6 } v^2 \e^2 \left(\partial^2_\e s_e\right ) \chi\;,\end{aligned}\ ] ] where @xmath145 and i have suppressed the dependence of the various quantities on spatial position and photon energy .
Abstract:

Answer: i study the formation of comptonization spectra in spherically symmetric , fast moving media in a flat spacetime . i analyze the mathematical character of the moments of the transfer equation in the system - frame and describe a numerical method that provides fast solutions of the time - independent radiative transfer problem that are accurate in both the diffusion and free - streaming regimes . i show that even if the flows are mildly relativistic ( @xmath0 , where @xmath1 is the electron bulk velocity in units of the speed of light ) , terms that are second - order in @xmath1 alter the emerging spectrum both quantitatively and qualitatively . in particular , terms that are second - order in @xmath1 produce power - law spectral tails , which are the dominant feature at high energies , and therefore can not be neglected . i further show that photons from a static source are upscattered by the bulk motion of the medium even if the velocity field does not converge . finally , i discuss these results in the context of radial accretion onto and outflows from compact objects .