Circumbinary Planet Stability With REBOUND

ISEF Category: Physics and Astronomy

Subcategory: Astronomy and Cosmology  ·  Difficulty: Advanced  ·  Setup: Home Setup  ·  Time: Full Year

The Hook

A planet can orbit two stars and still stay stable, but only in a narrow zone. That zone shifts when the stars are eccentric and the planet sits near a resonance. With laptop simulations, you can map where those orbits survive and where they fail. That gives you a real way to study habitability, not just draw pretty orbit diagrams.

What Is It?

Circumbinary planets orbit both stars in a binary system. Think of the two stars like dancers holding hands while a planet circles around them. The planet has to stay far enough away from the stars’ chaotic tug, or its orbit can get kicked out, stretched, or broken apart.

REBOUND is a software package for N-body simulations, which means it calculates how multiple bodies pull on each other over time. You can use it to test many fake planetary systems and see which ones stay stable. The key idea here is long-term stability, or whether an orbit survives for thousands or millions of orbits without becoming too wild.

A Mean-Motion Resonance happens when two objects have orbital periods that line up in a simple ratio, like two-to-one or three-to-one. Those resonances can create gaps or islands of stability. In this project, you can map a resonance-based habitability shore, which is a boundary where stable, planet-friendly orbits begin to appear around an eccentric binary.

Why This Is a Good Topic

This is a strong science fair topic because you can test a clear question with simulation data, not expensive lab gear. The system has real physics, real limits, and real astronomy data to compare against. You can study orbital stability, resonance structure, and binary eccentricity in a way that produces graphs, maps, and statistical comparisons. That gives you a project that is both easy to define and hard to fake.

Research Questions

• How does binary eccentricity change the inner edge of stable circumbinary orbits?
• What is the effect of planet starting distance on long-term orbital survival around an eccentric binary?
• Does adding a small initial orbital tilt change stability near a mean-motion resonance?
• To what extent do different resonance ratios predict stable versus unstable orbit zones?
• Which binary mass ratios create the widest stable circumbinary region?
• How does the simulated stability boundary compare with the orbits of known Kepler circumbinary planets?

Basic Materials

• Laptop or desktop computer with enough memory for repeated simulations.
• REBOUND software installed in Python.
• Python with NumPy, SciPy, and Matplotlib.
• Spreadsheet software or a CSV viewer for organizing outputs.
• Scientific calculator for quick checks of orbital ratios.
• Notebook for logging simulation settings and run IDs.
• Access to published orbital data for known Kepler circumbinary planets.

Advanced Materials

• Laptop or workstation with more cores and memory for batch runs.
• REBOUND and REBOUNDx for extended orbital physics tests.
• Python with pandas, NumPy, SciPy, Matplotlib, and seaborn.
• Jupyter Notebook for documenting simulation sweeps.
• ImageJ or other plotting tool for comparing resonance maps as images.
• Archived orbital elements from NASA Exoplanet Archive or published Kepler papers.
• Version control software such as Git for tracking code changes.

Software & Tools

• REBOUND: Simulates gravitational interactions among multiple bodies and tracks orbital stability over time.
• Python: Runs simulation scripts, processes output files, and automates parameter sweeps.
• Jupyter Notebook: Lets you document each run, inspect results, and keep code and notes together.
• Matplotlib: Makes stability maps, orbit plots, and resonance graphs from your simulation data.
• pandas: Organizes batches of run results into tables for filtering and comparison.

Experiment Steps

1. Define the exact binary and planet parameters you will vary, then keep the rest fixed.
2. Choose a stability metric, such as survival time, orbital crossing, or large eccentricity growth.
3. Build a small pilot set of runs to confirm your code, logging, and output format work.
4. Expand to a parameter sweep that tests distance, binary eccentricity, and resonance location.
5. Compare your simulated stability boundary with published circumbinary planets and their orbital ratios.
6. Analyze the pattern with plots and statistics so you can explain where the habitability shore appears.

Common Pitfalls

• Running too few simulation trials, which makes random orbit outcomes look like real trends.
• Changing several orbital parameters at once, which hides the role of resonance in the stability boundary.
• Using inconsistent units across scripts, which can break orbital periods and distance scales.
• Treating a short simulation as proof of long-term stability, which can miss slow chaotic drift.
• Comparing your results to Kepler planets without matching the same mass ratio, eccentricity, and period ratio context.

What Makes This Competitive

A stronger project will test a clear physical claim, not just make orbit animations. You can raise the level by comparing several binary eccentricities, checking multiple resonance bands, and using a formal stability metric across many trials. Good analysis matters here. If you add uncertainty estimates, compare against real Kepler systems, and explain where your model breaks, your work starts to look much more serious.

Project Variations

• Test how changing the binary mass ratio shifts the stable circumbinary zone for the same resonance pattern.
• Compare stability maps for planets with different initial inclinations to see when tilt destroys resonance protection.
• Focus on one known Kepler circumbinary system and ask how closely your simulation reproduces its observed orbit.

Learn More

• NASA Exoplanet Archive: Search for orbital data on known circumbinary planets and compare your model systems with observed systems.
• NASA ADS: Search published papers on circumbinary planet stability, resonance, and Kepler systems.
• The Astrophysical Journal: Search for review articles and case studies on circumbinary orbital dynamics.
• REBOUND documentation: Read the official simulation guide and example scripts for N-body modeling.
• MIT OpenCourseWare Classical Mechanics: Use the orbital mechanics material to review gravity, energy, and two-body motion before adding the binary system.