{"id":"46bd6e26-873a-4a4c-8783-d81719947cf0","arxiv_id":"2412.10397","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A boundary-condition adjustment that caps Alfvén speed and raises minimum plasma beta removes unphysical fast streams in COCONUT coronal simulations while preserving magnetic topology.","lead":"This paper changes how a global solar corona model sets its inner boundary, capping the Alfvén speed and raising the minimum plasma beta in strong-field regions. The change removes artificially fast, extremely hot outflow streams while leaving the magnetic field structure mostly unchanged.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Constraining the boundary with VA,max and beta_min removes the artifact, but the choice of constraint values and the claim that topology is unaffected need a quantitative check against the unconstrained solution.","rationale":"The reader's weakest assumption correctly identifies the transferability of Bifrost/literature ranges to strong-field active regions and the lack of an observational baseline for the removed stream. My proposed test directly checks whether the topology claim holds quantitatively and whether the artifact removal is sensitive to the specific constraint values. The paper itself acknowledges large uncertainties in the derived ranges (Section 2.2), and the 'not significantly affected' claim is only visually supported. However, these are addressable concerns, not fatal flaws; the CONDITIONAL verdict already captures this. The explicit self-critique in the paper ('large uncertainties may exist in them', Section 2.2) supports my concern. I do not see an internal inconsistency that would make the method unsound; the issue is external validation, which is the standard for a CONDITIONAL verdict.","tokens_in":9679,"tokens_out":1529,"duration_ms":12308,"concrete_test":"Quantitatively compare the magnetic topology between the unconstrained and constrained (VA,max = 1 Mm/s, beta_min = 0.01) solutions by computing a field-line connectivity metric, e.g., the fraction of field lines that map to the same footpoint polarity or the difference in the squashing factor (Q) over the boundary region. If the Q maps differ by less than 10% in the affected region, the topology claim holds; if they differ substantially, the claim overstates the impact. Additionally, run the constrained case with VA,max = 2 Mm/s and beta_min = 0.003 (the Bifrost-derived values) rather than the tuned values and check whether the stream still appears; if it does, the artifact removal depends on the chosen constraint values, not on physically justified limits.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the high-speed stream is a numerical artifact of the homogeneous boundary condition and that removing it via VA,max and beta_min constraints leaves the magnetic topology essentially unchanged. The load-bearing concern is that the constraint values are not independently justified for the strong-field active region at hand: the Bifrost simulation (ch024031_by200bz005) provides a maximum VA of about 2 Mm/s and minimum beta of about 0.003, but that simulation is not an active-region case; the paper applies these global ranges to a solar-maximum active region with |B| up to 50-100 G. If the true coronal conditions above that active region permit higher VA or lower beta, then the method is removing a real feature rather than an artifact. The paper declares the stream 'inexplicable' and 'unphysical' because it is hotter than 30 MK, but no observational baseline is given for temperatures in such streams. Also, the claim that topology is 'not significantly affected' is supported only by visual comparison of field lines (Figure 6) and a color plot; there is no quantitative metric for connectivity or field-line shape. If the constraints are chosen to exactly remove the artifact, the method is circular without an independent check that the constrained solution matches observations. Section 3 and Table 1 show that VA,max and beta_min together set the temperature at the boundary, so the selected values are partly tuned to produce a desired temperature. The weakest assumption is thus that the Bifrost/literature ranges transfer to strong-field active regions on the global Sun and that the removed stream is not a real coronal feature.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a boundary-condition adjustment for the COCONUT global coronal MHD model. Motivated by Bifrost simulations and literature values, the authors constrain the plasma beta and Alfvén speed at the inner boundary by modifying the prescribed pressure and density through smooth tanh-transition functions (Eqs. 1-4). They apply the method to the March 9, 2016 eclipse case (CR2174) and show that high-speed streams (speeds above 1 Mm/s, temperatures above 30 MK) are reduced or removed while the magnetic field lines are reported to remain broadly similar. The paper concludes that this technique is a practical, low-cost way to remove boundary-condition artifacts during solar-maximum conditions.","tokens_in":9962,"tokens_out":5489,"duration_ms":50113,"significance":"If the central claim is supported, the method is a useful and inexpensive fix for a known boundary-condition problem in global coronal MHD modeling: it uses physically motivated proxy constraints, is simple to implement, and avoids the computational cost of adding a resolved transition region. The use of independent Bifrost and literature ranges is a positive feature, as is the explicit smooth transition that avoids convergence issues. However, the demonstration currently rests on a single event, the artifact diagnosis is based on thresholds that are not validated against observations, the topology claim is qualitative, and part of the parameter selection is tuned to suppress the very feature that the paper wants to classify as an artifact. These issues leave the central claim plausible but not fully established.","major_comments":[{"comment":"The decision to label the stream as an artifact rests entirely on the statements in Section 3 that speeds above 1 Mm/s and temperatures above 30 MK are 'unexpected' and 'much hotter than realistic,' but no observational baseline or comparison is provided for these thresholds. In particular, the paper does not compare the constrained and unconstrained solutions against EUV or white-light observations of the 2016 eclipse or against any independent coronal diagnostic. Without such a check, the method could be removing a feature that is real or the thresholds could be stricter than necessary, so the central claim that the fix removes a numerical artifact is not yet established.","section":"Section 3, Figures 2-3"},{"comment":"The Bifrost run ch024031_by200bz005 is used to set VA,max ~2 Mm/s and beta_min ~0.003, but the manuscript does not establish that this simulation is representative of the strong-field active-region conditions (|B| ~50-100 G) present in the CR2174 magnetogram. Because these ranges justify the chosen constraint values and are then applied to an active region, the authors should either demonstrate that the Bifrost distribution covers that regime (e.g., plot beta and VA as a function of |Bz| for active-region pixels) or supplement with active-region-specific estimates. As written, the applicability of the Bifrost-derived limits to the target configuration is an assumption, not a demonstrated fact.","section":"Section 2.2, Figure 1"},{"comment":"The claim that the magnetic topology is 'not significantly affected' is supported only by visual inspection of field lines and a qualitative statement that the connectivity 'remain[s] very similar.' This is load-bearing because the practical value of the method is precisely that it removes the artifact without changing the coronal magnetic field. The authors should add a quantitative measure, such as footpoint displacements in degrees, a connectivity matrix overlap, or an open/closed flux area comparison between the unconstrained and constrained runs.","section":"Section 3, Figure 6"},{"comment":"The parameter scans and Table 1 show that VA,max and beta_min together determine the boundary temperature and that the 'most effective' values (beta_min between 1e-3 and 1e-2, VA,max ~1e6 m/s) are selected by their ability to suppress the artifact. This introduces a circular step: the constraints are tuned to remove the very feature whose physical status is the premise of the paper. The authors should break the circularity by applying the method to a second Carrington rotation or to a different magnetogram product without retuning, and by reporting the resulting speed and temperature to show that they emerge from the independently justified constraints rather than from the tuning.","section":"Section 3, Table 1"}],"minor_comments":[{"comment":"The transition width dtan is set to 10% of beta_min or VA,max based on numerical experiments; please state the range of dtan tested and whether the results are sensitive to this choice.","section":"Section 2.3, Eqs. (2) and (4)"},{"comment":"The color bars and panel labels in Figures 2-5 are very small and hard to read; please make them legible and explicitly state that the inner-surface colors in Figures 4 and 5 show prescribed boundary values rather than solution values.","section":"Figures 2-5"},{"comment":"The manuscript lacks a data and code availability statement; at minimum, the authors should state whether the COCONUT input files and the Bifrost data products used to produce Figures 1-6 are accessible to other researchers.","section":"General"},{"comment":"Section 2.2 reports a Bifrost maximum VA of about 2 Mm/s, but Section 3 uses VA,max = 1 Mm/s as the target with the phrase 'roughly corresponding'; one sentence explaining why 1 Mm/s rather than 2 Mm/s is chosen would clarify the comparison.","section":"Section 2.2 vs. Section 3"},{"comment":"The terms 'inexplicable' and 'unexpected' are used interchangeably for the fast stream; please define what is meant by these terms, for example whether the stream is inconsistent with observations, with the prescribed magnetic field, or with a reference model.","section":"Abstract and Section 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a plausible methods note and within the scope of the journal, but the central claim is currently under-supported by a single event and by qualitative criteria. I would encourage a revision that adds at least one more Carrington rotation and quantitative topology/validation metrics; without those, the paper reads more like a technical report than a fully supported A&A methods paper. No concerns about citation practice or novelty beyond the need for stronger validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Michaela and colleagues have done a clean, honest piece of engineering here: the inner boundary of COCONUT is given a smooth, magnetic-field-dependent pressure and density prescription via tanh transitions on plasma beta and Alfven speed. The method is simple, clearly stated, and directly addresses a known deficiency in the homogeneous boundary condition that leads to spurious fast streams over strong-field regions. The use of Bifrost and published ranges to set physically plausible limits is a reasonable choice, and the paper correctly notes the uncertainty in these ranges. The demonstration on the 2016 eclipse case shows the artifact can be suppressed while the large-scale field lines appear mostly unchanged; that is genuinely useful for operational space weather modeling.\n\nThe soft spots are real, but not disqualifying. First, the validation rests on one event, and the removal of the stream is judged visually and by temperature, not by comparison with any observational proxy for the plasma in that region. The assertion that the stream is 'unphysical' because it exceeds 30 MK is plausible but lacks a quantitative baseline. Second, the claim that magnetic topology is 'not significantly affected' is based on visual inspection of field lines; a quantitative measure (e.g., open/closed field fraction, footpoint connectivity) would make it much stronger. Third, the most effective values of beta_min and V_A,max are selected by their success at removing the artifact, so part of the demonstration is circular. The authors acknowledge this implicitly by testing a sweep, and the chosen values sit within the observationally motivated range, but an independent check against a second Carrington rotation or against observed streamer boundaries would reduce the circularity.\n\nNone of these issues undermines the central mechanical claim: the constraint formulation works as described, and it is a low-cost remedy. The paper is not the last word on the physics of inner boundary conditions, but it is a useful step for global coronal models used in space weather forecasting. I would send it to a competent referee; my own verdict would be minor revision, mostly asking for a quantitative topology metric and at least one additional test case.","headline":"A clean, practical fix for COCONUT's inner boundary artifact, but the validation is a single qualitative case study and the constraint values are partly chosen by their ability to suppress the artifact.","tokens_in":10541,"tokens_out":2361,"would_cite":true,"duration_ms":20345,"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":"Imposing smooth observation-based caps on boundary Alfvén speed and plasma beta removes inexplicable high-speed streams from the COCONUT global corona model without changing its magnetic topology.","keywords":["coronal magnetohydrodynamics","inner boundary condition","plasma beta","Alfvén speed","global coronal model","space weather forecasting","solar corona","COCONUT"],"falsifier":"An observational test would measure the actual plasma beta and Alfvén speed in the low corona over a strong active region at about 10 Mm, for example through radio or MHD-seismology loop oscillations or high-resolution spectropolarimetry, and then compare a COCONUT run with unconstrained and constrained boundaries. If real observations showed plasma beta below about 0.003 or Alfvén speeds above roughly 2 Mm/s at that height in the capped regions, or if the stream identified as an artifact were detected in coronal imaging or spectroscopy, the central claim would fail.","tokens_in":9494,"feed_emoji":"☀️","tokens_out":6357,"duration_ms":50159,"temperature":0.7,"pith_summary":"Global coronal magnetohydrodynamic models need thermodynamic conditions at their inner boundary, but direct observations there are sparse, so codes typically prescribe uniform pressure and density. This paper argues that in regions of strong magnetic field such uniform values produce plasma beta far below and Alfvén speeds far above what the low corona actually contains, and that these extremes drive spurious high-speed streams (above 1 Mm/s, over 30 MK) in the solution. The paper proposes replacing the uniform prescription, in the affected regions only, with a smooth constraint: pressure is raised so that plasma beta never falls below a chosen minimum, and density is raised so that Alfvén speed never exceeds a chosen maximum. In a 2016 solar-eclipse simulation for an active-region-rich Carrington rotation, this removes the inexplicable streams while leaving the magnetic topology essentially unchanged. The payoff is a low-cost boundary-condition fix that could make solar-maximum coronal modelling more reliable without global grid changes.","feed_headline":"Capping Alfvén speed and beta removes fake solar streams","feed_subtitle":"A smooth floor on plasma beta at the inner boundary erases 30 MK artifact streams without changing magnetic field lines.","key_machinery":"The carrying object is the double-sided hyperbolic tangent transition profile $\\zeta_{\\text{tan}} = \\tfrac12 + \\tfrac12\\tanh(\\pi\\Delta_{\\text{tan}})$, used twice. For pressure, $\\Delta_{\\text{tan}} = (\\beta_{\\min} - \\beta)/d_{\\text{tan}}$ so that where $\\beta$ falls below the floor the prescribed pressure is pushed toward $p_{\\text{mag}}\\beta_{\\min}$; for density, $\\Delta_{\\text{tan}} = (V_A - V_{A,\\max})/d_{\\text{tan}}$ so that where Alfvén speed exceeds the ceiling the prescribed density is pushed toward $|\\mathbf{B}|^2/(V_{A,\\max}^2 \\mu_0)$. The ghost-cell state is set as $p_g = 2p' - p_i$ and $\\rho_g = 2\\rho' - \\rho_i$ so the constrained values hold exactly on the boundary, and the transition width $d_{\\text{tan}}$ is chosen about 10% of the limit so convergence is preserved while the target values are reached within 0.2% at 10% away from the limit.","core_discovery":"On the paper's own terms, the central discovery is that the inner-boundary artifact is correctable by constraining two derived plasma parameters rather than trying to prescribe pressure and density from sparse observations. Using literature values and a high-resolution solar-atmosphere simulation, the paper sets observationally motivated bounds: maximum Alfvén speed on the order of 2 Mm/s and minimum plasma $\\beta$ around 0.003 at about 10 Mm, with limits chosen per experiment. The boundary state is then adjusted only where the default prescription violates these bounds, via formulas $p' = \\zeta_{\\text{tan}} p_{\\text{mag}}\\beta_{\\min} + p_0(1-\\zeta_{\\text{tan}})$ and $\\rho' = \\zeta_{\\text{tan}} |\\mathbf{B}|^2/(V_{A,\\max}^2 \\mu_0) + \\rho_0(1-\\zeta_{\\text{tan}})$, with a double-sided hyperbolic tangent transition factor that avoids discontinuous derivatives. In the CR2174 March 2016 eclipse test, capping $V_{A,\\max}$ at $10^6$ m/s and setting $\\beta_{\\min}$ between $10^{-3}$ and $10^{-2}$ removes the high-speed stream almost entirely, while setting $\\beta_{\\min}$ too high ($2\\times10^{-2}$, $5\\times10^{-2}$) deforms the flow and field lines. The combination of $V_{A,\\max}$ and $\\beta_{\\min}$ also gives a knob for the resulting boundary temperature, since $T \\propto p/\\rho$.","pith_inferences":["A natural stress test would be to run the same constrained boundary on multiple Carrington rotations and check whether the removed streams disappear consistently rather than moving elsewhere; the paper demonstrates one eclipse case only.","Because the constraints act as proxies for pressure and density, their success suggests that simple empirical maps of beta and Alfvén speed derived from magnetogram-based proxies might eventually replace homogeneous inner-boundary prescriptions in operational heliospheric forecast chains.","The same smooth-constraint machinery could in principle transplant to other global coronal MHD codes that share the homogeneous-boundary limitation, with the transition layer width tuned per solver.","If the technique becomes standard, then boundary-condition-induced artifacts may no longer be mistaken for coronal physics in solar-maximum simulations, which changes how model-data comparisons are interpreted."],"forward_implications":["Capping boundary Alfvén speed at about $10^6$ m/s removes the radial width of the inexplicable stream, and adding $\\beta_{\\min}$ in the range $10^{-3}$ to $10^{-2}$ removes it almost entirely.","For $\\beta_{\\min}$ values above about $2\\times10^{-2}$, the boundary forcing becomes nonphysical and deforms the flow and magnetic field lines.","Magnetic connectivity and overall field-line structure above constrained active regions stay essentially unchanged at physically justifiable constraint levels.","$V_{A,\\max}$ and $\\beta_{\\min}$ can be chosen jointly to target a specific boundary temperature, since temperature scales as $p/\\rho$.","The technique is a partial remedy until high-resolution observations from Solar Orbiter and Parker Solar Probe allow more direct coronal boundary prescription."],"supporting_citations":[{"why":"Supplies the COCONUT solver and its default uniform boundary conditions for pressure and density.","marker":"Perri et al. (2022)"},{"why":"Provides the expected plasma-beta range at about 10 Mm used to justify the beta floor.","marker":"Gary (2001)"},{"why":"Provides MHD-seismology estimates of coronal-loop Alfvén speeds used to justify the speed ceiling.","marker":"Anfinogentov & Nakariakov (2019)"},{"why":"Supplies the high-resolution solar-atmosphere solver whose simulation data set the resolved ranges of beta and Alfvén speed.","marker":"Gudiksen et al. (2011)"},{"why":"Supplies the specific simulation data probed at the inner-boundary height to extract maximum Alfvén speed and minimum plasma beta.","marker":"Finley, A. J. et al. (2022)"},{"why":"Demonstrates that nonphysical thermodynamic boundary conditions in COCONUT can create unexpected streams.","marker":"Kuzma et al. (2023)"},{"why":"Discusses the inaccuracy of homogeneous boundary conditions and their effect on COCONUT solutions.","marker":"Brchnelova et al. (2023)"},{"why":"Provides the HMI magnetogram data used to prescribe the magnetic field for the 2016 eclipse simulation.","marker":"Scherrer et al. (2012)"}],"fun_headline_variants":["Alfvén cap and beta floor erase phantom coronal streams","Boundary tweak kills spurious solar wind features","Plasma parameter limits fix coronal model artifacts","Removing fake streams: cap Alfvén, floor beta","Coronal boundary fix: constrain beta and Alfvén speed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the removed high-speed stream is a numerical artifact of the uniform boundary condition rather than a real coronal feature, and that the Alfvén-speed and plasma-beta ranges taken from existing atmosphere simulations and literature (about 2 Mm/s and 0.003 at about 10 Mm) are representative of active-region low-corona conditions on the global Sun.","fun_headline_variants_meta":{"raw":{"variants":["Alfvén cap and beta floor erase phantom coronal streams","Boundary tweak kills spurious solar wind features","Plasma parameter limits fix coronal model artifacts","Removing fake streams: cap Alfvén, floor beta","Coronal boundary fix: constrain beta and Alfvén speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000321,"raw_usage":{"total_tokens":1926,"prompt_tokens":1184,"completion_tokens":742,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":800,"completion_tokens_details":{"reasoning_tokens":661}},"tokens_in":800,"tokens_out":742,"duration_ms":6838,"temperature":1.0,"reasoning_tokens":661,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:17:36.814408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An observational test would measure the actual plasma beta and Alfvén speed in the low corona over a strong active region at about 10 Mm, for example through radio or MHD-seismology loop oscillations or high-resolution spectropolarimetry, and then compare a COCONUT run with unconstrained and constrained boundaries. If real observations showed plasma beta below about 0.003 or Alfvén speeds above roughly 2 Mm/s at that height in the capped regions, or if the stream identified as an artifact were detected in coronal imaging or spectroscopy, the central claim would fail.","supporting_citations":[{"cited_title":"2022, ApJ, 936","cited_arxiv_id":null,"evidence_quote":"Supplies the COCONUT solver and its default uniform boundary conditions for pressure and density."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the expected plasma-beta range at about 10 Mm used to justify the beta floor."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides MHD-seismology estimates of coronal-loop Alfvén speeds used to justify the speed ceiling."},{"cited_title":"2023, A&A, 676, A83","cited_arxiv_id":null,"evidence_quote":"Discusses the inaccuracy of homogeneous boundary conditions and their effect on COCONUT solutions."}],"review_version":1}