Čtyři hvězdy v krabičce od Merkuru

Why the Four-Body Problem Has No Solution

Imagine one planet orbiting one star. We can solve that exactly — Kepler did it in the 17th century. An ellipse. Done.

Add a second star. Now you have three bodies, and a problem that troubled Newton, Lagrange, Poincare, and hundreds of others. Three bodies interacting through gravity have no general analytical solution. Orbits are often chaotic — a small change in initial conditions leads to a completely different outcome.

Now add a fourth body.

The four-body problem is worse than three-body. More interactions, more resonances, more ways the system can fall apart. Most quadruple star systems that nature creates are unstable — they last a geologically short time before one star escapes or two collide.

And yet: TIC 120362137 exists. Four stars, compact as Mercury's box, pulsating like a metronome.

Architecture: 3+1

In 2026, a team of astrophysicists published the discovery of the most compact known 3+1 quadruple system in Nature Communications.

Architecture:

Inner binary (Aa+Ab): two stars orbiting each other with a period of 3.28 days. Both more massive and hotter than the Sun. They form an eclipsing binary — from our perspective, they regularly overlap.

Third star (B): orbits the inner pair with period P_mid = 51.3 days. Also more massive than the Sun. This entire trio — three stars, three orbits — fits into a space smaller than Mercury's orbit around our Sun. Mercury orbits at 0.39 AU. The three stars are closer together.

Fourth star (C): orbits the entire inner system with period P_out = 1,046 days (2.86 years). A Sun-like star. More distant, but still gravitationally bound.

Spectral lines of all four stars were detected — this is rare. Usually in multiple systems, we see only one or two components. Here we see all four.

Future of the system: modeling shows the entire system will likely end up as a pair of white dwarfs. Stars will evolve, inflate, lose mass — but the gravitational bond will survive.

What We Hear in the Data

We analyzed 46,953 brightness measurements from the TESS satellite, spread across 9 sectors and 5 years (2019-2024). We performed Lomb-Scargle analysis per-sector and combined.

The spectrum shows three distinct groups:

1. Harmonic series from the modulation frequency:
f0 = 0.098 c/d (period 10.2 days) with clean harmonics 2f0, 3f0, 4f0, 5f0. Amplitudes decrease monotonically (0.222, 0.159, 0.106, 0.047, 0.016). Key relation: f0 = 5/P_mid = 5/51.3 = 0.0975 c/d. The entire harmonic series consists of multiples of the orbital frequency of third star B (10th, 15th, 20th, 25th harmonics of P_mid). The 3rd harmonic 3f0 = 0.290 c/d (P = 3.45 d) is close to, but not identical with, the inner binary orbital period P_in = 3.28 d (Borkovits et al. 2026). The dominant signal in the light curve is the modulation envelope from the EB-third star interaction, not the eclipsing period itself.

2. p-mode cluster around 23 c/d:
Frequencies 22.9-23.3 c/d (period ~1 hour). Short-period pulsations of one component — likely delta Scuti oscillations. A pulsating star in an eclipsing system: rare and valuable, because mass and radius are independently known from eclipses.

3. Secondary peaks around 7.5 and 30.5 c/d:
Likely combination frequencies, aliases, or additional pulsation modes. Present mainly in sectors with longer cadence (1800s).

Key finding: the fundamental frequency is present in all 9 sectors with zero deviation. A four-body system behaving like clockwork.

Five harmonics of P_mid (5/P_mid through 25/P_mid), stable for 5 years. P_in = 3.28 d near 3f0 = 3.45 d. Modulation from third star dominates the spectrum.

Entropy 0.49 — Order in Chaos

The four-body gravitational problem has no analytical solution and is generally chaotic. Yet TIC 120362137 exhibits spectral entropy of 0.486 — deep in the ordered region.

For comparison on our tuning fork-chaos scale:

— Przybylski's Star (scroll-L1-027): low entropy, one dominant signal — tuning fork
— TIC 120362137: entropy 0.49, harmonic series — guitar string
— VY Canis Majoris (scroll-L1-028): entropy 0.87, no stable signal — falling drum

Energy concentration: 60.5% of all spectral energy sits in just 10 frequencies. Four gravitationally interacting stars manage to distribute energy into predictable modes instead of scattering it into noise.

Why? Hierarchy. A 3+1 system is stable because it is hierarchical: the inner pair is tight, the third star orbits the pair on a wider orbit, the fourth orbits the entire trio even further out. Each hierarchy level is separated in time — periods differ by orders of magnitude (days vs 51 days vs 1,046 days). Systems that violate this (periods close together) fall apart. Those that respect it survive billions of years.

TIC 120362137 survived. And it sings.

Why Feynman Explained to Children

Richard Feynman used to say that if you can't explain something to a six-year-old, you don't understand it yourself.

So:

Imagine two dancers holding hands and spinning around. That's the inner binary. Now add a third who circles around them in a bigger loop. And a fourth who circles around all three even further out.

Four dancers, each on a different path, each at a different speed. And yet they never collide. Ever. We watched them for five years and they never faltered.

How do they do it? They have a rule: never run at the same speed as anyone else. The inner pair spins every 3.28 days. The third dancer orbits in 51 days. The fourth in three years. They are so far apart in speed that they don't interfere.

And all of them fit inside a box that could hold Mercury's orbit around the Sun. Four stars in Mercury's box.

And the most beautiful thing? One of those stars also pulsates. Once an hour, it inflates a little and shrinks again. Like a heart beating inside a dancing machine. Four bodies, one heart, zero chaos.

Fundamental Frequency Stability

TESS observed TIC 120362137 in nine independent sectors from 2019 to 2024. Cadence improved progressively: 1800s (2019-2020), 600s (2021-2022), 200s (2024).

Fundamental frequency f0 = 0.098 c/d appeared in all nine sectors with identical values. RMS variability of 49,000-68,000 ppm (5-7%) confirms the signal is real and strong.

The p-mode cluster (22.9-23.3 c/d) is stable in sectors with sufficient sampling (cadence 600s and shorter). In sectors 14 and 26 (1800s cadence), these frequencies are near the Nyquist limit.

Total data span: 1,850 days (5.06 years). Frequency resolution of the combined periodogram: 0.00054 c/d. This allows resolving frequencies with precision ~0.6 uHz.

Technical Parameters

System TIC 120362137 (= TYC 2651-975-1). Coordinates: RA 18h59m41s, Dec +36d33m27s (Lyra constellation).

Data: TESS sectors 14, 26, 40, 41, 53, 54, 74, 80, 81. Pipeline TESS-SPOC (sectors 14-74) and QLP (sectors 80-81). Total 46,953 quality data points. Total span 1,850 days.

Analysis: Lomb-Scargle periodogram, streaming per-sector processing (3,000 frequency points, range 0.05-36 c/d), combined periodogram as weighted average of sectors. Peak detection with scipy.signal.find_peaks.

Where We Meet Standard Physics

Hierarchical stability of 3+1 systems requires sufficiently separated periods (Mardling & Aarseth 2001). TIC 120362137 satisfies this condition. But why does a four-body system that should be chaotic exhibit a perfect harmonic series over 5 years of observation? Celestial mechanics describes stability. It does not describe why it is so perfect.