dynamical friction is the reaction force that a perturber feels from the gravitational wake it generates in the ambient medium .
this fundamental interaction has therefore had a wide range of applications in the formation of astronomical structures .
applications to collisionless systems include mass segregation in star clusters ( e.g. kim and morris 2003 , portegies et al 2004 ) , galaxy distribution in hierarchical cluster formation ( e.g. frenk et al 1996 , bullock and johnston 2005 ) , the decay of satellites in galactic halos and galaxy mergers ( binney and tremaine 2008 and references therein ) .
dynamical friction also occurs in collisional particle systems and was applied to the study of the growth of planetesimals ( hornung et al 1985 , stewart and whetherill 1988 ) , the eccentricity excitation of planetary embryo orbits ( ida 1990 , kokubo and ida 1996 , namouni et al 1996 ) , and the confinement of planetary rings ( goldreich and tremaine 1982 ) . for gaseous systems , applications include the orbital decay of compact stars orbiting supermassive black holes ( narayan 2000 ) , and the heating of intercluster gas by decelerating supersonic galaxies ( el - zant et al .
2004 , faltenbacher et al 2005 , kim et al 2005 ) .    following the pioneering work of chandrasekhar ( 1943 ) on the gravitational drag experienced by a perturber in collisionless star systems , the derivations of dynamical friction
assume that the perturber has a constant velocity in the ambient medium .
the assumption of uniform motion is made for collisionless systems ( binney and tremaine 2008 ) , collisional particle systems ( hornung et al 1985 , stewart and whetherill 1988 , ida 1990 , goldreich and tremaine 1982 ) as well as gaseous media ( kim and kim 2007 , kim et al 2008 , ostriker 1999 ) . the resulting dynamical friction force which depends on the perturber s velocity
is subsequently used to model the evolution of the perturber s motion whose velocity varies with time .
only the numerical simulations of the nonlinear fluid equations are free from this assumption ( snchez - salcedo and brandenburg 1999 , 2001 ) .
the neglect of the perturber s acceleration in dynamical friction is justified by the small amplitude of the force that is proportional to the square of the pertuber s mass .
neglecting acceleration is independent of the type of motion and is mainly motivated by the lack of an analytical framework where the perturber s accelerated motion could be taken into account . to overcome this obstacle
, we derive in appendix a a new expression of dynamical friction that applies to a general type of motion and depends on the perturber s motion and the boundary conditions but not the density perturbation .
in particular , dynamical friction is shown to be non - local as it depends on the whole trajectory of the perturber . although non - locality does not appear if the perturber s velocity is constant , its origin is easy to identify : as the perturber s wake propagates in the medium , a velocity change because of acceleration affects the density enhancement near the perturber and offsets the force from its local value as the far parts of the wake were generated by a different velocity . therefore when acceleration is taken into account , the amplitude of dynamical friction may change with respect to the value derived assuming a constant velocity because in modelling the force on an accelerated perturber as that on a constant velocity perturber
, one assumes that the gravitational wake produced by accelerated motion is similar to that of uniform motion .
we shall see that this difference between the two wakes makes dynamical friction dependent on the initial velocity of the perturber .
another difference between the wakes generated by constant velocity motion and accelerated motion is the absence of a perturbed region with front - back symmetry with respect to a subsonic accelerating perturber . in effect
a constant velocity perturber raises a gravitational wake in the ambient medium that is symmetric with respect to it and truncated by the sonic shockwave .
this truncation defines two regions , one that lacks front - back symmetry located behind the perturber , and one with front - back symmetry around the perturber that becomes the largest region inside the wavefront for velocities smaller than half the sound speed .
dynamical friction arises from the former region as shown by ostriker ( 1999 ) .
when acceleration is taken into account , front - back symmetry is lost everywhere .
this would intuitively imply that dynamical friction on an accelerated perturber should be much larger than the value derived assuming uniform motion for subsonic perturbers .
curiously , we shall show that the force on an accelerating subsonic perturber initially at rest is similar to the uniform friction force .
lastly , accelerated motion friction would differ from uniform motion friction because the former never reaches steady state . in effect , dynamical friction on a constant subsonic velocity perturber is known to be time - independent ( ostriker 1999 ) .
it is the object of this paper to quantify how taking into account acceleration modifies dynamical friction on a perturber travelling in an infinite homogeneous gaseous medium .
the perturber s motion is assumed to have a constant acceleration , a choice that makes the analytical developments simpler .
physically such constant acceleration approximates the deceleration felt by the perturber because of dynamical friction as the force amplitude is small and the corresponding instantaneous velocity change is equally small .
constant acceleration also models the motion of a young star with an asymmetric jet system moving in interstellar gas ( hirth et al 1994 , hartigan and morse 2007 , namouni 2005 , 2007 ) . in what follows
, we characterize the density perturbation associated with accelerated motion in section 2 .
we derive the force on an accelerated perturber in section 3 using the general expression of dynamical friction ( derived in appendix a ) that does not depend on the density enhancement caused by the perturber but only on the perturber s motion and the medium s boundary conditions . in particular we show that including acceleration makes dynamical friction dependent on the initial velocity of the perturber .
section 4 contains concluding remarks . in appendix
b , the new expression of dynamical friction derived in appendix a is validated by applying it to uniform rectilinear motion and recovering the known results on dynamical friction .
consider a point - like perturber of mass @xmath0 on an accelerated rectilinear trajectory , of initial velocity @xmath1 and acceleration @xmath2 , that travels inside an inviscid isothermal gaseous medium of density @xmath3 and sound speed @xmath4 .
assuming the density perturbation associated with the perturber s wake is small compared to the medium s unperturbed density @xmath5 , the corresponding density enhancement @xmath6 defined as @xmath7 satisfies the standard forced sound wave equation ( appendix a , equation [ wav2a ] ) : @xmath8 { \cal h}(t ) , \label{waveeq}\ ] ] where @xmath9 is the perturber s gravitational potential , @xmath10 is the heaviside function , @xmath11 is the gravitational constant , @xmath12 is the perturber s trajectory chosen along the @xmath13direction and @xmath14 is dirac s delta function .
the solution of this equation is the linard - wiechert potential which may be written in integral form using the retarded green function as ( [ dens2a ] ) : @xmath15 { \cal h}(u ) } { |{\bm x}-{\bm \xi}(u)|}\ { \rm d}u .
\label{density}\end{aligned}\ ] ] the acceleration @xmath2 introduces two natural time and length scales given by : @xmath16 in these units , the perturber s trajectory may be written as @xmath17 where @xmath18 is the initial mach number .
the mach number is a function of time and reads : @xmath19 .
the scaling time @xmath20 is the duration to increase the mach number by unity and the scaling length @xmath21 is the corresponding travelled distance .
applied to the density enhancement these scalings yield : @xmath22^{1/2}\right ) { \cal h}(u ) } { \left[r^2+\left(z-{\cal m}_0u- u^{2}/2\right)^2\right]^{1/2}}\ { \rm d}u , \label{density2}\ ] ] where @xmath23 is the cylindrical radius ; unless it is stated otherwise , all length and time variables and constants are normalized to ( [ scalings ] ) . in order to derive the density profile we use the property of the @xmath24-function @xmath25=\sum_i \alpha(u_i ) { \delta(u - u_i)}/{|\beta^\prime(u_i)|}$ ] where @xmath26 are the roots of @xmath27 , and @xmath28 is the derivative of @xmath29 .
the roots of the argument of the @xmath24function in equation ( [ density2 ] ) satisfy the quartic equation : @xmath30 where @xmath31 is the spherical radius .
this equation is the counterpart of the quadratic equation derived by ostriker ( 1999 ) for uniform motion .
the former may be reduced to the latter once time and length dimensions are restored .
equation ( [ roots ] ) , however , can not be solved analytically because its varying coefficients do not permit simple and unique expressions of the roots .
instead we determine the roots numerically and substitute them into the density expression obtained by applying the @xmath24function property cited above .
the density expression reduces to : @xmath32 fig .
( 1 ) follows the level curves of @xmath6 for perturbers initially at rest ( @xmath33 ) as they evolve from @xmath34 to mach 6 . the front - back symmetry that characterizes the density enhancement of subsonic constant velocity motion in a uniform medium is absent for both subsonic and supersonic motion .
the perturber reaches mach 1 at @xmath35 but remains inside the sonic wavefront of radius @xmath36 as the latter has been propagating faster than the perturber s velocity .
once the velocity is supersonic , a mach cone is formed inside the sonic wavefront and a shallow density depression develops and trails the perturber ( @xmath37 ) .
eventually as mach 2 is reached the perturber catches up with the sonic shockwave and exits the sonic region to form an external mach cone .
( 2 ) shows the level curves of @xmath6 associated with decelerating perturbers with @xmath38 as a function time until the perturber s velocity vanishes .
the density enhancement of subsonic perturbers has a structure similar to that of uniform motion ( ostriker 1999 ) but lacks its front - back symmetry with respect to the perturber .
we show in the following section that under certain circumstances , dynamical friction may be not significantly affected by the loss of front - back symmetry even as closed level curves around the perturber contribute to the force unlike those of uniform motion ( see introduction ) .
the deceleration of supersonic perturbers modifies the structure of the density enhancement and reduces the extent of the mach cone but does not erase it as the velocity vanishes ( bottom row , right panel ) indicating the presence of a finite friction force .
the reaction force exerted on a point - like perturber in a homogeneous gaseous medium by its gravitational wake is given by equation ( [ for3a ] ) and is written in dimensional variables as : @xmath39 , r\leq ct}r\sin\theta{\rm d}r\ , { \rm d}\theta\ , { \rm d}\varphi   \ \frac{{\bm y}-{\bm
\delta}}{|{\bm y}-{\bm \delta}|^3},\label{force}\ ] ] where @xmath40 , @xmath41 and @xmath42 are the spherical coordinates of the vector @xmath43 defined as the relative position within the medium with respect to the perturber s position at the retarded time @xmath44
. it may be written as @xmath45 where @xmath46 is the position of fluid element in the medium .
this definition of @xmath43 influences the force integral through the boundary conditions of the medium denoted by @xmath47 .
the phase @xmath48 , the upper radial cutoff at the retarded sonic shockwave @xmath49 , and the boundary conditions are all that is necessary to compute the friction force . in particular
, the absence of an explicit expression for the density enhancement is irrelevant for the derivation of the force .
the derivation of expression ( [ force ] ) assumes that the perturber is assembled at @xmath34 a fact that introduces the size of the retarded sonic shockwave as an upper cutoff .
this assumption is not restrictive as it allows us to study the time - dependent evolution of dynamical friction as well approximate steady state behaviour ( see section 3.2 ) .
it also describes turning on a numerical simulation of dynamical friction ( snchez - salcedo and brandenburg 1999 , 2001 ) .    for rectilinear motion along the @xmath13direction in an infinite medium ,
the force is equally along @xmath50 and its expression can be reduced further by integrating over the longitude @xmath51 and the variable @xmath52 as : @xmath53 , r\leq t}\ \frac{{\rm d}r } { \delta^2}\ \left[\frac { r-\tau\delta}{\left|1+(\delta / r)^2 - 2(\delta / r)\tau\right|^{1/2}}\right]^{\tau_{\rm max}}_{\tau_{\rm min}},\label{forcex}\ ] ] where the scalings ( [ scalings ] ) have been applied ; the perturber s trajectory @xmath54 , the time , @xmath36 , and the position @xmath46 and its spherical components are dimensionless hereafter .
the expressions of @xmath55 and @xmath56 as a function of the perturber s trajectory and the variable @xmath40 are determined from the boundary conditions . for uniform rectilinear motion , @xmath57 is constant and the phase @xmath58 .
calling @xmath59 the gravitational potential s small distance cutoff scaled to @xmath21 ( [ scalings ] ) , the integral that appears in the expression ( [ forcex ] ) may be written as ( appendix b , equation [ ostria ] ; ostriker 1999 ) : @xmath60\left[1+{\rm sign}(1-{\cal m})\right]-\log \left|\frac{1+{\cal m}}{1-{\cal m}}\right| -\left[1-{\rm sign } ( 1-{\cal m})\right ] \ \log\left|\frac{(1-{\cal m})t}{r_0}\right|\label{ostriker}\right).\end{aligned}\ ] ] for accelerated motion , the phase @xmath61 . in order to examine how acceleration modifies the dynamical friction force , we consider the three cases of perturbers initially at rest ( @xmath33 ) , accelerating perturbers ( @xmath62 ) and decelerating perturbers ( @xmath63 ) .    before deriving dynamical friction for accelerated motion
, a choice of the small distance cutoff @xmath59 must be made .
this cutoff is usually taken as the perturber s accretion radius or the system s typical size .
the former dimensional cutoff reads @xmath64 and is a function of the mach number whereas the latter is a constant that does not depend on the flow s variables . in this work ,
the perturber is treated as a compact object whose gravitational potential is truncated at small distances by the perturber s size . to satisfy such a condition and be consistent with a point - like potential
, the perturber s size needs to be much smaller than the accretion radius @xmath64 .
this choice of cutoff is in agreement with the numerical simulation of the fluid equations in the context of bondi - hoyle - lyttleton accretion ( edgar 2005 ) . for
larger size or extended objects , equation ( [ force ] ) no longer applies and dynamical friction may be obtained from the general expression derived for an arbitrary perturbing potential ( appendix a , equation [ forc2a ] ) .      in an infinite medium , the only relevant physical boundary condition is that of the small distance cutoff @xmath65 .
, @xmath66 and @xmath59 remain scaled to @xmath21 ( [ scalings ] ) .
] for a perturber initially at rest , @xmath67 , @xmath68 and the mach number is @xmath69 .
the variables @xmath40 and @xmath70 therefore need to satisfy the three conditions : ( i ) @xmath71 , ( ii ) @xmath72/(2t - r)\equiv f(r , t , r_0)$ ] , and ( iii ) @xmath73 .
the last condition implies that @xmath74 .
examination of the three conditions shows @xmath75 , and that @xmath56 and the remaining bounds on @xmath40 in the integral ( [ forcex ] ) depend on the value of @xmath36 as the velocity increases .
there are five time intervals that need to be considered : @xmath76 case i is that of transitory motion where the density perturbation is in the perturber s vicinity . prior to @xmath77 ,
the force is zero as the three conditions above can not be satisfied .
the integration domain of @xmath40 and @xmath56 are given as : @xmath78^{1/2}. \end{aligned}\ ] ] calling @xmath79 the integral in equation ( [ forcex ] ) and using these boundaries for @xmath40 , we find : @xmath80 _ { r_2^{f}}^{t } , \label{i0}\\ i_{\rm i}&=&-\frac{r_0(2+r^f_2)-(2+r_0+r^f_2)t -r^f_2t^2+t^3}{2r_0t^2}-\frac{2t - r_0}{2t^3}\ , \log\left|\frac{2t - r_0}{tr^f_2-r_0}\right|.\label{case1}\end{aligned}\ ] ] examination of the variation of @xmath81 as a function of time shows that it decays smoothly to zero as @xmath77 is approached .
the expression of the force in this time interval is the equivalent of that derived for uniform motion in appendix b ( equation [ transient ] ) when @xmath82 and the wake has not left the perturber s vicinity . in that case , the force read @xmath83+{\cal m}-{r_0}/{2ct}+(1-{\cal m}^2)ct/2r_0\}$ ] ; it helped explain the transient phase observed before saturation in the numerical simulation of the nonlinear fluid equations ( snchez - salcedo and brandenburg 1999 ) .
case ii covers subsonic motion and the corresponding integration boundaries are : @xmath84^{1/2 } , \label{dom21}\\ r^c_1\leq r\leq   t&\mbox{with } & \tau_{\rm max}=1.\label{dom22}\end{aligned}\ ] ] the latter domain ( [ dom22 ] ) does not contribute to the force .
the corresponding integrand is proportional to @xmath85 and the sign argument is positive for @xmath86 .
the expression of @xmath79 for domain ( [ dom21 ] ) may be obtained from equation ( [ i0 ] ) by substituting @xmath87 for @xmath36 .
the result is : @xmath88 in the limit @xmath89 and @xmath90 , @xmath91 .
this is exactly the friction force acting on subsonic perturbers with uniform motion given by equation ( [ ostriker ] ) .
it is interesting that although the density enhancements of subsonic accelerated motion and subsonic uniform motion differ significantly , the corresponding friction forces are identical when the limit @xmath92 is taken for the accelerated motion friction . )
is independent of @xmath59 and that no assumptions on @xmath59 are made to derive its expression . ]
this result is the more surprising as the density enhancement associated with accelerated motion does not involve a region with front - back symmetry as explained in the introduction leading one to guess wrongly that friction would be stronger for accelerated motion than it is for uniform motion .
the explanation of this agreement comes from the friction force s dependence on the phase @xmath93 ( dimensional variables ) . as the contribution of the force for subsonic motion comes mostly from the vicinity of the perturber s retarded position ( @xmath94 ) , we have @xmath95 and the phase becomes @xmath96 .
the latter is exactly the phase of uniform motion with a mach number @xmath97 regardless of the type of motion @xmath98 as long as , ( 1 ) @xmath89 which implies in dimensional variables that @xmath99 , ( 2 ) the velocity remains subsonic and ( 3 ) the motion does not have a history with supersonic episodes ( see section 3.3 ) .
case iii concerns motion between mach 1 and mach 2 where the perturber travels at supersonic speed but has not yet crossed the sonic shock wave .
the integration domains are given as : @xmath100^{1/2 } , \label{dom31}\\ r^c_2\leq r\leq   r_3^c&\mbox{with } &   \tau_{\rm max}=1,\ \
\mbox{where}\ \ \ r_3^{c}=t-1+[(t-1)^2 - 2r_0]^{1/2},\label{dom32}\\ r^c_3\leq r\leq r^c_1&\mbox{with } &   \tau_{\rm max}=f(r , t , r_0),\label{dom33}\\ r^c_1\leq r\leq   t&\mbox{with } & \tau_{\rm max}=1.\label{dom34}\end{aligned}\ ] ] the integration radii satisfy @xmath101 .
this shows that domain ( [ dom34 ] ) does not contribute a finite term to the force in the same way as domain ( [ dom22 ] ) .
domains ( [ dom31 ] ) and ( [ dom33 ] ) contribute terms similar to that of equation ( [ i0 ] ) except that for the former @xmath36 is replaced with @xmath102 and for the latter @xmath36 is replaced with @xmath87 , and @xmath103 with @xmath104 .
the second domain ( [ dom32 ] ) contributes the term @xmath105 given as : @xmath106_{r_2^{c}}^{r_3^{c } } , \label{i2}\end{aligned}\ ] ] which when combined with the contributions of domains ( [ dom31 ] ) and ( [ dom33 ] ) yields : @xmath107 case iv defines a small time interval near mach 2 where the perturber exits the sonic wavefront .
the corresponding integration domains are : @xmath108 the expression of @xmath79 is obtained by summing up the terms associated with domains ( [ dom31 ] ) and ( [ dom32 ] ) of case iii , and substituting @xmath36 for @xmath109 in domain ( [ dom33 ] ) . the result is : @xmath110 case v concerns supersonic motion that exited the sonic wavefront @xmath111 .
the integration domains are given as : @xmath112 the integral @xmath79 is obtained as previously by combining the domains of case iv and making the necessary boundary substitutions . the result is : @xmath113 for @xmath89 and @xmath114 , the scaled force reads : @xmath115.\label{case5l}\ ] ] this expression differs from that of equation ( [ ostriker ] ) through the first term in the bracket @xmath116 implying that the force on accelerated motion of a perturber initially at rest , is larger than that on uniform motion .
( 3 ) shows the friction force ( [ case1 ] ) , ( [ case2]),([case3 ] ) , ( [ case4 ] ) , ( [ case5 ] ) as a function of time , or equivalently the mach number , for the three values of the small distance cutoff @xmath117 , @xmath118 , and @xmath119 .
the force obtained by using the uniform motion expression , @xmath120 , ( [ ostriker ] ) with @xmath69 is also shown for the same values of @xmath59 . for subsonic motion , the agreement between the two forces is best when @xmath59 is smallest .
the force @xmath120 overestimates maximal drag by about 20% ; it underestimates the drag on supersonic motion by the factor @xmath121 .
lastly , we remark that the derivation of dynamical friction in the five time intervals shows that there are no force singularities associated with accelerated motion such as that at mach 1 for uniform motion ( [ ostriker ] ) .      to an accelerating perturber of trajectory ,
@xmath122 with @xmath62 , corresponds the phase @xmath123 where @xmath19 is the mach number .
it could therefore seem sufficient to take the force expressions of perturbers initially at rest and substitute @xmath57 for @xmath36 .
there is however a key difference between the two cases : the condition from the retarded sonic wavefront , @xmath73 , refers to time regardless of the type of motion .
it should therefore be replaced with @xmath124 that we refer to as condition ( iv ) .
this difference makes the force dependent on the initial velocity in contrast to the expressions derived assuming a constant velocity ( [ ostriker ] ) .
the physical reason for such a dependence was discussed in the introduction : as the velocity varies , the far part of the wake that was launched in the medium corresponds to an earlier velocity . combining the force from the far part of the wake to that near the perturber results in a force that has a memory of the earlier velocity
.    inspection of the conditions ( i , ii ) of the previous section along with ( iv ) shows that the two cases @xmath125 and @xmath126 may be considered separately . in the following
we do not derive the force expressions for the counterparts of cases i and iv as their space extensions vanish when @xmath59 tends to zero .
for an accelerated perturber with an initial supersonic speed ( @xmath125 ) , the boundary conditions ( i , ii , iv ) away from the transient phase ( case i ) are identical to those of case v ( [ dom - v ] ) except that in @xmath103 and @xmath102 , @xmath36 is replaced with @xmath57 , and the outer integration domain is @xmath127 .
provided that @xmath128^{1/2}$ ] that we call case vi , the friction force is given as : @xmath129{\cal m } + [ r^c_2+r^f_2+{\cal m}_0(r^f_2-r^c_2)]{\cal m}^2+(r^f_2 -r^c_2){\cal m}^3}{2 r_0 { \cal m}^2({\cal m}+{\cal m}_0 ) } \nonumber \\ & & -\frac{2{\cal m}+r_0}{2{\cal m}^3}\,\log\left|\frac{2{\cal m}+r_0}{r_2^{c}{\cal m}-r_0}\right| - \frac{2{\cal m}-r_0}{2{\cal m}^3}\,\log\left|\frac{2{\cal m}-r_0}{r_2^{f}{\cal m}-r_0}\right| -\frac{2}{{\cal m}^2}\log\left|\frac{{\cal m}-{\cal m}_0}{{\cal m}+{\cal m}_0}\right|.\label{case6}\end{aligned}\ ] ] new terms depending the initial mach number appear in the expression of the force . as @xmath18
is set to zero the force reduces to ( [ case5 ] ) with @xmath130 .
this shows that including acceleration in the perturber s motion makes the medium store the information about the perturber s initial state which is in turn reflected in the reaction force felt by the perturber .
taking the limit @xmath89 , the force becomes : @xmath131 three extra terms appear with respect to the uniform velocity expression ( [ ostriker ] ) .
as @xmath132 , the force tends to that of case v ( [ case5l ] ) of perturbers initially at rest as the density perturbation of a larger mach number dominates the initial density perturbation .    for an accelerated perturber with an initial subsonic speed ,
the boundary conditions ( i , ii , iv ) show that the force expression is given as : @xmath133^{1/2}\leq { \cal m}\leq1+(2r_0)^{1/2 } : & & i = i_{\rm ii},\label{case61}\\ 1+(2r_0)^{1/2}\leq { \cal m}\leq 1+\left[(1-{\cal m}_0)^2 - 2r_0\right]^{1/2 } : & & i = i_{\rm iii},\label{case62}\\ 1+\left[(1-{\cal m}_0)^2 + 2r_0\right]^{1/2}\leq { \cal m } : & & i = i_{\rm vi},\label{case63}\end{aligned}\ ] ] where in the expressions of @xmath134 and @xmath135 , @xmath36 is replaced with @xmath57 .
we recover the general statement we made about the force on subsonic motion in case ii .
in particular , the friction force on an initially subsonic perturber depends on the initial velocity as soon as the velocity is supersonic and becomes larger than @xmath136 . in fig .
( 4 ) we show the friction force with @xmath137 for an initial subsonic velocity with @xmath138 using equations ( [ case61 ] ) , ( [ case62 ] ) and ( [ case63 ] ) , and an initial supersonic velocity with @xmath139 using ( [ case6 ] ) .
we remark that the dependence of dynamical friction on the initial velocity of the perturber is not an artifact of choosing to assemble the perturber at @xmath34 .
this starting condition can be replaced with @xmath140 and @xmath141 in turn can be made to tend to infinity thus approximating a steady state .
examination of the general force expression shows that @xmath142 is equivalent to the condition @xmath34 where the limits @xmath143 and @xmath144 are taken for accelerating and decelerating perturbers respectively .
for positively accelerating perturbers ( case iv ) , the force is independent of the initial mach number as @xmath145 ( [ case60 ] ) . in the next section ,
we show that for decelerating perturbers the dependence on @xmath18 remains as the velocity decays to zero .      for decelerating perturbers ( @xmath63 ) , we seek the dynamical friction force from @xmath34 until the velocity vanishes at @xmath146 .
the choice of a negative mach number instead of a negative acceleration makes the derivations easier but results in a positive friction force ( opposed to the direction of motion ) .
the phase @xmath123 is now negative and the boundary conditions that determine the integration domains of ( [ forcex ] ) are ( i ) , ( vi ) and ( v ) @xmath147 .
the latter condition shows that @xmath148 .
three cases need to be considered : @xmath149 , @xmath150 and @xmath151 . for subsonic motion ( @xmath149 ) and away from the transient regime ( i.e. @xmath152^{1/2}\leq { \cal m}$ ] ) , the integration domains are given as : @xmath153 as for case ii ( equations [ dom21 ] and [ dom22 ] ) , the latter domain ( [ decesub2 ] ) does not contribute to the force .
the corresponding integrand is proportional to @xmath154 and the sign argument is positive for @xmath155 .
the contribution of domain ( [ decesub1 ] ) is identical to that of case ii and @xmath156 ( [ case2 ] ) with @xmath36 replaced by @xmath57 . dynamical friction on a decelerating subsonic perturber is therefore independent of its initial velocity .
this confirms the general statement we made about subsonic motion in an infinite medium in case ii .    for supersonic perturbers with an initial mach number @xmath150
, we need to consider three cases depending on the value of the mach number @xmath57 . the first which we call case vii corresponds to @xmath157^{1/2}\leq { \cal m}\leq -1+(2r_0)^{1/2}$ ] whose boundaries are : @xmath158^{1/2 } , \label{decesup1}\\ r^f_1\leq r\leq { \cal m}-{\cal m}_0 , & \mbox{with } & \tau_{\rm min}=-1 . \label{decesup2}\end{aligned}\ ] ] domain ( [ decesup1 ] ) contributes a term similar to ( [ i0 ] ) with @xmath87 and @xmath159 substituted for @xmath36 and @xmath103 respectively . domain ( [ decesup2 ] ) contributes a term similar to ( [ i2 ] ) with @xmath159 and @xmath160 substituted for @xmath104 and @xmath102 respectively .
the force is then given as : @xmath161 { \cal m}+[r^c_1   + r^f_1 + { \cal m}_0 ( r^f_1-   r^c_1 ) ] { \cal m}^2 + ( r^f_1-r^c_1 ) { \cal m}^3}{2r_0{\cal m}^2({\cal m}+{\cal m}_0)}+\nonumber\\ & & + \frac{2{\cal m}+r_0}{2 { \cal m}^3}\,\log\left|\frac{2{\cal m}+r_0}{r_1^{f}{\cal m}+r_0}\right| + \frac{2{\cal m}-r_0}{2{\cal m}^3}\,\log\left|\frac{2{\cal m}-r_0}{r_1^{c}{\cal m}+r_0}\right| + \frac{2}{{\cal m}^2}\log\left|\frac{{\cal m}-{\cal m}_0}{{\cal m}+{\cal m}_0}\right| \label{case7}.\end{aligned}\ ] ] the second case which we call case viii concerns @xmath162^{1/2}$ ] .
the corresponding integration domains are : @xmath163^{1/2},\label{dom83}\\ r^f_1\leq r\leq { \cal m}-{\cal m}_0 , & \mbox{with } & \tau_{\rm min}=-1.\label{dom84}\end{aligned}\ ] ] the integration radii satisfy @xmath164 .
in contrast to case iii , this ordering shows that the second domain ( [ dom82 ] ) does not contribute a finite term to the force . domain ( [ dom81 ] ) contributes a term equal to @xmath134 , and domain ( [ dom83 ] ) a term similar to ( [ i2 ] ) where @xmath165 and @xmath159 are substituted for @xmath36 and @xmath103 respectively . the last domain ( [ dom84 ] ) contributes a term similar to that of ( [ decesup2 ] ) in case vii .
summing up all terms , we find : @xmath166{\cal m}}{2r_0{\cal m}^2({\cal m}+{\cal m}_0)}+\nonumber\\ & & + \frac{[r^c_1   + r^f_1 + r^f_2 - r^f_3 + { \cal m}_0 ( r^f_1 + r^f_2 - r^f_3 -   r^c_1 ) ] { \cal m}^2   + ( r^f_1 + r^f_2 - r^f_3-r^c_1){\cal m}^3}{2r_0{\cal m}^2({\cal m}+{\cal m}_0)}+\nonumber\\ & & + \frac{2{\cal m}+r_0}{2 { \cal m}^3}\,\log\left|\frac{2{\cal m}+r_0}{r_1^{f}{\cal m}+r_0}\right| + \frac{2{\cal m}-r_0}{2{\cal m}^3}\,\log\left|\frac{2{\cal m}-r_0}{r_3^{f}{\cal m}-r_0}\right| + \frac{2}{{\cal m}^2}\log\left|\frac{{\cal m}-{\cal m}_0}{{\cal m}+{\cal m}_0}\right|-\frac{2{\cal m}-r_0}{2{\cal m}^3}\,\log\left|\frac{r_1^c{\cal m}+r_0}{r_2^f{\cal m}-r_0}\right| . \label{case8}\end{aligned}\ ] ] the third case concerns @xmath152^{1/2}\leq { \cal m}\leq 0 $ ] whose boundaries turn out to be those of subsonic motion @xmath167 .
the corresponding force is therefore given by ( [ case2 ] ) where @xmath36 is replaced @xmath57 .
( 5 , left panel ) shows the force on a perturber with an initial velocity @xmath168 using equations ( [ case7 ] ) , ( [ case8 ] ) and ( [ case2 ] ) alongside the uniform motion friction force ( [ ostriker ] ) . in the supersonic regime
the uniform motion force overestimates the drag by about 10% at its maximum . in the subsonic regime
the force on accelerated motion departs significantly from the uniform velocity force before joining it at @xmath169 to assume the expression ( [ case2 ] ) . in particular , when the velocity vanishes so does the force .
lastly , for a supersonic perturber with an initial velocity @xmath151 , we find that : @xmath170^{1/2}\leq { \cal m}\leq -1+(2r_0)^{1/2 } : & & i = i_{\rm vii},\label{case8x } \\
-1+(2r_0)^{1/2}\leq { \cal m}\leq 0 : & & i = i_{\rm viii}.\label{case8x2}\end{aligned}\ ] ] fig . ( 5 , right panel ) shows the force on a perturber with an initial velocity @xmath171 using equations ( [ case8x ] ) and ( [ case8x2 ] ) alongside the uniform velocity friction force ( [ ostriker ] ) . in the supersonic regime
the uniform motion force overestimates the drag by about 20% at its maximum . in the subsonic regime
the force on accelerated motion departs significantly from that on uniform motion and is finite at @xmath172 . taking the limit where @xmath59 is small
, we find @xmath173 and : @xmath174 at @xmath172 , @xmath175 in contrast to the uniform motion force which vanishes with the perturber s velocity .
the knowledge of the history of the perturber s motion is therefore required to evaluate dynamical friction .
simply applying uniform motion dynamical friction force ( [ ostriker ] ) may yield incorrect results regarding the perturber s evolution in the ambient medium as evidenced in fig .
( 6 ) where we plot the relative error between the accelerated motion and uniform motion forces for decelerating supersonic perturbers .
a similar conclusion was previously reached by snchez - salcedo and brandenburg ( 1999 ) who numerically tested the accuracy of the uniform motion friction formula ( [ ostriker ] ) using the full nonlinear fluid equations . following the evolution of a decelerating initial perturber with @xmath176 up to @xmath177
, they found that equation ( [ ostriker ] ) underestimated the friction force as @xmath57 became smaller than unity ( see their fig .
5 ) which corresponds to our cases vii and viii .
the authors attributed this discrepancy to the ` memory of the wake ' .
in this paper , we studied the influence of acceleration on dynamical friction in a gaseous medium .
our main conclusions are stated in the abstract .
this work also constitutes an illustration of the use of the general expression of dynamical friction derived in appendix a. in particular , applying that formulation to a perturber with accelerated motion allowed us to illustrate how dynamical friction is not a local force and may depend the perturber s initial state .
the meaning of the initial state depends on the problem s configuration .
for instance , it may correspond to the time the perturber enters the gaseous medium or to the initial time of a numerical simulation of dynamical friction .
we note that non - locality is not peculiar to rectilinear motion and should be true of the dynamical friction associated with rotational ( circular or eccentric ) motion .
the main technical result is the set of expressions @xmath81 to @xmath178 that describe dynamical friction as a function of the initial and current velocities as the perturber is accelerated ( positively or negatively ) in the ambient medium .
our derivations have assumed a constant acceleration because the calculations are simpler to carry out in this case .
dynamical friction for constantly accelerating perturbers applies readily to the motion of young stars with asymmetric jet systems .
we believe that the general features we found will hold for a general type of acceleration as long as the instant velocity change remains small . including the deceleration generated by dynamical friction instead of a constant negative acceleration is likely to further reduce the friction force because the acceleration s amplitude increases with time .
the dependence of dynamical friction on the perturber s size , @xmath59 , is most important for supersonic motion . in the force expressions ,
the dimensionless radius , @xmath59 , stands in for the ratio @xmath179 of dimensional parameters .
the reaction force therefore depends explicitly on acceleration even in the limit @xmath180 as in case vi .
dynamical friction on subsonic motion is independent of @xmath59 when @xmath180 .
using the force expression requires the knowledge of the perturber s acceleration and the medium s sound speed . in fig .
( 3 ) to ( 5 ) , the value @xmath137 was used so as to make accelerated motion friction closer to uniform motion friction ( see fig .
the assumption of constant velocity in the derivation of dynamical friction is also made for collisionless star systems .
although they have different magnitudes , the uniform motion friction forces in a gaseous medium and a collisionless star system behave similarly ( as @xmath181 for supersonic motion and @xmath182 for subsonic motion ; see ostriker 1999 , snchez - salcedo and brandenburg 2001 ) .
this similarity suggests that the conclusions we found are likely to hold for collisionless systems .
the author thanks an anonymous referee and scott tremaine for useful comments .
in this appendix , we derive a general expression of dynamical friction on a perturber of gravitational potential @xmath183 with a general trajectory @xmath54 that travels inside an inviscid homogeneous gaseous medium .
the perturber is set in motion at time @xmath34 .
the medium s response to the perturber s presence is described by the standard fluid equations : @xmath184 where @xmath185 and @xmath186 are the perturbed velocity and density fields , @xmath187 is the pressure , and @xmath10 is the heaviside function .
prior to the perturber s motion , the gaseous medium is assumed to be at rest @xmath188 $ ] , homogeneous of density @xmath5 , pressure @xmath189 and isothermal with a sound speed @xmath190 .
the linearization of medium s equation in the vicinity of the rest state yields : @xmath191 where we used @xmath7 and @xmath192 .
the combined equations in turn yield the forced sound wave equation satisfied by the density perturbation @xmath6 as : @xmath193 the solution of this equation is known in electromagnetism theory as the linard - wiechert potential ( jackson 1998 , chapter 14 ) and represents the electric potential of a moving electric charge distribution .
it is written in integral form using the retarded green function as : @xmath194 \
\nabla^2\phi_p[{\bm x}^\prime - { \bm \xi}(u)]\ { \cal h}(u ) } { |{\bm x}-{\bm x}^\prime|}. \label{densa}\end{aligned}\ ] ] in reaction to the perturber s excitation of the wake , the gaseous medium exerts a force on the perturber given as : @xmath195\   \nabla^2\phi_p[{\bm x}^\prime - { \bm \xi}(u)]\ { \bm \nabla } \phi_p[{\bm x } - { \bm \xi}(t)]\ { \cal h}(u ) } { |{\bm x}-{\bm x}^\prime|}\label{forc1a}\ ] ] the symbol @xmath196 denotes the boundaries of the spacial integration domain defined by the medium s physical boundary as well as the lower distance cutoff in the perturber s vicinity such as its size or its corresponding accretion radius .
the boundaries @xmath196 are used to truncate the perturber s potential at small and large distances to avoid divergence . applying the variable changes , @xmath197 and @xmath198 and integrating the force over @xmath199 yields the general expression of dynamical friction as : @xmath200 , |{\bm y}|\leq ct}{\rm d}^3x\frac{{\rm d}^3y}{|{\bm y}| }   \   \ \rho_p({\bm x}-{\bm y}+{\bm \delta})\ { \bm \nabla } \phi_p({\bm x } ) , \label{forc2a } \end{aligned}\ ] ] where @xmath201 is the perturber s density distribution , @xmath202 , and @xmath203
. for a given perturbing potential , dynamical friction therefore depends only on the phase @xmath204 as well as the boundary conditions specified by @xmath205 and @xmath206 .
the dependence of @xmath204 on the whole trajectory @xmath54 shows that dynamical friction is manifestly non - local . deriving dynamical friction through equation ( [ forc2a ] )
is in contrast with the previous two - step methods of estimating dynamical friction for uniform motion by first calculating explicitly ( [ densa ] ) then applying an analytical or numerical calculation to the first integral in equation ( [ forc1a ] ) to derive the force ( ostriker 1999 , kim and kim 2007 , kim et al 2008 ) . as a result
such methods obscure the true dependence of the force on the perturber s motion and the boundary conditions . for a point - like perturber , @xmath207 and @xmath208 , where @xmath0 is the perturber s mass .
the forced sound wave equation , density enhancement and friction force are reduced to : @xmath209 { \cal h}(t ) , \label{wav2a}\\ \varrho({\bm x},t ) & = & \frac{gm}{c^2}\int_{-\infty}^{+\infty } \frac{\delta \left[u - t+|{\bm x}-{\bm \xi}(u)|/c\right ] { \cal h}(u ) } { |{\bm x}-{\bm \xi}(u)|}\ { \rm d}u , \label{dens2a}\\ { \bm f}&=&\frac{{\cal h}(t)(gm)^2\rho_0}{c^2}\int_{\partial v[{\bm y}+{\bm \xi}(t - r / c ) ] , r\leq ct}r\sin\theta{\rm d}r\ , { \rm d}\theta\ , { \rm d}\varphi   \ \frac{{\bm y}-{\bm \delta}}{|{\bm y}-{\bm \delta}|^3},\label{for3a }   \end{aligned}\ ] ] where @xmath45 .
the condition on the time of perturbation ( @xmath210 ) represented by the heaviside function in equation ( [ wava ] ) and the expression ( [ forc1a ] ) now implies that the spacial domain that enters the expression of the force is equal to the size of the wavefront , @xmath211 .
however , the spacial domain refers to the relative position with respect to the perturber not at time @xmath36 but at the retarded time @xmath44 .
this explains why the force ( [ for3a ] ) describes subsonic as well as supersonic motion that perturbs the medium outside the sonic wavefront of radius @xmath211 . in the force expression ,
the boundaries @xmath212 $ ] are treated as truncation boundaries . in the next appendix ,
the force expression is applied to a point - like perturber in uniform rectilinear motion to recover the corresponding known results .
we validate the expressions ( [ forc2a ] ) and ( [ for3a ] ) by deriving dynamical friction for a perturber set in uniform rectilinear motion in an infinite homogeneous medium .
the perturber s trajectory is written as @xmath213 where @xmath50 is the unit vector along @xmath13 , the phase @xmath214 where @xmath57 is the mach number .
the friction force is equally along @xmath50 and its expression is simplified as follows : @xmath215 , r\leq ct}r\sin\theta{\rm d}r\ , { \rm d}\theta\ , { \rm d}\varphi   \
\frac{r\cos \theta-{\delta}}{|r^2+\delta^2 - 2r\delta\cos\theta|^{3/2}},\\ f&=&\frac{2\pi { \cal h}(t)(gm)^2\rho_0}{v^2}\int_{\partial v[{\bm y}+{\bm \xi}(t - r / c ) ] , r\leq ct } \ \frac{{\rm d}r } { r } \left[\frac{1-{\cal m}\tau}{({1+{\cal m}^2 - 2 { \cal m}\tau})^{1/2}}\right]_{\tau_{\rm min}}^{\tau_{\rm max } } , \label{forc3a}\end{aligned}\ ] ] where in the last expression we performed the integrals over @xmath42 and @xmath52 . the extremal values @xmath216 and @xmath56 are determined from the boundary conditions
. calling @xmath59 the perturber s size , the force is truncated at small distances by @xmath217 . expressing these conditions in terms of the variables @xmath40 and @xmath70 that appear in equation ( [ forc3a ] ) shows that we need to find the domains that satisfy : ( i ) @xmath218 , ( ii ) @xmath219 , and ( iii ) @xmath49 .
the first two conditions yield the following integration domains : @xmath220 the inner domain ( [ dom1a ] ) corresponds to the vicinity of the perturber s retarded position .
the radius @xmath221 defines a region around the perturber that does not contribute to dynamical friction as the first two conditions are not satisfied .
the third condition depends on the distance @xmath211 travelled by the retarded wavefront .
there are two force regimes : the local perturbation regime ( @xmath222 ) where the perturbation remains in the perturber s vicinity and the distant perturbation regime ( @xmath223 ) where the perturbation has expanded beyond the perturber s vicinity .
calling @xmath79 the integral that appears in ( [ forc3a ] ) , the integration for the distant perturbation regime extends over the two domains ( [ dom1a ] ) and ( [ dom2a ] ) with @xmath224 and @xmath79 is given as : @xmath225\ { \rm d}r-\left[1-{\rm sign } ( 1-{\cal m})\right]\int_{\frac{r_0}{|1-{\cal m}|}}^{ct}\frac{{\rm d}r}{r},\\i&=&-\log \left|\frac{1+{\cal m}}{1-{\cal m}}\right|+\frac{1}{2}\left[1+{\cal m}-|1-{\cal m}|\right]\left[1+{\rm sign}(1-{\cal m})\right ] -\left[1-{\rm sign } ( 1-{\cal m})\right ] \
\log\left|\frac{(1-{\cal m})ct}{r_0}\right| .
\label{ostria}\end{aligned}\ ] ] this expression agrees with those derived by ostriker ( 1999 ) for subsonic and supersonic perturbers moving in an infinite gaseous medium using explicit expressions of the density perturbation . in particular , the force , @xmath226 , is time - independent and behaves as @xmath227 for small mach numbers .
we note the interesting result that only the vicinity of a subsonic perturber s retarded position ( [ dom1a ] ) contributes to dynamical friction .
we show in section 3 that this property implies that the force on subsonic perturbers with a general type of motion ( not necessarily rectilinear or uniform ) is similar to ( [ ostria ] ) provided that the perturber s motion history does not involve supersonic episodes and that the perturber s acceleration @xmath228 is much smaller than @xmath229 .    for the local perturbation regime ( @xmath222 )
, the force is zero before @xmath230 and only the domain ( [ dom1a ] ) is relevant to the integral @xmath79 which is written as : @xmath231\ { \rm d}r= -\log\left[\frac{ct(1+{\cal m})}{r_0}\right]+{\cal m}-\frac{r_0}{2ct}+\frac{(1-{\cal m}^2)ct}{2r_0}. \label{transient}\end{aligned}\ ] ] in fig .
( 7 ) , we show the friction force as a function of @xmath36 for subsonic and supersonic perturbers as @xmath36 increases smoothly from 0 to a value in the distant perturbation regime .
the time dependence of the force agrees with that obtained from the numerical simulation of the nonlinear fluid equations ( snchez - salcedo and brandenburg 1999 , figure 2 ) .
however , the attempted fit of the force in the local perturbation regime from numerical simulations ( snchez - salcedo and brandenburg 2001 , equation 22 ) does not reproduce the correct dependence of ( [ transient ] ) on @xmath57 , @xmath59 and time .                                                           from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ] from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ] +   from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ] from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ] +   from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ]   from a perturber initially at rest .
the level curves are determined for values of @xmath232 .
top left panel : subsonic motion , ( @xmath233 , @xmath234 ) .
the density increases nearer the perturber .
note how the level curves do not possess front - back symmetry with respect to the perturber .
top right panel : motion at mach 1 .
middle left panel : supersonic motion inside the sonic wavefront .
the solid level curves correspond to ( @xmath235 , @xmath234 ) and ( @xmath236 , @xmath237 ) near the perturber s mach cone .
the dashed level curves correspond to a shallow depression around @xmath238 .
middle right panel : the perturber reaches the sonic shockwave .
the solid level curves correspond to ( @xmath235 , @xmath234 ) .
the dashed level curves correspond to shallow depressions around @xmath239 inside the sonic wavefront and @xmath240 inside the mach cone .
bottom left panel : supersonic motion at @xmath241 .
the solid level curves correspond to ( @xmath242 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath243 .
bottom right panel : supersonic motion at @xmath244 .
the solid level curves correspond to ( @xmath245 , @xmath234 ) .
the dashed level curves correspond to a shallow depression around @xmath246 . for supersonic motion , larger values of @xmath247
are confined to the sides of the mach cone where the discontinuity with the unperturbed medium is largest .
, title="fig:",width=245 ]     from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity .
, title="fig:",width=245 ] from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity .
, title="fig:",width=245 ] +   from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity .
, title="fig:",width=245 ] from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity .
, title="fig:",width=245 ] +   from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity .
, title="fig:",width=245 ]   from decelerating perturbers .
the level curves are shown for a subsonic perturber with @xmath248 ( left column ) and a supersonic perturber with @xmath249 ( right column ) .
note how front - back symmetry is lost everywhere for the subsonic perturber .
this however does not significantly modify dynamical friction from its value obtained using a constant velocity . , title="fig:",width=245 ]    , @xmath118 , and @xmath119 .
the solid curves correspond to the friction force that takes acceleration into account given by equations ( [ case1 ] ) , ( [ case2 ] ) , ( [ case3 ] ) , ( [ case4 ] ) and ( [ case5 ] ) , and the dashed curves to the uniform motion friction force expression ( [ ostriker]).,width=321 ]    ) and an initial supersonic velocity ( right panel , @xmath139 ) .
the solid curves correspond to the friction force that takes acceleration into account given by equations ( [ case6]),([case61 ] ) , ( [ case62 ] ) and ( [ case63 ] ) , and the dashed curves to the uniform motion friction force ( [ ostriker ] ) . the small cutoff radius is @xmath137 .
the dotted curve in the left panel is the friction force on a perturber initially at rest.,title="fig:",width=321 ] ) and an initial supersonic velocity ( right panel , @xmath139 ) .
the solid curves correspond to the friction force that takes acceleration into account given by equations ( [ case6]),([case61 ] ) , ( [ case62 ] ) and ( [ case63 ] ) , and the dashed curves to the uniform motion friction force ( [ ostriker ] ) . the small cutoff radius is @xmath137 .
the dotted curve in the left panel is the friction force on a perturber initially at rest.,title="fig:",width=321 ]     ( left panel ) , @xmath250 ( right panel ) .
the solid curves correspond to the friction force that takes acceleration into account given by equations ( [ case2 ] ) , ( [ case7 ] ) and ( [ case8 ] ) for @xmath251 and by equations ( [ case8x ] ) and ( [ case8x2 ] ) for @xmath250 .
the dashed curves correspond to the uniform motion friction force ( [ ostriker ] ) .
the small distance cutoff is @xmath137.,title="fig:",width=321 ] ( left panel ) , @xmath250 ( right panel ) .
the solid curves correspond to the friction force that takes acceleration into account given by equations ( [ case2 ] ) , ( [ case7 ] ) and ( [ case8 ] ) for @xmath251 and by equations ( [ case8x ] ) and ( [ case8x2 ] ) for @xmath250 .
the dashed curves correspond to the uniform motion friction force ( [ ostriker ] ) .
the small distance cutoff is @xmath137.,title="fig:",width=321 ]     from equations ( [ ostriker ] ) and ( [ case7 ] ) , on the friction force when using the uniform motion expression instead of the accelerated motion expressions for decelerating perturbers with initial supersonic speeds .
the small distance cutoff is @xmath137 .
the curves correspond to different initial velocities.,width=321 ]    , is given by ( [ forc3a ] ) , ( [ ostria ] ) , and ( [ transient ] ) and scaled to @xmath252 ; time is scaled to @xmath253 .
the perturber s motion is initiated at @xmath34 .
the mach number is indicated to the left of each curve .
the time dependence is similar to that obtained from the simulation of the non - linear fluid equations in figure 2 of snchez - salcedo and brandenburg ( 1999).,width=283 ]