REVIEW 2 major objections 6 minor 9 references
A 3D thermophysical model for binary asteroid systems: Application to the BYORP effect on (175706) 1996 FG3
T0 review · 2 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (2)
- [§2.4, Eq. (7); §4, Eqs. (10)-(11)] 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.
- [§2.4] 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.
minor comments (6)
- [§3.3, paragraph on Figure 7] 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.
- [Abstract and §5] 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).
- [Table 1] 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.
- [§2.1.2 and §3.1] 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.
- [Figure 7 and Table 2] 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.
- [Eq. (7)] 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.
Circularity Check
No circularity: the BTM's thermal corrections are genuine model outputs, not fitted parameters or imports from the authors' own theory.
full rationale
The paper's central derivation is the Binary Thermophysical Model: it solves the 1D heat diffusion equation (Eq. 1) with a surface energy balance, ray-traced insolation, view factors, eclipses, and mutual radiation, and then computes the BYORP coefficient B from the simulated facet forces via Eq. 7. No parameter is fitted to the target quantity B: thermal inertias are varied over a range chosen from the literature, and albedo, emissivity, shape, and orbital parameters are taken from published observations or stated assumptions. The validation step reproduces the analytical BYORP theory of McMahon and Scheeres (2010b) in the limiting case of zero thermal inertia, no eclipses, and no mutual radiation, matching Bo to better than 0.01%. Although two coauthors are authors of that theory, the theory is an independently derived analytical benchmark whose assumptions explicitly exclude the effects under study; using it as a numerical consistency check does not make the new eclipse, thermal-inertia, or mutual-radiation corrections equivalent to the theory's inputs. The headline 7% eclipse effect is obtained by toggling eclipses on and off in otherwise identical model runs, not by fitting or by renaming a prior result. The perihelion-only evaluation is a scientific caveat about secular averaging, not a circularity. The paper is therefore self-contained against its benchmark and the new claims have independent content, so the circularity score is 0.
Assumptions & free parameters
free parameters (4)
- Thermal inertia Gamma =
0, 40, 80, 160, 320, 700, 2500 tiu
- Bond albedo A_B =
0.011
- Emissivity epsilon =
0.95
- Mutual orbit eccentricity and inclination =
e = 0, i = 0
assumptions (6)
- domain assumption 1D heat conduction is sufficient for modeled facets
- domain assumption Lambertian diffuse reflection and emission with tangential forces canceling
- domain assumption Single scattering suffices for view factors
- domain assumption Absorption of radiation has no secular effect on BYORP
- domain assumption Zeroth-order BYORP coefficient is independent of heliocentric distance
- ad hoc to paper Circular, zero-inclination mutual orbit
Cite this review
Pith. "Pith review of A 3D thermophysical model for binary asteroid systems: Application to the BYORP effect on (175706) 1996 FG3." pith.science (2026). https://pith.science/paper/KJ7MEQ4V
@misc{pith2026250818568,
author = {Pith},
title = {Pith review of: A 3D thermophysical model for binary asteroid systems: Application to the BYORP effect on (175706) 1996 FG3},
year = {2026},
howpublished = {\url{https://pith.science/paper/KJ7MEQ4V}},
note = {Machine review of arXiv:2508.18568}
}
read the original abstract
Differential heating and radiation on asymmetric asteroids can cause measurable changes in their rotation rates and spin axes, known as the YORP effect. In binary systems, such radiation-driven torques can change the mutual asteroid orbits, termed the binary YORP or BYORP effect. To study how binary asteroid shapes and thermophysical properties affect surface temperatures and BYORP, we developed a new 3D thermophysical model which balances insolation, 1D conduction, visible light reflection, and mutual heating through scattered infrared radiation. Using 3D ray tracing, we include eclipses, shadowing from horizons and topography, and mutual radiation exchange between the primary and secondary asteroids. We perform global modeling of the binary asteroid (175706) 1996 FG3, a Janus mission target. At perihelion, we find that the 1996 FG3 system experiences temperatures between 100 and 475 K. We find that eclipses and thermal inertia can alter secondary surface temperatures by up to 14%, with a mean difference due to radiation from the primary of just over 1%. We also present a model for calculating the BYORP effect using binary thermophysical model results. This model compares well to analytical approximations of the BYORP coefficient B, and suggests that thermal effects like eclipses and thermal inertia can reduce torque in the 1996 FG3 system and alter the BYORP coefficient by up to several percent. For 1996 FG3, eclipses alter B by approximately 7%, resulting in a lower torque on the secondary. Though small, in the absence of tidal effects this would reduce the contraction of the semimajor axis by about 20 meters over 10,000 years. Our findings suggest that thermal effects can alter temperatures and BYORP calculations sufficiently that they should be included when modeling binaries. The relative importance of each effect is predicted to vary with the properties of the studied system.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1950]
https://doi.org/10.1086/145272. J.-P. Williams, B. T. Greenhagen, D. A. Paige, N. Schorghofer, E. Sefton-Nash, P. O. Hayne, P. G. Lucey, M. A. Siegler, and K. Michael Aye. Seasonal Polar Temperatures on the Moon. Journal of Geophysical Research: Planets, 124(10):2505–2521, 10 2019. ISSN 2169-9100. https://doi.org/10.1029/2019JE006028. URL https://onlineli...
-
[1992]
Surface Properties of Asteroids from Mid-Infrared Observations and Thermophysical Modeling
ISSN 0004-637X. https://doi.org/10.1086/170906. Harold F. Levison, Catherine B. Olkin, Keith S. Noll, Simone Marchi, James F. Bell, Edward Bierhaus, Richard Binzel, William Bottke, Dan Britt, Michael Brown, Marc Buie, Phil Christensen, Joshua Emery, Will Grundy, Victoria E. Hamilton, Carly Howett, Stefano Mottola, Martin P¨ atzold, Dennis Reuter, John Spe...
-
[2002]
https://doi.org/10.1006/ICAR.2002.6856
ISSN 0019-1035. https://doi.org/10.1006/ICAR.2002.6856. Andrew F. Cheng, Andrew S. Rivkin, Patrick Michel, Justin Atchison, Olivier Barnouin, Lance Benner, Nancy L. Chabot, Carolyn Ernst, Eugene G. Fahnestock, Michael Kueppers, Petr Pravec, Emma Rainey, Derek C. Richardson, Angela M. Stickle, and Cristina Thomas. AIDA DART asteroid deflection test: Planet...
arXiv 2002
-
[2008]
https://doi.org/10.1038/nature07078
ISSN 1476-4687. https://doi.org/10.1038/nature07078. URL https://www.nature.com/articles/nature07078. 26 Fred L. Whipple. A comet model. I. The acceleration of Comet Encke. Astrophysical Journal, 111:375 – 394,
-
[2009]
https://doi.org/10.1016/J.ICARUS.2009.03.003
ISSN 0019-1035. https://doi.org/10.1016/J.ICARUS.2009.03.003. Cristina A. Thomas, Shantanu P. Naidu, Peter Scheirich, Nicholas A. Moskovitz, Petr Pravec, Steven R. Chesley, Andrew S. Rivkin, David J. Osip, Tim A. Lister, Lance A.M. Benner, Marina Brozovi´ c, Carlos Contreras, Nidia Morrell, Agata Ro˙ zek, Peter Kuˇ snir´ ak, Kamil Hornoch, Declan Mages, P...
-
[2021]
https://doi.org/10.1016/J.PSS.2021.105268
ISSN 0032-0633. https://doi.org/10.1016/J.PSS.2021.105268. B. Rozitis and S. F. Green. Directional characteristics of thermal–infrared beaming from atmosphereless planetary surfaces – a new thermophysical model. Monthly Notices of the Royal Astronomical Society , 415(3):2042–2062, 8
-
[2022]
https://doi.org/10.1016/j.icarus.2021.114826
ISSN 10902643. https://doi.org/10.1016/j.icarus.2021.114826. L Rezac and Y Zhao. Astronomy Astrophysics Accuracy of view factor calculations for digital terrain models of comets and asteroids. A&A, 642:167, 2020. https://doi.org/10.1051/0004-6361/202038462. URL https://doi.org/10.1051/0004-6361/202038462. J. H. Roberts, O. S. Barnouin, M. G. Daly, K. J. W...
arXiv 2021
-
[2023]
https://doi.org/10.1038/s41586-023-05878-z
ISSN 1476-4687. https://doi.org/10.1038/s41586-023-05878-z. URL https://www.nature.com/articles/s41586-023-05878-z. Rachel H Cueva, Jay W. McMahon, Alex J Meyer, Daniel J Scheeres, Masatoshi Hirabayashi, Sabina D Raducan, Seth A Jacobson, and Colby C Merrill. The Secular Dynamical Evolution of Binary Asteroid System (65803) Didymos Post-DART. Planetary Sc...
Show all 9 references
-
[2024]
standard
ISSN 26323338. https://doi.org/10.3847/PSJ/ad2173. URL https://doi.org/10.3847/PSJ/ad2173. Matija ´Cuk and Joseph A. Burns. Effects of thermal radiation on the dynamics of binary NEAs. Icarus, 176 (2):418–431, 8 2005. ISSN 0019-1035. https://doi.org/10.1016/J.ICARUS.2005.02.00...
2005
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.