Recognition: unknown
Life is But a Stream: The Distribution of Planetary Systems Along Stellar Streams and their Properties
Pith reviewed 2026-05-08 15:41 UTC · model grok-4.3
The pith
Stars near the edges of stellar streams retain unperturbed planetary systems while those near the centers experience more disruptions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Direct N-body simulations demonstrate that stars with early cluster escape times retain all their planets as they spend most of their time in the cluster's low-density outskirts. Stars with later escape times experience a range of survival fractions due to different local densities and encounter types. In the stellar stream formed by the cluster's dissolution, stars near the edge are more likely to have unperturbed planetary systems, while stars near the centre have a higher chance of having planets pushed to eccentric orbits, inclined orbits, or stripped from the system entirely. The suite of simulations provides an estimate of the probability that a star will host a planet with a given the
What carries the argument
Stellar stream position Delta phi as a proxy for a star's escape time and cumulative exposure to gravitational perturbations during cluster dissolution.
If this is right
- Early-escaping stars retain all their planets.
- Late-escaping stars show a wide range of planet survival fractions.
- Edge stars in the stream are likely to have unperturbed planetary systems.
- Central stars in the stream are more likely to have eccentric, inclined, or stripped planets.
- Probability estimates link initial planet semi-major axis to stream location for hosting likelihood.
Where Pith is reading between the lines
- Observations of exoplanets in known stellar streams could be prioritized by position to increase detection rates of intact systems.
- This approach offers a new way to probe the effects of cluster environments on planet formation and evolution without being limited to bound clusters.
- Accounting for stream position may help reconcile the difference in exoplanet occurrence rates between cluster and field stars.
- Extending the simulations to include additional physics like galactic tides could refine the position-dependent probabilities.
Load-bearing premise
The N-body simulations fully capture the relevant dynamics without needing to include effects like gas drag, stellar evolution, or external galactic tides that might change escape times or interaction outcomes.
What would settle it
Surveying planetary systems in a specific stellar stream and finding that the fraction of planets with certain semi-major axes does not vary with position along the stream as predicted by the simulations.
Figures
read the original abstract
The majority of discovered exoplanets have been observed orbiting field stars as opposed to within a star cluster. To determine whether the lack of observed exoplanets in star clusters is due to gravitational perturbations or observational limitations, we consider the possibility of studying exoplanets in stellar streams. We present the results of direct $N$-body simulations of planetary systems around stars that orbit within a star cluster. Our simulations demonstrate that stars with early cluster escape times tend to retain all their planets as they spend most of their time orbiting in the cluster's low-density outskirts. Alternatively, stars with later escape times can have a wide range of survival fractions as they are subjected to a range of local densities and encounter types. With respect to the stellar stream that forms as the result of the cluster's dissolution, stars near the edge of the stream are therefore more likely to have unperturbed planetary systems. Conversely, stars near the centre of the stream have a higher chance of having planets pushed to eccentric orbits, inclined orbits, or stripped from the system entirely. From our suite of simulations, we provide an estimate of the probability that a star will host a planet with a given initial semi-major axis $a_0$ based on the star's location along a stellar stream $\Delta \phi$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses direct N-body simulations of planetary systems embedded in a dissolving star cluster to argue that a star's escape time from the cluster determines the degree of gravitational perturbation its planets experience. Early escapers, which spend more time in low-density regions, retain unperturbed systems, while later escapers encounter higher densities and more disruptive interactions. This leads to the central claim that, once the cluster dissolves into a stellar stream, stars near the stream edges (small |Δϕ|) are more likely to host unperturbed planets, whereas stars near the stream center (larger |Δϕ|) have elevated probabilities of eccentric or inclined orbits or complete planet stripping. The work supplies numerical estimates of the survival probability for a planet with given initial semi-major axis a0 as a function of the star's stream position Δϕ.
Significance. If the mapping holds, the result supplies a direct, observationally testable link between stellar-stream morphology and exoplanet demographics, potentially explaining the paucity of detected planets in open clusters and offering a new diagnostic for the dynamical history of field stars. The forward-simulation approach (no fitted parameters) is a methodological strength that avoids circularity and allows falsifiable predictions for specific a0 values.
major comments (3)
- [§2] §2 (N-body simulation setup): The integrations omit the external galactic tidal field. This directly affects escape times, the rate of cluster dissolution, and the differential orbital evolution that sets the observed Δϕ distribution along the stream. Because the central claim equates Δϕ with perturbation level via escape time, the probability tables P(survive | Δϕ, a0) are uncontrolled and could shift once tides are included.
- [Abstract and §3] Abstract and §3 (results): No details are given on initial conditions (cluster density profile, number of stars, planetary architectures), numerical resolution, integrator, or validation against standard cluster-dissolution benchmarks. Without these, the reported survival fractions and probability estimates lack the supporting evidence needed to assess numerical convergence or robustness.
- [§4] §4 (discussion): No uncertainty quantification, sensitivity tests, or error bars accompany the probability estimates. This is load-bearing for the claim that stream position can be used to predict planet survival, as small changes in encounter statistics could alter the quoted probabilities substantially.
minor comments (2)
- The symbol Δϕ is introduced in the abstract without an explicit definition or reference to a figure showing its measurement along the stream; a short equation or diagram would improve clarity.
- The manuscript would benefit from a brief comparison table of survival fractions versus escape time to make the mapping from early/late escapers to stream edges/center more quantitative.
Simulated Author's Rebuttal
We are grateful to the referee for their thorough review and positive evaluation of the paper's significance. We have carefully considered each major comment and will make the necessary revisions to strengthen the manuscript.
read point-by-point responses
-
Referee: [§2] §2 (N-body simulation setup): The integrations omit the external galactic tidal field. This directly affects escape times, the rate of cluster dissolution, and the differential orbital evolution that sets the observed Δϕ distribution along the stream. Because the central claim equates Δϕ with perturbation level via escape time, the probability tables P(survive | Δϕ, a0) are uncontrolled and could shift once tides are included.
Authors: We agree that the external galactic tidal field plays an important role in determining escape times and the structure of the resulting stellar stream. Our simulations were designed to isolate the impact of close stellar encounters within the cluster on planetary systems. To address this concern, we will incorporate a galactic tidal field in follow-up simulations and update the reported probabilities P(survive | Δϕ, a0) accordingly. This revision will confirm whether the qualitative trend between stream position and planet survival persists under more complete physical conditions. revision: yes
-
Referee: [Abstract and §3] Abstract and §3 (results): No details are given on initial conditions (cluster density profile, number of stars, planetary architectures), numerical resolution, integrator, or validation against standard cluster-dissolution benchmarks. Without these, the reported survival fractions and probability estimates lack the supporting evidence needed to assess numerical convergence or robustness.
Authors: We will revise the manuscript to provide a comprehensive description of the simulation setup in §3, including the cluster initial conditions, the number of stars, the planetary system architectures considered, the numerical integrator and resolution parameters, as well as comparisons to established benchmarks for cluster dissolution. These additions will enable a proper assessment of the numerical robustness of our results. revision: yes
-
Referee: [§4] §4 (discussion): No uncertainty quantification, sensitivity tests, or error bars accompany the probability estimates. This is load-bearing for the claim that stream position can be used to predict planet survival, as small changes in encounter statistics could alter the quoted probabilities substantially.
Authors: We acknowledge the importance of uncertainty quantification. In the revised §4, we will add error bars to the probability estimates, derived from the ensemble of simulations, and include sensitivity tests to variations in initial conditions. This will provide a more rigorous basis for the predictive use of stream position. revision: yes
Circularity Check
No circularity; probabilities are direct outputs of forward N-body integrations
full rationale
The paper's central results consist of probability estimates for planetary survival as a function of stream position Δϕ and initial semi-major axis a0. These estimates are generated by running direct N-body simulations of cluster dissolution plus planetary systems and then tabulating outcomes; the mapping from escape time to Δϕ to perturbation level is therefore an emergent computational result rather than an input definition, fitted parameter, or self-referential relation. No load-bearing self-citations, uniqueness theorems, or ansatzes are invoked to close the argument, and the derivation does not rename or re-express prior empirical patterns as new theoretical content. The chain is self-contained as a forward-modeling experiment.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Newtonian gravity and point-mass interactions govern stellar and planetary motions
- domain assumption Planetary systems begin stable and circular before cluster encounters
Reference graph
Works this paper leans on
-
[1]
Aarseth, S. J. 2003, Gravitational N-Body Simulations (Cambridge University Press)
2003
-
[2]
Adams, F. C. 2010, ARA&A, 48, 47, doi: 10.1146/annurev-astro-081309-130830
-
[3]
2018, A&A, 619, A130, doi: 10.1051/0004-6361/201834285
Adibekyan, V., de Laverny, P., Recio-Blanco, A., et al. 2018, A&A, 619, A130, doi: 10.1051/0004-6361/201834285
-
[4]
Arzoumanian, D., Arakawa, S., Kobayashi, M. I. N., et al. 2023, ApJL, 947, L29, doi: 10.3847/2041-8213/acc849
-
[5]
, year = 1993, month = oct, volume =
Backer, D. C., Foster, R. S., & Sallmen, S. 1993, Nature, 365, 817, doi: 10.1038/365817a0
-
[6]
Batista, S. F. A., Adibekyan, V. Z., Sousa, S. G., et al. 2014, A&A, 564, A43, doi: 10.1051/0004-6361/201423645
-
[7]
Batista, S. F. A., & Fernandes, J. 2012, NewA, 17, 514, doi: 10.1016/j.newast.2011.12.001
-
[8]
Bonaca, A., & Price-Whelan, A. M. 2025, NewAR, 100, 101713, doi: 10.1016/j.newar.2024.101713
-
[9]
The Astrophysical Journal Supplement Series , author =
Bovy, J. 2015, ApJS, 216, 29, doi: 10.1088/0067-0049/216/2/29
-
[10]
Brown, A. G. A., Portegies Zwart, S. F., & Bean, J. 2010, MNRAS, 407, 458, doi: 10.1111/j.1365-2966.2010.16921.x
-
[11]
Cai, M. X., Kouwenhoven, M. B. N., Portegies Zwart, S. F., & Spurzem, R. 2017, MNRAS, 470, 4337, doi: 10.1093/mnras/stx1464
-
[12]
X., Meiron, Y., Kouwenhoven, M
Cai, M. X., Meiron, Y., Kouwenhoven, M. B. N., Assmann, P., & Spurzem, R. 2015, ApJS, 219, 31, doi: 10.1088/0067-0049/219/2/31
-
[13]
2019, MNRAS, 489, 4311, doi: 10.1093/mnras/stz2467
Spurzem, R. 2019, MNRAS, 489, 4311, doi: 10.1093/mnras/stz2467
-
[14]
Cai, M. X., Portegies Zwart, S., & van Elteren, A. 2018, MNRAS, 474, 5114, doi: 10.1093/mnras/stx3064
-
[15]
Cassan, A., Kubas, D., Beaulieu, J. P., et al. 2012, Nature, 481, 167, doi: 10.1038/nature10684
-
[16]
Curtis, J. L., Ag¨ ueros, M. A., Mamajek, E. E., Wright, J. T., & Cummings, J. D. 2019, AJ, 158, 77, doi: 10.3847/1538-3881/ab2899
-
[17]
Daffern-Powell, E. C., Parker, R. J., & Quanz, S. P. 2022, MNRAS, 514, 920, doi: 10.1093/mnras/stac1392
-
[18]
2023, AJ, 166, 219, doi: 10.3847/1538-3881/acff67
Dai, Y.-Z., Liu, H.-G., Yang, J.-Y., & Zhou, J.-L. 2023, AJ, 166, 219, doi: 10.3847/1538-3881/acff67
-
[19]
2016, MNRAS, 461, 1734, doi: 10.1093/mnras/stw1094
Drass, H., Haas, M., Chini, R., et al. 2016, MNRAS, 461, 1734, doi: 10.1093/mnras/stw1094
-
[20]
Grondin, S. M., Webb, J. J., Lane, J. M. M., Speagle, J. S., & Leigh, N. W. C. 2024, MNRAS, 528, 5189, doi: 10.1093/mnras/stae203 6 https://tess.mit.edu/science/
-
[21]
Monthly Notices of the Royal Astronomical Society , keywords =
Hanse, J., J´ ılkov´ a, L., Portegies Zwart, S. F., & Pelupessy, F. I. 2018, MNRAS, 473, 5432, doi: 10.1093/mnras/stx2721
-
[22]
2003, The Gravitational Million-Body Problem: A Multidisciplinary Approach to Star Cluster Dynamics
Heggie, D., & Hut, P. 2003, The Gravitational Million-Body Problem: A Multidisciplinary Approach to Star Cluster Dynamics
2003
-
[23]
Kroupa, P. 2001, MNRAS, 322, 231, doi: 10.1046/j.1365-8711.2001.04022.x
-
[24]
Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57, doi: 10.1146/annurev.astro.41.011802.094844
-
[25]
Luhman, K. L. 2007, ApJS, 173, 104, doi: 10.1086/520114 Mart´ ın, E. L.,ˇZerjal, M., Bouy, H., et al. 2025, A&A, 697, A7, doi: 10.1051/0004-6361/202450793
-
[26]
Newton, E. R., Mann, A. W., Kraus, A. L., et al. 2021, AJ, 161, 65, doi: 10.3847/1538-3881/abccc6
-
[27]
2021, ApJ, 912, 162, doi: 10.3847/1538-4357/abeaac
Pang, X., Li, Y., Yu, Z., et al. 2021, ApJ, 912, 162, doi: 10.3847/1538-4357/abeaac
-
[28]
Pelupessy, F. I., van Elteren, A., de Vries, N., et al. 2013, A&A, 557, A84, doi: 10.1051/0004-6361/201321252
-
[29]
2020, ApJ, 897, 60, doi: 10.3847/1538-4357/ab9533
Pfalzner, S., & Vincke, K. 2020, ApJ, 897, 60, doi: 10.3847/1538-4357/ab9533
-
[30]
Piatti, A. E., & Carballo-Bello, J. A. 2020, A&A, 637, L2, doi: 10.1051/0004-6361/202037994
-
[31]
Plummer, H. C. 1911, MNRAS, 71, 460, doi: 10.1093/mnras/71.5.460 Portegies Zwart, S. 2021, A&A, 647, A136, doi: 10.1051/0004-6361/202038888 Portegies Zwart, S., McMillan, S. L. W., van Elteren, E.,
-
[32]
Computer Physics Communications , keywords =
Pelupessy, I., & de Vries, N. 2013, Computer Physics Communications, 184, 456, doi: 10.1016/j.cpc.2012.09.024 Portegies Zwart, S., McMillan, S., Harfst, S., et al. 2009, NewA, 14, 369, doi: 10.1016/j.newast.2008.10.006 Portegies Zwart, S. F. 2009, ApJL, 696, L13, doi: 10.1088/0004-637X/696/1/L13 Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 201...
-
[33]
Rein, H., & Liu, S. F. 2012, A&A, 537, A128, doi: 10.1051/0004-6361/201118085
-
[34]
Rein, H., & Spiegel, D. S. 2015, MNRAS, 446, 1424, doi: 10.1093/mnras/stu2164
-
[35]
Scholz, A., Jayawardhana, R., Muzic, K., et al. 2012, ApJ, 756, 24, doi: 10.1088/0004-637X/756/1/24
-
[36]
1993, in Astronomical Society of the Pacific Conference Series, Vol
Sigurdsson, S. 1993, in Astronomical Society of the Pacific Conference Series, Vol. 36, Planets Around Pulsars, ed. J. A. Phillips, S. E. Thorsett, & S. R. Kulkarni, 173–179
1993
-
[37]
1958, ApJ, 127, 17, doi: 10.1086/146435
Spitzer, Lyman, J. 1958, ApJ, 127, 17, doi: 10.1086/146435
-
[38]
Spitzer, L. 1987, Dynamical evolution of globular clusters 12 van Elteren, A., Portegies Zwart, S., Pelupessy, I., Cai, M. X., & McMillan, S. L. W. 2019, A&A, 624, A120, doi: 10.1051/0004-6361/201834641
-
[39]
2015, MNRAS, 450, 4070, doi: 10.1093/mnras/stv817
Wang, L., Spurzem, R., Aarseth, S., et al. 2015, MNRAS, 450, 4070, doi: 10.1093/mnras/stv817
-
[40]
Weatherford, N. C., & Bonaca, A. 2026, ApJ, 997, 90, doi: 10.3847/1538-4357/ae21e0
-
[41]
Webb, J. J., & Bovy, J. 2022, MNRAS, 510, 774, doi: 10.1093/mnras/stab3451
-
[42]
J., Price-Jones, N., Bovy, J., et al
Webb, J. J., Price-Jones, N., Bovy, J., et al. 2020, MNRAS, 494, 2268, doi: 10.1093/mnras/staa788
-
[43]
Winn, J. N., & Fabrycky, D. C. 2015, ARA&A, 53, 409, doi: 10.1146/annurev-astro-082214-122246
-
[44]
2021, AJ, 162, 171, doi: 10.3847/1538-3881/ac1f1f
Ye, X., Zhao, J., Zhang, J., Yang, Y., & Zhao, G. 2021, AJ, 162, 171, doi: 10.3847/1538-3881/ac1f1f
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.