{"id":"03d710b6-ffe1-4d82-95c5-cfc6bd2ee032","arxiv_id":"2509.07048","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A hydrogen-atom variational calculation shows cold planets reach a maximum radius near Jupiter's radius at a mass of about 3.6 Jupiter masses.","lead":"The authors derive a simple formula, modeled on the hydrogen atom, for the largest size a cold planet can reach before gravity shrinks it. The formula predicts a maximum radius near Jupiter's and a mass of a few Jupiters, matching more detailed calculations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted coefficient in Eq. (14) depends on an unconstrained lattice-geometry parameter, so the 'parameter-free' quantitative claim is not robust.","rationale":"The reader identified the toy-model energy as the weakest assumption; this concern is a concrete instance of that, namely the unconstrained geometric relation between the variational length a and the volume per atom. The coefficient in Eq. (14) is not determined by fundamental constants alone but depends on the arbitrary β from Eq. (5). However, the qualitative conclusion — that cold planets have a maximum radius of order N*^{1/3} a0, which is about Jupiter's radius — survives for any O(1) β. The paper is explicitly a toy-model estimate, and the authors acknowledge its crudeness. Therefore, the concern does not overturn the reader's ACCEPT verdict; it would be a reviewer comment asking for a sensitivity check or a more careful statement that the numerical coefficient is model-dependent.","tokens_in":3219,"tokens_out":14907,"duration_ms":156497,"concrete_test":"Take Eq. (5) and replace with a = β(V/N)^{1/3}. Re-derive the maximization to obtain Rmax(β) and Mmax(β) as above. Then evaluate for β ∈ {0.25, 0.37, 0.5, 0.62, 1.0} and tabulate against the detailed-model values (1.165 RJ, 3.31 MJ). If the toy-model predictions bracket the detailed values but change by more than a factor of two across this physically plausible range, the specific coefficient in Eq. (14) is an artifact of the chosen lattice geometry rather than a fundamental prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The scale R ∝ N*^{1/3} a0 is robust, but the numerical coefficient in Eq. (14) is fixed only by the ad hoc relation Eq. (5), which sets a = (1/2)(V/N)^{1/3} (cubic lattice with spacing 2a). Replacing Eq. (5) with a general geometric relation a = β(V/N)^{1/3} changes Eqs. (14)–(15) to Rmax = (1/2)√(5/(4π)) β^{-3/2} N*^{1/3} a0 and Mmax = (5γ0/(3β))^{3/2} N* mH, where γ0 = (3/(4π))^{1/3}. For β = 0.25, 0.5, 0.62, 1.0 the predicted Rmax is roughly 2.1 RJ, 0.75 RJ, 0.53 RJ, and 0.26 RJ respectively; the detailed-model value 1.165 RJ corresponds to β ≈ 0.37. Thus the stated 'no free parameters' (or quantitative precision) is not justified unless the geometric factor is independently derived. This is a real soft spot, though the paper explicitly labels itself a toy model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a toy variational model for the maximum radius of cold hydrogen planets. A planet is approximated as N non-interacting hydrogen atoms on a cubic lattice of spacing 2a, with total energy N times the single-atom hydrogen expectation value evaluated at a rescaled length a, plus the uniform-density gravitational energy. Minimizing this energy with respect to a gives a 'modified Bohr radius' and hence R(N); maximizing R(N) over N yields Rmax = sqrt(5/(2pi)) N_star^(1/3) a0 ~ 0.75 RJ and Mmax = (5/3) sqrt(10/pi) N_star mH ~ 3.6 MJ. These are compared with detailed cold-sphere models (Rmax ~ 1.165 RJ, Mmax ~ 3.31 MJ) and the agreement is described as surprisingly good.","tokens_in":3535,"tokens_out":11716,"duration_ms":108353,"significance":"The qualitative insight that the Jupiter-scale maximum radius follows from a balance between electrostatic, gravitational, and quantum kinetic energies is appealing and pedagogically valuable. The variational algebra is internally consistent, and the scaling Rmax ~ N_star^(1/3) a0 is robust; the comparison with detailed models is made after the derivation, so the logic is not circular. The paper has no fitted parameters, and the presentation is transparent. However, the numerical coefficient of the scaling law is fixed by an unstated geometric assumption, and the many-body energy model is asserted rather than derived, so the 'quantitative' claim is weaker than the abstract implies.","major_comments":[{"comment":"The numerical coefficients in Eqs. (14) and (15) depend on the arbitrary relation a = (1/2)(V/N)^(1/3). Replacing it with a general relation a = beta (V/N)^(1/3) changes the results to Rmax = (1/2) sqrt(5/(4pi)) beta^(-3/2) N_star^(1/3) a0 and Mmax = ((5 gamma0)/(3 beta))^(3/2) N_star mH with gamma0 = (3/(4pi))^(1/3). For beta = 0.25, 0.5, 0.62, 1.0, Rmax is approximately 2.1, 0.75, 0.53, and 0.26 RJ, respectively. Thus the 'no free parameters' or quantitative precision claim is not justified unless the geometric factor is independently derived. The paper's toy-model framing mitigates this, but the text should explicitly state that the coefficient carries an O(1) geometric uncertainty and that the comparison with 1.165 RJ is order-of-magnitude, not a precise prediction.","section":"Section II, Eq. (10)"},{"comment":"The total energy E = N [hbar^2/(2ma^2) - tilde(e)^2/(4pi epsilon_0 a)] is assumed, not derived. The paper does not specify a many-body Hamiltonian for the collection of atoms, so the variational calculation is not an upper bound on the true ground-state energy of the planet; it is a heuristic energy model. This is a legitimate simplification for an order-of-magnitude paper, but the wording 'variational principle very similar to hydrogen atom' can overstate the rigor. Please add an explicit sentence clarifying that Eq. (10) is a model energy, not the expectation value of an actual many-body Hamiltonian.","section":"Section II, Eq. (10)"}],"minor_comments":[{"comment":"The integral notation in Eq. (1) is cramped; please define dv as the volume element and add spaces for readability.","section":"Section II, Eq. (1)"},{"comment":"The line 'H psi_0(r) = E_i psi_0(r)' should likely be 'H psi_0(r) = E_0 psi_0(r)' or the subscript i should be defined.","section":"Section II, text after Eq. (2)"},{"comment":"The text calls tilde(e) a 'modified charge,' but Eq. (11) defines tilde(e)^2. Please say 'modified squared charge' to avoid confusion.","section":"Section II, Eq. (11)"},{"comment":"The concluding paragraph says the approach applies to white dwarfs and iron cores, but the derivation uses hydrogen atoms with a specific mH and N_star. The scaling may generalize, but the numerical maximum mass is a Jupiter mass, not a white-dwarf mass; clarify the intended scope.","section":"Conclusions"},{"comment":"The phrase 'the size starts to decrease instead of increase' would be clearer as 'the size starts to decrease instead of continuing to increase.'","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is a short pedagogical contribution. The central concern is the arbitrary lattice-geometry factor beta; this is fixable by softening the quantitative claim and explicitly stating the O(1) uncertainty. If the target journal is astro-ph.IM, consider whether the paper might be better suited to an education-oriented venue. The authors are clearly aware that this is a toy model, but the abstract and Eqs. (14)-(15) currently overstate the precision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a nice little paper that does what it says—gives a variational argument for why cold planets top out around Jupiter's radius, using a hydrogen-atom analogy. The derivation is clean, the assumptions are stated up front, and it's honest about being a toy model. Worth a look if you teach intro astro or quantum mechanics.\n\nWhat's actually new: the closed-form expressions for Rmax and Mmax, Eqs. (14)-(15), and the recasting of the problem as a 'modified Bohr radius' with an effective charge. These aren't in the cited Zapolsky-Salpeter or Fortney et al. work. The result that Rmax scales as N*^{1/3} a0 and Mmax scales as N* mH is a nice scaling insight, and it reproduces known numbers to order of magnitude.\n\nWhere it's soft: the numerical coefficient in Eq. (14) is not as 'no free parameters' as it looks. The paper sets the lattice spacing a = (1/2)(V/N)^{1/3} for a cubic lattice, but that's an assumption. If you use a general geometric relation a = β(V/N)^{1/3}, the coefficients change: Rmax ∝ β^{-3/2}. For β ranging from 0.25 to 1.0, Rmax goes from 2.1 RJ to 0.26 RJ. The detailed-model value 1.165 RJ corresponds to β≈0.37. So the 'parameter-free' claim is only as good as that geometric guess. The paper could be more careful here; the scaling is robust, but the coefficient is not. That said, the paper explicitly calls itself a toy model and only claims order-of-magnitude agreement, so this is a moderate caveat, not a fatal flaw.\n\nThe energy model itself (N times a hydrogen atom with rescaled length, plus uniform-density gravity) is admittedly ad hoc, but the authors don't oversell it. For a teaching paper, that's fine. The comparison with detailed models is fair and they note the composition issue.\n\nOverall: the central argument holds up for what it claims—an accessible explanation of the Jupiter-scale radius limit. The coefficient sensitivity is real but secondary. I'd bring this to a reading group interested in pedagogical derivations, and I'd cite it as an example of a simple variational estimate if the occasion arises. It deserves peer review—send it to a journal like AJP or Eur. J. Phys., where the audience will appreciate the didactic value. Accept.","headline":"A clean toy-model derivation that explains Jupiter's radius cap with a hydrogen-atom analogy, but the advertised 'no free parameters' coefficient actually depends on an unstated lattice-geometry assumption.","tokens_in":3986,"tokens_out":2605,"would_cite":false,"duration_ms":25853,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Cold planets top out at Jupiter size by a hydrogen-atom argument","keywords":["cold planets","maximum radius","Jupiter radius","variational principle","hydrogen atom","gravitational compression","mass-radius relation","white dwarfs"],"falsifier":"Compute the zero-temperature mass-radius relation for pure hydrogen using a realistic dense-matter equation of state. If the maximum radius is not within tens of percent of 0.75 RJ, or occurs at a mass far from 3.6 MJ, then the single-rescaled-Bohr-radius toy model is not the mechanism behind Jupiter's maximum radius.","tokens_in":3127,"feed_emoji":"🪐","tokens_out":3803,"duration_ms":40313,"temperature":0.7,"pith_summary":"This paper tries to show that the observed fact that planets have radii no larger than Jupiter's is not a coincidence of planetary formation but a consequence of the same variational principle that fixes the size of a hydrogen atom. Treating a cold pure-hydrogen planet as N non-interacting hydrogen atoms squeezed by gravity, the authors write the total energy as a function of a single scale a, balance the electrostatic, gravitational, and quantum kinetic terms, and derive a closed-form maximum radius Rmax = sqrt(5/(2π)) N*^(1/3) a0 ≈ 0.75 RJ and corresponding mass Mmax ≈ 3.6 MJ. The gravity-to-electric ratio enters only through the dimensionless number N* = (e²/(4πϵ0 G mH²))^{3/2}, which is fixed by fundamental constants. A sympathetic reader would care because the result explains a Jupiter-scale upper limit without free parameters and is within tens of percent of more detailed cold-sphere models. The paper is explicitly a pedagogical order-of-magnitude estimate, not a full equation-of-state calculation.","feed_headline":"A hydrogen-atom trick caps cold planet radius at 0.75 Jupiter","feed_subtitle":"A variational energy balance with no free parameters gives maximum radius ~0.75 Jupiter and mass ~3.6 Jupiter masses.","key_machinery":"The central object is the dimensionless ratio N* = (e²/(4πϵ0GmH²))^{3/2}, whose 2/3 power is the ratio of electric repulsion to gravitational attraction between two protons. It carries the argument by replacing the single charge e in the hydrogen-atom energy with an effective charge e-tilde that grows with N; the Bohr-radius minimization then becomes a radius maximization, and the maximum radius and mass are algebraic functions of N*, a0, and mH. The companion object is the trial wavefunction e^(−r/a): it lets compression be represented by rescaling one length scale a.","core_discovery":"On the paper's own terms, the discovery is that the maximum radius and mass of a cold planet follow from a variational calculation nearly identical to the hydrogen atom. The trial wavefunction e^(−r/a) gives a single-atom energy E(a); packing N such atoms in a cubic lattice of spacing 2a and adding the uniform-density gravitational energy −3GM²/(5R) changes the effective charge e² in the Coulomb term to e-tilde² = e²[1 + (3/5)(πN²/(6N*²))^{1/3}]. Minimizing the total energy over a yields a modified Bohr radius, and maximizing R(N) = a0(6N/π)^{1/3}/[1 + (3/5)(πN²/(6N*²))^{1/3}] gives Eqs. (14)–(15). Numerically, Rmax ≈ 0.75 RJ and Mmax ≈ 3.6 MJ, which the authors compare with more sophisticat","pith_inferences":["Including helium, the actual gas-giant composition, would change the effective particle mass and mean charge, shifting the numeric prefactor but leaving the N* structure intact; one test is whether a helium-fraction-corrected version lands closer to the observed ~1.2 RJ.","The framework suggests that any cold object made of atoms with Bohr radius a0 and particle mass mH has a maximum radius proportional to a0 times a fundamental-constant ratio; comparing objects made of heavier atoms would probe this scaling.","The maximum-radius mass may show up as a characteristic turnover in the mass–radius diagram of cool exoplanets, so a survey of transiting planets could test whether the turnover sits near a few Jupiter masses."],"forward_implications":["Below N*, radius grows as N^{1/3} and volume scales linearly with particle number; above N*, radius shrinks as N^{−1/3}, so adding mass compresses the planet.","The maximum mass is independent of ħ: only the radius scale depends on quantum mechanics through the Bohr radius a0.","The same energy balance applies to cold degenerate objects without fusion, such as cool white dwarfs and iron cores, giving a qualitative upper radius for those objects too.","The derivation produces a parameter-free estimate of the Jupiter scale—0.75 RJ and 3.6 MJ for pure hydrogen—with more detailed models giving 1.165 RJ and 3.31 MJ.","The toy model is designed to be usable in undergraduate courses as an accessible derivation of why gas giants share a common radius scale."],"supporting_citations":[{"why":"Supplies the detailed cold-sphere mass-radius relation for low-mass hydrogen spheres that the toy model's 0.75 RJ and 3.6 MJ results are compared against.","marker":"1"},{"why":"Supplies the more detailed planetary radius models across mass and stellar insolation used as the comparison for gas-giant composition.","marker":"2"}],"fun_headline_variants":["Hydrogen-atom trick explains planet size cap","Cold planets stop growing at 0.75 Jupiter radii","Variational math sets cold planet max radius","Quantum analogy predicts planet maximum size","Cold planet max radius pinned at 0.75 Jupiter"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The argument rests on assuming a cold pure-hydrogen planet can be described as N non-interacting hydrogen atoms whose only response to compression is a rescaling of the single Bohr radius a, with uniform density for gravity; if dense hydrogen's real equation of state makes interactions or degeneracy dominate, the derived maximum radius and mass lose their grounding.","fun_headline_variants_meta":{"raw":{"variants":["Hydrogen-atom trick explains planet size cap","Cold planets stop growing at 0.75 Jupiter radii","Variational math sets cold planet max radius","Quantum analogy predicts planet maximum size","Cold planet max radius pinned at 0.75 Jupiter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1393,"prompt_tokens":640,"completion_tokens":753,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":384,"completion_tokens_details":{"reasoning_tokens":696}},"tokens_in":384,"tokens_out":753,"duration_ms":7756,"temperature":1.0,"reasoning_tokens":696,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T23:06:12.751309+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the zero-temperature mass-radius relation for pure hydrogen using a realistic dense-matter equation of state. If the maximum radius is not within tens of percent of 0.75 RJ, or occurs at a mass far from 3.6 MJ, then the single-rescaled-Bohr-radius toy model is not the mechanism behind Jupiter's maximum radius.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the detailed cold-sphere mass-radius relation for low-mass hydrogen spheres that the toy model's 0.75 RJ and 3.6 MJ results are compared against."}],"review_version":1}