No Arabic abstract
We consider the gravitational force exerted on a point-like perturber of mass $M$ travelling within a uniform gaseous, opaque medium at constant velocity $V$. The perturber irradiates the surrounding gas with luminosity $L$. The diffusion of the heat released is modelled with a uniform thermal diffusivity $chi$. Using linear perturbation theory, we show that the force exerted by the perturbed gas on the perturber differs from the force without radiation (or standard dynamical friction). Hot, underdense gas trails the mass, which gives rise to a new force component, the heating force, with direction $+V$, thus opposed to the standard dynamical friction. In the limit of low Mach numbers, the heating force has expression $F_mathrm{heat}=gamma(gamma-1)GML/(2chi c_s^2)$, $c_s$ being the sound speed and $gamma$ the ratio of specific heats. In the limit of large Mach numbers, $F_mathrm{heat}=(gamma-1)GML/(chi V^2)f(r_mathrm{min}V/4chi)$, where $f$ is a function that diverges logarithmically as $r_mathrm{min}$ tends to zero. Remarkably, the force in the low Mach number limit does not depend on the velocity. The equilibrium speed, when it exists, is set by the cancellation of the standard dynamical friction and heating force. In the low Mach number limit, it scales with the luminosity to mass ratio of the perturber. Using the above results suggests that Mars- to Earth-sized planetary embryos heated by accretion in a gaseous protoplanetary disc should have eccentricities and inclinations that amount to a sizeable fraction of the discs aspect ratio, for conditions thought to prevail at a few astronomical units.
We present a rigorous calculation of the dynamical friction force exerted on a spherical massive perturber moving through an infinite homogenous system of field stars. By calculating the shape and mass of the polarization cloud induced by the perturber in the background system, which decelerates the motion of the perturber, we recover Chandrasekhars drag force law with a modified Coulomb logarithm. As concrete examples we calculate the drag force exerted on a Plummer sphere or a sphere with the density distribution of a Hernquist profile. It is shown that the shape of the perturber affects only the exact form of the Coulomb logarithm. The latter converges on small scales, because encounters of the test and field stars with impact parameters less than the size of the massive perturber become inefficient. We confirm this way earlier results based on the impulse approximation of small angle scatterings.
In this paper, we explore the impact of a galactic bar on the inspiral time-scale of a massive perturber (MP) within a Milky Way-like galaxy. We integrate the orbit of MPs in a multi-component galaxy model via a semi-analytical approach including an accurate treatment for dynamical friction generalized to rotationally supported backgrounds. We compare the MP evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. We find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time, implying that dynamical friction alone is greatly insufficient to fully characterize the orbital decay. The effect of the bar is more prominent for initially in-plane, prograde MPs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MPs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated to the non-symmetric mass distribution is thus not advisable even within our idealized, smooth Milky Way model, and it should be avoided when dealing with more complex and realistic galaxy systems. This has important implications for the orbital decay of massive black holes in late-type spirals, the natural candidate sources to be detected with the Laser Interferometer Space Antenna (LISA).
We compute the dynamical friction on a small perturber moving through an inviscid fluid, i.e., a superfluid. Crucially, we account for the tachyonic gravitational mass for sound waves, reminiscent of the Jeans instability of the fluid, which results in non-zero dynamical friction even for subsonic velocities. Moreover, we illustrate that the standard leading order effective theory in the derivative expansion is in general inadequate for analysing supersonic processes. We show this in two ways: (i) with a fluid treatment, where we solve the linearized hydrodynamical equations coupled to Newtonian gravity; and (ii) with a quasiparticle description, where we study the energy dissipation of a moving perturber due to phonon radiation. Ordinarily a subsonic perturber moving through a superfluid is kinematically prohibited from losing energy, however the Jeans instability modifies the dispersion relation of the fluid which can result in a small but non-vanishing dynamical friction force. We also analyse the soft phonon bremsstrahlung by a subsonic perturber scattered off an external field.
In many astrophysical situations, as in the coalescence of supermassive black hole pairs at gas rich galactic nuclei, the dynamical friction experienced by an object is a combination of its own wake as well as the wakes of its companions. Using a semi-analytic approach, we investigate the composite wake due to, and the resulting drag forces on, double perturbers that are placed at the opposite sides of the orbital center and move on a circular orbit in a uniform gaseous medium. The circular orbit makes the wake of each perturber asymmetric, creating an overdense tail at the trailing side. The tail not only drags the perturber backward but it also exerts a positive torque on the companion. For equal-mass perturbers, the positive torque created by the companion wake is, on average, a fraction ~40-50% of the negative torque created by its own wake, but this fraction may be even larger for perturbers moving subsonically. This suggests that the orbital decay of a perturber in a double system, especially in the subsonic regime, can take considerably longer than in isolation. We provide the fitting formulae for the forces due to the companion wake and discuss our results in light of recent numerical simulations for mergers of binary black holes.
We use three-dimensional hydrodynamic simulations to investigate the nonlinear gravitational responses of gas to, and the resulting drag forces on, very massive perturbers moving on circular orbits. This work extends our previous studies that explored the cases of low-mass perturbers on circular orbits and massive perturbers on straight-line trajectories. The background medium is assumed to be non-rotating, adiabatic with index 5/3, and uniform with density rho0 and sound speed a0. We model the gravitating perturber using a Plummer sphere with mass Mp and softening radius rs in a uniform circular motion at speed Vp and orbital radius Rp, and run various models with differing R=rs/Rp, Mach=Vp/a0, and B=G*Mp/(a0^2*Rp). A quasi-steady density wake of a supersonic model consists of a hydrostatic envelope surrounding the perturber, an upstream bow shock, and a trailing low-density region. The continuous change in the direction of the perturber motion makes the detached shock distance reduced compared to the linear-trajectory cases, while the orbit-averaged gravity of the perturber gathers the gas toward the center of the orbit, modifying the background preshock density to rho1=(1 + 0.46B)*rho0 depending weakly on Mach. For sufficiently massive perturbers, the presence of a hydrostatic envelope makes the drag force smaller than the prediction of the linear perturbation theory, resulting in F = 4*pi*rho1*(G*Mp/Vp)^2 * (0.7/etaB) for etaB = B/(Mach^2 -1) > 0.1; the drag force for low-mass perturbers with etaB < 0.1 agrees well with the linear prediction. The nonlinear drag force becomes independent of R as long as R < etaB/2, which places an upper limit on the perturber size for accurate evaluation of the drag force in numerical simulations.