{"id":"27d0c286-6296-4f19-8650-644eb8796d7a","arxiv_id":"2412.02455","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A Puiseux-series analysis of the Kottler metric produces a taxonomy of massive-particle orbits and a galaxy-rotation-based upper bound on a negative cosmological constant.","lead":"This paper classifies the possible orbits of massive particles around a central mass when a cosmological constant is added to Einstein's equations, using power-series expansions to organize the trajectory types. It also derives an upper bound on a negative cosmological constant from galaxy rotation speeds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The asserted root r=rg in Sec. 4 is absent from the claimed finite-motion configurations: p(rg)>0 there, so the elliptic-integral orbit formulas and the trajectory classification are unsupported.","rationale":"The reader's weakest-assumption call is correct, and the problem is sharper than a single counterexample. The algebra of p(r_g) is independent of the approximate Puiseux table; it uses only the definitions of K and the fact that A/B/C configurations have j > r_g. Since p(r_g) > 0, the 'root at r_g' used in Sec. 4 cannot exist. I therefore agree with the reader's REJECT verdict and with the identification of the load-bearing flaw. The §2 rotational-curve bound is an additional independent error: for Λ = 0 the claimed formula does not reduce to the Keplerian circular speed, so the bound (2.5) is not trustworthy. But the root-at-rg error attacks the advertised classification itself, so it is the most load-bearing concern. No refinement of the rotation-curve estimate would repair the elliptic-integral derivation. I set verdict_should_be to UNCHANGED because this stress-test does not move the reader's verdict; it reinforces it.","tokens_in":9693,"tokens_out":13934,"duration_ms":135288,"concrete_test":"Take Λ = -10^{-8}, r_g = 1, j = 10, K = 0.005, a parameter set satisfying the configuration C (Λ<0, K>0) conditions in Table 2. Evaluate Eq. (4.1) at r = r_g: p(1) ≈ 1.005, not zero. Then numerically count and order the positive roots of p(r) for these parameters with a Sturm sequence and compare with the factorization used in Eq. (4.3). If p(r_g) ≠ 0 and the root structure differs, the Sec. 4 reduction to elliptic integrals is not valid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Substituting r = r_g into Eq. (4.1) for Λ = -|Λ| gives p(r_g) = r_g^3 (k^2 + (|Λ|/3)(j^2 - r_g^2)), because K = k^2 - 1 + |Λ|j^2/3. In the finite-motion configurations A, B, C of Table 2 the stated conditions force j > r_g: A and C list j > r_g explicitly, and B implies it from j^2 > K r_g and j^2 K > r_g. Therefore p(r_g) > 0 in every such configuration. The Sec. 4 assertion that 'the root r = r_g is present only once in each of those configurations' is false. Consequently the integrand of Eq. (4.3) contains no (r - r_g) factor, so the expansion (4.4) and the elliptic-integral formulas (4.5)-(4.7) do not follow from the trajectory equation. The classification's completeness also rests on leading-order Puiseux terms in Table 2 rather than a demonstration of exact root counts and ordering, but the root-at-rg contradiction is enough by itself to break the central construction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies massive-particle orbits in the Kottler metric with a cosmological constant Λ. Starting from the Hamilton–Jacobi equation and the Binet equation, it derives a trajectory equation whose radicand is the polynomial F(r) in Eq. (1.11), equivalently the quintic p(r)=F(r)/r in Eq. (4.1). For both signs of Λ it proposes a classification of non-circular trajectories based on leading-order Puiseux expansions of the roots of p(r) (Table 2), singles out three configurations A, B, C that permit finite motion, and writes the corresponding orbits as combinations of incomplete elliptic integrals. For negative Λ it also proposes a circular-velocity formula (2.2) and an upper bound |Λ| ≲ 10^-13 v_min^6 / M^2 from galaxy rotation curves.","tokens_in":9950,"tokens_out":16770,"duration_ms":171278,"significance":"The topic is of continuing interest, and the paper has some useful ingredients: the reduction to a single quintic is standard and correctly carried out, the Newton-polygon/Puiseux viewpoint is a legitimate heuristic for organizing parameter regimes, and the explicit elliptic-integral formulas, if valid, would be convenient. The authors also provide a promising heatmap tool for scanning the (j, K) parameter plane. However, the main results as stated are not reliable: the factorization underlying the non-circular orbit integrals assumes a root at r = r_g that is absent in the very configurations to which it is applied, the rotational-curve formula has a wrong Λ = 0 limit, and the claimed completeness of the classification is not demonstrated. These are load-bearing defects, not presentation issues.","major_comments":[{"comment":"The assertion after Eq. (4.3) that “the root r = r_g is present only once in each of those configurations and it’s the smallest positive root” is false for the Λ < 0 finite-motion configurations to which the paragraph refers. Substituting r = r_g into Eq. (4.1) with Λ = -|Λ| and K = k^2 - 1 + |Λ| j^2 / 3 gives p(r_g) = r_g^3 [k^2 + (|Λ|/3)(j^2 - r_g^2)]. Configurations A and C explicitly require j > r_g, so p(r_g) > 0 there; configuration B, while not requiring j > r_g explicitly, lists r_g as a root in Table 2 without any condition that would make p(r_g) vanish. Therefore r_g is not a root of p(r) in the cases where the text asserts it is, the factor (r - r_g)^(-1/2) in Eq. (4.4) is absent, and the elliptic-integral formulas (4.5)–(4.7) do not follow from the trajectory equation.","section":"Sec. 4, Eqs. (4.3)–(4.7) and Table 2"},{"comment":"Eq. (2.2) cannot be the classical circular-speed formula. For Λ = 0 it reduces to v^2/c^2 = (1/2)(r_g/r)^2 - (3/2)(r_g/r), which is negative for every r > r_g, rather than the Keplerian circular speed. The derivation of Eq. (2.2) from Eq. (2.1) is not shown, and the subsequent minimum-speed expressions (2.3)–(2.4) and the bound (2.5) inherit this problem. Moreover, substituting r_g = 2GM/c^2 into the preceding formula in Eq. (2.5) gives a prefactor 32/9 with a factor 1/(G^2 M^2 c^2), not 1/(G^2 c^4), and the simplified form “≈10^-13 v_min^6/M^2” is dimensionally inconsistent unless additional factors are intended.","section":"Sec. 2, Eqs. (2.2)–(2.5)"},{"comment":"The completeness of the proposed classification is not established. Table 2 lists only the first term of each Puiseux expansion, and the conditions in the table are Newton-polygon inequalities for term dominance. These determine the leading asymptotic form of roots, not the exact number, reality, or ordering of the positive roots of the quintic (4.1) for all parameter values. The Appendix states that the method is intended for “extracting valuable information about all possible sets of roots” and not for precise formulas, yet the paper’s central claim is a complete classification. No theorem or exhaustive numerical verification is supplied to close this gap. Because the finite-motion regions and the orbit formulas are read off from the root configurations, this gap is load-bearing.","section":"Sec. 3, Table 2, and Appendix"}],"minor_comments":[{"comment":"The horizon equation is written as 1 - r_g/r - Λ r^2 = 0, while the metric component (1.10) contains Λ r^2 / 3; the relation between r_g, the mass parameter, and the actual event horizon should be stated consistently.","section":"Eq. (1.12)"},{"comment":"The layout is extremely dense because conditions and roots for four sign combinations share one row; separating the four cases into distinct columns or tables and marking which listed roots are real and positive would make the classification checkable.","section":"Table 2"},{"comment":"The heatmap lacks axis labels, numerical scales, and the fixed values of Λ and r_g; without these the claimed parameter-scanning tool cannot be reproduced.","section":"Fig. 7"},{"comment":"The symbols X_s, R_n, V_n, α, α_1, and the elliptic-integral arguments are introduced abruptly; please define every symbol and state which roots correspond to a, b, c, d.","section":"Eqs. (4.6)–(4.7)"},{"comment":"The manuscript contains several typographical and grammatical errors (e.g., “rotatonal curves”, inconsistent capitalization, duplicated Russian/English abstracts); a thorough language edit is needed.","section":"Throughout"}],"recommendation":"reject","confidential_remarks":"For the editor: the paper is not ready for publication in its present form. The Sec. 4 root-at-r_g error is a definitive algebraic contradiction, not merely a missing proof, and the Sec. 2 rotation-curve derivation appears inconsistent in the Λ→0 limit. Even with these local fixes, the completeness claim would require either a rigorous root-counting argument or an exhaustive numerical verification of Table 2. I recommend rejection rather than major revision, because the central construction would need to be reworked."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper is not ready for publication. Its central classification of timelike orbits in the Kottler metric relies on a root that is not actually a root, and the galaxy-rotation bound on Lambda is dimensionally inconsistent. That said, the Puiseux-expansion approach to the trajectory quintic is genuinely new to me, and the closed-orbit formula for Lambda<0, rg=0 is worth a second look.\n\nWhat the paper does well: it sets up the trajectory equation cleanly from the Hamilton-Jacobi equation, writes down the quintic p(r), and uses Newton polygons to organize parameter space into root configurations. The claim of a new stable circular root for large angular momentum in the anti-de Sitter case is plausible and not present in the cited references. Equation (4.10), for rg=0, gives an explicit closed orbit in negative Lambda with K>0, which looks like a real result. The authors also admit some limitations and offer to share their script.\n\nThe soft spots are serious, not minor. First, the rotation-curve section: equation (2.2) has no derivation, and for Lambda=0 it gives v^2/c^2 = (1/2)(rg/r)^2 - (3/2)(rg/r), which is negative at large r and does not reduce to the Keplerian or Schwarzschild circular speed. The bound (2.5) has units that do not work out to 1/m^2; the numerical prefactor is meaningless without a proper derivation. Second, and more damaging, Section 4 claims that r=rg is a root of p(r) in configurations A, B, and C. From the paper's own equation (4.1), p(rg)=rg^3(k^2+|Lambda|(j^2-rg^2)/3), and the listed conditions force j>rg in those configurations, so p(rg)>0. The expansion (4.4) and the elliptic-integral formulas (4.5)-(4.7) therefore do not follow. Third, Table 2's root conditions come from leading-order Puiseux terms; there is no proof that these give the exact number and ordering of positive roots, so the completeness claim is not established. The paper also does not compare with the existing exact geodesic classifications of Schwarzschild-(anti-)de Sitter spacetimes, which weakens the novelty.\n\nWho this is for: a specialist in exact geodesics might find the Puiseux taxonomy useful as a heuristic, and the rg=0 orbit formula could be a nice exercise. But as a contribution to the literature it needs major revision: drop or fix the rotation-curve bound, remove the false root-at-rg expansion, and compare with known classifications.\n\nMy recommendation: I would not publish this version. If it crossed my desk, I would send it to a referee mainly to check whether the Puiseux method and the rg=0 result can be salvaged; the current form should not be accepted. I would not cite it in its present state.","headline":"A novel Puiseux-root approach to Kottler orbits is undercut by a false root-at-rg expansion and a dimensionally inconsistent rotation-curve bound; the rg=0 closed orbit is worth a look, but the paper is not publishable as is.","tokens_in":10443,"tokens_out":7319,"would_cite":false,"duration_ms":75985,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.10.Eg","98.80.Es"],"model":"deepseek-v4-flash","headline":"The paper claims a complete classification of massive-particle trajectories in the Kottler metric for both signs of the cosmological constant, obtained by Puiseux-series root expansions, plus an upper bound on negative Lambda from galaxy…","keywords":["Kottler metric","cosmological constant","massive particle orbits","Puiseux series","Newton polygon","elliptic integrals","galaxy rotation curves","general relativity"],"falsifier":"Evaluate $p(r)$ for the parameter set $\\Lambda=-1$, $r_g=1$, $j=2$, $K=4/3$ in the paper's units: $p(1)=2\\neq 0$, so the root $r=r_g$ is absent there, contradicting the assumption that $r_g$ is one of the positive roots in the classification. A denser check is to compute the positive roots of $p(r)$ over a grid of the Table 2 conditions and compare their count and ordering with the claimed configurations.","tokens_in":9488,"feed_emoji":"🌀","tokens_out":14318,"duration_ms":140801,"temperature":0.7,"pith_summary":"This paper tries to establish how the cosmological constant $\\Lambda$ changes every possible trajectory of a massive test particle orbiting a spherically symmetric central body. Working in the Kottler metric, the authors reduce the motion to a quintic polynomial $p(r)$ whose positive roots are turning points, and claim that expanding the associated algebraic curve in Puiseux series—with Newton-polygon inequalities—gives a complete classification of orbits for both signs of $\\Lambda$. The classification assigns each parameter set to one of eight root configurations and determines which configurations allow finite, bound, spiral, or unbound motion, with explicit trajectory formulas built from incomplete elliptic integrals. It also turns galaxy rotation-curve data into an upper bound $|\\Lambda|\\lesssim 10^{-13}\\, v_{\\min}^6/M^2$ (SI units) for negative $\\Lambda$. If the classification is right, orbit observations become a direct probe of the sign and magnitude of the cosmological constant.","feed_headline":"All massive-particle orbits in the Kottler metric, classified","feed_subtitle":"A Puiseux-root scheme maps every parameter set to a finite, spiral, or unbound orbit, and bounds negative Λ from galaxy data.","key_machinery":"The object that carries the argument is the trajectory polynomial $p(r)$: its positive roots are the turning points of the orbit, and the number, reality, and ordering of those roots decide whether the orbit falls inward, escapes, spirals, or oscillates between two radii. The paper treats $p(r)$ as an algebraic curve in its coefficients and analyzes it with the Newton-Puiseux method: each monomial contributes a point to a Newton polygon, the slopes of the polygon edges give the leading Puiseux exponents of the roots, and the inequalities defining each convex hull become the parameter conditions in Table 2. The same expansion isolates the factor $(r-r_g)^{-1/2}$, which is expanded as $1/\\sqrt{r}+(r_g/2)r^{-3/2}+\\dots$ so that the orbit integral reduces to a sum of incomplete elliptic integrals of the first, second, and third kinds.","core_discovery":"The central claim is that the quintic $p(r)=(\\Lambda/3)r^5+K r^3+r_g r^2-j^2 r+j^2 r_g$ (with $\\Lambda$ replaced by $-|\\Lambda|$ when it is negative) determines all non-circular orbits of massive particles in the Kottler metric, and that Puiseux expansion of this polynomial's algebraic curve yields every possible configuration of its positive roots. For $\\Lambda<0$ and $\\Lambda>0$ the paper lists five numbered and three lettered configurations (Table 2), with inequalities that decide how many roots are real, positive, and ordered. Motion is possible only where $p(r)>0$; in configurations A, B, and C a particle is confined between two turning radii, and the trajectory integral is evaluated as a linear combination of incomplete elliptic integrals, with explicit formulas (4.6) and (4.7). In the limit $r_g\\to 0$ the paper derives a closed precessing orbit (4.10)–(4.12), and for $\\Lambda<0$ it derives from the rotation curve the bound $|\\Lambda|\\lesssim 10^{-13}\\, v_{\\min}^6/M^2$ in SI units, consistent with earlier estimates for a $10^5\\,M_\\odot$ black hole.","pith_inferences":["Beyond the paper, the Table 2 conditions could be stress-tested by a random numerical scan of the positive roots of $p(r)$; any mismatch between the computed roots and the claimed configuration would reveal a missing or mislabeled case.","Beyond the paper, the bound implies a mass–velocity relation $v_{\\min}\\propto M^{1/3}$ for galaxies if $\\Lambda$ is universal, so the lower envelope of a $v_{\\min}$ versus $M$ plot from real rotation-curve catalogs is a direct observational test.","Beyond the paper, the Newton-polygon approach is not tied to this particular metric and could classify orbits for other spherically symmetric line elements whose trajectory polynomial has higher degree."],"forward_implications":["For negative $\\Lambda$, every rotation curve has a minimum circular speed $v_{\\min}$ at a radius near $(3r_g/(4|\\Lambda|))^{1/3}$, and observed speeds above $v_{\\min}$ translate into the bound $|\\Lambda|\\lesssim 10^{-13}\\, v_{\\min}^6/M^2$ (SI units).","For negative $\\Lambda$, a stable circular orbit switches from the Schwarzschild-like radius $r\\approx 2j^2/r_g$ to a cosmological-constant-dominated radius $r\\approx (3j^2/|\\Lambda|)^{1/4}$ once the dimensionless angular momentum exceeds a critical value.","For positive $\\Lambda$, unbound trajectories outside $r_{\\max}=\\sqrt{3|K|/\\Lambda}$ are modified hyperbolic spirals when $K>0$ and modified hyperbolas when $K<0$, both covered by a single formula.","In the confined configurations A, B, and C, trajectories between $r_{\\min}$ and $r_{\\max}$ are explicit linear combinations of incomplete elliptic integrals, so any concrete parameter set can be integrated without numerical ODE solving.","With $r_g=0$, negative $\\Lambda$, and positive $K$, the motion is a closed bound orbit; a nonzero $r_g$ adds perihelion precession whose magnitude is given by an elliptic integral."],"supporting_citations":[{"why":"Introduces the Kottler metric, the spacetime whose massive-particle motion the paper classifies.","marker":"[2]"},{"why":"Supplies galaxy rotation-curve data and an earlier mass/Lambda estimate that the upper bound in Sec. 2 is compared with.","marker":"[5]"},{"why":"Provides Descartes' rule of signs used to count the possible positive roots of p(r) for each sign of Lambda.","marker":"[18]"},{"why":"Gives the handbook formulas for incomplete elliptic integrals used in the trajectory solutions (4.6), (4.7), and (4.12).","marker":"[19]"},{"why":"Gives the computational procedure for Puiseux expansion that the paper adapts into a root-finding method.","marker":"[20]"},{"why":"States the Newton-Puiseux theorem for algebraic curves, the theoretical basis for expanding the roots of p(r).","marker":"[21]"},{"why":"Provides the Binet equation used to derive the effective potential and the rotation-curve relation.","marker":"[22]"}],"fun_headline_variants":["Every massive Kottler orbit, mapped via Puiseux roots","One quintic polynomial classifies all Kottler orbits","Negative Λ bounded by galaxy curves in Kottler study","Full orbit classification in Kottler metric with Λ"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classification hangs on the unproved premise that the leading Puiseux-root approximations and the Newton-polygon inequalities in Table 2 always give the true number, reality, and ordering of the positive roots of the quintic—and that one of those roots is exactly $r=r_g$ so the $(r-r_g)^{-1/2}$ expansion can be used.","fun_headline_variants_meta":{"raw":{"variants":["Every massive Kottler orbit, mapped via Puiseux roots","One quintic polynomial classifies all Kottler orbits","Negative Λ bounded by galaxy curves in Kottler study","Full orbit classification in Kottler metric with Λ"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000451,"raw_usage":{"total_tokens":2258,"prompt_tokens":919,"completion_tokens":1339,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":1271}},"tokens_in":535,"tokens_out":1339,"duration_ms":15044,"temperature":1.0,"reasoning_tokens":1271,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:26:49.694753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate $p(r)$ for the parameter set $\\Lambda=-1$, $r_g=1$, $j=2$, $K=4/3$ in the paper's units: $p(1)=2\\neq 0$, so the root $r=r_g$ is absent there, contradicting the assumption that $r_g$ is one of the positive roots in the classification. A denser check is to compute the positive roots of $p(r)$ over a grid of the Table 2 conditions and compare their count and ordering with the claimed configurations.","supporting_citations":[{"cited_title":"Über die physikalischen Grundlagen der Einsteinschen Gravitationstheorie.Annalen der Physik , 1918, vol","cited_arxiv_id":null,"evidence_quote":"Introduces the Kottler metric, the spacetime whose massive-particle motion the paper classifies."},{"cited_title":"Estimation of mass and cosmological constant of nearby spiral galaxies using galaxy rotation curve","cited_arxiv_id":"1307.1824","evidence_quote":"Supplies galaxy rotation-curve data and an earlier mass/Lambda estimate that the upper bound in Sec. 2 is compared with."},{"cited_title":"Handbook of Elliptic Integrals for Engineers and Scientists","cited_arxiv_id":null,"evidence_quote":"Gives the handbook formulas for incomplete elliptic integrals used in the trajectory solutions (4.6), (4.7), and (4.12)."},{"cited_title":"How to Compute a Puiseux Expansion","cited_arxiv_id":"0807.4674","evidence_quote":"Gives the computational procedure for Puiseux expansion that the paper adapts into a root-finding method."},{"cited_title":"Classical Mechanics","cited_arxiv_id":null,"evidence_quote":"Provides the Binet equation used to derive the effective potential and the rotation-curve relation."}],"review_version":1}