{"id":"77d71d63-e373-466a-a2f8-d875595c3771","arxiv_id":"2411.12895","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Hill-sphere Rossby number set by flow speed versus orbital shearing predicts whether an exoplanet outflow forms a spherical bubble or a tidal stream and the sign of its transit velocity gradient.","lead":"Close-in planets lose their atmospheres in shapes that range from round expanding bubbles to long thin streams, and this paper finds that a single dimensionless number, the Rossby number at the planet's Hill sphere, controls which shape appears. The result gives transit observers a direct way to read outflow temperature and mass-loss geometry from the Doppler shift of escaping gas.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central 'Rossby number alone' collapse in Figure 5 is not tested against variations in mass-loss rate or stellar wind strength, so the claimed universality is broader than the simulation suite supports.","rationale":"The reader identified the population-level temperature modeling as the weakest assumption, which is a real downstream sensitivity. My concern is more upstream: the central relation between Ro,H and ξ_los is calibrated on a simulation suite that does not vary mass-loss rates or stellar wind strength, even though Section 5.4 says those quantities shape outflows. This makes the abstract claim 'Rossby number alone is sufficient' stronger than the evidence. The concern does not invalidate the dimensionless organization of the 27 reported runs, nor does it undermine the useful qualitative distinction between bubbles and streams. It does, however, justify a CONDITIONAL verdict rather than full acceptance, because the key predictive relation needs a broader parameter scan before it can be applied to the heterogeneous exoplanet population. The reader's temperature sensitivity is related but distinct: both are examples of model dependence, but the mass-loss and stellar-wind coverage is internal to the paper's own simulation design and therefore more directly load-bearing for the central claim. My concrete test would settle whether the Figure 5 collapse is universal or an artifact of fixed auxiliary parameters, and it requires only modest additional computational effort beyond the existing suite.","tokens_in":1000,"tokens_out":848,"duration_ms":42581,"concrete_test":"Run a set of Athena++ models at fixed Ro,H ≈ 1 (for example, a = 0.025 au, Rp = 1 RJ, λp = 4) while varying the planetary mass-loss rate over 10^10, 10^11, and 10^12 g/s and the stellar wind escape parameter over λ* = 5, 15, 30, then recompute ξ_los from the same line-of-sight integrals as Equations (7)–(11). If ξ_los changes by more than about 0.1 at fixed Ro,H, the Figure 5 collapse is not unique and the 'Rossby number alone' claim must be weakened to 'Rossby number given fixed mass-loss and stellar-wind parameters.'","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's strongest claim, stated in the abstract and conclusions, is that the Hill-sphere Rossby number (Equation 6) is sufficient to predict outflow morphology and the transit kinematic gradient ξ_los. The evidence is Figure 5, where 27 simulations appear to collapse onto a single relation ξ_los ≈ ln Ro,H. However, the simulation grid varies only semi-major axis a, planetary radius Rp, and planetary escape parameter λp (Section 3.2). The planetary mass-loss rate and stellar wind escape parameter are held fixed at 10^11 g/s and λ* = 15, respectively. Section 5.4 explicitly states that stellar winds redirect initially bubble-like or stream-like flows into cometary tails, and that the outcome depends on the ratio of stellar to planetary mass-loss rate. Thus a parameter explicitly excluded from Ro,H is acknowledged to alter outflow morphology and, presumably, ξ_los. Without runs varying Mdot_p, Mdot_*, or λ*, the collapse in Figure 5 could be a property of a one-dimensional slice of parameter space rather than a universal organizing principle. The population-level prediction in Section 5.3 inherits this issue: even if outflow temperatures are correctly estimated by the sunset/sunbather models, systems at the same Ro,H could still have different observable morphologies if their mass-loss rates or stellar wind environments differ. This is a correctness risk for the headline claim, not merely for the temperature calibration, and it is addressable by targeted simulations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates the morphology of hydrodynamic outflows from close-in exoplanets. Using a suite of 27 Athena++ simulations in which semi-major axis, planetary radius, and planetary escape parameter are varied, the authors identify two morphologies: nearly isotropic 'bubbles' and thin tidal 'streams.' They define a Hill-sphere Rossby number, Ro_H = c_s/(Omega r_H) (Eq. 6), and show that simulation diagnostics collapse approximately onto a relation between Ro_H and the dimensionless transit line-of-sight velocity gradient xi_los (Fig. 5). They then use photoevaporation-model sound speeds to predict Ro_H for the known exoplanet population (Fig. 6) and discuss observing strategies.","tokens_in":16709,"tokens_out":6166,"duration_ms":58239,"significance":"The clean dimensional analysis in Section 2 and the release of the simulation source code are strengths. The proposed Rossby-number framework is simple and testable: if it holds, measuring transit velocity gradients would constrain outflow sound speeds and hence outflow temperatures, independent of spectral-line thermal broadening. The paper also makes a falsifiable population-level prediction. The main caveats are that the simulation grid fixes stellar wind and mass-loss parameters and that the population prediction inherits model-dependent temperatures; both need to be addressed before the headline claim is as broad as stated.","major_comments":[{"comment":"Section 3.2 and Table 1 restrict the simulation grid to variations of a, R_p, and lambda_p, fixing Mdot_p = Mdot_* = 10^11 g/s and lambda_* = 15. Section 5.4, however, states that stellar winds redirect initially stream- or bubble-like flows into cometary tails whose dynamics depend on the ratio of stellar to planetary mass-loss rates. Since Ro_H (Eq. 6) does not contain Mdot_p, Mdot_*, or lambda_*, the Figure 5 collapse and the abstract's claim that 'Rossby number alone is sufficient' are not supported outside this one-dimensional slice of parameter space. I recommend either adding simulations with varied Mdot_p, Mdot_*, and lambda_* to test whether the collapse persists, or explicitly restricting the claim to the weak-stellar-wind regime and revising the abstract and conclusions accordingly.","section":"§3.2, Table 1, §5.4, Fig. 5"},{"comment":"The predicted population distribution of outflow morphologies is based on sound speeds from the sunset and sunbather photoevaporation models (Linssen et al. 2024a,b), which assume a stellar spectral energy distribution, solar metallicity, and a one-dimensional radiative-transfer treatment. The authors acknowledge in Section 5.3 that HAT-P-32b and HAT-P-67b require outflows cooler than these predictions by a factor of roughly two (Nail et al. 2024a). Because Ro_H is proportional to c_s, such a temperature offset propagates directly into the predicted Ro_H and hence the assigned morphology. Please add a sensitivity test (e.g., halving and doubling the adopted sound speeds) and state in the text that Figure 6 is a model-dependent prediction rather than a direct inference from observations.","section":"§5.3, Eq. (6), Fig. 6"},{"comment":"The claim that Ro_H predicts 'kinematic gradients across transit' is based on the mass-weighted, optically-thin line-of-sight velocity xi_los. Section 4.4 correctly notes that this quantity is not directly observable because spectral line formation depends on opacity and the line-forming region; the abstract and conclusions should carry this caveat or soften 'kinematic gradients' to 'kinematic gradients in the optically thin limit.' As written, the headline overstates the direct observability of the modeled diagnostic.","section":"Abstract, §4.4, Eqs. (8)-(11)"}],"minor_comments":[{"comment":"Introduction: 'eg.' should be 'e.g.' (e.g., 'metal, hydrogen (eg. Lyman α) ...').","section":"§1"},{"comment":"Figure 2 caption: 'the the only factor' contains a duplicated article; remove the second 'the.'","section":"Fig. 2 caption"},{"comment":"Section 5.3: 'absorbtion' should be 'absorption' in 'He 1083 nm absorbtion' (the same typo appears elsewhere in the text).","section":"§5.3"},{"comment":"Section 5.1: 'it may be able to use measured constraints' is ungrammatical; 'it may be possible to use' is intended.","section":"§5.1"},{"comment":"Figure 4 caption: 'displaces' should be 'displays' in 'Figure 4 displaces density slices.'","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal well and the core simulations are reproducible. My main concern is the breadth of the headline claim relative to the fixed mass-loss grid; a major revision with either additional runs or careful re-scoping would make the paper publishable. I do not see a need to question the authors' integrity or novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper deserves a serious referee. The new thing is the Hill-sphere Rossby number as an organizing parameter for tidal shaping of exoplanet outflows, backed by a 27-model 3D hydrodynamic suite that shows a clean collapse of the transit velocity gradient xi_los onto ln(Ro_H). That is a real addition; McCann et al. and others had the pieces, but nobody had made this the organizing variable or demonstrated the collapse over a systematic grid. The dimensional analysis in Section 2 is clean, the simulations are described in enough detail to reproduce, and the energy-angular momentum phase-space treatment in Section 4.3 is a nice extra. Credit also where due: Section 5.3 openly flags that the sunset/sunbather outflow temperatures are uncertain, and Section 5.4 admits that stellar winds redirect flows into cometary tails.\n\nThe soft spot is exactly what the stress-test note says. The abstract's 'Rossby number alone is sufficient' is stronger than the evidence. The grid varies a, Rp, and lambda_p only; stellar mass, planet mass, both mass-loss rates, and lambda* are fixed. Section 5.4 acknowledges that increasing the stellar wind relative to the planetary wind redirects both bubbles and streams into cometary tails, with the outcome depending on the ratio of mass-loss rates. So the Figure 5 collapse is established on a one-slice grid; it could shift or broaden when winds vary. The fix is cheap: soften the 'alone' language to 'in the absence of strong stellar wind shaping,' or run a few additional suites varying Mdot_p and Mdot_*. The population prediction in Figure 6 inherits this and the temperature-model dependence, but because the paper itself flags the temperature issue, I weigh that as minor rather than central.\n\nOne more minor point: the xi_los ~ ln(Ro_H) relation has no quoted scatter or uncertainty, so it is hard to know how predictive it is. Not fatal, but it would help.\n\nBottom line: a solid, honest paper with a useful framework. It should go to peer review, with minor-to-major revisions: qualify the headline claim, add a paragraph on expected sensitivity to winds (or a couple of runs), and report the scatter. The central argument about tidal shaping holds; the 'alone' is the misstep.","headline":"A genuinely useful Rossby-number framework for tidal shaping of exoplanet outflows, but the 'alone' claim overreaches the fixed-parameter simulation grid.","tokens_in":17093,"tokens_out":2722,"would_cite":true,"duration_ms":28230,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single dimensionless number, the Hill-sphere Rossby number, predicts whether a planet's escaping atmosphere forms a spherical bubble or a thin tidal stream.","keywords":["exoplanet atmospheric escape","tidal streams","Rossby number","transit spectroscopy","hydrodynamic simulations","Hill sphere","exoplanet outflows","phase-resolved kinematics"],"falsifier":"A decisive test would target a transiting evaporating planet whose predicted Rossby number is high ($\\mathrm{Ro}_H\\gtrsim 2$) and measure the phase-resolved velocity gradient of its metastable helium or Lyman-$\\alpha$ absorption with radiative-transfer modeling. The bubble prediction is $\\xi_\\mathrm{los}\\approx 1$, with outflow kinematics tracking the planet; observing instead a shallow or inverted gradient, or excess absorption that demands a geometrically thin stream, would falsify the single-parameter mapping proposed here.","tokens_in":16017,"feed_emoji":"🪐","tokens_out":11261,"duration_ms":100671,"temperature":0.7,"pith_summary":"Most escaping exoplanet atmospheres are not symmetric winds: depending on the system, the gas can leave as a roughly spherical bubble or be stretched into a thin stream that trails along the orbit. This paper argues that the shape is controlled by a single dimensionless number, the Hill-sphere Rossby number $\\mathrm{Ro}_H = c_s/(\\Omega r_H)$, which compares the outflow speed to the orbital shear across the planet's Hill sphere. Using three-dimensional gas-dynamic simulations, the authors show that high-$\\mathrm{Ro}_H$ outflows form bubbles while low-$\\mathrm{Ro}_H$ outflows form streams, and that this same number predicts the velocity gradient seen across transit. If the claim holds, transit spectra can be used to measure outflow temperatures directly, and the known exoplanet population should contain a mix of bubble-like and stream-like mass loss that observers must account for when interpreting absorption signals.","feed_headline":"Rossby number decides: bubble or stream for escaping planets","feed_subtitle":"The same number predicts the velocity gradient seen in transit, letting observers read outflow temperature from spectra.","key_machinery":"The central object is the Hill-sphere Rossby number, $\\mathrm{Ro}_H = c_s/(\\Omega r_H) = \\sqrt{3/\\lambda_H}$, where $c_s$ is the outflow sound speed, $\\Omega$ is the planet's orbital angular velocity, $r_H=(M_p/3M_\\ast)^{1/3}a$ is the Hill radius, and $\\lambda_H=GM_p/(c_s^2 r_H)$ is the Hill-sphere escape parameter. The number quantifies how much orbital shear and Coriolis acceleration divert the wind before it crosses the Hill sphere: when $\\mathrm{Ro}_H\\gg 1$ the wind wins and forms a bubble, when $\\mathrm{Ro}_H\\ll 1$ the tide wins and channels gas into streams. The companion diagnostic is the dimensionless transit velocity gradient $\\xi_\\mathrm{los}$, the mass-weighted line-of-sight velocity difference between egress and ingress normalized to the planet's own velocity gradient; the simulations show it tracks $\\ln \\mathrm{Ro}_H$, saturating at $\\xi_\\mathrm{los}\\to1$ for bubbles and falling toward zero or negative values for streams.","core_discovery":"On the paper's own terms, the central discovery is that the morphology and observable kinematics of a planetary outflow are determined, to good approximation, by the Hill-sphere Rossby number $\\mathrm{Ro}_H$. Defined as $c_s / (\\Omega r_H)$ with $c_s$ the outflow sound speed, $\\Omega$ the planet's orbital angular velocity, and $r_H$ the Hill radius, it is also $\\sqrt{3/\\lambda_H}$ in terms of the Hill-sphere escape parameter. In 27 hydrodynamic simulations spanning factors of four in orbital distance, planetary radius, and escape parameter, the authors find that flows with $\\mathrm{Ro}_H\\gtrsim$ a few are quasi-spherical and bounded by bow shocks, while flows with $\\mathrm{Ro}_H\\lesssim 1$ are channeled through the inner and outer Lagrange points into thin, dense streams. The dimensionless transit velocity gradient $\\xi_\\mathrm{los}$ collapses onto a single curve when plotted against $\\ln \\mathrm{Ro}_H$, saturating at $\\xi_\\mathrm{los}\\to1$ in the bubble limit and approaching zero or negative values in the stream limit. Applied to the known exoplanet population, the model predicts a continuum of shapes, with detected evaporating systems falling in stream-like parts of the diagram.","pith_inferences":["The same Rossby-number criterion should organize mass loss from any body embedded in a stronger tide, such as circumplanetary disks or donor stars in close binaries, provided the sound-speed estimate is replaced by the relevant outflow speed.","If real outflow temperatures are systematically cooler than the photo-evaporation models assume, the population diagram shifts many systems into the stream regime, and the observed sample may be biased toward bubbles simply because they are easier to detect near the planet.","A survey of a dozen evaporating planets spanning the predicted $\\mathrm{Ro}_H$ range would test the claimed collapse more strongly than the two stream systems available so far; the paper's own optically-thin caveat means synthetic spectra from the simulated densities are the right next check.","The opposite signs of $\\xi_\\mathrm{los}$ in the bubble and stream limits make a single well-measured transit of an evaporating planet a powerful morphology discriminator, suggesting a targeted observing program rather than a statistical one."],"forward_implications":["Phase-resolved transit spectra become a thermometer for escaping atmospheres: the sign and slope of $\\xi_\\mathrm{los}$ across transit constrain the outflow sound speed without assuming line widths are purely thermal.","Stream-like outflows create significant excess absorption outside optical transit, so simple in-transit versus out-of-transit subtraction biases measured line depths and mass-loss rates; observing baselines should be chosen from the predicted Rossby number.","The known exoplanet population should show a continuum of outflow geometries rather than one spherical-wind template, with bubble-like and stream-like systems both common.","In the stream limit the outflow is kinematically cold in energy–angular momentum phase space despite being spatially extended, so its material moves nearly on the planet's own orbit and can be treated as ballistic streams.","Detected stream-like systems fall where the diagram predicts streams, and a system with mostly in-transit absorption falls where the diagram predicts a confined outflow, giving initial support to the population forecast."],"supporting_citations":[{"why":"Supplies the numerical groundwork and tidal/orbital terms that this paper generalizes into the Rossby-number classification.","marker":"McCann et al. (2019)"},{"why":"Sets the methodology and the stellar-wind shaping baseline that the new bubble/stream framework sits alongside.","marker":"MacLeod & Oklopčić (2022)"},{"why":"Provides the observed stream-like HAT-P-32b/HAT-P-67b geometry that motivates cooler outflows and anchors the stream morphology.","marker":"Nail et al. (2024a)"},{"why":"Provides the catalog of photoevaporative outflow temperatures used to assign Rossby numbers to the known exoplanet population.","marker":"Linssen et al. (2024a)"},{"why":"Underlies the one-dimensional radiative-transfer modeling choices that set the predicted outflow temperatures and hence Rossby numbers.","marker":"Linssen et al. (2024b)"},{"why":"Supplies the numerical hydrodynamics solver used for the three-dimensional outflow simulations.","marker":"Stone et al. (2020)"},{"why":"Explains the stream-like flow through Lagrange points used to interpret sonic crossing near the Hill sphere.","marker":"Lubow & Shu (1975)"}],"fun_headline_variants":["Rossby number predicts planet outflow shape and spectra","Bubble or stream: one number decides planetary outflows","Hill-scale Rossby number controls escaping planet morphology","Outflow shape and transit gradient collapse to RosH curve","Single dimensionless number sets planet escape morphology"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise, adopted in Section 5.3, is that photo-evaporation models give each planet's outflow sound speed correctly; if those model temperatures run hot by the factor of about two that the observed stream systems already suggest, the predicted Rossby numbers and the claimed population split shift accordingly.","fun_headline_variants_meta":{"raw":{"variants":["Rossby number predicts planet outflow shape and spectra","Bubble or stream: one number decides planetary outflows","Hill-scale Rossby number controls escaping planet morphology","Outflow shape and transit gradient collapse to RosH curve","Single dimensionless number sets planet escape morphology"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000497,"raw_usage":{"total_tokens":2460,"prompt_tokens":993,"completion_tokens":1467,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":1394}},"tokens_in":609,"tokens_out":1467,"duration_ms":9863,"temperature":1.0,"reasoning_tokens":1394,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:03:12.681034+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would target a transiting evaporating planet whose predicted Rossby number is high ($\\mathrm{Ro}_H\\gtrsim 2$) and measure the phase-resolved velocity gradient of its metastable helium or Lyman-$\\alpha$ absorption with radiative-transfer modeling. The bubble prediction is $\\xi_\\mathrm{los}\\approx 1$, with outflow kinematics tracking the planet; observing instead a shallow or inverted gradient, or excess absorption that demands a geometrically thin stream, would falsify the single-parameter mapping proposed here.","supporting_citations":[],"review_version":1}