1 · The restricted N-body framework
Toomre & Toomre's 1972 demonstration that bridges and tails are tidal relics of close encounters was a restricted calculation: each galaxy was treated as a single point mass, and each disk was a swarm of non-self-gravitating test particles on initially circular orbits. The decoupling is the entire trick — the disk does not pull on itself, the disk does not pull on the companion, and the two galaxy centers fall along an exact two-body trajectory. Test-particle equations of motion run in parallel; cost is linear in particle count rather than the N² of a true N-body code. The morphology that emerges — the long ratty tail, the broad facing-side bridge, the prograde-retrograde asymmetry — is set by the relative orbit and the disk inclinations alone, and it survives essentially unchanged when self-gravity is reintroduced in full simulations.
This entry generalizes the 1972 setup in three ways: (i) the point-mass galaxies are replaced by analytic Hernquist 1990 halos plus thin Miyamoto-Nagai 1975 disks, so the disk rotation curve is realistic at small radii; (ii) integration is done with a kick-drift-kick leapfrog at 0.5 Myr step, symplectic in the canonical variables of each test particle, with the galaxy centers integrated in lockstep so the orbit is consistent at every drift; (iii) the relative orbit is selectable across parabolic, hyperbolic, and eccentric-bound regimes; (iv) a Chandrasekhar 1943 dynamical-friction prescription drags the center orbit through each rigid halo's density field (lnΛ = 2 fixed, Maxwellian F(X), σ ≈ Vc/√2), so bound passages decay and the centers coalesce; below 1 kpc separation the two potentials are locked to their common center of mass and the test particles phase-mix in the combined potential — the closest a restricted calculation can come to a remnant. What that remnant can and cannot be trusted for is in the validity panel below.
Equations of motion
2 · Configurations available
Five presets, four of which reproduce iconic interacting pairs at the qualitative level the 1972 paper claimed. The fifth (retrograde-retrograde) is the control: same orbit, same masses, opposite disk spins, no tails. Click a preset to load it, then change disk inclinations or pericenter and click Recompute.
| Preset | q (kpc) | e | MA (1011) | MB (1011) | iA, ωA | iB, ωB | Comment |
|---|---|---|---|---|---|---|---|
| Antennae (NGC 4038/9) | 19 | 1.00 | 1.0 | 1.0 | 60°, 30° | 60°, 60° | parabolic; Toomre 1972 best-fit geometry |
| Mice (NGC 4676) | 16 | 1.00 | 0.7 | 0.7 | 15°, 0° | −10°, 90° | parabolic; Barnes 2004 disposition |
| M51 (NGC 5194/5) | 22 | 1.00 | 1.2 | 0.4 | 20°, 170° | 30°, 0° | single-passage (Salo & Laurikainen 2000) |
| Bound coalescing (toy) | 8 | 0.50 | 1.0 | 1.0 | 45°, 0° | 45°, 90° | elliptic, T ≈ 420 Myr; friction decays the orbit through repeated passages → coalescence + spheroidal remnant |
| Retrograde control | 19 | 1.00 | 1.0 | 1.0 | 120°, 30° | 120°, 60° | parabolic, both disks retrograde — no tails |
Orbital geometry & disk orientation
3 · Boundary of validity mandatory
Where rigid-potential restricted N-body is trustworthy — and where it isn't
Two things are removed from the equations: (i) self-gravity of the disks and (ii) all gas physics — pressure forces, shocks, cooling, star formation, and feedback. A third, dynamical friction, is no longer absent but is included only as the Chandrasekhar 1943 analytic prescription on rigid Hernquist profiles with a fixed lnΛ = 2 — a calibrated stand-in for the emergent wake, not the wake itself. The failure modes do not appear at the same time in the merger sequence.
What stays right
The morphology and the timing of tidal tail and bridge formation at first pericenter. Tail length grows linearly with time after first pericenter at ≈ vtail ≈ vesc(q), reproducing Hibbard & Mihos 1995 and Barnes & Hibbard 2009 to within ~10%.
Prograde / retrograde asymmetry
Toomre & Toomre's central result that prograde disks form long thin tails and retrograde disks barely respond is purely kinematic. It is reproduced exactly by the restricted calculation — the retrograde control preset is the clean demonstration.
Orbit decay: prescribed, not emergent
The orbital separation now shrinks between passages via the Chandrasekhar 1943 friction term, reproducing the qualitative Karl+ 2010 Antennae sequence (second pericenter a few hundred Myr after the first as friction lowers the apocenter). But lnΛ is a fixed constant and the host density rigid, so the rate of decay carries a factor-of-~2 uncertainty — coalescence times are indicative, not predictive (Boylan-Kolchin+ 2008).
Disk self-gravity → bar/spiral feedback
Real disks respond to a perturbation by raising bars and spiral arms that amplify the tidal response. The restricted calculation gives only the linear tidal kick. Tail mass and width are accurate; inner disk structure (bars, ovals, m=2 spirals) is not reproduced. M51-like grand-design arms in the secondary are partly tidal and partly self-gravitating; see Dobbs+ 2010 for the breakdown.
No gas, no starbursts, no AGN fueling
Mergers funnel gas inward and trigger 10–100 Myr starbursts (Mihos & Hernquist 1996) and obscured AGN (Hopkins+ 2008). None of that is in this calculation. The tool produces a stellar-tracer prediction, not a luminosity or color prediction.
Final coalescence and remnant
Coalescence now occurs: below 1 kpc the centers are locked to their COM and the disks phase-mix into a spheroidal swarm in the combined potential. But the real remnant is shaped by violent relaxation in a time-varying self-consistent potential (Barnes 1992, Hopkins+ 2009) — here the potentials stay frozen, so the de Vaucouleurs profile, kinematic misalignments, and boxy/disky isophotes are not reproduced. The remnant shown is a morphological cartoon: right shape class, wrong profile.
Bottom line. The restricted N-body model is a kinematic predictor of the first 300–500 Myr after pericenter: tail length, bridge geometry, sense of asymmetry, dependence on disk inclination. With the friction prescription it now also gives a plausible decay-and-coalescence sequence and a spheroidal end state — but merger timescales are indicative only, and remnant profile, kinematics, and every baryonic observable remain out of scope. The crossover to full SPH/GADGET/FIRE is well-mapped in the literature cited below.
References
- Toomre & Toomre 1972. Galactic Bridges and Tails. ApJ 178, 623. ADS — the founding paper; restricted three-body code, NGC 4038/9 + 4676 + M51 + Arp 295.
- Hernquist 1990. An analytical model for spherical galaxies and bulges. ApJ 356, 359. ADS — Φ(r) = −GM/(r+a); halo potential used here.
- Miyamoto & Nagai 1975. Three-dimensional models for the distribution of mass in galaxies. PASJ 27, 533. ADS — thin-disk potential added on top of the Hernquist halo.
- Barnes & Hernquist 1992. Dynamics of interacting galaxies. ARA&A 30, 705. ADS — canonical review of the field; benchmark against full N-body.
- Barnes & Hernquist 1992. Formation of dwarf galaxies in tidal tails. Nature 360, 715. Nature — tidal-dwarf clumps that need self-gravity.
- Hibbard & Mihos 1995. Dynamical modeling of NGC 7252. AJ 110, 140. ADS — tail-length test for restricted models.
- Mihos & Hernquist 1996. Gasdynamics and starbursts in major mergers. ApJ 464, 641. ADS — the gas piece this model omits.
- Salo & Laurikainen 2000. N-body model for M51. MNRAS 319, 377. MNRAS — single-passage vs. bound encounter for the M51/5195 system.
- Barnes 2004. Shock-induced star formation in a model of the Mice. MNRAS 350, 798. ADS — parabolic encounter NGC 4676; first passage 170 Myr ago.
- Karl, Naab, Johansson+ 2010. One moment in time — modeling star formation in the Antennae. ApJ 715, L88. IOP — orbit fit places observation 40 Myr after 2nd pericenter.
- Dobbs, Theis, Pringle, Bate 2010. Simulations of the grand design galaxy M51. MNRAS 403, 625. MNRAS — tidal vs. self-gravitating components of the arms.
- Renaud, Bournaud, Emsellem+ 2015. A parsec-resolution simulation of the Antennae. MNRAS 446, 2038. MNRAS — high-resolution AMR comparison.
- Barnes & Hibbard 2009. Identikit 1: a modeling tool for interacting disk galaxies. AJ 137, 3071. ADS — restricted-N-body fitting framework descendant of TT72.
- Chandrasekhar 1943. Dynamical Friction. I. ApJ 97, 255. ADS — the decay term, applied here as an analytic prescription on the center orbit.
- Boylan-Kolchin, Ma & Quataert 2008. Dynamical friction and galaxy merging time-scales. MNRAS 383, 93. ADS — calibration of merger times against full N-body; the factor-of-~2 caveat on Chandrasekhar estimates.
- Hopkins, Hernquist, Cox, Kereš 2008. A cosmological framework for the co-evolution of QSOs, galaxies, and their halos. ApJS 175, 356. ADS — merger-driven AGN fueling.