{"id":"55d18fab-868a-425a-a11e-599c043304d7","arxiv_id":"2604.01028","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A non-spherical source surface, defined as an isosurface of magnetic field strength, is used in a finite-element potential field extrapolation to produce more open solar magnetic flux and better match PSP IMF observations.","lead":"This paper introduces a new model of the Sun's outer magnetic field that replaces the usual spherical boundary with a wavy, non-spherical surface shaped by the field strength itself, which lets more magnetic field lines open into space. It matters because the standard model has long underestimated the amount of open magnetic flux, and this new model claims to fix that while still matching observed coronal shapes and spacecraft data.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-step NSSS extraction is not self-consistent; the open-flux increase may be an artifact of the initial PFSS guess.","rationale":"The reader's weakest assumption identifies the same load-bearing premise: the validity of an isosurface of |B| from an initial PFSS solution as a physical source surface after recalculating the field. The paper's own §2.2.3 admits the iterative self-consistency is not demonstrated, and the optimal configuration is selected by matching observed open flux, which weakens the claim of prediction. My concern sharpens this by pointing to the specific consequence: the NSSS dipping to 1.05 R_sun, a height at which plasma beta is typically low, so the assumed opening mechanism is physically questionable. This does not invalidate the model as a useful constructive tool, especially since the code is publicly available and the FEM solver is a genuine technical contribution. However, the abstract's strong claim that the model 'successfully reproduces' observed topologies, IMF properties, and source mappings is premature without demonstrating that the NSSS is at least approximately self-consistent or that the results are robust to the one-shot extraction. The reader's CONDITIONAL verdict appropriately captures this; my analysis does not move it.","tokens_in":14982,"tokens_out":5864,"duration_ms":55688,"concrete_test":"Run one additional iteration of the NSSS extraction for R_ini=2.2: take the NSPF solution, extract the |B| isosurface at the same isovalue used originally, recompute the potential field with that new NSSS, and compare open flux, loop-height distribution, and PSP metrics to the single-step result. If the second iteration yields substantially different open flux or a shrinking NSSS (as the authors report), the single-step NSSS is an artifact of the initial guess, not a converged solution. Additionally, report the relative variation (std/mean) of |B| over the final NSSS; a large variation (e.g., >10%) would quantitatively confirm that the final field is not an isosurface.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The NSPF model's central innovation is the NSSS defined as an isosurface of |B| from an initial PFSS solution (§2.2.2), followed by a potential-field solve with u=0 on that surface (§2.2.3). The headline claim—more open flux while preserving realistic loop heights—depends on this NSSS being a valid physical source surface. However, the authors explicitly concede in §2.2.3 that they cannot demonstrate convergence to a fully self-consistent NSSS that is simultaneously an isosurface of both |B| and the magnetic potential, and that 'the NSSS tends to shrink during the iterative process.' Because the surface is extracted only once, the final field is generally not one in which the NSSS is an isosurface of |B|; the boundary is inherited from a different solution. The open-flux enhancement (13.01 Gs·R_sun^2 for R_ini=2.2) is achieved by allowing the NSSS to dip to 1.05 R_sun (Table 1), far below typical streamer cusps and the expected Alfvén surface, where plasma beta is usually <<1; the physical mechanism for opening field at such heights is not established. Additionally, the optimal R_ini=2.2 is selected by matching observed open flux and PSP IMF (§2.2.3, §3.3), so part of the agreement is imposed rather than predicted. Therefore the key result is conditional on an unvalidated, non-self-consistent surface definition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the Non-Spherical Potential Field (NSPF) model, which replaces the spherical source surface of the classic PFSS model with a non-spherical source surface (NSSS) defined as an isosurface of |B| extracted from an initial PFSS solution. The potential field is then recomputed with u=0 on the NSSS, and the field is extended through a Schatten-type current-sheet layer to 10 R_sun and then as a Parker spiral. For Carrington Rotation 2282, the authors compare NSPF fields (for initial radii 2.2, 2.5, 3.0 R_sun) with PFSS+PFCS models, using EUV and white-light images, total open flux, and PSP in-situ IMF measurements. They report that the NSPF model with R_ini=2.2 R_sun yields an open flux of 13.01 Gs·R_sun^2 (close to the observed value), produces more complex open-field regions consistent with EUV observations, and better predicts IMF polarity reversals than PFSS+PFCS with the same flux. The paper emphasizes that the NSSS dips to 1.05 R_sun, preserving loop heights up to ~1.08 R_sun, thereby solving the open-flux problem without uniformly lowering the source surface.","tokens_in":15369,"tokens_out":5078,"duration_ms":49410,"significance":"If the approach is valid, the NSPF model offers a promising, practical alternative to the common practice of lowering the PFSS source surface to fix the open-flux problem, and it could improve solar wind source mapping. The paper is honest about its main limitation—the lack of demonstrated self-consistency of the NSSS—and provides reproducible open-source code (FEM-PFS based on FEniCS/DOLFINx), which is a strength. The quantitative open-flux comparison in Table 1 is useful, and the idea of extracting a non-spherical surface from an isosurface of |B| is a concrete step forward. However, the central claims that the model 'successfully reproduces' coronal topology, IMF properties, and solar wind source regions are weakened by two load-bearing issues: (i) the NSSS is not shown to be a self-consistent source surface for the recomputed field, and (ii) the optimal configuration is selected using the same observed open flux and PSP data that are later used for validation. The physical mechanism for opening field at R_SS,min=1.05 R_sun is not quantitatively supported. These issues may be addressable in a revision, but they are central to the paper's headline conclusions.","major_comments":[{"comment":"The NSSS is extracted from an initial PFSS solution and then used as a Dirichlet boundary (u=0) for a new potential-field solve. The final field is not guaranteed to have |B| nearly constant on the NSSS; the authors explicitly state they cannot demonstrate convergence to a self-consistent NSSS and that the surface 'tends to shrink during the iterative process.' This is load-bearing because the physical justification in §2.2.2 assumes the NSSS is simultaneously a |B| isosurface and a zero-potential open-field boundary. If the final NSSS is not an isosurface of the recomputed field, the increased open flux (Table 1) may be an artifact of imposing an irregular boundary rather than a physically self-consistent representation of the source surface. Please quantify the deviation (e.g., a histogram or map of |B| on the final NSSS), test at least one additional iteration step, and report how the","section":"§2.2.3"},{"comment":"The optimal NSPF configuration (R_ini=2.2 R_sun) is selected by matching the observed open flux and PSP IMF measurements: the text states 'varying the initial spherical source-surface radius, and then use observational constraints to select the optimal NSSS configuration.' The resulting A=1.0 metric in Fig. 5a is therefore a fit to the same data that are used to claim successful reproduction of IMF properties. This circularity is central to the paper's validation. Please provide an out-of-sample test—for example, fix R_ini based on a different Carrington rotation or use independent constraints (such as coronal hole areas or white-light streamer positions) before comparing to PSP—or explicitly reframe the R_ini dependence as calibration rather than prediction. Without this, the headline agreement is not a falsifiable prediction.","section":"§2.2.3 and §3.3"},{"comment":"The NSSS reaches a minimum radius of 1.05 R_sun in the preferred configuration. At these heights, the coronal plasma beta and dynamic pressure are generally expected to be much less than unity, and the Alfvén surface is typically well above 1.05 R_sun. The physical argument in §2.2.2—that field lines open where ram pressure exceeds magnetic tension—is not quantified for this regime, and the NSSS is defined via |B|, not via beta or the stated pressure balance. The claim that the NSPF model solves the open-flux problem 'physically' requires evidence that the field can indeed open at such low heights (e.g., comparison with MHD simulations, observed streamer cusp heights, or a quantitative beta/ram-pressure analysis along the NSSS). Please add such a check or temper the physical interpretation accordingly.","section":"§3.2 and Table 1"}],"minor_comments":[{"comment":"The caption reads 'interstellar magnetic field'; this should be 'interplanetary magnetic field.'","section":"Fig. 2 caption"},{"comment":"The reference 'Neukrich, T.' appears to be a typo for 'Neukirch, T.' (in Zhu et al. 2022). Please check all author names.","section":"References"},{"comment":"The comparison to EUV and white-light images is qualitative ('provide the best match to the observed coronal rays'). For a model that claims to 'successfully reproduce' coronal topology, a quantitative metric (e.g., overlap of open-field regions with coronal holes, location of streamer rays) would strengthen the conclusion. As written, the topology comparison is illustrative rather than definitive.","section":"§3.1"},{"comment":"The linear form L(v) is written with 'brvds' in the text; this appears to be a LaTeX rendering issue for b_r v ds. Please ensure the symbol is defined (b_n is used in Eq. 1).","section":"Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The paper is an honest and technically interesting contribution, and the open-source FEM solver is a useful community resource. However, the two load-bearing concerns—lack of demonstrated self-consistency of the NSSS and circular selection/validation with PSP data—prevent me from recommending acceptance in its current form. A revision that quantifies the |B| variation on the final NSSS, tests iterative convergence, and provides an out-of-sample validation would substantially strengthen the claims. The qualitative topology comparison could also be quantified. If the authors can demonstrate that the open-flux enhancement is robust to the single-step extraction and not purely a boundary artifact, the paper could become a valuable addition to the coronal modeling literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a genuine step forward in making non-spherical source surfaces tractable. The FEM implementation is clean, the code is public, and the authors are unusually candid about where the method falls short. The headline result — more open flux than PFSS while keeping loop heights realistic — is worth taking seriously. The comparison with EUV and LASCO is qualitative but suggestive, and the source-surface extraction as an |B| isosurface is a practical, clever proxy.\n\nThe soft spots are exactly where the reader puts them. The surface is extracted once from an initial PFSS solution, and after solving Laplace with u=0 on that surface, the final field is not generally one where that surface is an isosurface of |B|. The authors say so explicitly in §2.2.3, and also note the NSSS shrinks under iteration. That means the model is not self-consistent, and the open-flux increase is partly an artifact of where the initial surface happens to sit. The dip to 1.05 R_sun is far below typical streamer cusps; the physical mechanism for opening field there is not established. Also, the choice of R_ini=2.2 is made by matching the observed open flux and PSP IMF, so A≈1.0 is in part a fit, not a prediction. One Carrington rotation and no error bars makes the validation thin, and Kruse et al. (2020) is never quantitatively benchmarked despite being the closest prior work.\n\nNone of this kills the paper. The authors frame the model as a first practical implementation, not a fully converged theory, and they flag the limitations themselves. The fact that they get a plausible topology and a reasonable flux budget without artificially lowering the whole source surface is more than prior conceptual work achieved. The method is reproducible and the honest treatment of the self-consistency problem is a point in its favor.\n\nThis deserves a serious referee. A good referee would push for iterative convergence tests (even if they fail, report how), a multi-rotation validation with magnetograms from cycle minimum and maximum, quantitative comparison against Kruse et al., and a clearer separation between fitted parameters and genuinely predicted outputs. If those are addressed, the model could be widely used as a practical alternative to PFSS for solar wind connectivity studies.","headline":"A practical but not yet self-consistent non-spherical source-surface model; the open-flux match is partly fitted, partly real.","tokens_in":15892,"tokens_out":2477,"would_cite":true,"duration_ms":24241,"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":"This paper argues that swapping the standard spherical source surface for a non-spherical surface of constant magnetic-field strength restores the Sun's missing open magnetic flux without shrinking coronal loops.","keywords":["coronal magnetic field","open flux problem","source surface","potential field extrapolation","solar wind mapping","heliospheric magnetic field","helmet streamers","non-spherical source surface"],"falsifier":"Run the NSSS extraction on the magnetic field of an MHD coronal simulation where the true open/closed boundary is known, then test whether the constant-|B| surface coincides with that boundary and reproduces the observed open flux without choosing an initial radius; if the match fails, the geometrical mechanism is not the cause of the extra flux.","tokens_in":14833,"feed_emoji":"☀️","tokens_out":8313,"duration_ms":72073,"temperature":0.7,"pith_summary":"The paper develops a coronal magnetic-field extrapolation in which the surface where field lines are declared open — the source surface — is no longer a sphere but a surface of constant field strength extracted from an initial potential-field solution. Because this surface dips downward beneath helmet-streamer cusps, the model opens field lines preferentially near open–closed boundaries, producing about three times as much open magnetic flux (13.0 versus 4.6 in standard units) while keeping loop heights in the observed range from 0.05 to 1.08 solar radii. The authors show the modeled field matches extreme-ultraviolet and coronagraph images and reproduces near-Sun interplanetary field polarity and magnitude better than the standard model under the same open-flux constraint. The paper states that the surface is not iterated to full self-consistency and that the best variant is selected by matching observed open flux, so part of the flux increase is imposed by construction rather than freely predicted.","feed_headline":"Curved source surface nearly triples open solar flux","feed_subtitle":"Dipping the field-opening boundary, not lowering it everywhere, restores missing flux while keeping loops realistic.","key_machinery":"The central object is the Non-Spherical Source Surface (NSSS), a constant-|B| isosurface extracted from an initial potential-field solution and then used as the zero-potential upper boundary of a finite-element Laplace solver. The load-bearing property is that it is not a level surface in radius: it dips toward the Sun beneath current sheets at the cusps of helmet streamers, so the number of field lines that reach it and become open is larger locally than a sphere at the same mean height would allow. The surrounding structure is computed with a current-sheet layer, an exit sphere at 10 solar radii, and a spiral mapping to interplanetary space.","core_discovery":"The central claim is that the long-standing open-flux problem — models systematically predicting less open magnetic flux than spacecraft observe — can be substantially reduced by changing the shape of the source surface rather than by lowering it everywhere. The authors construct a Non-Spherical Potential Field model whose upper boundary is the isosurface of |B| from an initial potential-field source-surface solution. The isosurface automatically forms concave pockets under external current sheets at the bases of helmet streamers, so the potential-field layer is thinner there; lower loops open near separatrix boundaries while taller loops remain closed elsewhere. For a solar-maximum interval","pith_inferences":["If the constant-|B| source surface is as representative as the results suggest, future models could treat the source surface as an output of the field solution rather than an adjustable parameter, removing the residual tuning against observed open flux.","The compactness of modeled solar-wind source regions implies that the cross-hemisphere footpoint drift seen in spherical models may be an artifact of an over-high uniform boundary; slow-wind source mapping for space-weather connections could become more reliable.","The same finite-element extraction could be applied to MHD coronal solutions or to magnetograms of other stars, making the open-flux geometry testable beyond the Sun."],"forward_implications":["With one magnetogram as input, the same workflow yields coronal topology, current-sheet structure, and interplanetary field predictions, so source-surface height no longer has to be tuned separately for each purpose.","Open flux can be raised to observed levels without forcing all loops below a small height, removing a known artifact of simply lowering the spherical source surface.","Solar-wind source footpoints become more compact and localized than in models with a spherical surface, concentrating the missing flux near open–closed boundaries.","The resulting field can serve as a more realistic initial condition for global magnetohydrodynamic simulations of the corona and heliosphere."],"fun_headline_variants":["Curved boundary nearly triples open solar flux","Non-spherical model boosts open solar flux","Concave solar surface opens more magnetic flux","Reshaped source surface fixes open-flux problem","Bent boundary frees solar magnetic flux"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The model rests on assuming that a constant-|B| surface extracted from an initial spherical solution is a good stand-in for the physical surface where field lines open; the paper admits it cannot yet prove this surface would be self-consistently both a field-strength isosurface and a zero-potential surface, and the best-case variant is chosen to match observed open flux.","fun_headline_variants_meta":{"raw":{"variants":["Curved boundary nearly triples open solar flux","Non-spherical model boosts open solar flux","Concave solar surface opens more magnetic flux","Reshaped source surface fixes open-flux problem","Bent boundary frees solar magnetic flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000861,"raw_usage":{"total_tokens":3558,"prompt_tokens":717,"completion_tokens":2841,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":2772}},"tokens_in":461,"tokens_out":2841,"duration_ms":19148,"temperature":1.0,"reasoning_tokens":2772,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T16:55:48.337799+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the NSSS extraction on the magnetic field of an MHD coronal simulation where the true open/closed boundary is known, then test whether the constant-|B| surface coincides with that boundary and reproduces the observed open flux without choosing an initial radius; if the match fails, the geometrical mechanism is not the cause of the extra flux.","supporting_citations":[],"review_version":1}