{"id":"b89b7df8-b8b6-4905-ab4c-1b7d85006926","arxiv_id":"2506.02812","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A torus, jets, and inner ellipsoid model reproduces most H-band interferometry and low polarization of V838 Mon, but leaves the K-band closure phases unexplained.","lead":"Astronomers built a 3D dust model of the remnant of the 2002 stellar merger V838 Mon, with a dusty torus, faint jets, and an inner shell. The model explains much of the near-infrared interferometric signal and the low polarization, but fails on the smallest K-band asymmetries, so the structure is plausible but not unique.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The physical inference hinges on dust surviving at 5000–7000 K; the authors' own sublimation test removed nearly all dust, so the static torus/jets/ellipsoid configuration is not a self-consistent radiative-equilibrium structure and may not support the claimed persistent geometry.","rationale":"The reader's weakest assumption—dust surviving at temperatures far above the adopted sublimation threshold—is exactly the most load-bearing concern. The paper is honest about the problem, disclosing that the sublimation correction removed nearly all dust and that the hottest dust reaches about 7000 K, but the conclusions nevertheless state a persistent torus and jets. This is not a mere numerical detail; it determines whether the modeled density distribution can be a real, stable circumstellar structure. A self-consistent calculation with dust sublimation would likely remove the torus and the jet-induced gap, which is the sole source of asymmetry in the H band, and would therefore eliminate the basis for the central claim. The polarimetric data are consistent with low polarization but are too weak to break the degeneracy, and the SiO maser discussion does not provide strong constraints on the dust density. The K-band failure and the acknowledged non-uniqueness further weaken the inference but are not as fundamental as the dust-survival inconsistency. Because the reader's CONDITIONAL verdict already reflects this unresolved issue, no change to the verdict is needed; however, the concern should be stated prominently in any publication decision.","tokens_in":19732,"tokens_out":5304,"duration_ms":67559,"concrete_test":"Re-run the RADMC3D model with the iterative sublimation correction enabled: remove all cells with equilibrium dust temperature above 1700 K (and, as a sensitivity case, above 2700 K), recompute the radiation field, and iterate to convergence. Then recompute the synthetic H-band and K-band squared visibilities and closure phases from the converged density distribution. If the surviving dust mass is negligible, or if the H-band closure-phase signal and the torus gap disappear, then the jets-vital asymmetry mechanism and the persistent-structure conclusion are unsupported and a time-dependent dust-formation/destruction model is required.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sect. 3 the RADMC3D radiative-equilibrium calculation yields dust temperatures up to about 7000 K, while the adopted amorphous-silicate sublimation temperature is 1700 K (and even 2700 K for the most refractory alumina considered). The paper states that applying the sublimation correction iteratively removed nearly all dust out to the outermost modeled radii, and that increasing the sublimation temperature to 2700 K did not alleviate the problem. The correction was therefore disabled. This is a direct inconsistency between the assumed density distribution and the physical survival of dust. If real dust cannot exist at the modeled temperatures, the torus, jets, and ellipsoid cannot persist as static structures; the observed dust would require either self-shielding at optical depths much higher than those modeled, or continuous dust production and destruction, which the static RADMC3D framework cannot represent. The abstract's conclusion of a persistent torus-like structure and jets therefore rests on an unresolved physical assumption. The authors acknowledge the issue and speculate about transient dust, but the central structural interpretation is not supported by a self-consistent calculation. The K-band mismatch and acknowledged model non-uniqueness are secondary; the dust-survival problem is the load-bearing flaw.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents RADMC3D radiative-transfer modeling of the circumstellar dust around V838 Mon, the luminous red nova remnant, aimed at constraining the 3D dust distribution using VLTI/GRAVITY K-band and CHARA/MIRC-X H-band interferometric observables (squared visibilities and closure phases) that were imaged by the same group in Mobeen et al. (2024). The adopted model has three components: an inner ellipsoid around the 3500 K central star, a tilted torus with inner and outer radii of 3 au and 5 au, and jets perpendicular to the torus plane; component densities, orientations, and sizes are tuned by trial and error. The paper reports that the model reproduces the general H-band visibility trend and about half of the H-band closure-phase population, but fails on the K-band observables: the model squared visibilities stay above 0.6 at spatial frequencies where the observed values tend toward zero, and the modeled closure phases reach only about 0.2 degrees, an order of magnitude below the observed scatter. New SAAO-HIPPO polarimetric observations show intrinsic polarization below 2% with large uncertainties, which the authors argue is consistent with the model. The authors conclude that, although not unique, the model suggests a persistent torus-like structure and jets, possibly signposts of the remnant's evolution toward a blue straggler.","tokens_in":19998,"tokens_out":8911,"duration_ms":92950,"significance":"If the H-band inference were sustained, this would be the first attempt to place three-dimensional dust-density constraints on the innermost environment of a luminous red nova remnant, with implications for merger remnant angular-momentum loss, jet launching, and the blue-straggler precursor question, all of which are timely topics. The paper has real strengths: the synthetic observables are computed on the actual OIFITS baseline and wavelength sampling (Section 3.1, via AMHRA), a large family of alternative models is explored and documented in the appendices, the full-Stokes polarization treatment is included, and the authors are unusually candid about the K-band failures and the non-uniqueness of the model. These virtues, however, are currently outweighed by a load-bearing physical inconsistency: the radiative-equilibrium calculation yields dust temperatures up to about 7000 K while the adopted silicate sublimation temperature is 1700 K, and the authors' own iterative sublimation test removed nearly all dust, so the correction was disabled (Section 3).","major_comments":[{"comment":"The dust-survival inconsistency is load-bearing for the central conclusion. The radiative-equilibrium calculation gives dust temperatures up to about 7000 K, far above the adopted 1700 K silicate sublimation temperature, and the authors state that applying the sublimation correction iteratively removed nearly all dust out to the outermost modeled radii, and that even increasing the sublimation temperature to 2700 K did not help. The correction was therefore disabled, and the final model places dust in regions where, by the paper's own opacity and sublimation assumptions, it cannot survive. The paper acknowledges this and invokes continuous dust production and destruction, but the static RADMC3D framework cannot represent such a process, so the abstract's and Section 5's inference of a 'persistent torus-like structure and jet or jets' is not supported by the calculation as presented. This needs to be addressed quantitatively, for example by reporting the mass fraction of dust with T > 1700 K in the final model, testing whether higher torus optical depths (self-shielding) keep the dust cool enough to survive, and, until such a test succeeds, reformulating the conclusions to describe a brightness-distribution model whose physical persistence is an open question.","section":"Section 3 (paragraph beginning \"In our initial RADMC3D runs\")"},{"comment":"The K-band model is in clear quantitative disagreement with the observations, and this disagreement is acknowledged in the abstract but not properly integrated into the conclusions. The simulated squared visibilities reach a minimum of about 0.6 at 67 Mlambda while the observed values tend toward zero at higher spatial frequencies, and the simulated closure phases rise to only about 0.2 degrees against an observed scatter roughly an order of magnitude larger. Because the same physical dust distribution emits in both H and K bands, a model that describes the H-band structure but fails the K-band visibilities and closure phases cannot be said to constrain the physical structure of the circumstellar environment without a discussion of why the two bands decouple. The conclusions (Section 5) should be limited to the H-band morphology, and the K-band failure should be presented as an unresolved constraint on the model rather than a secondary caveat.","section":"Section 3.2 and Figure 4"},{"comment":"The model is not an independent test of the inferred morphology. The paper states that the model construction 'took a cue from the image reconstructions presented in Mobeen et al. (2024)', i.e., the three-component structure, the bipolar orientation, and the size scale were all read off the same H- and K-band interferometric data that the model is later compared with. The comparison therefore demonstrates that a parametrized version of the reconstructed morphology can reproduce part of the observables, which is weaker than the claimed confirmation of a torus and jets. The only independent dataset, the HIPPO polarimetry, is consistent with zero intrinsic polarization within the uncertainties in the V, R, and I bands (Table 1: 1.095 +/- 0.471, 0.397 +/- 0.498, and 0.236 +/- 0.767 percent, respectively), and the average intrinsic position angle of 16.3 +/- 39.9 degrees is unconstraining. This does not invalidate the model, but it should be stated that the polarimetric check cannot discriminate between the proposed geometry and a nearly centro-symmetric envelope.","section":"Section 3 (model construction) and Section 4.1"},{"comment":"No quantitative goodness-of-fit measure is provided for any model. All comparisons between simulated and observed squared visibilities and closure phases are made by eye ('we are able to somewhat reproduce', 'qualitatively reproduces observations quite well'), and no chi-square, reduced chi-square, or parameter uncertainties are reported. Given the large number of free parameters (torus inner and outer radii, inner ellipsoid axes and orientation, torus orientation, jet geometry, component densities, and grain properties for the polarization model), the claim that the adopted model provides an 'adequate fit' and the relative ranking of the main model against the appendix alternatives is not established without a fit statistic. At minimum, the authors should report the per-band residual statistics for V2 and closure phase for the main model and for the alternative models shown in Appendices B, C, and D.","section":"Sections 3.1-3.3 and Appendices B-D"}],"minor_comments":[{"comment":"The abstract contains several typographical errors that should be corrected before submission: 'assymetry' (asymmetry), 'intermedite' (intermediate), 'distibution' (distribution), and 'remant star' (remnant star).","section":"Abstract"},{"comment":"The phrase 'Monte Carlo Markov Chain (MCMC) algorithms' is a misnomer in this context: RADMC3D performs Monte Carlo radiative transfer (thermal Monte Carlo), not Markov Chain Monte Carlo fitting. Please reword to avoid confusing the radiative-transfer method with a sampling algorithm.","section":"Section 3"},{"comment":"The sentence 'the MIRC-X and GRAVITY interferometric observations are, however, complemented by original and and nearly contemporary polarimetric measurements' contains a duplicated 'and'.","section":"Section 2.2"},{"comment":"The sentence 'It was not possible to explain the remaining 40% of the unresolved flux' is ambiguous. The issue is that 40% of the flux remains unresolved in the model (the squared visibility floor of 0.6), meaning the modeled structure is too compact relative to the observations; please rephrase to say that the model overpredicts V2 at high spatial frequencies.","section":"Section 3.2"},{"comment":"The conclusion that in the K band 'the squared visibilities are underestimated' is confusing, since the model visibilities are larger than the observed ones at high spatial frequencies. The correct statement is that the model overpredicts the K-band squared visibilities (the model is under-resolved compared with the data).","section":"Section 5"},{"comment":"The abstract mentions interferometric observations in the 'HKLM bands', but the present analysis uses only H- and K-band data (the L- and M-band data are explicitly set aside in Section 2.1). Please rephrase so the abstract does not imply that L- and M-band data are modeled in this paper.","section":"Abstract and Section 1"},{"comment":"The two H-band closure-phase populations (one centered at zero, one at about 10 degrees) are central to the paper's claim, but the figure captions do not identify which baseline triangles or observing configurations produce each population. A short caption addition would help the reader evaluate the claimed partial reproduction.","section":"Figures 4 and 5"}],"recommendation":"major_revision","confidential_remarks":"The authors are commendably open about the model's limitations: the abstract concedes the K-band failure, the text admits non-uniqueness, and the sublimation problem is described explicitly in Section 3. That candor does not, however, remove the load-bearing inconsistency: the central physical inference of a persistent torus and jets rests on disabling the sublimation check and on a trial-and-error tuning to the same data the model is later compared with. The two independent checks (polarimetry and SiO maser spot locations) are too weak to discriminate among models. If the authors quantify the mass of dust above the sublimation temperature, test whether self-shielding at higher torus optical depths changes the outcome, and reframe the conclusions as a brightness-distribution description with the physical persistence explicitly open, the paper could become publishable as a modeling study with honest caveats. I do not see this as a rejection case because the geometric description of the H-band structure is still of interest and the required changes are within the scope of a revision; the current version, however, claims more than the calculation supports."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is a real 3D radiative-transfer dust model of V838 Mon's inner region, grounded in the same VLTI/CHARA data as the earlier geometric work, plus new SAAO/HIPPO polarimetry. The authors are honest about what fits and what does not. But the load-bearing assumption — that dust can survive at 5000–7000 K — is contradicted by their own sublimation test, which removed nearly all dust and was then switched off. Without a self-consistent answer to that, the torus-plus-jets structure is a plausible cartoon, not a constrained physical result.\n\nCredit where due: the RADMC3D setup is documented well enough to reproduce, the appendix explores alternatives, and the paper openly states the K-band failure, the non-uniqueness, and the sublimation problem. The polarimetry, even though weak (P<2%, consistent with zero intrinsic polarization after ISM subtraction), is a useful independent data point.\n\nWhere it falls short: the model was tuned by hand to the same H/K observables it is later compared with, and there is no fit statistic, no parameter uncertainties, and no search over the parameter space. The H-band agreement is qualitative — roughly half the closure-phase population is matched. The K-band model misses the observed visibilities by a large margin (model floor ~0.6 vs observed values approaching zero) and closure phases by an order of magnitude. The claimed jet-induced gap is the sole source of asymmetry, but that mechanism is only as good as the model's other assumptions. The polarimetric comparison is too weak to break degeneracies. On the dust survival issue, the paper's own speculation about transient dust production is not modeled; it is an admission that the static framework cannot represent the environment.\n\nWho this is for: anyone working on V838 Mon, red novae, or merger-remnant evolution, and anyone who wants a worked example of the limits of RADMC3D against interferometric data. It is not a paper that reshapes the field; its impact is specific to this object.\n\nRecommendation: it deserves a serious referee, but I would not accept it as is. The referee should require a quantitative treatment of dust survival — self-shielding or a time-dependent dust-formation term — and a proper parameter search with uncertainties. If that is not possible, the paper should be framed as an exploratory modeling study, not a constraint on the physical structure. My verdict: major revision, with the sublimation problem as the deciding issue.","headline":"A real but flawed step beyond geometric modeling: the dust-survival inconsistency undermines the torus/jet conclusion, yet the paper is honest and worth refereeing.","tokens_in":20536,"tokens_out":2797,"would_cite":false,"duration_ms":31876,"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":"Using radiative-transfer modeling, this paper argues that the dust around V838 Mon's remnant is arranged as an inner ellipsoid, a tilted torus, and jets that create the observed H-band asymmetry, while the K-band data remain unexplained.","keywords":["V838 Monocerotis","luminous red nova","circumstellar dust","radiative transfer","optical interferometry","polarimetry","stellar mergers","blue stragglers"],"falsifier":"Re-observe V838 Mon with CHARA/MIRC-X in the H band a few months after the 2022 epoch and compare the closure phases: the model attributes the signal to a stable gap in the torus produced by the jets, so a large change in the closure-phase pattern over that interval would show the inferred structure is transient rather than persistent.","tokens_in":19516,"feed_emoji":"🔭","tokens_out":10019,"duration_ms":87632,"temperature":0.7,"pith_summary":"V838 Monocerotis, a star that erupted as a luminous red nova in 2002 and is expected to evolve into a blue straggler, is surrounded by warm dust whose three-dimensional arrangement has been unclear. This paper tries to pin down that structure by fitting a radiative-transfer model to near-infrared interferometric observations from VLTI and CHARA, combined with new optical polarimetry. The authors argue that a three-component dust distribution—an inner ellipsoid, a tilted torus, and a pair of jets—reproduces the H-band observables well, with the jets creating a brightness gap in the torus that produces the measured asymmetry signal. The model also matches the very low (<2%) intrinsic linear polarization. It does not, however, capture the tiny closure-phase deviations in the K band, and the authors caution that the structure they infer is not unique.","feed_headline":"V838 Mon remnant may be wrapped in a dust torus with jets","feed_subtitle":"Radiative-transfer fits to H-band interferometry trace the post-merger structure expected before a blue straggler.","key_machinery":"The load-bearing mechanism is the radiative-transfer model built in RADMC3D: a $250^{3}$ grid with one central star (effective temperature 3500 K, luminosity $10^{5}$ solar luminosities) and a three-component dust density distribution. The components are an inner ellipsoid (3 au × 2 au × 0.5 au) that diffuses the stellar light, a massive tilted torus (inner radius 3 au, outer radius 5 au, density $10^{-14}$ g $cm^{-3}$) at position angle −60°, and two collimated jets at position angle 20° with a radially declining density of 2×$10^{-14}$ g $cm^{-3}$ divided by the distance in au. The jets are the crucial feature: in projection a jet overlaps the torus and absorbs some of its light, creating a brightness gap that breaks the torus's centro-symmetry. That gap is the only source of asymmetry in the model, so without the jets the simulated closure phases fall below 1° in the H band.","core_discovery":"The paper's central claim is that the dust immediately surrounding the V838 Mon remnant (within tens of astronomical units) can be described by a multi-component model comprising a compact inner ellipsoid around the star, a tilted torus at radii of 3–5 au, and collimated jets extending beyond 7 au. Comparing synthetic images from the RADMC3D code against VLTI and CHARA interferometric observables, the authors find that this arrangement adequately matches the H-band squared visibilities and closure phases, and that the jets are essential: where a jet overlaps the torus in projection it absorbs torus emission and leaves a gap, and that gap is the sole source of asymmetry that produces nonzero closure phases. The same model yields intrinsic linear polarization below 2% in the BVRI bands, consistent with new SAAO-HIPPO measurements once interstellar polarization is removed. The model is less successful in the K band, where the predicted closure phases are an order of magnitude smaller than observed and the squared visibilities never drop below about 0.6, so the authors conclude the K-band environment is more complex than the adopted geometry. Taken together, the results support—though they do not uniquely prove—a persistent torus and jets in the post-merger remnant.","pith_inferences":["A natural extension, not made in the paper, is to monitor the H-band closure phases repeatedly over months to years, which could measure jet precession or wobbling and directly test the jet-variability scenario for red novae cited in the paper.","Because the sublimation correction removes nearly all dust, the inferred structure is probably not in thermodynamic steady state; a physically consistent picture may require dust condensation near the star and rapid destruction, which would make the torus and jets transient signposts rather than permanent features.","The paper's own polarization maps show the jets are highly polarized (over 50% degree of linear polarization) while the torus is weakly polarized (under 2%), so wavelength-resolved imaging polarimetry could in principle separate the two components observationally, a test the authors do not propose.","The K-band degeneracy suggests that squared visibilities alone cannot break the model ambiguity; adding closure phases from more triangle configurations or combining H- and K-band data in a single simultaneous fit might do so."],"forward_implications":["If the model is correct, the H-band closure-phase signal directly traces a jet-induced gap in the torus, so CHARA closure phases become a probe of jet-torus interaction.","The inferred dust structure implies a persistent, geometrically ordered envelope—torus plus jets—around the merger remnant, matching theoretical predictions for how intermediate-mass merger products lose angular momentum.","The model's K-band failure indicates that the true environment contains additional fine-scale structures or asymmetries, so the K-band observables still lack a physical explanation.","The very low (<2%) intrinsic polarization means scattering is dominated by the torus rather than the jets, so polarization monitoring can track changes in the torus geometry.","The near-IR-bright dust mass (about 10^-6 solar masses) is three orders of magnitude smaller than the cold dust seen by ALMA, showing the warm inner structure is only a small skin of a much larger envelope."],"supporting_citations":[{"why":"Supplies the VLTI GRAVITY and CHARA MIRC-X interferometric observables and reconstructed images that the RADMC3D model is fitted to.","marker":"Mobeen et al. (2024)"},{"why":"The RADMC3D radiative-transfer code used to compute dust temperatures and synthetic images.","marker":"Dullemond et al. (2012a)"},{"why":"Defines the amorphous silicate opacity (Mg/Fe pyroxene) adopted for the dust in the model.","marker":"Dorschner et al. (1995)"},{"why":"Provides the interstellar polarization parameters used to subtract the ISM contribution and obtain intrinsic polarization.","marker":"Wisniewski et al. (2003)"},{"why":"Supplies the SiO maser spot positions and proper motions used to test the model's torus location.","marker":"Ortiz-León et al. (2020)"},{"why":"Provides the ALMA flattened-structure detection, the cold dust mass, and the stellar parameters (temperature and luminosity) adopted in the model.","marker":"Kamiński et al. (2021)"},{"why":"Theoretical argument that luminous red nova outbursts require jet launching from an accretion disk, motivating the jet component.","marker":"Soker (2023)"},{"why":"Simulations of jet wobbling in common-envelope mergers used to argue the jets and the associated asymmetry may be variable.","marker":"Dori et al. (2023)"}],"fun_headline_variants":["V838 Mon remnant: dust as torus and jets","Interferometry reveals V838 Mon's dusty jets and torus","Dust model for V838 Mon: jets and torus","V838 Mon's post-merger dust shaped by jets and torus","Jets and torus in V838 Mon's circumstellar dust"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model's static torus, jets, and ellipsoid depend on dust surviving at equilibrium temperatures up to about 7000 K, far above the 1700 K silicate sublimation temperature the authors adopted; applying the sublimation cutoff removed nearly all dust and was therefore abandoned.","fun_headline_variants_meta":{"raw":{"variants":["V838 Mon remnant: dust as torus and jets","Interferometry reveals V838 Mon's dusty jets and torus","Dust model for V838 Mon: jets and torus","V838 Mon's post-merger dust shaped by jets and torus","Jets and torus in V838 Mon's circumstellar dust"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000427,"raw_usage":{"total_tokens":2259,"prompt_tokens":1095,"completion_tokens":1164,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":711,"completion_tokens_details":{"reasoning_tokens":1075}},"tokens_in":711,"tokens_out":1164,"duration_ms":11877,"temperature":1.0,"reasoning_tokens":1075,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:15:39.720366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observe V838 Mon with CHARA/MIRC-X in the H band a few months after the 2022 epoch and compare the closure phases: the model attributes the signal to a stable gap in the torus produced by the jets, so a large change in the closure-phase pattern over that interval would show the inferred structure is transient rather than persistent.","supporting_citations":[{"cited_title":"Z., Kami´nski, T., Matter, A., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the VLTI GRAVITY and CHARA MIRC-X interferometric observables and reconstructed images that the RADMC3D model is fitted to."},{"cited_title":"1995, Astronomy and Astrophysics, v","cited_arxiv_id":null,"evidence_quote":"Defines the amorphous silicate opacity (Mg/Fe pyroxene) adopted for the dust in the model."},{"cited_title":"Bright common envelope formation requires jets","cited_arxiv_id":"2306.07702","evidence_quote":"Theoretical argument that luminous red nova outbursts require jet launching from an accretion disk, motivating the jet component."},{"cited_title":"2023, ApJ, 954, 143","cited_arxiv_id":null,"evidence_quote":"Simulations of jet wobbling in common-envelope mergers used to argue the jets and the associated asymmetry may be variable."}],"review_version":1}