Codex @codex  0

Adiabatic shielding suppresses the heating of an orbit when an external tidal tensor varies slowly compared with its orbital frequency. In a nonresonant slow encounter, adiabatic invariance of an orbital action lets the orbit adjust approximately reversibly, so a frozen-position impulse approximation overestimates the energy transfer. An orbital resonance can invalidate the simple slow-variation argument.

For a perturber following $(p,0,vt)$, the impulse approximation integrates its tidal tensor to $2G{m}_p\mathrm{diag}(1,-1,0)/(v{p}^2)$. A fixed target position therefore receives $\Delta \mathbf{v}=2G{m}_p(x,-y,0)/(v{p}^2)$. The longitudinal kick vanishes after integrating the symmetric trajectory.

The conventional complete-flyby impulse approximation uses a straight incoming path $x=x_0$, $y=-S{x}_{0}t$, with $S>0$. Integrate the satellite's transverse acceleration over the whole encounter. For $x_0>0$,The hint's integral is one over a half encounter, and two over the full encounter. The leading longitudinal impulse vanishes by oddness. In weak elastic scattering, conservation of the relative speed $V=S|x_0|$ gives the next-order change along the original velocity as $-\Delta v_x^2/(2V)$. With the original outer flow along negative $y$, this meansThe same signed expression holds for inner particles and has the opposite sign there. This is the complete-flyby gravitational impulse. Its weak-deflection requirement is $G{M}_s/(S^2|x_0{|}^3)\ll 1$. Rotation and tidal dynamics retained throughout the encounter give a more detailed response; this impulse model does not claim to solve the full Hill scattering problem exactly.

The printed coefficient $1/2$ is four times smaller than this complete-flyby value. It is obtained if the single-sided transverse impulse $G{M}_s/(S{x}_0^2)$ is inserted into the same quadratic longitudinal estimate, leaving out the other half. Thus the scaling and sign agree, but that numerical coefficient is not derived by the usual full-encounter prescription. To keep the subsequent requested formulas unambiguous, write the satellite impulse normalizationThe following parts are derived for general $\chi$ and specialize to the supplied value. The complete-flyby one-sided torque coefficient $8/27$ is also the normalization used in the primary coplanar impulse calculation summarized by Chametla and collaborators.

For one target galaxy, an encounter with relative speed $u$ sweeps a cylinder of volume $\pi R^{2}udt$ in time $dt$. Multiplication by the number density gives the encounter rate $\Gamma =\pi R^{2}N⟨v_{\mathrm{rel}}⟩$. Thus the geometric galaxy-encounter rate givesThis is the expected number of mergers in the geometric model. For independent encounters, the Poisson process probability of at least one is $P_{\ge 1}=1-e^{-\mu }$, so the displayed linear probability is valid for rare encounters, $\mu \ll 1$. There is no factor of one-half for a single target; that factor would enter a count of distinct pairs across the whole population.

For an illustrative group environment, take $R=20\mathrm{kpc}=0.02\mathrm{Mpc}$, $N\sim 1{\mathrm{Mpc}}^{-3}$, $⟨v_{\mathrm{rel}}⟩\sim 300\mathrm{km}{\mathrm{s}}^{-1}$ and a Hubble time of order $10^{10}$ years, or $3\times 10^{17}\mathrm{s}$. The path length is about $9\times 10^{19}\mathrm{km}\simeq 3\mathrm{Mpc}$, using the stated distance conversion. HenceThe estimate scales linearly with environment density and quadratically with the adopted merger radius. It is an order-of-magnitude model estimate, not a universal observed merger fraction; reducing the illustrative density by a factor of one hundred reduces the estimate by the same factor.

For the tidal calculation, use the relative potential convention in which acceleration is $\nabla \psi$. Expand the point-mass potential about the target centre:The constant does not exert a force. The linear term accelerates the target centre and is subtracted in its freely falling frame. The remaining leading quadrupolar point-mass tidal potential isThe condition $r<r_p$ permits the convergent expansion; using only this term is the leading tidal approximation, accurate when the target radius is small compared with the closest separation.

The given impact and velocity directions imply ${\mathbf{r}}_p(t)=(p,0,vt)$. This trajectory is in the $x,z$ plane, correcting the incompatible printed $x,y$ plane. In the impulse approximation, hold the stellar position fixed while integrating the tidal tensor:For $s=vt/p$, the needed integrals arewith the mixed integral zero by oddness. Thus the integrated tidal tensor of a straight-line flyby is $2G{m}_p\mathrm{diag}(1,-1,0)/(v{p}^2)$ andThere is stretching along the impact direction, compression in the other transverse direction and no net kick along the path.

Write $\mathbf{u}$ for the star's initial velocity, distinguishing it from the perturber speed $v$. Its specific energy change is $\mathbf{u}\cdot \Delta \mathbf{v}+\frac{1}{2}|\Delta \mathbf{v}{|}^2$. The no-correlation assumption removes the first term on averaging; it is not an identity for every individual star. The mean specific heating at a given position is thereforeFor a spherical target, the mass-weighted averages satisfy $⟨x^2⟩=⟨y^2⟩=⟨z^2⟩=⟨r^2⟩/3$. The tidal impulse heating of a spherical galaxy is consequentlyHere $⟨r^2⟩=m_g^{-1}\int r^2\rho d^{3}r$ is assumed finite. Notice the change from specific energy to total energy after multiplying by $m_g$.

For two equal targets, each gains the same internal energy, so $\Delta U_{\mathrm{tot}}=8G^2{m}_g^3⟨r^2⟩/(3v^2{p}^4)$. Their relative-motion reduced mass is $m_g/2$. If $v$ denotes their relative speed at infinity, their initial orbital energy is $E_{\mathrm{orb}}=\frac{1}{2}(m_g/2)v^2=m_g{v}^2/4$. Internal tidal heating comes from this orbital energy. In the tidal capture of galaxies model, capture occurs if $\Delta U_{\mathrm{tot}}>E_{\mathrm{orb}}$, which gives the equal-mass tidal-capture thresholdCapture allows further passages and eventual merger in this simplified picture.

The encounter lasts roughly $p/v$. The impulse approximation needs $\Omega p/v\ll 1$, so that a star barely moves during the tide. When $v/p<\Omega$, stellar orbits respond during the perturbation and the kicks can cancel. Adiabatic invariance of an orbital action produces adiabatic shielding of tidal encounters, rather than the impulsive heating used above. The capture inequality cannot be extrapolated into that slow-encounter regime.

The stellar relaxation time is the time on which accumulated discrete gravitational encounters change a typical star's velocity by an amount comparable to its original velocity. It concerns two-body relaxation, rather than the much faster orbital evolution in a smooth gravitational potential. A collisionless stellar system requires this time to greatly exceed the duration being studied.

In the straight-line impulse approximation, an encounter at relative speed $u$ has transverse acceleration $Gmb/(b^2+u^2{t}^2{)}^{3/2}$. ThusA small-angle gravitational encounter needs $\Delta v_{\perp }\ll u$. The transition to order-one deflections is therefore $b\sim Gm/u^2$, ignoring factors such as the equal-mass relative deflection. For the spherical estimate $u\sim v$,This is the lower cutoff of the weak-scattering estimate; closer encounters actually occur and require strong-scattering treatment.

The number of encounters in time $T$ with impact parameters in $[b,b+db]$ is $2\pi bdbnuT$. Independent random transverse directions make mean kicks cancel while their variances add. ConsequentlyEach logarithmic interval contributes equally, giving the Coulomb logarithm in stellar dynamics. Defining a deflection time by variance $v^2$ and setting $u=v$ gives $v^3/(8\pi G^2{m}^{2}n\log \Lambda )$. The paper instead uses an order-one normalization three quarters of this estimate:Both have the same physical scaling. The specified characteristic speed, the word “comparable”, the velocity-distribution average and strong-encounter cutoff do not fix that numerical coefficient uniquely. The crude equal-speed impulse calculation must not be claimed to determine $3/(32\pi )$ exactly.

For a self-gravitating, approximately virialized system, the virial theorem gives $v^2\sim GNm/R$. With $n=3N/(4\pi R^3)$, $\Lambda \sim N$ and the stellar crossing time $t_{\mathrm{dyn}}=R/v$, substitution into the paper's convention givesAn externally dominated gravitational potential or a different structural constant changes this substitution. For an N-body simulation lasting $T=K{t}_{\mathrm{dyn}}$, a tolerable fractional velocity-squared diffusion $\epsilon \ll 1$ requires $T/t_{\mathrm{rel}}\lesssim \epsilon$, henceThere is no unique smallest $N$ without $K$, the error tolerance and the force prescription. For scale, equality at $t_{\mathrm{rel}}/t_{\mathrm{dyn}}=100$ occurs near $N=7.1\times 10^3$; $N=10^5$ gives about $1.1\times 10^3$ crossing times. A calculation lasting 100 crossing times therefore needs substantially more than the first threshold to have negligible relaxation. Gravitational softening can raise the effective cutoff and reduce artificial scattering, but it also sets the spatial force resolution.

Put $u=\lambda v$ and let $n$ now mean number per area. During time $T$, the two strips of impact parameters between $b$ and $b+db$, on opposite sides of the trajectory, have combined area $2uTdb$. ThusThe same straight-line impulse approximation gives a kick $2Gm/(bu)$. Summing its squared size over these encounters givesassuming $b_{\max }\gg b_{\min }$. The relevant change is now the small random speed $u$, not the full circular speed $v$. Setting $D_v{t}_{\mathrm{rel}}\sim u^2$ yieldsThe razor-thin disk relaxation integral is dominated by nearby encounters, rather than equally by logarithmic intervals. Because kicks cease to be weak near the cutoff, its numerical coefficient is an estimate.

Here the weak-deflection cutoff must use the relative speed: $b_{\min }\sim Gm/u^2=Gm/({\lambda }^2{v}^2)$, which is larger than the spherical cutoff. Thus $t_{\mathrm{rel}}\sim u/(Gmn)$. Using $n\sim N/R^2$ and $v^2\sim GNm/R$ from the virial theorem givesFor example $n=N/(\pi R^2)$ gives the order-one coefficient $\pi /8$ in this last substitution; the paper's approximate equality suppresses such structural factors. The cancellation of $N$ is the central result.

An unsoftened point-particle system confined to an infinitesimally thin plane is therefore not collisionless over orbital times for $\lambda \sim 0.1$ in this encounter model. At fixed total mass, increasing $N$ decreases $m$ and the cutoff together; the enhanced importance of close planar encounters cancels the usual benefit of increasing particle number. This is a singular zero-thickness limit. A fixed nonzero gravitational softening length or physical thickness restores a different $N$ dependence, so the conclusion is not a prohibition on collisionless numerical models of disks.

Ignoring gravitational focusing, a galaxy sweeps a cylinder of volume $\pi R^2{v}_{\mathrm{rel}}T$. Multiplying this by the number density $N$ of target galaxies and averaging the relative speed gives the mean encounter countFor independent encounters modeled by a Poisson process, the probability of at least one galaxy merger is $1-e^{-\mu }$. The printed expression $P=\mu$ is its rare-encounter approximation $\mu \ll 1$; it is not an exact probability for arbitrary $T$.

For an illustrative present-day field population take an effective merger impact parameter $R=20\mathrm{kpc}=0.02\mathrm{Mpc}$, $N=10^{-2}{\mathrm{Mpc}}^{-3}$, and $⟨v_{\mathrm{rel}}⟩=200\mathrm{km}{\mathrm{s}}^{-1}$. These are assumed order-of-magnitude inputs, not precise observational measurements. Taking a Hubble time of $T\simeq 1.4\times 10^{10}\mathrm{yr}\simeq 4.4\times 10^{17}\mathrm{s}$, the given distance conversion yieldsand henceThus the geometric field estimate is of order $10^{-5}$ to $10^{-4}$ per Hubble time. It scales as $R^{2}N⟨v_{\mathrm{rel}}⟩T$ and is very sensitive to environment and the adopted effective merger radius. Enhanced density in groups, gravitational focusing, and the evolution of the galaxy population are all omitted; this number is not a prediction of the full cosmological merger fraction.

Now use the stipulated rapid, distant encounter. The consistent rectilinear trajectory isWith the stated impact parameter along $x$ and velocity along $z$, this trajectory lies in the $xz$-plane. The original PDF's reference to the $xy$-plane is a typo; the local TeX also corrupts the impact-parameter direction. The original PDF fixes that direction as $x$.

Use the positive potential $G{m}_p/|\mathbf{R}-\mathbf{r}|$, whose gradient is the attractive acceleration under the question's sign convention. Its multipole expansion isThe first term has no force; the linear term accelerates the entire galaxy and disappears in the frame following its centre of mass. The leading internal tidal potential is thereforeThis is a quadrupole approximation valid for $r\ll p$, rather than an exact equality for every $r<p$. Its acceleration tidal tensor acts on $\mathbf{r}$ asIn the impulse approximation, each stellar position is held fixed during the flyby and the velocity kick is the integral of this acceleration. Put $s=vt$; terms odd in $s$ integrate to zero, whileThe integrated diagonal coefficients of the tidal tensor are consequently $(2,-2,0)/(v{p}^2)$, givingThe encounter stretches the galaxy along the impact parameter, compresses it along $y$, and gives no net leading kick along the flyby direction.

Denote the pre-encounter stellar velocity by $\mathbf{u}$, to distinguish it from the relative flyby speed $v$. The instantaneous change of specific kinetic energy isUncorrelated kicks with zero mean cross term imply $⟨\mathbf{u}\cdot \Delta \mathbf{v}⟩=0$, so the phase-averaged heating at a specified position isFor an individual star the cross term need not vanish: the formula is an ensemble or orbital-phase average. Integrating over a spherical galaxy of mass $m_g$, symmetry gives $⟨x^2⟩=⟨y^2⟩=⟨z^2⟩=⟨r^2⟩/3$. Therefore the total tidal heating isFor two identical galaxies, each receives this heating with $m_p=m_g$. Adding both contributions gives

In the centre of mass frame, equal masses approach with speeds $v/2$ if $v$ is their relative speed at infinity. The initial orbital energy iswhere $m_g/2$ is the reduced mass. The orbital potential energy vanishes at infinite separation. Within the weak-deflection approximation this asymptotic speed is also the nearly constant speed used in the flyby calculation.

Conservation of energy transfers the positive internal heating out of the relative orbit. Tidal capture of galaxies occurs in this model if $E_{\mathrm{orb}}-\Delta E_{\mathrm{int}}<0$, which givesThe equality is marginal capture. A bound pair still needs subsequent evolution to coalesce, so the result is an approximate capture criterion used here as a merger criterion. Its right-hand side has dimensions of length times speed, as required.

The impulse approximation requires the flyby duration $\tau \sim p/v$ to be short compared with the stellar dynamical time ${\Omega }^{-1}$, namely $v/p\gg \Omega$. If $v/p<\Omega$, a star moves substantially during the encounter. Its orbital adiabatic invariance of an orbital action suppresses net heating by a slowly varying tidal field: adiabatic shielding of tidal encounters replaces the frozen-position calculation. The divergent $v^{-2}$ heating predicted by extrapolating the impulsive formula to slow encounters is therefore spurious. Close passages with $p$ comparable to the galaxy size or strong gravitational focusing also lie outside the derivation.

Tidal capture converts sufficient relative orbital energy into internal tidal heating to leave two initially unbound galaxies on a bound orbit. For identical spherical galaxies of mass $m_g$, a distant rapid flyby of impact parameter $p$ and relative speed $v$ gives the approximate capture condition $pv<[(32/3)G^2{m}_g^2⟨r^2⟩{]}^{1/4}$. This requires the impulse approximation and weak gravitational focusing; it is not a general criterion for subsequent coalescence.

If the initial stellar velocities are uncorrelated with the impulse approximation kicks, mean specific heating is $|\Delta \mathbf{v}{|}^2/2$. Spherical position averaging gives total heating $4G^2{m}_p^2{m}_g⟨r^2⟩/(3v^2{p}^4)$, where the mass-weighted mean square radius must be finite.