{"id":"67dea1d3-8a99-4a59-82dd-4f8637629a3b","arxiv_id":"2411.17561","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Three hydrodynamic mechanisms, shock levitation, Rayleigh-Taylor plumes, and convective rebound, lift cool molecular gas above evolved stars and can explain dimming events like Betelgeuse's Great Dimming.","lead":"Using 3D radiation-hydrodynamics simulations, this paper identifies three mechanisms that lift cool molecular gas above the surfaces of evolved stars, creating dark patches that dim the star. The results offer a physical explanation for events like the 2019/2020 Great Dimming of Betelgeuse.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Betelgeuse link rests on extrapolation from 5-8 Msun models that are acknowledged to have deficient interiors; no model at M>8 Msun tests the mechanism.","rationale":"The central claim has good support within its own domain: the simulations are global 3D RHD, the dark patches are emergent, the opacity/molecular interpretation is checked with post-processed spectra, and the authors explicitly flag the model deficiencies in Sect. 6.6. My concern is not that the mechanisms are impossible; it is that the move from 'these processes occur in 1-8 Msun models' to 'they can account for a 15-20 Msun RSG like Betelgeuse' is a parameter-space extrapolation across an unmodeled gap. The missing radiation pressure in the interior is not a small detail for M>10 Msun: it changes the equation of state, the pulsation eigenfunctions, and the convective fluxes, i.e., precisely the drivers of all three dimming scenarios. The paper honestly says the resulting amplitudes are reduced and the atmosphere is too compact, so the model is not a faithful Betelgeuse even before mass extrapolation. This could cut either way: a real Betelgeuse with stronger pulsations might produce more frequent or deeper dimming (favoring the claim), or the more extended, lower-density atmosphere might suppress molecule formation (weakening it). Because the sign and magnitude are unknown, the paper's conclusion should remain conditional until a Betelgeuse-like model or a quantitative 1D pulsation test is run. The rebound scenario being based on a single event is a secondary concern; the RT-plume mechanism provides an independent route, and one robust mechanism is sufficient for the claim 'can account for'.","tokens_in":26363,"tokens_out":7779,"duration_ms":80343,"concrete_test":"Create a Betelgeuse-like model by first building a 1D stellar envelope (M=15 Msun, L~1e5 Lsun, Teff~3600 K, solar abundances, including radiation pressure in the EOS), extracting the time-dependent radial velocity history from its pulsation, and feeding it into a high-resolution 3D CO5BOLD atmospheric box. Run for at least 10 pulsation periods and check whether cool, optically thick molecular clouds (T<2500 K at tau_Ross~1) form. If they do not, the extrapolation in Sect. 6.6 fails; if they do, quantify their size, contrast, and occurrence rate and compare with the 5/8 Msun models and with Betelgeuse observations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is not the mechanism itself but its transfer to Betelgeuse. The 5 and 8 Msun models are the only members of the grid with 3000-3400 K atmospheres, and Sect. 6.6 states that CO5BOLD models of higher-mass RSGs suffer from an insufficient description of the stellar interior at higher temperatures (likely missing radiation pressure, possibly missing self-gravity, and a contribution from the damping inner boundary condition), which reduces pulsation amplitudes and large-scale convective flows and presumably makes the atmosphere too compact. Since all three scenarios (Sects. 3-5) depend on the interaction of convective flows with pulsationally driven shocks, a systematic suppression of either ingredient could change the existence, frequency, and contrast of the dark patches. The grid also has no model between 8 Msun and Betelgeuse, so the statement that the mechanisms 'will also work' in massive RSGs is an assumption, not a result; the conclusion itself says 'We assume' and 'We speculate'. If the missing physics changes pulsation-convection coupling in the opposite direction, the model could either under- or over-produce dark patches; either way, the specific quantitative claim for Betelgeuse is not yet established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a set of global 3D radiation-hydrodynamical 'star-in-a-box' simulations of evolved stars (M⋆ = 1.0, 1.5, 5.0, and 8.0 M⊙) computed with the CO5BOLD code, which include self-excited pulsations, convection, and shock-dominated atmospheric dynamics. It identifies three mechanisms that lift cool, dense gas above the top of the convection zone and produce dark patches in bolometric intensity maps: (i) Rayleigh-Taylor plumes triggered by the combination of large convective upflows and strong radial pulsations (Sect. 3, demonstrated in the 5 M⊙ model st35gm04n045); (ii) a rare 'convective rebound' event following the merging of deep downdrafts, which produced a single prominent dimming episode in the 5 M⊙ model st35gm04n048 (Sect. 4); and (iii) nearly global obscuration by cool post-shock gas in the lower-mass, cooler AGB models (Sect. 5). The authors combine intensity maps, radial profiles, ray-based spectrum synthesis, and light curves to argue that the associated cooling, molecular-opacity increase, and TiO/H2O absorption can account for the dark patches and strong dimming events of red supergiants such as Betelgeuse, while acknowledging (Sect. 6.6) that the application to Betelgeuse requires an extrapolation beyond the mass range of the models. The final section explicitly frames the transfer to massive RSGs as an assumption and the dust-forming corollary as speculation.","tokens_in":26630,"tokens_out":12748,"duration_ms":171930,"significance":"If the mechanisms hold up, the paper provides a physically concrete, three-channel classification of dimming phenomena in evolved stars that can be confronted with interferometric imaging, photometric light curves, and time-resolved spectroscopy of AGB stars and red supergiants. The main strengths are that the dark patches are an emergent outcome of self-consistent 3D RHD models rather than a prescribed boundary condition; that the analysis connects the hydrodynamics (dynamical pressure, Rayleigh-Taylor instabilities, shocks) to observable signatures (molecular opacity feedback, TiO/H2O bands, wavelength-dependent darkening); and that the paper is unusually candid about its own limitations, including the manual selection of snapshots (Sect. 2.3), the single-event basis of the rebound scenario (Sect. 4), and the admitted interior deficiencies of CO5BOLD models for higher-mass RSGs (Sect. 6.6). The steps that currently carry the largest uncertainty are the transfer of mechanisms demonstrated at 5-8 M⊙ to a >10 M⊙ red supergiant, and the translation of local dark patches into disk-integrated dimming amplitudes; both are flagged as assumptions or extrapolations in the text.","major_comments":[{"comment":"The paper asserts that the local darkening produces 'a strong dimming of the disc-integrated light' with V-band fluxes 'decreasing by orders of magnitude' and that the processes have 'little effect on the shape of disk-integrated spectra', but the manuscript presents only ray-based spectra (Fig. 9) and bolometric observer-integrated normalized intensities (Figs. 5, 12, 14); no synthetic disk-integrated spectrum, broadband (e.g., V-band) light curve, or spectroscopically derived effective temperature is computed for any dimming epoch. Since these wavelength-dependent, disk-integrated quantities are precisely the observables used in Sect. 6.6 to connect the model to the Great Dimming of Betelgeuse, this is a load-bearing gap; the authors should compute synthetic disk-integrated photometry/spectra for the dimming events in st35gm04n045, st35gm04n048, and st28gm05n056, or explicitly downgrade the corresponding claims.","section":"§1, §3.5, §6.6"},{"comment":"The transfer of the mechanisms to Betelgeuse rests on an acknowledged extrapolation: Sect. 6.6 states that CO5BOLD models of higher-mass RSGs suffer from an insufficient description of the stellar interior (likely missing radiation pressure, possibly missing self-gravity, and a contribution from the damping inner boundary condition), which reduces pulsation amplitudes and large-scale convective flows and presumably makes the atmosphere too compact, and the grid contains no model between 8 M⊙ and Betelgeuse's mass. The conclusion that the mechanisms 'will also work in the slightly hotter atmospheres of massive RSGs' is labeled 'We assume' (Sect. 7). Because both the plume and rebound scenarios depend on the coupling between pulsations and convection, this is a correctness-risk rather than a presentation issue; I ask for a concrete scaling argument (for example, using the Pdyn/P and Mach-number profiles in Fig. 15 together with Betelgeuse's observed pulsation amplitudes and surface gravity) that shows the expected direction and magnitude of the bias, or a clear statement of which observations would discriminate the scenarios.","section":"§6.6, §7"},{"comment":"The claimed frequencies are not quantified. The abstract states that AGB stars show dark patches 'regularly during intermediate phases' while more massive stars show them 'only intermittently', and that the rebound process is 'the third and rarest scenario', but Sect. 2.3 states that snapshots were manually selected, and Sect. 4 reports only one prominent rebound event in 'around a hundred models' without specifying the search criterion, the number of models inspected, or the definition of 'prominent'. Either provide a systematic (even simple) detection statistic for dark patches and light-curve minima across the available model runs, or explicitly qualify these frequency statements as preliminary, since the scenario classification by rarity is part of the paper's claims.","section":"§2.3, §4, Abstract"}],"minor_comments":[{"comment":"The top-row caption lists 'st28gm06n056' for the model with Teff≈2848 K and refers to it as the opposite side of the model in Fig. 13, but Table 1 and Fig. 13 identify that model as st28gm05n056; the caption should be corrected.","section":"Fig. 1 caption"},{"comment":"The final sentence of the caption contains a typo: 'The timer for this model does not start at zero' should read 'The time for this model does not start at zero'.","section":"Fig. 14 caption"},{"comment":"The ray-based spectra in Fig. 9 are computed at resolution R=200 with a microturbulent velocity of 2.5 km/s and without velocity-field effects; given that the dark patch material moves at tens of km/s relative to the photosphere, a brief statement of how Doppler shifts would affect the molecular-band depths would strengthen the quantitative interpretation.","section":"§3.5"},{"comment":"The caption of Fig. 1 states that the snapshots of the 5 and 8 M⊙ models were 'selected such that large-scale — comparably rare — dark patches are visible'; it would be useful to state in the text how many snapshots per model were inspected and how the displayed instances were chosen, to give the reader a sense of the selection bias.","section":"§2.3"},{"comment":"The symbol 'Maradial' in the caption appears to be the radial Mach number with a subscript that has not rendered properly; please fix the notation.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is likely to be of interest to A&A readers, and the emergent-model nature of the result is a genuine strength. The main risk is that the abstract and conclusions state the Betelgeuse connection more strongly than the presented analysis supports; the authors themselves hedge in Sects. 6.6 and 7. I recommend major revision focused on the disk-integrated photometry/spectra computation and on a scaling argument for the extrapolation, rather than rejection, because the mechanisms in the models are well documented. The single-event rebound example and the manual snapshot selection should also be addressed in the revision; if a systematic statistic cannot be provided, the frequency claims in the abstract should be explicitly downgraded."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the mechanism taxonomy, not for the Betelgeuse verdict. The paper does something real: it uses 3D RHD simulations of evolved stars (1–8 Msun) to separate three distinct ways cool molecular gas gets lifted above the photosphere and causes dark patches – shock-driven levitation (familiar from 1D AGB models), Rayleigh-Taylor plumes seeded by convective upflows plus radial pulsations, and a rarer convective rebound after downdraft merging. The Rayleigh-Taylor plume analysis, with its time sequences, radial profiles, and spectra showing deep TiO/H2O absorption, is the most convincing part. The models are legitimate 'star-in-a-box' CO5BOLD runs; the code and opacity tables are independently validated in earlier work, so the self-citation here is not a circularity problem.\n\nThe soft spots are real but mostly transparent. The rebound scenario rests on one strong event in one model; the snapshots are manually selected; and the Betelgeuse application is an extrapolation. Section 6.6 says outright that CO5BOLD models of higher-mass RSGs suffer from an insufficient interior description (likely missing radiation pressure, possibly missing self-gravity, and a damping inner boundary condition), which reduces pulsation amplitudes and convective flows. There is no model between 8 Msun and Betelgeuse's mass. So the claim that these mechanisms 'will also work' in massive RSGs is, as the authors say, an assumption. They label it as such, which is honest. But it means the paper's quantitative relevance to the 2019/2020 dimming is not yet established.\n\nThat said, the core result – the identification and characterization of the three mechanisms in the simulated regime – stands on its own. The extrapolation problem is a scope limitation, not an internal contradiction. The paper is carefully argued, the analysis is detailed, and the distinction between mechanisms is a useful organizing frame for future work and for interpreting interferometric and photometric observations.\n\nI'd send this to a serious referee. A good referee should push on the higher-mass extrapolation and on the statistical weight of the rebound event, but neither issue warrants a desk reject. People working on evolved-star variability and on the Betelgeuse dimming will want to cite this.","headline":"A solid mechanism-taxonomy paper: three distinct ways cool molecular gas dims evolved stars, with an honest but real extrapolation gap to Betelgeuse.","tokens_in":27116,"tokens_out":2103,"would_cite":true,"duration_ms":19982,"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":"Episodic levitation of dense gas clumps above the stellar surface, cooling into molecules, explains dark patches and strong dimming events in evolved stars such as Betelgeuse.","keywords":["red supergiants","asymptotic giant branch stars","radiation-hydrodynamics","convection","stellar pulsations","shock waves","molecular opacities","Betelgeuse dimming"],"falsifier":"Observe a resolved dimming event on a Betelgeuse-class star with simultaneous imaging and spectroscopy: the molecular-gas model predicts a dark patch with strong TiO and H2O absorption from gas cooler than about 2600 K, outward motion within the patch, the strongest flux drop in the visual band, and little spectral change outside the patch. A patch that instead shows gray extinction without molecular bands, mid-infrared thermal dust emission, or a dark feature moving across the disk as a foreground cloud would falsify the model for that event.","tokens_in":26196,"feed_emoji":"🔭","tokens_out":8676,"duration_ms":71358,"temperature":0.7,"pith_summary":"This paper aims to establish that the dramatic dimming seen in evolved cool stars—most prominently Betelgeuse's 2019/2020 Great Dimming—can be caused by clouds of cool molecular gas levitated above the stellar surface, with dust playing at most a secondary role. The authors analyze global 3D radiation-hydrodynamics models of stars from 1 to 8 solar masses and identify three distinct mechanisms that lift dense gas into the lower atmosphere: material carried outward by strong shocks, Rayleigh-Taylor plumes created when a large convective upflow combines with strong radial pulsation, and rebounding deep downdrafts after an enhanced convective flow. In every case the lifted gas cools to roughly 2000 K, molecules such as TiO and water form, the opacity jumps sharply, and the gas appears as a dark patch against the bright stellar disk; a sufficiently large patch dims the disk-integrated light. A sympathetic reader would care because this connects two already observed features of red supergiants, convection and pulsation, to the resolved dark patches and light-curve minima, and it explains the dimming without requiring a foreground dust cloud.","feed_headline":"Cool gas clouds, not dust, drive Betelgeuse-style dimming","feed_subtitle":"3D simulations show three ways convection and pulsation lift dense gas that darkens and blocks starlight.","key_machinery":"The central object is the 'star-in-a-box' global 3D radiation-hydrodynamics model computed with the CO5BOLD code: a Cartesian box containing a whole star, with a spherically symmetric gravitational potential, nonlocal radiative transfer, gray Rosseland opacities extended by low-temperature molecular opacity tables, a tabulated equation of state that includes H2 formation, and a small energy-source core. The machinery couples large-scale convection in the stellar interior, self-excited radial and non-radial pulsations, and strong radiative shocks in the atmosphere; from its output, bolometric intensity maps, radial density-temperature-opacity profiles, and synthetic spectra are used to show where, when, and why a dark patch forms.","core_discovery":"The paper's central claim is that the episodic levitation of dense gas clumps above the stellar surface, followed by the formation of complex molecules in the cooling gas and possibly dust grains at a later stage, can account for the dark patches and strong dimming events of supergiant stars such as Betelgeuse. In the simulations the dimming is caused by gas alone: dense material is pushed above the top of the convection zone by one of three mechanisms, expands and radiatively cools, and once its temperature drops below roughly 2000 to 2600 K the steeply rising molecular opacities, TiO in the optical and H2O in the near-infrared, make the clump optically thick and dark against the photosphere. The same physics produces nearly global obscuration in low-mass AGB stars during pulsation minima and rarer, localized patches in 5 to 8 solar-mass models, and the authors propose the Rayleigh-Taylor plume and convective-rebound scenarios as the plausible mechanisms for Betelgeuse.","pith_inferences":["Extension: if improved higher-mass models produced stronger pulsations and less compact atmospheres, plumes and rebounds would likely become more frequent, so the current 1 to 8 solar-mass results may underestimate how often Betelgeuse-class stars dim.","Extension: the three scenarios are hard to distinguish in a single image but have different anchors—plumes form above global convective upflows, rebounds above downdrafts, and post-shock clumps are not tied to either—so time-resolved Doppler maps of a forming patch could identify the mechanism in a real event.","Extension: a testable prediction follows: during a dimming minimum, spectroscopy of the dark patch should show strong TiO and H2O bands from gas cooler than about 2600 K with outward line-of-sight velocities, and no coincident thermal dust emission; if those signatures are absent, a molecular-gas cloud is not the cause."],"forward_implications":["Dark patches appear in all of the models whenever dense gas is levitated above the convection zone, cools below roughly 2000 to 2600 K, and forms molecules; the patches are never uniform and show substructure.","In low-mass AGB models the same mechanism produces almost global obscuration at pulsation minima, while in 5 to 8 solar-mass models only localized, intermittent patches form, mostly near luminosity minima.","For Betelgeuse, the Rayleigh-Taylor plume and convective-rebound scenarios are the viable candidates, and the dimming should be tied to the pulsation cycle and stay in one region of the disk, as the resolved images show.","Because the obscuring material is molecular gas, the visual flux can drop by orders of magnitude while the spectral energy distribution outside the patch changes little; dust is not required for the dimming itself."],"supporting_citations":[{"why":"Supplies the CO5BOLD hydrodynamics and radiative-transfer solver used for all models.","marker":"Freytag et al. 2012"},{"why":"Establishes the global AGB-star models with self-excited pulsations and realistic convective surface structure that this paper extends.","marker":"Freytag et al. 2017"},{"why":"Provides the molecular opacity tables responsible for the low-temperature opacity rise and the synthetic spectra.","marker":"Aringer et al. 2016"},{"why":"Supplies the pulsation-period and radial-velocity diagnostics used to characterize the model dynamics.","marker":"Ahmad et al. 2023"},{"why":"Shows how density and temperature inhomogeneities in the gas imprint on dust, grounding the dust-related discussion.","marker":"Höfner & Freytag 2019"},{"why":"Demonstrates dust formation above a dark patch in AGB models, supporting the claim that dust can form after the molecular gas.","marker":"Freytag & Höfner 2023"},{"why":"Provides resolved images of Betelgeuse showing the darkening appearing and disappearing in the same disk region, the key observation to explain.","marker":"Montargès et al. 2021"},{"why":"Gives a two-component gas model of Betelgeuse's dimming, a benchmark the simulations support without requiring dust.","marker":"Harper et al. 2020b"}],"fun_headline_variants":["Gas clouds, not dust, drive Betelgeuse-style dimming","Three gas-cloud scenarios explain Betelgeuse's dimming","Simulations reveal how gas clouds dim evolved stars","Betelgeuse dimming caused by levitated gas clouds","Not dust: gas clouds darken supergiants like Betelgeuse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the mechanisms seen in models of 1 to 8 solar-mass stars also operate in Betelgeuse-class supergiants above 10 solar masses, even though the models themselves are acknowledged to have too-compact atmospheres and reduced pulsation amplitudes at higher masses.","fun_headline_variants_meta":{"raw":{"variants":["Gas clouds, not dust, drive Betelgeuse-style dimming","Three gas-cloud scenarios explain Betelgeuse's dimming","Simulations reveal how gas clouds dim evolved stars","Betelgeuse dimming caused by levitated gas clouds","Not dust: gas clouds darken supergiants like Betelgeuse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":1969,"prompt_tokens":1099,"completion_tokens":870,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":715,"completion_tokens_details":{"reasoning_tokens":782}},"tokens_in":715,"tokens_out":870,"duration_ms":7981,"temperature":1.0,"reasoning_tokens":782,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:58:34.914378+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a resolved dimming event on a Betelgeuse-class star with simultaneous imaging and spectroscopy: the molecular-gas model predicts a dark patch with strong TiO and H2O absorption from gas cooler than about 2600 K, outward motion within the patch, the strongest flux drop in the visual band, and little spectral change outside the patch. A patch that instead shows gray extinction without molecular bands, mid-infrared thermal dust emission, or a dark feature moving across the disk as a foreground cloud would falsify the model for that event.","supporting_citations":[{"cited_title":"2012, J.Comput.Phys., 231, 919","cited_arxiv_id":null,"evidence_quote":"Supplies the CO5BOLD hydrodynamics and radiative-transfer solver used for all models."},{"cited_title":"2016, MNRAS, 457, 3611","cited_arxiv_id":null,"evidence_quote":"Provides the molecular opacity tables responsible for the low-temperature opacity rise and the synthetic spectra."}],"review_version":1}