{"id":"49c0b460-fb31-435c-bc9b-f573523e23d5","arxiv_id":"2412.10633","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In 3D simulations, partial ionization makes emerging solar magnetic fields arch-like with little twist and produces faster ejective eruptions than the fully ionized case.","lead":"Using 3D computer simulations, this paper compares how solar magnetic fields rise through the Sun's atmosphere when the plasma is only partially ionized versus fully ionized. It finds that partial ionization changes the shape of the emerging magnetic field and makes the resulting solar eruptions faster, which matters for predicting space weather.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PI vs FI comparison confounds partial ionization with a 2.5x difference in initial field strength; no control run separates the two.","rationale":"The reader's weakest assumption correctly identifies the confound. I agree. The paper's same-beta design is not an adequate control because beta is a ratio; a factor-2.5 change in B0 changes the magnetic energy density by a factor of ~6 and shifts the plasma beta at every atmospheric height, which is known to affect flux emergence dynamics. For example, the PI run's earlier and faster emergence could simply reflect greater magnetic buoyancy from the stronger field. The paper presents no argument that B0 differences are dynamically irrelevant, and the existing diagnostics (field-line topology, height-time profiles, energy traces) cannot distinguish the two causes. A matched-B0 FI control is the minimal experiment that would settle this. I also note a secondary inconsistency in §3.5 where 'larger value of νz in the FI case' is used to explain why FI kinetic energy is lower; the surrounding argument suggests 'smaller' was intended, but this typo does not change the primary concern. The verdict stays CONDITIONAL pending the control run.","tokens_in":15035,"tokens_out":7572,"duration_ms":72447,"concrete_test":"Run a fully ionized (FI) simulation with B0 = 7882 G, keeping the FI EOS and all other parameters (grid, atmosphere, twist, boundary conditions) identical to the original FI run. Compare the resulting emergence structure, flux-rope formation height, and eruption kinematics against the original FI (B0=3150 G) and PI (B0=7882 G) runs. If the high-B0 FI run exhibits arch-like fieldlines, delayed/coronal flux-rope formation, and faster eruptions similar to the PI run, the attributed partial-ionization effects are instead field-strength effects. For a stronger test, run both EOSs at both B0 values in a 2x2 matrix.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central attribution—that partial ionization changes emergence topology, flux-rope formation height, and eruption speed—rests on a single PI/FI pair that is not controlled for initial magnetic field strength. Section 2 states both flux tubes begin with the same plasma beta, but B0_PI = 7882 G while B0_FI = 3150 G. Because the runs differ simultaneously in both the ionization treatment and B0, the observed differences (arch-like vs twisted field, low-atmosphere vs coronal flux-rope formation, faster PI eruptions) could be produced by the stronger initial field in the PI run, which lowers beta at every height and increases magnetic pressure, rather than by partial ionization itself. The paper does not report a run varying B0 independent of the EOS, and its own diagnostics (e.g., normalized B2, kinetic energy) are not sufficient to separate these effects. A control simulation is needed before the claim 'partial ionization causes...' is warranted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents 3D MHD simulations of magnetic flux emergence from the solar interior into an initially unmagnetized corona, comparing a partially ionized (PI) plasma with a fully ionized (FI) plasma. The authors report that, in the PI case, the emerging field consists of arch-like field lines with little twist, a new flux rope forms at coronal heights via reconnection of J-shaped field lines, the first eruption carries neutrals into the corona, and eruptions are overall faster. In the FI case, a flux rope forms in the low atmosphere and eruptions involve denser, more strongly magnetized plasma. The paper continues the authors' previous work (Paper I) and includes diagnostics of field-line topology, height-time profiles, energy and axial-flux evolution, and neutral fractions.","tokens_in":15209,"tokens_out":6321,"duration_ms":54741,"significance":"If the reported differences are robust, this would be the first 3D demonstration that partial ionization qualitatively changes the location of flux-rope formation and the speed of eruptions in flux-emergence simulations, with implications for understanding eruption onset and the role of neutrals in the corona. The paper has several strengths: it presents multiple independent diagnostics (field-line topology, height-time profiles, energy and axial-flux time series, neutral fraction maps), it checks the ambipolar drift speed against published warnings about the EOS (Section 3.7), and it explicitly compares with previous fully ionized simulations. These checks show scientific care. However, the central PI-versus-FI comparison is not cleanly controlled, and the quantitative support for the 'faster eruptions' claim is internally inconsistent. The significance of the work for the field would be high if the control issue were resolved, but as it stands the causal attribution to partial ionization is not yet established.","major_comments":[{"comment":"The PI and FI runs are not controlled for initial magnetic field strength. The manuscript states that both flux tubes begin with the same plasma beta but that B0_PI = 7882 G and B0_FI = 3150 G. Consequently, the two models differ simultaneously in both the ionization treatment and the initial field strength. Because the magnetic field strength determines the Lorentz force, the plasma beta distribution, and the dynamics of emergence and eruption, the observed differences (arch-like vs twisted field, coronal vs low-atmosphere flux-rope formation, faster PI eruptions) cannot be unambiguously attributed to partial ionization. A control run that varies B0 while keeping the EOS fixed, or vice versa, is required to separate these effects. Without such a control, the central causal claim of the paper is not supported by the evidence presented.","section":"Section 2"},{"comment":"There is a clear internal inconsistency regarding vertical velocities and kinetic energy. In Section 3.5 the text says, 'the plasma acceleration during the eruptions in the FI case is less compared to the PI case. This is due to the larger value of νz in the FI case.' A larger νz in the FI case would tend to produce higher, not lower, kinetic energy, so this statement is logically contradictory. The Conclusions then state that 'the PI eruptions ... are faster by a factor of two,' which appears to contradict the claim that FI has the larger νz. The manuscript needs to reconcile these statements by presenting direct measurements of eruption speed (e.g., from the height-time profiles in Figure 4) or by clearly defining what quantity is being compared. As written, the evidence for the central 'faster eruptions in PI' claim is unclear and internally inconsistent.","section":"Section 3.5 and Conclusions"},{"comment":"The claim that PI eruptions are 'faster by a factor of two' is not supported by the diagnostics shown. The kinetic energy in Figure 7 is integrated over a fixed volume at 40 Mm height, but the normalization and units are not specified, and a factor-of-two difference in kinetic energy does not directly imply a factor-of-two difference in velocity unless the density is identical. The height-time profiles (Figure 4) are described qualitatively (slow rise, fast rise, different starting heights) but no quantitative speeds are extracted. The manuscript should either present a direct speed measurement (e.g., fitting the fast-rise phase and reporting the e-folding time or terminal velocity) or moderate the claim to a qualitative statement about the kinetic energy peaks.","section":"Section 3.5 and Conclusions"}],"minor_comments":[{"comment":"The phrase 'the axis of the twisted emerging flux tube in FI rises above the photosphere and it stays within the photosphere' is ambiguous; please clarify whether the axis remains above or within the photosphere.","section":"Section 3.1"},{"comment":"The sentence 'the actual emergence to the corona proceeds a bit earlier in the PI case, since the formation of the new FR that eventually erupts in the corona, it occurs earlier in the FI case' is grammatically confusing and appears to contain a logical error; please rephrase.","section":"Section 3.1"},{"comment":"The normalization of the kinetic energy E_kin(t) is not defined; state whether it is normalized and by what quantity, or provide units.","section":"Section 3.5"},{"comment":"The text uses 'plasma acceleration' but seems to mean kinetic energy or velocity; please correct the terminology.","section":"Section 3.5"},{"comment":"There are several typos, including 'ellaborate' (Section 3.1), 'substancially' (Section 3.2), 'occured' (multiple places), 'procees' (Section 3.6), and 'fiedlines' (Section 3.2). A thorough proofread is recommended.","section":"Section 3.6"},{"comment":"The conclusion states 'the normalized axial flux, which remains below the photosphere is less in the PI case', but earlier text says the opposite (PI has about 40% below, FI just below 30%); please verify and correct.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The central problem is the uncontrolled initial field strength between the PI and FI runs. This is a standard experimental-design issue and should be fixable by adding a control simulation (e.g., a PI run with B0 = 3150 G or an FI run with B0 = 7882 G). If the authors can show that the qualitative differences persist when B0 is held fixed, the paper would be a solid contribution. As it stands, the causal attribution is not yet justified. The internal inconsistency about velocity/kinetic energy also needs careful correction. I recommend major revision rather than rejection because the underlying scientific question is important and the authors have demonstrated the capability to perform the needed control runs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper has a real new 3D result: partial ionization still permits ejective eruptions in flux emergence, and the first flux rope forms by coronal-height reconnection of J-like fieldlines, not low-atmosphere reconnection as in the FI case. That is a useful counterpoint to the 2.5D Leake & Linton (2013) suggestion. The catch is that the PI and FI runs are not controlled for initial field strength: same plasma beta, but B0 = 7882 G in PI and 3150 G in FI. Partial ionization and field strength vary together, so the differences in topology, density, and eruption speed cannot be unambiguously attributed to PI.\n\nWhat is good: the diagnostics are varied—fieldline topology, height-time profiles, energy and axial-flux time series, neutral fraction maps. The authors check the ambipolar drift speed and engage with the more realistic EOS of Nóbrega-Siverio et al., which is a good-faith limitation. The EOS discussion in Sec. 3.7 is honest.\n\nSoft spots: the B0 confound is the biggest. It is not minor: a factor of 2.5 in field strength changes magnetic pressure, buoyancy, and Alfvén speeds at every height. No control run separates PI from B0. There is also one run per case and no resolution study, so numerical convergence is unknown. The kinetic-energy discussion in Sec. 3.5 is muddled: it says acceleration is lower in FI despite larger v_z, then explains higher PI kinetic energy via coronal Alfvén speed. The 'faster by factor of two' conclusion needs clearer support. And the axial-flux plots are normalized to each run's own initial B0, which can hide absolute differences.\n\nNone of this makes the central result implausible. The topology and neutral-fraction evolution are coherent, and the authors are careful not to overclaim the eruption mechanism. But the causal language—partial ionization changes the eruption site and speed—is ahead of the evidence as presented.\n\nThis is a paper for the solar MHD community. It deserves a serious referee, with a request for a B0-matched control run and cleanup of the energy diagnostics. I would not desk-reject.","headline":"Useful first 3D result on PI flux emergence, but the PI/FI comparison is confounded by different initial field strengths, so the central causal claims need a control run.","tokens_in":15756,"tokens_out":4338,"would_cite":true,"duration_ms":39019,"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":"Partially ionized plasma shifts where solar eruptions begin and makes them faster, 3D MHD simulations show.","keywords":["solar eruptions","magnetic flux emergence","partial ionization","ambipolar diffusion","flux rope formation","MHD simulations","solar corona","neutral atoms"],"falsifier":"A matched-field-strength simulation would settle it: rerun the PI case with the same initial field strength as the FI case ($B_0 = 3150\\,\\mathrm{G}$) instead of the same plasma $\\beta$. If the arch-like emergence, J-shaped reconnection at coronal heights, and faster eruptions persist, partial ionization is the cause; if they vanish or invert, the difference in initial field strength is the actual driver.","tokens_in":14815,"feed_emoji":"🌞","tokens_out":8458,"duration_ms":67755,"temperature":0.7,"pith_summary":"This paper argues that partial ionization of the low solar atmosphere changes the entire eruption sequence, not just the details of flux emergence. In the fully ionized comparison run, a new flux rope forms low in the atmosphere by reconnection of sheared field lines and then erupts. In the partially ionized run, the axis of the rising tube stays below the photosphere, the emerged field is an untwisted arch, and the flux rope is built instead by reconnection of J-shaped field lines at coronal heights. The paper also reports that the first partially ionized eruption carries neutral material into the corona and that these eruptions proceed faster than their fully ionized counterparts. If correct, models of solar eruptions must include neutral–ion coupling to predict where eruptions begin and how fast they go.","feed_headline":"Solar eruptions run faster when plasma is partly ionized","feed_subtitle":"3D simulations show the erupting flux rope is born in the corona, not the low atmosphere, and carries neutrals upward.","key_machinery":"The load-bearing machinery is the single-fluid MHD treatment of partial ionization: an ambipolar diffusion term in the induction equation plus a hydrogen equation of state that includes ionization and recombination. Ambipolar diffusion lets neutrals and ions drift relative to each other, which the paper argues changes the stratification of the emerging tube, suppresses the rise of the tube axis, and keeps the emerged field arch-like with very little twist. The comparison is made visible by running the same emergence experiment twice, once with this PI physics and once with a fully ionized plasma; the PI run's J-shaped field-line reconnection at coronal heights is the specific new mechanism that builds the erupting flux rope.","core_discovery":"On the paper's own terms, the discovery is that partial ionization changes the whole pathway from emergence to eruption. In the fully ionized run, the axis of the twisted rising tube emerges above the photosphere, sheared field lines reconnect low in the atmosphere, and the new flux rope is already present before the eruption begins. In the partially ionized run, the tube axis stays below the photosphere, so the emerged field is an untwisted arch; the erupting flux rope forms later, at coronal heights, when J-shaped field lines from the two flanks reconnect at a current sheet around $z \\approx 11.5$ Mm. The paper reports that the first PI eruption lifts a column of neutral material into the corona for a short time, that the normalized axial flux above the photosphere peaks near 0.85 before the first eruption in both runs, and that the PI eruptions have a higher $\\rho v_z^2$ signature and are faster by about a factor of two, even though the FI eruptions carry stronger normalized magnetic energy.","pith_inferences":["Because the same-beta setup also changes the initial field strength (7882 G versus 3150 G), the paper does not fully isolate partial ionization from field-strength effects; a matched-field-strength pair of runs would be the clean test.","The coronal J-shaped reconnection mechanism may operate in any emergence where the flux-tube axis fails to reach the photosphere, so variations in stratification or helicity could reproduce the same eruption pathway without invoking partial ionization.","If neutral material reaches coronal heights in real eruptions, spectral diagnostics of erupting prominences could be used to identify partially ionized flux ropes observationally."],"forward_implications":["If partial ionization keeps the tube axis below the photosphere, the coronal field that emerges is an arcade with very little twist, so twist-driven instabilities are not the trigger that creates the erupting flux rope.","A new flux rope can be assembled at coronal heights by reconnection of J-shaped field lines even when no twisted tube axis emerges, so the formation site of eruptive structures depends on the ionization state of the lower atmosphere.","Neutral atoms can be carried upward with the erupting flux rope for a short time, meaning partially ionized material can appear at coronal heights during an eruption.","Eruptions in the partially ionized run are faster, with a kinetic-energy proxy roughly a factor of two higher, because the reconnection outflows that drive them form where the Alfvén speed is higher.","The amount of axial flux above the photosphere reaches similar peaks in both runs, so the subsequent eruptions carry comparable axial flux even though the subphotospheric flux reservoirs differ."],"supporting_citations":[{"why":"Supplies the numerical setup and the earlier result that the flux-tube axis stays below the photosphere in the partially ionized case.","marker":"Paper I (Chouliaras et al. 2023)"},{"why":"A 2.5D partial-ionization emergence experiment whose subphotospheric-axis result the PI run reproduces.","marker":"Leake & Linton (2013)"},{"why":"Fully ionized emergence experiment that is the template for the FI run and for recurrent eruptions.","marker":"Syntelis et al. (2017)"},{"why":"Establishes the fully ionized mechanism of flux-rope formation by reconnection of sheared field lines.","marker":"Archontis & Török (2008)"},{"why":"Provides the equation of state used for the partially ionized plasma.","marker":"Leake et al. (2013)"},{"why":"Gives the ion-neutral drift velocity formula used to check that the ambipolar diffusion approximation is valid.","marker":"Nóbrega-Siverio et al. (2020b)"},{"why":"Provides the normalized axial flux comparison that the paper uses to contextualize its flux results.","marker":"Moreno-Insertis & Galsgaard (2013)"}],"fun_headline_variants":["Solar eruptions faster with partial ionization","Partial ionization speeds solar eruptions and changes flux rope origin","Neutrals ride erupting flux rope in partially ionized solar plasma","Untwisted arches emerge when solar plasma is partly ionized","Flux rope born in corona in partial ionization solar runs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison assumes that matching the plasma beta (the ratio of gas pressure to magnetic pressure) between the two runs is a fair control, even though that choice forces the partially ionized tube to start with a magnetic field of 7882 G and the fully ionized tube with 3150 G.","fun_headline_variants_meta":{"raw":{"variants":["Solar eruptions faster with partial ionization","Partial ionization speeds solar eruptions and changes flux rope origin","Neutrals ride erupting flux rope in partially ionized solar plasma","Untwisted arches emerge when solar plasma is partly ionized","Flux rope born in corona in partial ionization solar runs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1473,"prompt_tokens":1006,"completion_tokens":467,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":388}},"tokens_in":622,"tokens_out":467,"duration_ms":4994,"temperature":1.0,"reasoning_tokens":388,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:45:52.531063+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A matched-field-strength simulation would settle it: rerun the PI case with the same initial field strength as the FI case ($B_0 = 3150\\,\\mathrm{G}$) instead of the same plasma $\\beta$. If the arch-like emergence, J-shaped reconnection at coronal heights, and faster eruptions persist, partial ionization is the cause; if they vanish or invert, the difference in initial field strength is the actual driver.","supporting_citations":[{"cited_title":"2017, ApJ, 850, 95, doi: 10.3847/1538-4357/aa9612","cited_arxiv_id":null,"evidence_quote":"Fully ionized emergence experiment that is the template for the FI run and for recurrent eruptions."}],"review_version":1}