{"id":"ff4932bc-0aad-4314-b845-76fa00a5b986","arxiv_id":"2508.11013","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"3D RMHD simulations with MURaM and RH1.5D radiative synthesis show plasmoid-mediated turbulent reconnection in small-scale events can produce coexisting hot and cool plasma that explains UV bursts connected with Ellerman bombs.","lead":"This paper uses a 3D radiation-MHD simulation of the Sun's lower atmosphere to show that magnetic reconnection in small-scale events forms plasmoids and flux ropes, creating mixed hot and cool plasma. The authors synthesize H-alpha and Si IV emissions and argue this explains observations where Ellerman bombs and UV bursts occur together.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on identifying plasmoids at ~15 km grid resolution with no convergence check; if the structures are numerical dissipation artifacts of MURaM's grid-scale diffusion, the EB-UV burst explanation is not established.","rationale":"The reader's conditional verdict is appropriate. The simulation uses the well-established MURaM code and RH1.5D for synthesis, and the synthetic Hα wing enhancements and double-peaked Si IV profiles are concrete, falsifiable predictions. The weakest link is causal: the paper claims plasmoid-mediated turbulent reconnection is responsible for the multi-thermal structure, yet the plasmoid identification is purely morphological (Figs. 2–4) and the grid spacing leaves current sheets only a few cells wide. The authors themselves note in Sec. 2 that the resolution is 'still much lower than' their 2D AMR runs, and Sec. 4 lists only the coronal boundary limitation, not a resolution or convergence check. Since MURaM's dissipation is numerical hyperdiffusion, the effective Lundquist number and hence the onset and wavelength of plasmoids are resolution-dependent quantities. Without a convergence study, the central claim that plasmoid instability operates in 3D RMHD small-scale events and explains EB-UV burst coexistence is not established beyond reasonable doubt. No allegation of error is intended; this is a missing verification step. The verdict therefore remains conditional, unchanged from the reader's assessment.","tokens_in":20200,"tokens_out":5265,"duration_ms":62351,"concrete_test":"Perform a resolution study: re-run the identical flux-emergence setup with uniform grid spacing reduced by a factor of 2 (Δx ≈ 11.7 km, Δy = Δz ≈ 7.8 km) over a subdomain containing the ROI for the time interval t = 4301–4489 s, and compare plasmoid count, flux-rope diameter, current-sheet thickness, peak temperature, and synthesized Hα and Si IV intensities. If these quantities change by more than about 20% between the two resolutions, the plasmoid structures and the resulting EB/UV burst emissions are not numerically converged, and the central claim is not supported by the present run.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim—that turbulent plasmoid-mediated reconnection in 3D RMHD produces the coexisting cool/hot plasmas underlying EBs and UV bursts—depends on the physical reality of the plasmoid-like structures shown in Figs. 2–4. The only support is morphological: blob-like isosurfaces, twisted field-line bundles, and fragmented hot regions. MURaM has no explicit resistivity; magnetic diffusion is numerical hyperdiffusion (Rempel 2017; Przybylski et al. 2022). With grid spacings of 23.4 km in x and 15.63 km in y and z, the thin current sheets that in 2D AMR runs (Ni et al. 2021; Cheng et al. 2024) produced plasmoids are resolved here by only a few cells, and the effective Lundquist number is set by grid-scale dissipation rather than by chromospheric physics. The text (Sec. 2) acknowledges the resolution is much lower than in the prior 2D AMR simulations, yet no convergence study, resolution comparison, or quantitative plasmoid metric is provided to establish that the filamentation is a genuine tearing-mode cascade rather than a numerical noise pattern. If the plasmoids are artifacts, the predicted co-spatial and co-temporal EB and UV burst emissions, and the reported height inversion, would not follow. This is the load-bearing weak point; it is not a disagreement with prior results, but a missing verification that the newly claimed 3D phenomenon is converged.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a 3D radiation-MHD simulation of magnetic flux emergence in the solar lower atmosphere using the MURaM code, with a domain extending from the upper convection zone to 8 Mm in the atmosphere. A flat magnetic flux sheet inserted in the convection zone produces repeated small-scale magnetic reconnection events; the authors focus on one region of interest (ROI) and identify thin elongated current sheets with plasmoid-like structures and twisted flux ropes, with hot (>20,000 K) and cool (<10,000 K) plasma coexisting in a fragmented pattern. Using RH1.5D for Halpha and Si IV profiles and optically thin synthesis for Si IV 139.4 nm and EUI 17.4 nm images, they find wing-enhanced Halpha signatures characteristic of Ellerman bombs and double-peaked Si IV profiles characteristic of UV bursts, sometimes at the same reconnection region. They also report a cool EB-like blob above hot plasma at heights >2 Mm and another event with hot UV-burst plasma at about 0.7 Mm. They conclude that 3D turbulent reconnection mediated by the plasmoid instability can explain temporally and spatially connected EBs and UV bursts and the multi-thermal emission structure.","tokens_in":20391,"tokens_out":6123,"duration_ms":67223,"significance":"If the results hold, this is a notable step: it is, to my knowledge, the first 3D RMHD simulation that identifies plasmoid-mediated reconnection in the cool lower atmosphere of the Sun in the context of EBs and UV bursts, extending earlier 2D/2.5D models (Ni et al. 2021; Cheng et al. 2024) to a more realistic 3D geometry. The use of a well-tested code (MURaM), the standard RH1.5D radiative transfer code, optically thin CHIANTI-based synthesis, and explicit comparison with observational Halpha and Si IV line characteristics are clear strengths. The predicted coexistence and height inversion of cool and hot plasmas provide a falsifiable observational diagnostic. However, the case-study nature and the lack of resolution verification limit the strength of the general conclusions.","major_comments":[{"comment":"The central claim that plasmoid-mediated reconnection produces the multi-thermal structure rests on the physical reality of the small-scale filamentation, but the paper offers no quantitative resolution verification. The grid spacing is 23.4 km in x and 15.63 km in y and z, MURaM's magnetic diffusion is numerical hyperdiffusion, and Section 2 itself notes the resolution is much lower than in the previous 2D AMR runs. Without a resolution study, a measured current-sheet width, a local Lundquist-number estimate, or a plasmoid size distribution, the blob-like isosurfaces in Figs. 2-4 could plausibly be grid-scale numerical dissipation artifacts rather than a converged tearing-mode cascade. Since the height inversion and the co-spatial EB/UV emission claims follow directly from this filamentation, please add either a convergence test at an intermediate resolution or a quantitative demonstration that the plasmoid sizes are well above the grid scale and that the heating is not governed by the numerical dissipation.","section":"Section 2 and Figs. 2-4"},{"comment":"The text states that 'plasmoid instability appears in most of these reconnection events' and that 'most small-scale reconnection lead to the formation of twisted magnetic flux ropes,' but only one event (the ROI) is analyzed in detail, and a second event is shown in Section 3.3. No event census, selection criterion, or statistical measure is provided to support 'most.' Please either quantify the fraction of events showing plasmoids and flux ropes with clear selection criteria, or weaken the claim to a case study, because the general conclusion as written exceeds the presented evidence.","section":"Section 3.1 and Section 4 (conclusions)"},{"comment":"The claim that cool EB-like plasma is located above hot plasma at heights greater than 2 Mm is based on a single cool blob and a single synthesized Halpha profile. The wing enhancement in Fig. 8b is shown after subtracting a nearby background profile, but no analysis is given of how sensitive this result is to the size or location of the background region, nor is the statistical significance of the enhancement estimated. Since this height inversion is one of the paper's headline results, please provide a quantitative characterization of the blob (temperature, density, line-center optical depth, and background-subtraction sensitivity) or show the same behavior in additional events.","section":"Section 3.3 and Fig. 8"}],"minor_comments":[{"comment":"The abstract and Section 1 contain typos such as 'comaprison' and 'Figsures'; the manuscript would benefit from a careful proofread.","section":"Abstract and Section 1"},{"comment":"The color bars and axis labels in these figures appear garbled in the version I reviewed, with repeated tick values and duplicated labels; please check the figure production pipeline and provide clean, readable color scales.","section":"Figs. 2, 3, 4"},{"comment":"The text refers to the code as 'RH1.5' in the first sentence of Section 3.2 but as 'RH1.5D' elsewhere; please use the acronym consistently.","section":"Section 3.2"},{"comment":"The Cheng et al. (2024) references appear in two forms (arXiv:2402.07175 and ApJ 966, L29); please unify the citation entries.","section":"References"},{"comment":"The limitation concerning depletion of coronal material near the upper boundary is acknowledged, but the paper should state whether any of the analyzed events (e.g., the current sheet extending to about 6 Mm in Fig. 3) are affected by this boundary-induced deficit.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is essentially a case study that consolidates the authors' own 2D model in a 3D RMHD context; the main novelty is the 3D demonstration rather than a new physical mechanism. The key risk is the lack of any resolution or dissipation sensitivity test for the plasmoid structures, and I would ask for a strong quantitative response on that point before acceptance. The manuscript is within the scope of A&A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious referee. The paper's main claim—that plasmoid-mediated turbulent reconnection in 3D RMHD explains the coexistence and height inversion of EBs and UV bursts—is plausible but rests on plasmoid identification at roughly 15 km resolution with no convergence check.\n\nWhat is new: this appears to be the first 3D RMHD simulation with realistic radiative transfer to show plasmoid-like structures in small-scale reconnection events in the lower atmosphere. The authors use MURaM, a well-tested code, and synthesize Hα and Si IV with RH1.5D. The synthetic Hα wing profiles show the classic double-wing enhancement, the Si IV profiles are double-peaked, and the spatial arrangement—cool EB-like plasma above hot UV-burst plasma, hot plasma as low as 0.7 Mm—is a concrete, testable pattern. They also flag their own limitations, including the upper boundary and the intermittency of the plasmoid instability. The figures are informative.\n\nSoft spots, in order. First and main: the plasmoid identification is purely morphological. MURaM has no explicit resistivity; magnetic diffusion is numerical hyperdiffusion. With grid spacings of 15–23 km, the current sheets are resolved by only a few cells. The text asserts this resolution is adequate but gives no convergence study, no comparison with a higher-resolution run, and no quantitative plasmoid metric. If the blobs are numerical dissipation artifacts, the claimed causal link to the EB–UV burst emission does not follow. That is a load-bearing missing verification.\n\nSecond: the analysis is a case study. One target event at one time step, plus one additional event, is thin support for the statement that plasmoid instability appears in 'most' small-scale reconnection events. It may be true, but the evidence is not statistical.\n\nThird: the interpretation leans heavily on the authors' prior 2D simulations for both setup and spectral interpretation. That is not circular in the fitting sense, but it does mean the confirmation of their earlier model is partly based on the same assumptions.\n\nMinor: the writing has typos ('comaprison', 'resulted', 'Figsures') and some figures are hard to parse without the movies. These are fixable.\n\nIn sum, the paper is a plausible and honest new simulation result, useful for anyone working on small-scale reconnection or EB–UV burst connections. It deserves peer review. I'd ask the authors for a convergence check or, failing that, a clear statement of what the effective Lundquist number is; a quantitative plasmoid criterion; and some event statistics.","headline":"A plausible first 3D RMHD look at plasmoid-mediated reconnection in the quiet lower atmosphere, but the plasmoid identification has no convergence backing.","tokens_in":21056,"tokens_out":3326,"would_cite":true,"duration_ms":33646,"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":"Plasmoid-mediated reconnection produces Ellerman bombs and UV bursts from the same current sheet.","keywords":["magnetic reconnection","plasmoid instability","solar chromosphere","Ellerman bombs","ultraviolet bursts","flux emergence","radiation magnetohydrodynamics","synthetic spectroscopy"],"falsifier":"Repeat the simulation at significantly higher resolution or with a physical resistivity model and check whether plasmoids, their sizes, and the 20,000 to 90,000 K heating converge; if they do not, the central claim fails. Observationally, a well-resolved event where H-alpha wing and Si IV brightenings are never co-spatial or co-temporal, despite adequate resolution, would contradict the predicted coexistence.","tokens_in":19882,"feed_emoji":"☀️","tokens_out":6642,"duration_ms":63539,"temperature":0.7,"pith_summary":"The paper argues that the small-scale brightenings known as Ellerman bombs and ultraviolet bursts are not separate phenomena originating at different heights, but two views of one turbulent reconnection process. Using a three-dimensional radiation magnetohydrodynamic simulation of magnetic flux emerging into the Sun's lower atmosphere, the authors find thin current sheets that repeatedly break up into plasmoids, ejecting small twisted magnetic flux ropes. The simulated plasma is a patchwork in which material above 20,000 K sits next to material below 10,000 K, and synthetic spectral diagnostics reproduce both H-alpha wing brightenings (Ellerman bombs) and Si IV emission (UV bursts) in the same region. Cool EB-like plasma appears above hot UV-burst-like plasma, reaching altitudes beyond 2 Mm, while hot plasma descends as low as 0.7 Mm. The authors conclude that plasmoid-mediated turbulent reconnection can explain observed cases where UV bursts are temporally and spatially connected to Ellerman bombs.","feed_headline":"Plasmoids tie solar Ellerman bombs and UV bursts to one sheet","feed_subtitle":"A 3D simulation shows cool and hot plasma inside one reconnection region, explaining observed connections.","key_machinery":"The load-bearing mechanism is the plasmoid instability operating in three-dimensional current sheets: a thin, elongated reconnection layer becomes unstable and fragments into plasmoids, which in three dimensions are small twisted magnetic flux ropes ejected by bidirectional outflows. This instability converts a single smooth current sheet into a turbulent, multi-thermal structure where hot and cool plasma alternate in space. The supporting machinery is the simulation's radiation treatment and spectral synthesis, using radiative transfer for H-alpha and optically thin emission for Si IV, which allow the plasma state to be compared with observed Ellerman bomb and UV burst diagnostics.","core_discovery":"The central claim is that in the cool lower solar atmosphere, magnetic reconnection is routinely mediated by the plasmoid instability, and that this turbulence is what produces the coexistence of Ellerman bombs and UV bursts. In a three-dimensional RMHD simulation driven by flux emergence, reconnection forms thin, elongated current sheets in which plasmoid-like structures develop, becoming many small twisted magnetic flux ropes expelled along the sheet. The reconnection region becomes multi-thermal: hot plasma exceeding 20,000 K, with peaks near 90,000 K, and cool plasma below 10,000 K interleave in space. Synthetic H-alpha and Si IV observations made from the simulation show the characteristic signatures of EBs and UV bursts arising from the same current sheet, with cool EB plasma located above hot plasma at heights greater than 2 Mm and hot UV-burst plasma reaching down to about 0.7 Mm. The authors state this is the first time plasmoid instability has been shown to appear in most small-scale reconnection events relating to EBs and UV bursts in a three-dimensional RMHD simulation with radiation.","pith_inferences":["If the plasmoid instability is as widespread as the simulation suggests, some observed 'two-component' events may be projection effects of a single multi-thermal sheet, which could be tested by comparing line-of-sight velocities with the predicted bidirectional outflows.","The same plasmoid-mediated mechanism may apply to other small-scale brightenings and to events whose magnetic topology resembles flare-like configurations at smaller scales.","Because the effective resistivity in the simulation is numerical, a resolution study is the natural next test: current-sheet aspect ratios and plasmoid sizes should follow a power law if the instability is physical.","The simulation's limited coronal plasma above 8 Mm could alter how much reconnection-generated material is transported upward; extending the domain may change the high-altitude signatures."],"forward_implications":["Ellerman bombs and UV bursts can be two emission signatures of a single reconnection event, so joint observations in the two passbands should frequently find them co-located and synchronized.","Cool EB-like plasma can sit above hot UV-emitting plasma, so height alone does not separate the two phenomena; the 2 Mm and 0.7 Mm altitudes bracket the same current sheet.","UV bursts can occur in the lower chromosphere when reconnection fields are strong enough, not only at transition-region heights.","Twisted magnetic flux ropes ejected by the reconnection carry magnetic twist upward and may contribute to coupling the lower atmosphere to the corona.","Most small-scale reconnection in the simulated lower atmosphere is plasmoid-mediated, implying turbulent reconnection is a common heating channel there."],"supporting_citations":[{"why":"The MURaM code that carries out the 3D RMHD simulation; the basic numerical model this work extends with chromospheric and coronal physics.","marker":"Vögler et al. 2005"},{"why":"The specific MURaM version used, including coronal radiation cooling and time-step treatment for high Alfvén speeds.","marker":"Rempel 2017"},{"why":"Supplies the dissipation settings used in the presented simulations.","marker":"Przybylski et al. 2022"},{"why":"The radiative transfer code used to synthesize H-alpha and Si IV spectral line profiles from the simulation.","marker":"Pereira & Uitenbroek 2015"},{"why":"The CHIANTI atomic database used to compute optically thin Si IV 139.4 nm and EUI 17.4 nm synthetic images.","marker":"Dere et al. 1997"},{"why":"The 2.5D model showing plasmoid instability creates alternating hot and cool plasma; this paper validates and extends that model to 3D.","marker":"Ni et al. 2021"},{"why":"The earlier 2D study with radiative-transfer synthesis that established the EB-UV burst connection this work reproduces in 3D.","marker":"Cheng et al. 2024"},{"why":"The prior 3D RMHD work with UV bursts and EBs in separate, distant layers, the comparison case this paper's co-spatial result contrasts with.","marker":"Hansteen et al. 2019"},{"why":"Shows strong Si IV emission can form in the lower chromosphere for fields above about 500 G, supporting the low-altitude UV burst result.","marker":"Ni et al. 2022"}],"fun_headline_variants":["Plasmoid instability ties Ellerman bombs and UV bursts","One reconnection sheet unites cool and hot solar plasma","Plasmoid chaos explains Ellerman bomb–UV burst pairs","3D simulation links plasmoids to solar bomb and burst events"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulation grid spacing of roughly 15 to 23 km is fine enough that the plasmoids and the associated heating are real physics rather than numerical artifacts, even though the reconnection rate and heating depend on the code's numerical dissipation.","fun_headline_variants_meta":{"raw":{"variants":["Plasmoid instability ties Ellerman bombs and UV bursts","One reconnection sheet unites cool and hot solar plasma","Plasmoid chaos explains Ellerman bomb–UV burst pairs","3D simulation links plasmoids to solar bomb and burst events"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1704,"prompt_tokens":1090,"completion_tokens":614,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":545}},"tokens_in":706,"tokens_out":614,"duration_ms":6867,"temperature":1.0,"reasoning_tokens":545,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:28:57.301037+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the simulation at significantly higher resolution or with a physical resistivity model and check whether plasmoids, their sizes, and the 20,000 to 90,000 K heating converge; if they do not, the central claim fails. Observationally, a well-resolved event where H-alpha wing and Si IV brightenings are never co-spatial or co-temporal, despite adequate resolution, would contradict the predicted coexistence.","supporting_citations":[],"review_version":2}