{"id":"b5eb9abf-7678-49cf-a720-d63eccfc8176","arxiv_id":"2507.20383","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive, for the first time, Bošković's spherical trigonometric equations for solar rotation elements and confirm them against his 1777 sunspot data.","lead":"This paper derives the equations behind Bošković's 1785 spherical trigonometric method for computing the Sun's rotation axis and period from three sunspot positions. The authors show that these equations reproduce Bošković's results from his 1777 sunspot observations, providing a compact single-procedure alternative to his earlier methods.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Load-bearing concern: the derivation assumes P′C = P′C′ = P′C′′ (constant heliographic latitude, rigid rotation), yet the paper's own Table 6 shows i and Ω vary by 15° and 129° across the 20 triples; the 'no approximation' claim therefore depends on the same-parallel assumption and on the chosen…","rationale":"The reader's weakest assumption is correct and is the main threat to the central claim. The paper's reconstruction of Bošković's spherical method is algebraically coherent, and the agreement with Bošković's arcminute values for the chosen triple is real. But the claim that the method is complete and 'without any approximation' depends on the equal-distance identity, which is physically an idealization. The all-triples test is decisive because, under the model, any three positions of the same rigidly rotating spot should yield the same elements; the large scatter in Table 6 indicates that the selected triple, not the equations alone, is responsible for the agreement. This does not invalidate the historical reconstruction, but it requires the authors to state the assumption, propagate errors, and soften the precision and no-approximation claims. Hence the reader's CONDITIONAL verdict stands; no change is needed.","tokens_in":22756,"tokens_out":10798,"duration_ms":119450,"concrete_test":"Apply the paper's spherical-trigonometric equations (Eqs. 29, 30, 34, 36) to all 20 triples of the six positions in Table 1, not only the selected [136]. If the spread in i, Ω, and T′ exceeds the propagated arcminute-level measurement errors (roughly 0.01° and 0.01 d), the same-parallel assumption is not satisfied by the data and the 'single procedure, no approximation' claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations 26, 34, 43, 50–51 and the mirror symmetry of triangles △9 and △6 in §2.4.1 all rely on the equality P′C = P′C′ = P′C′′. This equality is not derived from the measurements; it is an assumption that the spot remains on a single heliographic parallel while the Sun rotates rigidly. The Conclusion's statement that the spherical solution is 'closed equations without any approximation' is therefore not supported. The assumption is quantitatively important: Table 6, computed with the VFM analogue, gives i between 3.512° and 18.696° and Ω between −41.197° and 87.811° for the 20 triples of the six observed positions, while the selected triple [136] gives i = 6.807° and Ω = 74.048°. Such scatter shows that the result is controlled by the choice of a triple that approximately satisfies the same-parallel condition (C1 ≈ C6, C3 at minimum latitude), not by a robust single-procedure inversion. Bošković's published arcminute values cannot validate the assumption, and the 10⁻⁵ arcsec precision quoted in §4.3 is inconsistent with arcminute input data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reconstructs the trigonometric spherical solution that Ruđer Bošković described verbally in his 1785 Opuscule II for determining the solar rotation elements (i, Ω, and the sidereal period T′) from three observed positions of the same sunspot. The authors derive a sequence of spherical-trigonometric equations for the full and the “short” solutions, apply them to Bošković's September 1777 sunspot positions (positions 1, 3, and 6), and obtain i = 6.80728° (≈6°48′26″), Ω = 74.04774° (≈74°03′), and T′ = 26.806 d (≈26.81 d), closely reproducing Bošković's published values from the planar-trigonometric solution. They also compare with a contemporary method (Roša et al. 2021) and with their own unpublished vector formalism method, and they examine the sensitivity of the result to the choice of triple via a 20-combination table.","tokens_in":22994,"tokens_out":8812,"duration_ms":89748,"significance":"If the reconstruction is correct, the paper closes a historical gap by turning Bošković's verbal description into a modern algorithm, and it appears to be the first complete equation set for that method. The derivation is self-contained and the worked example is internally consistent, with no free parameters, and it reproduces the arcminute-level published values—a concrete success that supports the historical claim. The main value is historical-methodological rather than a new solar-physics measurement. The paper's stronger claims (arcsecond-level precision, “without any approximation,” a single unambiguous procedure, and independent VFM confirmation) are not supported by the data or by the derivation as written, and these claims should be revised. With those revisions the paper would be a sound contribution to the history of solar rotation measurements.","major_comments":[{"comment":"The derivation is exact only under the assumption that the three positions lie on a single heliographic parallel, i.e. P′C = P′C′ = P′C′′, combined with rigid rotation. This equality is asserted from the mirror-symmetry of triangles △9 and △6 (§2.4.1) and is not derived from the measurements. The conclusion's statement that the spherical solution consists of “closed equations without any approximation” is therefore unsupported: the trigonometric algebra is exact, but the physical model contains the constant-latitude, rigid-rotation approximation. The sensitivity of this assumption is visible in Table 6, where the VFM calculation over the 20 triples gives i between 3.512° and 18.696° and Ω between −41.197° and 87.811°, while Bošković's selected triple [136] gives i=6.807° and Ω=74.048°. The selection criteria in §4.2 (C1≈C6, C3 at minimum latitude) are precisely a search for a triple that approximately satisfies the same-parallel condition, so the method is not a general three-position inversion without additional assumptions.","section":"§2.3, §2.4.1, Eqs. (26), (34), (43), (50)–(51); §4.2, Table 6"},{"comment":"The claimed precision of 10⁻⁵ arcseconds is not supported by the input data. Table 1 lists ecliptic longitudes and latitudes to arcminutes, and the historical measurement process is far coarser; nevertheless Table 4 reports i = 6.80728° = 6°48′26.20337″ and Ω = 74.04774° = 74°02′51.87646″. A 10⁻⁵ arcsecond digit is about nine orders of magnitude finer than the arcminute input. The statement in §4.3 that the results were “calculated using high precision and closed equations without loosing precision” conflates arithmetic precision with measurement precision. The paper should propagate the input uncertainties and quote i, Ω, and T′ to a justified number of significant digits; the agreement with Bošković's values is then at the arcminute level, which is sufficient for the historical claim but not for the advertised precision.","section":"§4.3, Tables 1, 2, and 4"},{"comment":"The VFM method used as a “contemporary method” check is unpublished (“paper is in preparation”). For the triple [136] it returns values identical to the spherical-trigonometric solution to eight significant figures (i=6.80727871°, Ω=74.0477436°, T′=26.8062322 d). If VFM is mathematically equivalent to the present derivation, this agreement is by construction and provides no independent confirmation; if it is not equivalent, the reader needs the VFM definitions to evaluate the comparison. The independent validation should rest on the comparison with Bošković's published values and with Roša et al. (2021); the VFM table should be moved to supplementary material or explicitly reported as an internal consistency check.","section":"§4.1, Table 4"},{"comment":"The equations use inverse trigonometric functions (arccos, arctan) without specifying branch selection. For example, Eq. (30) determines B′′−D from its cosine, which has two possible signs; Eqs. (42) and (55) have the same issue, and Eq. (7) assumes that P lies inside the angle CC′C′′ so that the two angles sum rather than subtract. A reader applying the method to arbitrary triples—as the Conclusion invites for exoplanets and stars—cannot reproduce the results without additional quadrant rules or case distinctions. The claim of a single closed-form procedure is therefore incomplete as stated.","section":"§2.2–2.4, Eqs. (7), (30), (42), (55)"}],"minor_comments":[{"comment":"The row for Eq. (6) appears twice with different values: first cos P C′′C′ = 0.107516 with angle 83.82785°, and later cos P C′′C′ = 0.108113 with angle 83.79345°. The second value is the one used in Eq. (28); the first appears to be a typographical error and should be harmonized.","section":"Table 2"},{"comment":"The numerical values quoted for C1 and C6 (20°27′ and 22°45′) do not match the entries in Table 1 as the table is currently formatted (20°37′ and 22°45′, if read literally); please reconcile the text with the table.","section":"§4.2"},{"comment":"The quantity A = 365.25 days should be defined explicitly as the Earth's orbital (or tropical/sidereal) period used in converting between synodic and sidereal rotation periods; the present text introduces it only in a footnote-like parenthetical.","section":"Eq. (37)"},{"comment":"The scanned historical page in Figure 5 is difficult to read in the arXiv version; a transcription of the relevant numbered paragraphs (№76–№81) would make the mapping from Bošković's text to the equations easier to verify.","section":"Figure 5"}],"recommendation":"major_revision","confidential_remarks":"The historical core of the paper is solid and the derivation appears internally consistent, but the authors should be urged to recalibrate the precision language and to present the VFM comparison transparently, since an unpublished method cannot serve as independent confirmation. The paper is more history-of-science than new solar-physics measurement, which is acceptable for Solar Physics if the editors see fit; no concerns about citation practice beyond the dependence on the in-preparation VFM paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it promises: it gives the first explicit spherical-trigonometric equations for Bošković's 1785 method, applied to his 1777 sunspot positions. The math is classical and the worked example is transparent; I checked the key numbers and they reproduce Bošković's published rounded values. That is a real contribution to a narrow historiographic puzzle.\n\nThe soft spots are mostly about framing. The conclusion's claim that the solution is \"closed equations without any approximation\" is too strong. The derivation assumes the three sunspot positions lie on the same heliographic parallel and that the spot rotates rigidly (P'C = P'C' = P'C''). That is a model assumption, not a consequence of the spherical geometry. The paper's own Table 6 shows how much the answer depends on triple choice: across the 20 triples, i varies from 3.5° to 18.7°, and Ω from -41.2° to 87.8°. The chosen triple [136] gives the historically nice result, but other triples are far off. The paper describes Bošković's selection criteria in Section 4.2, which is honest, but it does not reconcile that sensitivity with the \"no approximation\" claim.\n\nThe 10^-5 arcsec precision claim in Section 4.3 is also unsupported. The inputs are arcminute-level (Bošković's published positions are rounded to arcminutes), so the many-digit agreement between the full and short solutions is internal arithmetic consistency, not measurement precision. The authors should either propagate errors from the inputs or explicitly say the precision refers only to the consistency of the equations.\n\nUsing the unpublished VFM method as confirmation is weak. \"Paper in preparation\" is not a citable validation. If VFM is essential, the authors need to publish it or at least provide a preprint; otherwise it should be demoted to a sanity check or removed. Table 6 is fine as exploratory, but the labels \"near Carrington\" and \"near Spörer\" are doing a lot of work.\n\nThere are also a few minor typographical issues (e.g., the footnote about №79 confusing CC′ and C′C′′), but nothing that undermines the derivation.\n\nBottom line: if the authors soften the precision and approximation claims, add a real error analysis, and either make VFM citable or drop it as validation, this is a solid historical reconstruction. It deserves a serious referee and, with revisions, publication. I would not cite it in my own work, but for anyone working on 18th-century solar astronomy or the history of spherical trigonometry, it will be the reference for Bošković's method.","headline":"A genuinely useful first derivation of Bošković's spherical-trig method, but the 'no approximation' and 10^-5 arcsec claims need to be walked back.","tokens_in":23602,"tokens_out":2701,"would_cite":false,"duration_ms":31368,"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":"This paper derives, for the first time, the closed-form spherical-trigonometric equations behind Bošković's 1785 verbal method, and reproduces his 1777 solar rotation elements from three sunspot positions.","keywords":["Ruđer Bošković","sunspot positions","solar rotation elements","spherical trigonometry","solar equator inclination","ascending node longitude","sidereal rotation period","historical astronomy"],"falsifier":"Compute high-precision heliographic positions of a single sunspot at three times, apply these equations, and compare with an independent least-squares rotation-element fit over a longer track; if the three-position values differ from the fit by more than the measurement uncertainties, the rigid-parallel assumption is the source of the bias.","tokens_in":22557,"feed_emoji":"☀️","tokens_out":9619,"duration_ms":91012,"temperature":0.7,"pith_summary":"This paper claims that the spherical-trigonometric method Ruđer Bošković described in words in 1785 for finding solar rotation elements from three sunspot positions can be expressed as a complete, closed system of equations. The authors derive those equations and apply them to Bošković's own 1777 measurements, obtaining solar equator inclination $i = 6.80728^\\circ$, ascending-node longitude $\\Omega = 74.04774^\\circ$, sidereal rotation period $T' = 26.806232$ days, and synodic period $T'' = 28.929403$ days. These values agree with his published planar solution ($i = 6^\\circ49'$, $\\Omega = 74^\\circ03'$) and his sidereal period of 26.77 days at the arcminute level. The importance, if the derivation is right, is that all three rotation elements can now be extracted from just three spot positions in one procedure, with no iterative fitting.","feed_headline":"Bošković's 1785 solar-rotation recipe becomes equations","feed_subtitle":"Closed-form spherical trigonometry reproduces his 1777 Sun: tilt 6.81°, node 74.05°, spin period 26.81 days.","key_machinery":"The load-bearing object is the configuration of eight oblique spherical triangles (plus the short-solution triangles △10 and △11) built from the northern ecliptic pole $P$, the northern solar-equator pole $P'$, and three sunspot positions $C$, $C'$, $C''$ in ecliptic coordinates. Sides such as $CC'$ and $P'C''$ are obtained by repeated cosine and cotangent rules; midpoints $E$ and $E'$ create right triangles through which $P'E'$ is found; and mirror symmetry of triangles △9 and △6 sets $P'C' = P'C''$. The final elements come from the side $PP' = i$ via the cosine rule, the angle at $P$ giving the node $\\Omega$, and the angle at $P'$ combined with elapsed mean solar time giving $T'$.","core_discovery":"The paper's central claim is that Bošković's §VII, №76–№81 description corresponds to an exact chain of spherical triangles linking the northern ecliptic pole, the northern solar-equator pole, and three observed positions of one sunspot in ecliptic coordinates. Solving the chain with cosine and cotangent rules yields the arc between the poles as the inclination $i$, the angle at the ecliptic pole as the node $\\Omega$, and the angle at the equator pole divided by the elapsed mean solar time as the sidereal period $T'$. The reported results are $i = 6.80728^\\circ = 6^\\circ48'26.20337''$, $\\Omega = 74.04774^\\circ = 74^\\circ02'51.87646''$, and $T' = 26.806232$ days, closely matching Bošković's planar solution and his sidereal period. The same equations over all twenty triples of his six positions produce a wide spread in $i$ and $\\Omega$, so the choice of well-separated positions matters, and the paper concludes that the spherical solution is complete and gives all three elements without approximation.","pith_inferences":["The derivation's rigid-parallel assumption means that differential rotation will leak into the inferred elements; the three period estimates in Table 4 (26.77, 26.81, 26.85 days) already bracket the effect, so feeding the equations modern latitude-resolved spot data would quantify the bias.","A sharper historical test would re-reduce Bošković's apparent-disk measurements with modern ephemerides before applying the new equations; any shift in $\\Omega$ and $i$ away from his 1785 values would measure how much log-table rounding contributed to the original result.","The closed-form chain could be inverted symbolically to propagate measurement errors into $i$, $\\Omega$, and $T'$ analytically, giving uncertainty estimates without Monte Carlo sampling."],"forward_implications":["Three well-separated positions of one spot are enough to return $i$, $\\Omega$, and $T'$ simultaneously from closed-form trigonometry, so no least-squares fit or initial guess is required.","Applying the new equations to Bošković's six-day 1777 sequence reproduces his published values ($i \\approx 6^\\circ49'$, $\\Omega \\approx 74^\\circ03'$, $T' \\approx 26.8$ d), strengthening the case that his verbal procedure was numerically complete.","Relative to the Carrington and Spörer elements for 1777, the spherical solution's errors are about 2–3% in $\\Omega$, 2–7% in $i$, and under 1% in $T'$, making it a serviceable quick estimator for sparse historical data.","Because the same chain works for any three tracked surface features, the paper's stated generalization to other rotating bodies such as stars and exoplanets follows directly."],"supporting_citations":[{"why":"Supplies the verbal description (§VII, №76–№81) and the 1777 observation tables that the derived equations must reproduce.","marker":"Boscovich 1785b"},{"why":"The five-book compendium in which the sunspot procedure and its context were published.","marker":"Boscovich 1785a"},{"why":"Earlier reconstruction and tabulation of Bošković's rotation elements, providing the numerical baseline this work extends.","marker":"Husak et al. 2023"},{"why":"Repeated Bošković's original logarithmic calculations, giving the values the new spherical solution is compared against.","marker":"Husak, Brajša, and Špoljarić 2021a"},{"why":"Described the open problem of Bošković's missing equations, the gap this paper fills.","marker":"Husak, Brajša, and Špoljarić 2021b"},{"why":"Independent contemporary method using the same 1777 positions; its $i$, $\\Omega$, $T'$ values serve as comparison in Table 4.","marker":"Roša et al. 2021"},{"why":"Reference elements $\\Omega = 72.647620^\\circ$ and $i = 7.25^\\circ$ used for relative-error checks.","marker":"Carrington (1863)"},{"why":"Reference inclination $i = 6.97^\\circ$ used for relative-error checks.","marker":"Spörer (1874)"}],"fun_headline_variants":["Bošković's 1785 spherical trig becomes equations","Exact equations from Bošković's 1777 sunspot method","Spherical trig reproduces Bošković's solar rotation values","Bošković's solar-rotation recipe now closed-form equations","1777 sunspot data yield Bošković's trig solution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The chain assumes the spot stays at one solar latitude and the Sun rotates rigidly between the observations, so the pole-to-spot distance is the same for the first and third positions and the triangles on either side of the midpoint are mirror images.","fun_headline_variants_meta":{"raw":{"variants":["Bošković's 1785 spherical trig becomes equations","Exact equations from Bošković's 1777 sunspot method","Spherical trig reproduces Bošković's solar rotation values","Bošković's solar-rotation recipe now closed-form equations","1777 sunspot data yield Bošković's trig solution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000295,"raw_usage":{"total_tokens":1765,"prompt_tokens":1046,"completion_tokens":719,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":627}},"tokens_in":662,"tokens_out":719,"duration_ms":6742,"temperature":1.0,"reasoning_tokens":627,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:43:40.033230+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute high-precision heliographic positions of a single sunspot at three times, apply these equations, and compare with an independent least-squares rotation-element fit over a longer track; if the three-position values differ from the fit by more than the measurement uncertainties, the rigid-parallel assumption is the source of the bias.","supporting_citations":[{"cited_title":", Braj s a , R","cited_arxiv_id":null,"evidence_quote":"Earlier reconstruction and tabulation of Bošković's rotation elements, providing the numerical baseline this work extends."}],"review_version":1}