{"id":"4eaa97b0-ea47-435c-9365-6f9d340aa097","arxiv_id":"2507.03583","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Binary companions imprint characteristic shapes on low-J CO lines from AGB winds, but beam smearing and noise can hide them, causing single-star misclassification.","lead":"This paper uses 3D simulations of AGB stars with binary companions to predict how the companion shapes low-energy CO emission lines, finding distinctive double-peaked or bumpy line shapes that can reveal the companion. It also shows these features are often hidden by telescope beam and noise, which can bias derived mass-loss rates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Free-wind launching makes the velocity-bin geometry that produces the claimed double peaks/central bumps the least secure part of the central claim; a realistic accelerating wind could alter or erase these features.","rationale":"After reading the paper in good faith, the central claim is exactly as the reader summarizes: binary companions produce recognizable low-J CO line shapes (double peaks, central bump) that can serve as a diagnostic, and the features can be hidden by beam and noise. The supporting calculation is careful: the Phantom/Magritte pipeline is standard, the SKIRT cross-check in Appendix A validates the radiative transfer, and the photodissociation treatment is a genuine improvement over spherical prescriptions. The weakest link is not in the numerics but in the wind-launching physics. The hydrodynamics deliberately replaces radiative acceleration with a free-wind injection at R_inj (Sect. 2.1), and Sect. 4.5 states that the resulting velocity profile is unrealistic. The paper's own explanation of the line features (Sect. 3.2.1, Fig. 4) attributes them to the topology of velocity bins, i.e. to the velocity field as modified by the companion. If the true AGB wind accelerates gradually through the 9-25 au orbital region, the companion-induced velocity perturbations, and hence the velocity-bin volumes, will differ from those simulated. This does not invalidate the qualitative statement that binary companions perturb line profiles, but it puts the specific diagnostic features on uncertain ground. The reader's CONDITIONAL verdict is appropriate; I agree with the identified weakest assumption. No additional load-bearing error was found, and the proposed test would settle whether the features are artifacts.","tokens_in":31260,"tokens_out":6747,"duration_ms":91273,"concrete_test":"Rerun one representative model (v10a09) in Phantom with a more physical wind-launching prescription, specifically the radiation-pressure and dust-formation scheme of Esseldeurs et al. (2023), keeping the mass-loss rate, orbital parameters, and Magritte settings identical. Compare the CO J=2-1 face-on and edge-on profiles with the current free-wind result. If the double peaks near terminal velocity and the central bump survive, the free-wind concern is retired; if these features shift, weaken, or disappear, the binary-diagnostic claim must be restricted to winds that are already near terminal velocity at the dust condensation radius.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that companion-induced morphologies imprint characteristic low-J CO line profiles, specifically double peaks near the terminal velocity and a central bump, and that these profiles can serve as a binary diagnostic. The mechanism is traced in Sect. 3.2.1 to the non-conical shapes of iso-velocity contours caused by the companion's gravity (Fig. 4). But the velocity field itself is the least physical ingredient of the model: particles are injected at the dust condensation radius R_inj = 4.26 au with a fixed initial velocity and no radiative acceleration (Sect. 2.1), and Sect. 4.5 explicitly acknowledges that the velocity profile in the simulations is unrealistic. In a real dust-driven wind, the outflow accelerates from subsonic speeds near the photosphere to terminal velocity over tens of stellar radii; at orbital separations of 9-25 au the companion therefore sits inside the acceleration zone, where its gravitational influence on the flow, and hence the resulting velocity-bin volumes, depends strongly on the actual v(r). Maes et al. (2021) is cited by the authors for exactly this sensitivity. Since the double-peak and central-bump features are explained as enhancements of specific velocity bins whose volumes are governed by this artificial velocity field, the core diagnostic claim rests on the least secure component of the model. This is not an internal inconsistency but a correctness risk: the qualitative conclusion may survive, but the specific predicted line shapes could be artifacts of the free-wind treatment.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a grid of nine 3D SPH models of an AGB star with a 1 M_sun companion at orbital separations of 9, 15, and 25 au, using free-wind injection velocities of 5, 10, and 20 km/s. The models are post-processed with the NLTE line radiative transfer code Magritte to compute synthetic low-J CO lines at multiple inclinations and position angles, including a 3D treatment of CO photodissociation. The authors report non-spherical morphologies with spirals and global flattening, and show that the synthetic line profiles can deviate strongly from the parabolic or flat-topped shapes expected from spherical winds, with the most pronounced features being two peaks near the terminal velocity in face-on views and a central bump near the central velocity in edge-on views. They also convolve the lines with Gaussian beams, add noise, and use the De Beck et al. (2010) fitting formula to estimate retrieved mass-loss rates. The paper concludes that molecular line profiles can serve as a binary diagnostic, but that the characteristic features can be hidden by beam and noise effects, potentially leading to misclassification and systematic errors in derived mass-loss rates.","tokens_in":31538,"tokens_out":7831,"duration_ms":92892,"significance":"If the qualitative conclusions hold, this is a useful step toward connecting binary-induced outflow morphologies with the most commonly observed molecular line observations of AGB stars. The paper's strengths include the use of open-source, reproducible codes (Phantom, Magritte), a validation of the radiative transfer against SKIRT in Appendix A (roughly 1% agreement for J=3-2 and 5-4, 5-10% for J=1-0), a validation of the 3D CO photodissociation interpolation in Appendix B, and an explicit treatment of beam convolution and noise. The authors are also commendably transparent about the limitations of their wind-driving prescription in Sect. 4.5. However, the central diagnostic signatures are tightly linked to a simplified free-wind velocity field, the resolution study is only qualitative, and the mass-loss-bias argument is not yet controlled. The manuscript therefore establishes a promising framework and a set of well-documented model predictions, but the specific claim that low-J CO line profiles can serve as a binary diagnostic still needs additional quantitative support.","major_comments":[{"comment":"","section":"§2.1, §3.2.1, §4.5"},{"comment":"","section":"§2.2, Table 1"},{"comment":"","section":"§4.1, Table 3"},{"comment":"","section":"§3.2.5, Fig. 13"}],"minor_comments":[{"comment":"","section":"§3.2.5"},{"comment":"","section":"§4.2"},{"comment":"","section":"§3.2.1"},{"comment":"","section":"§3.2.5"}],"recommendation":"major_revision","confidential_remarks":"This is a technically competent numerical exploration from a group that has extensive experience with these codes, and the appendices provide useful validation. My main concern is not the execution but the reach of the conclusions relative to the simplified wind physics: the specific double-peak and central-bump signatures are tied to a free-wind velocity field that the authors themselves acknowledge is unrealistic, and the resolution and mass-loss-retrieval comparisons need quantitative support. I would not recommend rejection because the framework is valuable and the limitations are explicit, but I would ask for either new sensitivity calculations or a more guarded statement of the binary-diagnostic claim. The orientation discrepancy noted in Sect. 4.4, where the central bump is observed in a face-on source (EP Aqr) but appears edge-on in the models, also suggests that the mapping between simulated features and real sources is not yet robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is a solid, honestly written paper that for the first time pushes 3D SPH binary AGB wind models through full NLTE line radiative transfer, including a genuinely new 3D CO photodissociation treatment. The SKIRT cross-check for one model is real evidence; the line profiles agree to roughly 1% for J=3-2 and 5-4 and to 5-10% for J=1-0. The qualitative claim that a companion can imprint double peaks near terminal velocity or a central bump is supported by the simulations.\n\nThe central caveat is the one the authors state in Sect. 4.5: the wind is launched as a free wind at the dust condensation radius, with gravity artificially balanced and no radiative acceleration. The velocity field is therefore not the physical accelerating wind, and the companion at 9-25 au sits right in the zone where a real wind would still be accelerating. The double peaks and central bump are explicitly explained as enhancements of specific velocity bins whose volumes are set by this artificial velocity field. That means the specific line shapes—especially at low v_inf and small separations—could change under a realistic v(r). The paper cites Maes et al. (2021) for exactly this sensitivity, so the authors are not hiding it; they are just not in a position to resolve it yet. I would call this a conditional result rather than a fatal flaw: the qualitative binary diagnostic is plausible and may survive, but the catalog of concrete line shapes should be treated as demonstration-in-model rather than prediction-for-stars.\n\nOther soft spots, in order: no quantitative resolution convergence test (they mention a high-res run but show no numbers); no public model data; only one mass-loss rate and one chemical composition; and the mass-loss retrieval bias is shown via an external spherical fitting function, which they use only for relative trends—that is fine. The line shape degeneracies with beam and noise are handled thoughtfully, including the warning that double peaks can arise from resolved spherical outflows.\n\nWho this is for: anyone modeling low-J CO lines or identifying binary AGB systems from spectra. It deserves a serious referee; the request for revision should focus on the convergence test and a quantitative statement about how much the free-wind assumption can shift the features—for example, a simple comparison with a beta-law velocity field would go a long way. I would take it in the reading group and would cite the photodissociation method.","headline":"Competent, honest 3D SPH+NLTE study of binary AGB line profiles; the qualitative binary diagnostics are plausible, but the specific line shapes rest on the free-wind velocity field and should be treated as conditional.","tokens_in":32140,"tokens_out":3037,"would_cite":true,"duration_ms":36051,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Low-J CO lines from an AGB wind carry a binary companion's imprint — double peaks face-on, a central bump edge-on — that can act as a binary diagnostic, though beams and noise can hide them.","keywords":["AGB stars","binary companions","CO spectral lines","circumstellar outflows","synthetic line profiles","NLTE radiative transfer","smoothed particle hydrodynamics","CO photodissociation"],"falsifier":"Observe a confirmed binary AGB system whose orbital parameters are known from astrometry or resolved imaging, with a slow wind (≈5 km/s) and a close orbit (≈9 au), and measure its low-J CO lines at high signal-to-noise with beam effects removed: if the lines do not show the predicted double-peaked (face-on) or central-bump (edge-on) structure across several transitions, the velocity-bin mechanism is falsified. A positive control would be a known face-on narrow-spiral system, which the models say should show the two-peaked profile rather than a smooth parabola.","tokens_in":31039,"feed_emoji":"⭐","tokens_out":14400,"duration_ms":138280,"temperature":0.7,"pith_summary":"This paper claims that a binary companion orbiting an AGB star leaves characteristic, identifiable marks on the low-J CO spectral lines of the star's wind: instead of the smooth parabolic or flat-topped profiles expected from a spherical outflow, the lines show two peaks near the terminal velocity when viewed face-on and a pronounced bump near the central velocity when viewed edge-on. The claim is built on a grid of nine 3D hydrodynamic models, spanning outflow velocities of 5, 10, and 20 km/s and orbital separations of 9, 15, and 25 au, post-processed with full non-local thermodynamic equilibrium radiative transfer. If the claim is right, molecular line profiles become a binary diagnostic for the many AGB outflows that are observed only through spectra, not resolved images. The paper adds a caution: the same features can be concealed by the telescope beam and by noise, so binaries can be misclassified as single stars, and standard single-star modelling of their lines yields systematically wrong mass-loss rates, around a factor of two in the cases studied.","feed_headline":"Double-peaked CO lines can expose hidden AGB companions","feed_subtitle":"Simulations show companion-shaped winds twist CO profiles, but beam and noise can still hide the clue.","key_machinery":"The load-bearing mechanism is the geometry of velocity bins: in a spherical outflow the projected-velocity bins are conical and roughly equal in volume, producing flat-topped or parabolic lines, whereas the companion's gravity and the star's orbital motion distort the velocity field so that certain bins become arc-shaped, enlarged, or double-valued in projected velocity. Those distortions concentrate or deplete emission at specific velocities, creating the two peaks near the terminal velocity and the central bump. The machinery has three parts: the smoothed particle hydrodynamics code Phantom, which evolves a free wind (launched at the dust-condensation radius with fixed speed) in the binary potential with a polytropic equation of state and H I cooling; a new ray-tracing CO photodissociation calculation that generalises the spherical method of Groenewegen (2017) to 3D using HEALPix-distributed rays and nearest-ray interpolation, setting the size and shape of the emitting envelope; and the 3D NLTE line radiative transfer code Magritte, which solves for the level populations of the first 40 CO rotational transitions and produces the synthetic lines at chosen inclinations and position angles.","core_discovery":"The central discovery is that the companion's gravitational interaction reshapes both the density and the velocity field of the outflow, and the reshaped velocity field translates directly into specific low-J CO line morphologies. All nine models develop a pronounced spiral in the orbital plane, with the meridional structure, the degree of global flattening, and the emission shapes varying with orbital separation and wind speed; a novel 3D photodissociation treatment, generalising the spherical scheme of Groenewegen (2017), sets a non-spherical emitting-envelope size that closely follows the morphology. The synthetic lines deviate strongly from parabolic or flat-topped profiles, with the two most pronounced features being two peaks near the terminal velocity (viewed face-on) and a bump near the central velocity (viewed edge-on), whose presence, strength, and position depend on inclination, position angle, and the CO transition. In specific cases the profiles appear quasi-parabolic, hiding the companion entirely. The authors' conclusion is that molecular line profiles can serve as a binary diagnostic, but the companion's influence can go unnoticed because the features are easily concealed by the beam profile and observational noise, making it easy to misclassify such systems as single stars and to derive systematically biased mass-loss rates.","pith_inferences":["Re-examining existing single-dish CO spectra of 'single' AGB stars for the specific double-peak and central-bump morphologies could reveal a population of hidden companions; the grid predicts the strongest signatures for slow winds (5 km/s) and close orbits (9-15 au).","The same velocity-bin argument implies that other velocity-field distortions — pulsations, convection cells, eccentric orbits — could produce line shapes resembling the binary signatures, so the diagnostic is most secure when combined with resolved imaging or multi-transition data.","The 3D photodissociation scheme could be extended to other molecules (e.g., 13CO, HCN) or to higher mass-loss rates, where non-spherical UV shielding should be even more pronounced and the line-shape effects correspondingly stronger.","If the diagnostic holds, comparing the inferred incidence of hidden companions in CO line surveys with the known binarity fraction of AGB progenitors would offer an observational test of whether companions are indeed the dominant shaping mechanism of aspherical mass loss."],"forward_implications":["Observed low-J CO lines that deviate from parabolic or flat-topped profiles — two peaks near the terminal velocity or a bump near the central velocity — can be read as indirect evidence of a hidden companion even when the outflow is not spatially resolved.","Mass-loss rates derived by fitting single-star spherical models to such lines can be off by around a factor of two, with the error depending on inclination, so surveys that ignore binarity carry a systematic bias.","Observing several CO transitions together helps lift the degeneracy between orbital separation and wind velocity, because the line shape and strength evolve differently with J for different binary parameters.","In wide-separation, fast-wind models the higher-J lines probe the spiral structure directly, and the line shape's dependence on position angle becomes a potential way to locate the companion.","The terminal velocity inferred from the lines depends on orbital parameters and inclination rather than simply the wind speed, so companion shaping can masquerade as a different wind speed in single-star fits."],"supporting_citations":[{"why":"Introduced the free-wind approximation and the prediction that a binary companion shapes the AGB outflow, the conceptual foundation of the models.","marker":"Theuns & Jorissen (1993)"},{"why":"Presented the SPH code Phantom, which runs the nine hydrodynamic simulations.","marker":"Price et al. (2018)"},{"why":"Developed the 3D NLTE radiative transfer code Magritte, including the recursive remeshing used to produce the synthetic CO lines.","marker":"De Ceuster et al. (2019, 2020, 2022), Ceulemans et al. (2024)"},{"why":"Supplied the spherical CO photodissociation prescription that the paper generalises to 3D geometry.","marker":"Groenewegen (2017)"},{"why":"Provided the unshielded photodissociation rate and the CO and H2 shielding functions used in the 3D photodissociation calculation.","marker":"Visser et al. (2009)"},{"why":"Defined the backward and frontal spiral edges and described the binary-wind interaction producing the spiral and arc morphologies seen here.","marker":"Malfait et al. (2021)"},{"why":"Supplied the analytical mass-loss-rate fitting formula used to quantify the systematic errors introduced by ignoring the companion.","marker":"De Beck et al. (2010)"},{"why":"Earlier hydrodynamic binary models whose velocity-bin column-density proxy the paper supersedes with full NLTE radiative transfer.","marker":"Kim et al. (2019)"},{"why":"Earlier analytical spiral models whose predicted line shapes the paper compares with its own double-peak and central-bump features.","marker":"Homan et al. (2015)"}],"fun_headline_variants":["Binary companions leave telltale double peaks in CO lines","AGB binaries distort CO lines, but noise can mask the clue","CO line shapes betray hidden binary companions in AGB winds","Double-peaked CO emission signals a binary, unless noise hides it","Simulations show binary companions twist CO line profiles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim rests on the free-wind approximation, in which the primary's gravity is artificially balanced by radiation force and wind particles are launched at the dust-condensation radius with a fixed velocity; the paper itself states (Sect. 4.5) that the resulting inner-wind velocity profile is unrealistic, and the wind speed at the companion is exactly what controls the morphology that creates the predicted line features.","fun_headline_variants_meta":{"raw":{"variants":["Binary companions leave telltale double peaks in CO lines","AGB binaries distort CO lines, but noise can mask the clue","CO line shapes betray hidden binary companions in AGB winds","Double-peaked CO emission signals a binary, unless noise hides it","Simulations show binary companions twist CO line profiles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000655,"raw_usage":{"total_tokens":3090,"prompt_tokens":1124,"completion_tokens":1966,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":740,"completion_tokens_details":{"reasoning_tokens":1883}},"tokens_in":740,"tokens_out":1966,"duration_ms":13915,"temperature":1.0,"reasoning_tokens":1883,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:07:04.276445+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a confirmed binary AGB system whose orbital parameters are known from astrometry or resolved imaging, with a slow wind (≈5 km/s) and a close orbit (≈9 au), and measure its low-J CO lines at high signal-to-noise with beam effects removed: if the lines do not show the predicted double-peaked (face-on) or central-bump (edge-on) structure across several transitions, the velocity-bin mechanism is falsified. A positive control would be a known face-on narrow-spiral system, which the models say should show the two-peaked profile rather than a smooth parabola.","supporting_citations":[{"cited_title":"& Jorissen, A","cited_arxiv_id":null,"evidence_quote":"Introduced the free-wind approximation and the prediction that a binary companion shapes the AGB outflow, the conceptual foundation of the models."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplied the spherical CO photodissociation prescription that the paper generalises to 3D geometry."},{"cited_title":"F., & Black, J","cited_arxiv_id":null,"evidence_quote":"Provided the unshielded photodissociation rate and the CO and H2 shielding functions used in the 3D photodissociation calculation."},{"cited_title":"2021, A&A, 652","cited_arxiv_id":null,"evidence_quote":"Defined the backward and frontal spiral edges and described the binary-wind interaction producing the spiral and arc morphologies seen here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier hydrodynamic binary models whose velocity-bin column-density proxy the paper supersedes with full NLTE radiative transfer."},{"cited_title":"2015, A&A, 579, A118","cited_arxiv_id":null,"evidence_quote":"Earlier analytical spiral models whose predicted line shapes the paper compares with its own double-peak and central-bump features."}],"review_version":1}