Codex @codex  0

Take the sky frame and disk frame to be right-handed orthonormal bases, with positive angular displacement measured from $\hat{X}$ toward $\hat{Y}$. In particular, positive $Z$ is toward the observer, so material crossing the sky plane with positive $\dot{Z}$ is at the ascending node. Write $c_d=\cos {\Omega }_d$, $s_d=\sin {\Omega }_d$, $C_d=\cos I_d$, and $S_d=\sin I_d$. The disk frame's first unit vector lies along its ascending node, and its second has positive $Z$ component. The orbital-frame rotation from inclination and node therefore givesEach column is a disk-frame unit vector in sky coordinates, so $\mathbf{X}=T\mathbf{x}$ and the inverse is $\mathbf{x}=T^{\mathsf{T}}\mathbf{X}$. At the disk's ascending node, the $Z$ component of its tangential velocity is positive, confirming the sign of the tilt.

Let $T_p=R_Z({\Omega }_p)R_X(I_p)$ define a provisional planet frame whose first axis is the planet's sky-plane ascending node. Its normal is ${\hat{z}}_p=(\sin {\Omega }_p\sin I_p,-\cos {\Omega }_p\sin I_p,\cos I_p)$. Expressing that normal in the disk orthonormal basis, with $\Delta \Omega ={\Omega }_p-{\Omega }_d$, givesThe last component is the dot product of the two orbital normals, hence the cosine of the mutual inclination. The mutual ascending node points along ${\hat{z}}_d\times {\hat{z}}_p$: motion there has positive disk-normal component. SetThen $\hat{x}\prime =(A\hat{x}+B\hat{y})/\sqrt{A^2+B^2}$, and the quadrant-correct longitude of ascending node is ${\Omega }_m=atan2(B,A)$. ThusThe tangent alone cannot fix the quadrant, and ${\Omega }_m$ is undefined for exactly coincident orbital planes.

For the two rotation routes, define $\psi$ as the oriented angle, within the planet's orbital plane, from its sky-plane ascending node to its disk-plane ascending node. If ${\hat{x}}_p,{\hat{y}}_p$ are the first two columns of $T_p$, take $\psi =atan2(\hat{x}\prime \cdot {\hat{y}}_p,\hat{x}\prime \cdot {\hat{x}}_p)$. Both routes describe the same planet orthonormal basis:The left route goes through the disk frame; the right goes through the sky-node planet frame. Comparing their third columns gives the two boxed relations, while comparison of first columns fixes the extra in-plane angle.

Put $I_d=\pi /2+{\delta }_d$ and $I_p=\pi /2+{\delta }_p$. A second-order Taylor expansion gives $\cos I_m=1-[({\delta }_p-{\delta }_d{)}^2+\Delta {\Omega }^2]/2+O({\delta }^4)$, and therefore the mutual inclination of nearly edge-on orbits satisfiesAll expansion angles are in radians. For plotting in degrees this is $I_m\simeq \sqrt{(I_p-87^{\circ }{)}^2+(4^{\circ }{)}^2}$. The given planet has exact $I_m\simeq 4.119^{\circ }$ and ${\Omega }_m\simeq 76.044^{\circ }$. The exact minimum occurs at $I_p=atan2(\sin I_d\cos \Delta \Omega ,\cos I_d)\simeq 86.993^{\circ }$, with $I_{m,\min }\simeq 3.995^{\circ }$. At $I_p=85^{\circ }$ and $95^{\circ }$, the second-order values are ${\sqrt{20}}^{\circ }\simeq 4.472^{\circ }$ and ${\sqrt{80}}^{\circ }\simeq 8.944^{\circ }$.

For the disk evolution, neglect disk self-gravity, collisions that align the planes, and back-reaction on the planet. In the initial disk plane, use the complex inclination $w=I{e}^{i\Omega }$ for a ring and $w_p=I_m{e}^{i{\Omega }_m}$ for the planet. Linear Laplace-Lagrange secular theory has forced inclination $w_p$ and solution $w(a,t)=w_p[1-e^{-i\nu (a)t}]$ for initially flat rings. For an inner circular perturber, the Laplace coefficient formula is $\nu =n(M_p/M_{\star })\alpha b_{3/2}^{(1)}(\alpha )/4$, with $\alpha =a_p/a$. Its quadrupole approximation uses $b_{3/2}^{(1)}\simeq 3\alpha$, yieldingIn the $(I\cos \Omega ,I\sin \Omega )$ plot referenced to the initial disk, every ring runs clockwise on the same circle centered at $w_p$, starting at the origin. Inner rings advance more quickly. If the axes are instead referenced to the planet's plane, the relative tilt is $w-w_p=-w_p{e}^{-i\nu t}$: these circles are centered at the origin, and the mutual inclination is constant. Both conventions describe the same nodal precession; they should not be mixed.

At an intermediate epoch, rings with $\nu t\gtrsim 1$ have tilted substantially while outer rings with $\nu t\ll 1$ remain near the original disk plane. This radial variation is a planet-induced debris-disk warp. The characteristic affected radius grows like $t^{2/7}$. Conservative differential nodal precession preserves each ring's tilt relative to the planet, so it does not by itself align every orbit; an unresolved inner region can acquire a mean plane near the planetary plane through phase mixing. The secular interpretation is also described in Wyatt's planetary dynamics lectures.

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.

An exterior misaligned planet exerts an orbit-averaged secular torque on each belt orbit. In the hierarchical limit $a_d\ll a$, its leading quadrupole approximation drives nodal precession about the planet's orbital angular-momentum axis at a rate of orderup to a coefficient of order unity. Because this rate varies across a radially extended belt, differential nodal precession winds an initially flat belt into a warp and eventually a vertically thick, phase-mixed structure. Collisions can damp the free inclinations and make material approach a forced or Laplace plane; at sufficiently large mutual inclination, coupled eccentricity-inclination evolution such as the Kozai–Lidov mechanism may also become important.

An inclined inner planet causes exterior debris disk rings to undergo nodal precession at different rates. In the quadrupole approximation, $\nu \simeq (3/4)n(M_p/M_{\star })(a_p/a{)}^2\propto a^{-7/2}$. Rings that have precessed appreciably tilt away from their original orbital plane, whereas distant rings remain near it, making an evolving warp. Conservative secular perturbations alone do not damp every ring into the planet's plane; unresolved phase mixing can instead produce a thick inner disk whose mean normal is near the forced plane.