{"id":"19e2a965-f4ae-454c-a02e-8d461d36510c","arxiv_id":"2601.14886","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A new stellarator equilibrium, CIEMAT-pw1, comes close to piecewise omnigenity while meeting stability, transport, and divertor criteria in simulations.","lead":"Plasma physicists have designed a new stellarator, CIEMAT-pw1, whose magnetic field traps charged particles unusually well while staying stable. If the result holds, fusion reactor designers gain a new family of shapes beyond today's quasi-symmetric stellarators.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The piecewise maximum-J benefit hinges on Eq. (6)'s factorized B(s,θ,ζ)=f(s)BpwO(θ,ζ), but actual equilibria have B_mn~s^{m/2} so different harmonics decay at different rates; the paper never directly verifies J(s,α,E,μ) or ∂sJ<0 away from s=0.5, leaving the volume-integrated transport and stability","rationale":"This is a serious numerical design study, not a fraud or a trivial consistency failure. The use of DESC, SFINCS, stella, ASCOT, and COBRA with honest caveats about electrostatic turbulence, missing coils, and neutronics is appropriate. The reader's CONDITIONAL verdict is well calibrated: the simulations cannot be independently reproduced from the manuscript alone, and the strongest physics claim depends on an approximate radial extension. The weakest point is indeed Eq. (6), exactly as the reader identifies. I add two refinements: (i) the factorization is not merely approximate but formally ambiguous regarding how Bmin and Bmax are encoded in f(s); (ii) the direct verification that would settle the issue—computing J(s,α,E,μ) and its radial derivative—is absent, even though the paper has the machinery to do it. The paper's own statement that pwO quality degrades with distance from s=0.5 (Section I) makes this gap load-bearing rather than cosmetic. If a direct J computation confirmed region-wise α-independence and ∂sJ<0, the central claim would be substantially strengthened; until then, CONDITIONAL is the right verdict. The missing Zenodo link and the explicitly deferred coil/neutronics work reinforce the conditional status but are not the main physics risk.","tokens_in":14503,"tokens_out":10276,"duration_ms":118188,"concrete_test":"Using the equilibrium files promised in Appendix A (the Zenodo link currently reads '[?]'), compute J(s,α,E,μ)=2∫√(2(E−μB))dl between bounce points for representative trapped-particle classes (e.g., E/μ=Bmin+0.1(Bmax−Bmin) and E/μ=Bmax−0.1(Bmax−Bmin)) at s=0.06, 0.25, 0.5, 0.75, and 1.0 in the β=3% equilibrium. Verify (i) within each pwO region J is α-independent to <1% and (ii) ∂sJ<0. If either condition fails, Eq. (7) and the maximum-J-based turbulence and fast-ion arguments are unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that CIEMAT-pw1 simultaneously achieves reactor-level transport and piecewise maximum-J stability rests on the factorization in Eq. (6): B(s,θ,ζ)=f(s)BpwO(θ,ζ) with all pwO shape parameters constant except Bmin and Bmax. This is the only analytic bridge from the optimized surface at s=0.5 to the full plasma volume. It is insecure for two reasons. First, the text does not state how Bmin and Bmax are absorbed into f(s); if f(s) is a common radial multiplier, Bmin and Bmax must share one radial profile, which is not demonstrated, and if they vary independently, Eq. (6) is not a valid representation. Second, even if Bmin and Bmax were fixed, real equilibria have B_mn~s^{m/2}, so high-m Fourier components of the broad-spectrum BpwO decay faster than low-m components; B(s)/BpwO is therefore not separable and the B contours change shape with s. The paper acknowledges degradation and points to Fig. 6 and the Γ_c proxy, but never computes the actual second adiabatic invariant J(s,α,E,μ) or its radial derivative from the equilibrium. Since the turbulence stabilization and fast-ion confinement arguments are explicitly derived from Eq. (7) (piecewise maximum-J), the volume-integrated benefits remain an unverified extrapolation. The unresolved Zenodo placeholder in Appendix A and the acknowledged absence of coil/neutronics design are additional but secondary limitations.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents CIEMAT-pw1, a fixed-boundary stellarator equilibrium optimized with DESC to minimize the local deviation δB (Eq. 5) between the magnetic field strength at s=0.5 and the piecewise-omnigenous (pwO) model field of Eq. (2) with p=2. The authors report that the optimized equilibrium has δB below 1% relative to the p=2 surrogate (and ~10% relative to the strict p→∞ pwO limit), low neoclassical effective ripple and bootstrap coefficient compared with W7-X, Mercier and COBRA stability, a rotational-transform profile compatible with an island divertor, electrostatic gyrokinetic turbulent fluxes in the reactor-relevant range, and an alpha heating efficiency near 95% at β=3%. The central claim is that CIEMAT-pw1 satisfies the standard set of physics criteria for a viable reactor candidate, thereby demonstrating piecewise omnigenity as a practical stellarator design concept.","tokens_in":15019,"tokens_out":4225,"duration_ms":42862,"significance":"If the results hold, CIEMAT-pw1 would be the first MHD equilibrium explicitly designed for piecewise omnigenity that simultaneously exhibits low neoclassical and turbulent transport, small bootstrap current, fast-ion confinement, and MHD stability, substantially broadening the stellarator design space beyond quasisymmetric and quasi-isodynamic concepts. The paper draws on a credible set of established numerical tools (DESC, SFINCS, stella, ASCOT, COBRA) and provides quantitative comparisons with W7-X, which strengthens the empirical content. The main gap is that the volume-integrated transport and stability benefits are bridged from the optimized surface to the whole plasma by the approximate factorization Eq. (6); this bridge is not directly verified via the second adiabatic invariant. If the authors close that gap and resolve the reproducibility issue in Appendix A, the paper would be a strong contribution to the field.","major_comments":[{"comment":"The volume-integrated conclusions (turbulence mitigation, fast-ion confinement, maximum-J) rely on Eq. (6), but the factorization is not established. Eq. (6) assumes B(s,θ,ζ)=f(s)BpwO with all shape parameters except Bmin and Bmax constant. The text does not specify how Bmin and Bmax are absorbed into f(s); if f(s) is a common radial multiplier, their radial profiles must be identical, which is not demonstrated, and if they vary independently, Eq. (6) is not a valid representation. Moreover, in real equilibria the Fourier harmonics of B decay as B_mn ~ s^{m/2}, so the high-m content of BpwO decays faster than the low-m content and B(θ,ζ) changes shape with s. The paper acknowledges degradation but never computes J(s,α,E,μ) or ∂sJ from the equilibrium; Γc is only a phase-space average of |∂αJ/∂sJ|^2 and does not establish the sign of ∂sJ. Since Eq. (7) is the bridge from s=0.5 to the full","section":"§I, Eq. (6)–(7)"},{"comment":"The abstract claims 'unprecedented levels of piecewise omnigenity,' but the quantity actually targeted and reported is δB relative to the p=2 surrogate BpwO, not relative to the p→∞ piecewise-omnigenous limit. Fig. 6 states that the p→∞ relative difference is ~10% in specific surface regions and was not an optimization target. This is an order of magnitude larger than the reported <1% deviation from the p=2 surrogate and raises the question of whether the configuration demonstrates pwO at the level claimed, or only closeness to a smooth surrogate that shares transport properties with pwO. The distinction should be stated precisely in the abstract and the p→∞ deviation should be quantified as a function of s.","section":"§II, Fig. 6"},{"comment":"The Mercier criterion is reported as not satisfied near the core at the highest β values, with the text attributing this to simulations becoming 'unreliable.' The abstract's claim of 'robust MHD stability across a range of β values' depends on this being a numerical artifact. Please show convergence of D_mercier with resolution at β=3%, or otherwise state explicitly which β range is stable. COBRA negative growth rates for all cases are supportive but do not resolve the reported core Mercier inconsistency.","section":"§III, Fig. 9"},{"comment":"The maximum-J property is stated to follow from ∂sBmin > 0 and ∂sBmax > 0 up to s=0.5. Under Eq. (6), ∂s f > 0 requires both Bmin and Bmax to increase radially; Fig. 8 shows this only for 'sufficiently high β' and up to s=0.5. The paper does not demonstrate ∂s f > 0 across the whole volume, nor does it present the direct computation of J(s,α,E,μ) that would validate Eq. (7) away from the optimized surface. The alpha-particle and turbulence conclusions are drawn at β=3%, where ∂sB are positive, so this is a partial gap rather than a fatal flaw, but the radial and β range of the maximum-J property should be made explicit.","section":"§II, Fig. 8 and Eq. (7)"},{"comment":"The data availability statement is incomplete: the text says the configuration and scripts 'are available in the Zenodo link [?]' with a literal placeholder. For a computational design paper, the central claims cannot be independently checked without access to the equilibrium and optimization scripts. This must be resolved before publication.","section":"Appendix A"}],"minor_comments":[{"comment":"Typos and minor wording issues: 'abscence' in the abstract; 'collisionlity' near Fig. 7; 'estabilization' near Fig. 12; 'orbit-averated' in §I. Also, the notation Γ_c is introduced in §II without defining the subscript c or the exact normalization used in Fig. 8.","section":"Abstract and §II"},{"comment":"The bootstrap coefficient is labeled D31; it would be clearer to use D_31 with sub/superscript formatting consistent with the text and references.","section":"§II, Fig. 7"},{"comment":"Eq. (2) defines the pwO field with w2=π, while w1 depends on ι through Eq. (3). In §I, the text says 'all the parameters discussed below equation (2) (and ι) kept constant throughout the plasma volume' — but ι is a radial profile in the equilibrium, so this assumption needs an explicit caveat; otherwise the reader cannot tell whether the factorization refers to the local ι at each s or to a fixed value.","section":"§I, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"This is a potentially important contribution to stellarator optimization, and I would not reject on novelty grounds. The main risk is that the volume-integrated benefits are extrapolated from s=0.5 through an unverified factorization; this is fixable by directly computing J and its radial derivative or by explicitly mapping the region of validity of Eq. (7). The abstract should also be calibrated to the p→∞ deviation reported in Fig. 6. The unresolved Zenodo placeholder is a reproducibility issue that should be fixed before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a serious numerical design study, not a new physics principle. The pwO concept and the target ansatz come from the same group's earlier work; what is new here is the first deliberately optimized, fully MHD equilibrium that comes close to piecewise omnigenity, with a broad multi-criteria check. That is a real step forward for the stellarator design program.\n\nThe paper does several things well. It uses standard codes (DESC, SFINCS, stella, ASCOT, COBRA) and reports honest caveats. The results are internally consistent: the relative difference to the p=2 target is below 1% at s=0.5, effective ripple is below W7-X, bootstrap current is small, Mercier and ballooning are stable for the cases shown, and the alpha heating efficiency at beta=3% is around 95%. The authors are also candid about the 10% deviation from the sharp p->infinity pwO limit and about coils being future work.\n\nThe soft spots are moderate. The main one is the factorization assumption in Eq. (6): B(s,theta,zeta)=f(s)BpwO(theta,zeta) with only Bmin and Bmax varying radially. The paper acknowledges that harmonics decay like s^{m/2}, but it never directly computes J(s,alpha,E,mu) or its radial derivative to verify the piecewise maximum-J property away from s=0.5. The Gamma_c proxy is a step in that direction—it's a phase-space average of |partial_alpha J / partial_s J|^2—but it is not the same as checking partial_s J < 0 region by region. So the volume-integrated turbulence and alpha-confinement benefits are partly extrapolated. This doesn't sink the paper, but it should be tightened in revision.\n\nAlso, Appendix A still shows the Zenodo link as a placeholder '?'. For a numerical design paper, the scripts and inputs are part of the result; that has to be fixed before anyone can reproduce the configuration. Minor: the Mercier criterion is described as 'unreliable' near the core at the highest beta, and the 'reactor candidate' framing goes beyond what fixed-boundary physics can claim. The authors say that themselves, so it's a framing issue, not an error.\n\nBottom line: I'd send this to peer review. It deserves a serious referee. The missing scripts and the unverified max-J radial extension should be addressed, but the core demonstration is credible. I'll probably cite it when I next write about stellarator optimization.","headline":"CIEMAT-pw1 is a credible first demonstration of a deliberately optimized piecewise-omnigenous equilibrium, but the radial extension of the maximum-J benefits is argued, not shown.","tokens_in":15472,"tokens_out":3712,"would_cite":true,"duration_ms":40560,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.55.Hc","52.25.Fi","52.35.Kt"],"model":"deepseek-v4-flash","headline":"A new stellarator magnetic configuration satisfies the ideal MHD equilibrium equation while achieving unprecedented levels of piecewise omnigenity, combining low transport and stability in one reactor-scale design.","keywords":["piecewise omnigenity","stellarator reactor","MHD equilibrium","maximum-J","neoclassical transport","gyrokinetic turbulence","alpha particle confinement","island divertor"],"falsifier":"Compute the exact B(s, theta, zeta) from the equilibrium and compare it with f(s) B_pwO(theta, zeta) across the whole radius; if the residual is large enough that the derivative of J with respect to s changes sign, or a significant variation of J within a flux surface appears away from s = 0.5, the piecewise maximum-J claim fails. A direct test would be a full-volume neoclassical and gyrokinetic simulation at reactor collisionality checking that the radial energy flux stays below the reactor threshold at all radii, not just at the single optimized surface.","tokens_in":14464,"feed_emoji":"⚛️","tokens_out":3985,"duration_ms":39784,"temperature":0.7,"pith_summary":"The paper argues that piecewise omnigenity — a relaxed notion of omnigenity where the second adiabatic invariant is constant within regions of a flux surface rather than on the whole surface — is not just a theoretical construct but a practical design target. It presents a five-field-period stellarator equilibrium optimized so that the magnetic field strength at the middle flux surface differs from an ideal pwO form by less than about one percent. The authors claim that this configuration meets the standard set of reactor physics criteria: low neoclassical and turbulent transport, small bootstrap current, MHD stability, and roughly 95 percent alpha heating efficiency at beta equal to three percent. A sympathetic reader should care because pwO fields may allow a broader family of reactor candidates than exact omnigenous designs, with added robustness to changes in the rotational transform.","feed_headline":"Stellarator design passes reactor physics checklist","feed_subtitle":"A new equilibrium pairs near-perfect piecewise omnigenity with MHD stability, low transport, and ~95% alpha heating.","key_machinery":"The load-bearing object is the parametrized piecewise-omnigenous field B_pwO(theta, zeta): a smooth (p = 2) generalization of a parallelogram-shaped field-strength distribution whose width and slope parameters are fixed by the requirements of collisionless trapped-particle confinement and zero bootstrap current. A fixed-boundary equilibrium solver and optimizer is used to minimize the relative deviation delta_B between the equilibrium B and this B_pwO at the middle flux surface, while also targeting Mercier stability, rotational transform, elongation, and mirror ratio. The argument that pwO closeness pays off in turbulence and fast-ion confinement relies on the approximate factorization B(s,","core_discovery":"On the paper's own terms, the central discovery is that a fixed-boundary stellarator equilibrium can be optimized to be very close to piecewise omnigenity while still solving the ideal MHD force balance equation and passing the usual reactor-physics filters. The optimized design has magnetic field strength at s = 0.5 within about one percent of the target pwO form (and roughly ten percent of the strict p-to-infinity limit), with the approximation degrading only weakly across the plasma volume. This closeness yields small effective ripple, a bootstrap transport coefficient comparable to that of a proven optimized stellarator, maximum-J-like behavior at finite beta, Mercier and ballooning stab","pith_inferences":["The single-surface optimization at s = 0.5 may understate edge transport, since effective ripple rises toward the edge; a full reactor assessment should check whether finite-volume degradation remains acceptable when realistic profiles and electromagnetic effects are included.","The robustness to rotational-transform changes hints that pwO designs might tolerate larger bootstrap current or coil tolerances than quasi-isodynamic designs, a property worth testing with explicit coil-error studies.","The coil-feasibility estimate suggests a minimum plasma-to-coil distance of roughly 1.3 meters in the tightest region; whether this is compatible with a breeding blanket is an open question the paper leaves to future work.","Because the optimizer simultaneously drives many targets, it is unclear how much of the good performance is due to pwO closeness versus the other constraints; a useful control experiment would be to optimize with the pwO deviation excluded and compare reactor metrics."],"forward_implications":["If correct, piecewise omnigenity becomes a third viable design family alongside quasi-axisymmetry and quasi-isodynamicity.","pwO configurations can satisfy all standard reactor physics criteria simultaneously in a single equilibrium.","Robustness of transport properties to changes in the rotational transform — notably bootstrap-current-induced changes — suggests design flexibility for reactors.","Turbulence can be mitigated through the piecewise maximum-J property even without exact omnigenity.","The concept opens a wider configuration space for future coil and divertor optimization."],"fun_headline_variants":["Unprecedented piecewise omnigenity in a stable stellarator","Stellarator passes reactor checklist with new confinement design","Near-perfect piecewise confinement meets MHD equilibrium","Piecewise omnigenous stellarator: a viable reactor candidate","Stellarator achieves piecewise omnigenity without sacrificing stability"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The assumption that the magnetic field of the almost-pwO configuration can be written approximately as B = f(s) B_pwO with all shape parameters fixed except the minimum and maximum field strength — if the radial variation of different Fourier harmonics is too strong, the piecewise maximum-J property and its transport benefits do not extend away from the optimized surface.","fun_headline_variants_meta":{"raw":{"variants":["Unprecedented piecewise omnigenity in a stable stellarator","Stellarator passes reactor checklist with new confinement design","Near-perfect piecewise confinement meets MHD equilibrium","Piecewise omnigenous stellarator: a viable reactor candidate","Stellarator achieves piecewise omnigenity without sacrificing stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000267,"raw_usage":{"total_tokens":1445,"prompt_tokens":733,"completion_tokens":712,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":630}},"tokens_in":477,"tokens_out":712,"duration_ms":8214,"temperature":1.0,"reasoning_tokens":630,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:03:01.659693+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact B(s, theta, zeta) from the equilibrium and compare it with f(s) B_pwO(theta, zeta) across the whole radius; if the residual is large enough that the derivative of J with respect to s changes sign, or a significant variation of J within a flux surface appears away from s = 0.5, the piecewise maximum-J claim fails. A direct test would be a full-volume neoclassical and gyrokinetic simulation at reactor collisionality checking that the radial energy flux stays below the reactor threshold at all radii, not just at the single optimized surface.","supporting_citations":[],"review_version":1}