{"id":"a919e215-fa4d-437d-a584-7ddcd843131a","arxiv_id":"2508.18568","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new 3D thermophysical model for binary asteroids reports that eclipses reduce the BYORP torque on 1996 FG3 by about 7 percent, with smaller contributions from thermal inertia and mutual infrared heating.","lead":"The authors built a 3D thermal model of two asteroids that heat and shade each other, and applied it to the binary asteroid 1996 FG3. They find that eclipses and heat storage change the radiation torque that drives the pair's orbit, shifting the BYORP coefficient by about 7 percent.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Perihelion-only evaluation may make the quoted 7% eclipse effect a non-secular upper bound, not the secular BYORP change.","rationale":"The reader identified the circular-zero-inclination mutual orbit as the weakest assumption; that is an important input uncertainty, but it is explicitly bounded (e <= 0.07) or unconstrained and affects only the eclipse geometry. The perihelion-only evaluation is a methodological simplification that affects the headline 7% even under the assumed orbit, because the BYORP coefficient is defined as a secular quantity and Equation 7's normalization by P(R) is only exact when thermal emission scales instantaneously with insolation. The paper's own caveat about perihelion sensitivity is plausible but untested, and the dynamical conclusion uses the perihelion value. No circular fitting or internal inconsistency was found; the validation against the analytical B and Bennu temperatures provides independent support. The reader's CONDITIONAL verdict remains appropriate, perhaps even more strongly so, but does not need to change.","tokens_in":28161,"tokens_out":4950,"duration_ms":46476,"concrete_test":"Run the BTM eclipse/no-eclipse pair at aphelion (1.423 AU) and at the semimajor-axis distance (1.054 AU), keeping the same shape models, spin states, and mutual-orbit geometry with the Sun in the mutual orbit plane. Compute the eclipse-induced fractional change in B at each heliocentric distance. If the fractional change varies by more than about 1 percentage point, the perihelion-only 7% is not the secular change. To settle it fully, evaluate B at several true anomalies and compute the heliocentric-orbit-averaged value, weighting each mutual-orbit evaluation by the time spent at that heliocentric distance.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline result that eclipses alter the BYORP coefficient B by approximately 7% is computed at perihelion only (Section 2.4: 'Each BTM run returns ... the BYORP coefficient B, measured at perihelion'). Equation 7 removes the heliocentric distance dependence by dividing Fy by the instantaneous solar pressure P(R) and integrating over one mutual orbit. This yields the secular B only if Fy/P(R) is independent of heliocentric distance. That independence holds for the zero-thermal-inertia analytical theory, where emission is instantaneously proportional to insolation. It is not guaranteed once thermal inertia and eclipses are included, because the thermal emission term sigma*epsilon*T^4 in Eq. 6 depends on temperatures that lag and retain diurnal history, so the eclipse-induced deficit in Fy will not generally scale as 1/R^2. The authors acknowledge the possibility ('if the magnitude of BYORP is sensitive to thermophysical properties, we expect it to be most pronounced near perihelion') but do not quantify the variation. With the system's high heliocentric eccentricity (e=0.3497; insolation drops about 75% at aphelion), the perihelion-only 7% may overstate the secular change. The dynamical estimate of about 20 m over 10,000 years rests on this perihelion value. Thus the quantitative claim is not yet established for the real system.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a new 3D thermophysical model (\"BTM\") for binary asteroid systems, coupling 1D heat conduction with 3D ray tracing to include self-shadowing, eclipses, visible-light reflection, and mutual infrared radiation between the two bodies. The model is validated against analytical bowl-crater temperature solutions, against global temperature maps of Bennu (within about 1 K for the maximum), and against the analytical BYORP coefficient of McMahon and Scheeres (to better than 0.01%). Applying the model to the binary near-Earth asteroid (175706) 1996 FG3 at perihelion, the authors find surface temperatures of roughly 100 to 475 K, find that eclipses and thermal inertia can alter secondary surface temperatures by up to 14%, and find that mutual radiation is a small effect. Using a new force-integration method, they compute the BYORP coefficient B for different thermal inertias, with and without eclipses and mutual radiation. They report that eclipses change B by about 7% (making it less negative, i.e., lower torque), that thermal inertia alone changes B by up to about 0.2%, that the combination of eclipses and thermal inertia at Gamma = 160 tiu gives a change of about 6.3%, and that mutual radiation changes B by about 0.2%. They estimate that, absent tidal effects, the eclipse-induced change would reduce the semimajor-axis contraction by about 20 m over 10,000 years.","tokens_in":28368,"tokens_out":13375,"duration_ms":128267,"significance":"If the reported corrections are robust, this is the first quantitative treatment of second-order thermal effects on the BYORP coefficient for a real binary system, with direct relevance to interpreting measured orbital drift rates (e.g., the Scheirich et al. 2015 value for 1996 FG3) and to planning binary evolution models. The paper's strengths are the multiple validation benchmarks (analytical crater solution, Bennu temperatures within 1 K, and the analytical BYORP limit to 0.01%), the physical transparency of the model, and the fact that the new effects are genuine predictions based on stated physical parameters rather than fitted values. The BYORP comparison is a self-consistency check with the same theoretical framework (McMahon and Scheeres, two coauthors) rather than an independent validation of the new eclipse and thermal-inertia terms, but the temperature validation against Bennu provides independent support for the core thermal model.","major_comments":[{"comment":"The BYORP coefficient B used in the secular evolution estimates is computed from a single model run at perihelion (Section 2.4: 'Each BTM run returns ... the BYORP coefficient B, measured at perihelion'). In the zero-thermal-inertia limit the ratio F_y/P(R) in Eq. (7) is independent of heliocentric distance because emission is instantaneous and proportional to insolation, but for Gamma > 0 the sigma*epsilon*T^4 term in Eq. (6) carries memory of prior insolation, so F_y/P(R) can depend on R. Given the system's heliocentric eccentricity e = 0.3497 (insolation at aphelion is about 25% of that at perihelion), the perihelion-only 7% value is not automatically the secular B used in Eqs. (10)-(11). The paper's own caveat in Section 2.4 ('if the magnitude of BYORP is sensitive to thermophysical properties, we expect it to be most pronounced near perihelion') confirms the issue. Please compute B at multiple heliocentric distances and average over the heliocentric orbit, or demonstrate quantitatively that B varies by less than the claimed precision.","section":"§2.4, Eq. (7); §4, Eqs. (10)-(11)"},{"comment":"The model assumes a circular mutual orbit with zero inclination, although the observational bound is e <= 0.07 and the inclination is unconstrained. Eclipse duration, depth, and phase, which drive the dominant reported 7% change, depend directly on the assumed geometry. To make the 'approximately 7%' claim about 1996 FG3 rather than about the assumed geometry, the paper should include sensitivity runs over e in [0, 0.07] and over a plausible range of mutual inclinations. In addition, Table 1 notes that pole positions are given in ecliptic coordinates, but no pole values appear in the table or text; without the primary's pole orientation relative to the heliocentric orbit, the Sun's elevation above the mutual orbit plane, and hence the eclipse pattern, is not reproducible.","section":"§2.4"}],"minor_comments":[{"comment":"The sentence 'Holding the thermal inertia steady at 0 tiu and not including eclipses, the BYORP coefficient is altered by approximately 7%' appears to be a typo; it should read 'including eclipses' (or 'with eclipses enabled'), since the same paragraph attributes the 7% change to the inclusion of eclipses.","section":"§3.3, paragraph on Figure 7"},{"comment":"The 'approximately 7%' eclipse change is computed at perihelion and for zero thermal inertia; please state these conditions in the abstract and conclusions so that readers do not take it as the secular value for the realistic thermal inertia (the paper reports about 6.3% for eclipse plus thermal inertia at Gamma = 160 tiu).","section":"Abstract and §5"},{"comment":"The word 'Assummed' should be 'Assumed'; in addition, the table's note that 'Pole positions are given in ecliptic coordinates' is not supported by any pole values in the table, so either add the pole coordinates or remove the note.","section":"Table 1"},{"comment":"Two wording errors: 'shape models derived from from Earth-based radar data' has a duplicated 'from', and 'This causes a an irregular drop' has a doubled article.","section":"§2.1.2 and §3.1"},{"comment":"The reported B values and temperature metrics lack uncertainty or numerical precision estimates; given the authors' own sphere test showing numerical B errors at the ~10^-5 level (Section 4), a statement of the numerical precision of the quoted percent changes would help the reader judge their significance.","section":"Figure 7 and Table 2"},{"comment":"The normalization constant and the meaning of the integral (why dividing by 2*pi*R_RES^2 yields the A0(2) coefficient) could be explained in one sentence or with a specific reference to the McMahon-Scheeres definition, to make the paper more self-contained.","section":"Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The perihelion-only evaluation and the unconstrained mutual-orbit geometry are the core revision points; both are acknowledged in the text but not quantified. I do not see a fundamental flaw in the thermal model itself, and the validation against Bennu and the analytical crater case is convincing. The authors should also ensure that the code, shape models, and the mutual orbit pole orientation are explicitly documented, as the current manuscript mentions a link but does not provide the pole values needed to reproduce the eclipse geometry."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper builds a ray-traced binary thermophysical model, validates it sensibly, and shows that eclipses can change the BYORP coefficient B for 1996 FG3 by a few percent. That is new. Prior work on this system and on Didymos did not compute BYORP with full binary thermal effects. The model also gives useful temperature maps and quantifies mutual radiation, which is a nice bonus.\n\nThe validation is genuinely good for a new numerical model: bowl-crater analytical temperatures within 1 K, Bennu global temperatures within 1 K of an independent model, and the zero-thermal-inertia BYORP coefficient reproduced to 0.01%. The fact that two coauthors developed the analytical theory being matched is not a real problem here, because the BTM is an independent implementation and the match is at the level of a consistency check, not a fit. The parameter sweep over thermal inertia and the spherical-secondary control runs are also honest checks.\n\nNow the soft spots, and one of them is load-bearing. B is computed at perihelion only and then inserted into secular evolution equations. Equation 7 divides Fy by the instantaneous solar pressure P(R) and integrates over one mutual orbit. That gives a secular coefficient only if Fy/P(R) is independent of heliocentric distance. That holds for zero thermal inertia, but not once you have thermal lag and eclipse history. The paper explicitly says the zeroth-order Fourier coefficient is assumed independent of position along the heliocentric orbit, and then adds that thermophysical sensitivity should be most pronounced near perihelion. That is an admission that the assumption is suspect, and it is not quantified. With the system's eccentricity of 0.35, insolation at aphelion is about 25% of perihelion. The perihelion-only 7% could easily be an upper bound, not the secular value. The 20-meter drift over 10,000 years rests directly on that unquantified number.\n\nThe other soft spots are minor by comparison: the mutual orbit is assumed circular and zero-inclination because those are unconstrained, and the eclipse geometry drives the main result. No propagated uncertainties are given for B. No code or data are released, which hurts reproducibility but is not unusual at this stage.\n\nWho is this for? People modeling binary asteroid thermal environments or BYORP evolution, and mission teams needing temperature predictions for close binaries. The qualitative conclusions—eclipses reduce torque, thermal inertia has a small non-monotonic effect, mutual radiation is third-order—are credible and worth taking seriously. But the quantitative 7% claim needs either a multi-heliocentric-distance calculation or a clear justification for why perihelion dominates.\n\nRecommendation: send it to peer review. It deserves referee time. The authors should be pushed to test the heliocentric-orbit dependence and report uncertainties before the 7% number is quoted as the secular effect.","headline":"A credible first step toward including eclipses and thermal inertia in BYORP calculations, but the headline 7% number is computed at perihelion only and is not yet established as the secular effect.","tokens_in":29017,"tokens_out":1228,"would_cite":false,"duration_ms":14079,"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":"For the binary asteroid 1996 FG3, eclipses alter the BYORP coefficient B by about 7 percent, lowering the radiation torque on the secondary and changing its predicted orbital contraction.","keywords":["binary asteroids","BYORP effect","thermophysical model","thermal inertia","eclipses","mutual radiation","1996 FG3","YORP effect"],"falsifier":"Measure the semimajor-axis drift of 1996 FG3's mutual orbit precisely enough to separate BYORP from tidal forces, or observe the eclipse-induced temperature dips on the secondary with spacecraft infrared data; if the drift or temperature pattern matches the standard no-eclipse theory rather than the BTM prediction, the claimed 7 percent reduction in B is wrong. A shorter test: recompute B with the mutual orbit's true eccentricity and inclination once they are measured; if the approximately 7 percent eclipse signature disappears, the result is an artifact of the assumed geometry.","tokens_in":27891,"feed_emoji":"☄️","tokens_out":4504,"duration_ms":36255,"temperature":0.7,"pith_summary":"This paper develops a 3D thermophysical model for binary asteroid systems that combines heat conduction, insolation, visible-light reflection, ray-traced shadows and eclipses, and mutual infrared radiation between the two bodies, and applies it to the binary (175706) 1996 FG3. The authors aim to show that thermal effects that standard BYORP theory ignores, chiefly eclipses but also thermal inertia and mutual radiation, change the BYORP coefficient B in measurable ways. For 1996 FG3 they find eclipses alter B by about 7 percent, lowering the radiation torque on the secondary, while thermal inertia adds a smaller correction and mutual radiation a still smaller one. The point of the claim is that binary evolution models should include these thermal corrections, and that the coming wave of spacecraft infrared data on binaries will need such a model to be interpreted.","feed_headline":"Eclipses cut binary asteroid torque by 7 percent","feed_subtitle":"New model shows thermal effects slow 1996 FG3's orbit shrinkage by ~20 meters in 10,000 years.","key_machinery":"The central machinery is the Binary Thermophysical Model: a 1D heat-conduction solver coupled to a 3D ray tracer and view-factor calculator that tracks insolation, self-shadowing, eclipses, and mutual visible and infrared scattering between the two asteroid meshes. The load-bearing quantity is the BYORP coefficient B, defined as the normalized zeroth-order Fourier coefficient of the solar-radiation-pressure force parallel to the secondary’s orbital motion, computed by integrating the along-track force over one mutual orbit and dividing by the solar pressure and the secondary’s volume-equivalent radius squared. Eclipses enter through ray-traced shadows, thermal inertia through the subsurface conduction equation, and mutual radiation through single-scatter view factors using the M2 method. The model is validated against analytical crater temperatures and against Bennu thermophysical results before being applied to 1996 FG3 under the assumption of a circular, zero-inclination mutual orbit.","core_discovery":"On its own terms, the paper claims that the Binary Thermophysical Model (BTM) reproduces the analytical BYORP coefficient B to better than 0.01 percent when thermal effects are switched off, and that turning them on changes B in a systematic way. For 1996 FG3 at perihelion, including eclipses makes B about 7 percent less negative, meaning less torque on the secondary, by cooling facets on the leading and trailing sides whose thermal emission contributes to the along-track force. Thermal inertia produces a small non-monotonic correction: it first increases torque by introducing a phase lag, then decreases it by damping diurnal temperature amplitudes, with a tipping point between 160 and 320 tiu. Mutual radiation from the primary changes B by roughly 0.2 percent, an order of magnitude smaller. Over 10,000 years and in the absence of tides, the eclipse-induced reduction in B would shrink the semimajor axis by about 20 meters less than the standard theory predicts.","pith_inferences":["For binaries with larger secondaries or tighter separations than 1996 FG3, the eclipse-induced shift in B could be substantially larger than 7 percent, making thermal corrections mandatory rather than optional.","The non-monotonic dependence of B on thermal inertia suggests a testable prediction: a binary with very low thermal inertia should show a smaller eclipse-driven torque reduction than one with moderate inertia, a difference that spacecraft infrared light curves could detect.","The same ray-tracing machinery could be applied to single rubble-pile asteroids with strong topography, where self-shadowing plays the role of eclipses, to check whether YORP coefficients likewise shift with thermal inertia.","Because the secondary’s synchronously locked face is eclipsed every orbit, the model implies a persistent day-night asymmetry on that hemisphere, which may affect regolith mobility and observable photometric behavior."],"forward_implications":["BYORP coefficients computed from shape models alone will overestimate the radiation torque for eclipsing binaries by up to several percent.","Binary evolution models that ignore thermal inertia miss a small but nonzero torque that changes sign with thermal inertia, so the effect cannot be absorbed into a single B value.","Mutual radiation can safely be neglected for most binaries, but not for close, low-thermal-inertia systems where it patches a few percent of temperature and a fraction of a percent of B.","Spacecraft infrared observations of binaries can be interpreted with the BTM to constrain thermal inertia and regolith properties, not just shape.","The eclipse-driven reduction of B weakens the BYORP contraction of the mutual orbit, which shifts the inferred balance of tides and BYORP in tidal-BYORP equilibrium systems like 1996 FG3."],"supporting_citations":[{"why":"Supplies the BYORP coefficient B formalism and the analytical value the model is compared against.","marker":"McMahon and Scheeres (2010b)"},{"why":"Provides the YORP secular-averaging methodology that the BYORP calculation builds on.","marker":"Scheeres (2007)"},{"why":"Provides the 1996 FG3 mutual orbit parameters, shape models, and the observed semimajor-axis drift rate.","marker":"Scheirich et al. (2015)"},{"why":"The validated 1D heat-conduction solver that the BTM extends to 3D binary geometry.","marker":"Hayne et al. (2017)"},{"why":"The single-body thermophysical model used as the Bennu validation benchmark.","marker":"Rozitis and Green (2011)"},{"why":"Supplies the M2 view-factor method used for ray-traced mutual radiation.","marker":"Rezac and Zhao (2020)"},{"why":"Provides the tidal-BYORP equilibrium framework and the semimajor-axis/eccentricity evolution equations used to estimate the 20-meter effect.","marker":"Jacobson and Scheeres (2011)"},{"why":"Prior work showing eclipses can dominate binary Yarkovsky effects, motivating the eclipse treatment here.","marker":"Zhou et al. (2024)"}],"fun_headline_variants":["Eclipses reduce binary asteroid torque by 7%","Binary asteroid 1996 FG3: eclipses cut torque, shrink orbit 20m less","Thermal model: eclipses alter asteroid binary orbit by 20m over 10k yr","Eclipses slow 1996 FG3's orbit shrinkage by 20m via 7% torque cut","New thermophysical model: eclipses lower binary asteroid torque"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes a circular mutual orbit with zero inclination, and the eclipse geometry that drives the 7 percent shift in B depends directly on that assumption, while the real eccentricity is only bounded at e≤0.07 and the inclination is unconstrained.","fun_headline_variants_meta":{"raw":{"variants":["Eclipses reduce binary asteroid torque by 7%","Binary asteroid 1996 FG3: eclipses cut torque, shrink orbit 20m less","Thermal model: eclipses alter asteroid binary orbit by 20m over 10k yr","Eclipses slow 1996 FG3's orbit shrinkage by 20m via 7% torque cut","New thermophysical model: eclipses lower binary asteroid torque"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001538,"raw_usage":{"total_tokens":6243,"prompt_tokens":1123,"completion_tokens":5120,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":5012}},"tokens_in":739,"tokens_out":5120,"duration_ms":35156,"temperature":1.0,"reasoning_tokens":5012,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:55:35.632065+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the semimajor-axis drift of 1996 FG3's mutual orbit precisely enough to separate BYORP from tidal forces, or observe the eclipse-induced temperature dips on the secondary with spacecraft infrared data; if the drift or temperature pattern matches the standard no-eclipse theory rather than the BTM prediction, the claimed 7 percent reduction in B is wrong. A shorter test: recompute B with the mutual orbit's true eccentricity and inclination once they are measured; if the approximately 7 percent eclipse signature disappears, the result is an artifact of the assumed geometry.","supporting_citations":[],"review_version":2}