{"id":"c50b2c02-342d-49c3-bb16-44cd7a424e8a","arxiv_id":"2501.12468","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"First joint optimization of passive superconducting coil arrays and background fields produces low-error magnetic field solutions for four stellarators.","lead":"Researchers optimized arrays of passive superconducting coils that carry current induced by a background magnetic field, eliminating power supplies. They designed such arrays for four different stellarator fusion concepts and report accurate magnetic field shaping.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The feasibility claim depends on neglecting Meissner effects, but Appendix A's uniform-field Type-I argument does not cover the Type-II HTS reactor-scale regime where field gradients are large; this unresolved physics could change the optimized currents and magnetic surfaces.","rationale":"The reader's weakest assumption identifies the Meissner-effect neglect as the principal concern, and my reading agrees. The paper's own Appendix A limits its cancellation argument to a spatially uniform B0 on the scale of the coil (with a first-order correction scaling as R/L0) and to a Type-I circular wire. The optimized reactor-scale PSCs have radii up to 3.26 m and the TF fields vary strongly across them, so R/L0 is not guaranteed to be small; additionally, HTS REBCO tapes are Type-II with a superconductor layer thickness of order the penetration depth, where the surface-current approximation is acknowledged as unclear. Since the central claim is that PSC arrays can achieve B dot n errors of 1e-3 or below, a magnetic perturbation of that size from Meissner or vortex-magnetization currents would directly undermine the result. I therefore propose a finite-element test that quantifies this perturbation for the actual optimized geometry. The minor textual inconsistency in the quasi-axisymmetry error (Sec. III.A) is not load-bearing. The verdict should remain CONDITIONAL, as the authors could address the concern with the proposed calculation or with an experimental prototype.","tokens_in":12027,"tokens_out":12820,"duration_ms":136804,"concrete_test":"Model one optimized configuration (e.g., Sec. III.A) with a finite-element H-formulation (or T-A formulation) in which each PSC filament is replaced by a REBCO tape cross-section (width 4 mm, superconductor thickness 1 micron, London penetration depth 0.2 micron, critical-current density from e.g. [35]), with TF coil currents fixed at their optimized values. Compute the quasi-static field after one TF ramp and evaluate <B dot n>/<B> on the plasma surface. If this error exceeds roughly 2e-3 while the filament model gives 5.9e-4, the Meissner/magnetization contribution is not subdominant and the design is invalidated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central feasibility claim (Sec. III) rests on the flux-freezing model LI+Psi=0 (Eq. 2), which neglects Meissner screening entirely (Sec. I.A, Appendix A). Appendix A's subdominance argument assumes a Type-I superconductor, a circular wire, and a background field uniform on the scale R of the coil, with corrections suppressed by R/L0<<1. In the reactor-scale designs of Secs. III.A-C, PSC radii are 1.1-3.3 m while the TF field varies on scales comparable to the minor radius (1.7 m), so R/L0 is not small; and the proposed REBCO tapes are Type-II with lambda about 0.2 microns comparable to the 1-2 micron superconductor layer, a regime the appendix explicitly leaves unresolved. The paper itself states that the surface-current approximation is of unclear accuracy for such tapes. If Meissner or vortex-magnetization currents produce fields at the plasma surface comparable to the claimed B dot n errors (5.9e-4 to 2.4e-3), the presented Poincare surfaces would not be reproduced by a physical array.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper formulates and demonstrates the first large-scale joint optimization of passive superconducting coil (PSC) arrays with the background magnetic field for stellarators. The currents in the PSC array are determined by flux-freezing through the linear system LI + Ψ = 0 (Eq. 2), and the optimization variables include the PSC orientations, locations, and shapes, together with the active TF coil shapes and currents. The authors present solutions for four stellarators (Landreman–Paul QA and QH, Schuett–Henneberg QA, and CSX), reporting normalized B·n errors between 5.9e-4 and 2.4e-3, along with Poincaré plots indicating good magnetic surfaces. The derivatives are computed with an analytic chain rule and validated against finite differences. The paper explicitly restricts the physics to the flux-freezing model and argues in Appendix A that Meissner effects are subdominant for Type-I superconductors in a sufficiently uniform background field.","tokens_in":12281,"tokens_out":3783,"duration_ms":39445,"significance":"If the underlying physical model is accepted, this work provides a valuable new capability: a general and open-source (SIMSOPT-based) framework for designing passive superconducting coil arrays that shape stellarator fields without power supplies. The reported B·n errors are low, the Poincaré plots show good surfaces, and the derivative formulation in Eqs. (4)–(5) is an important technical contribution, validated numerically. The use of a sparse-regression-type coil removal and the discussion of PSC currents, forces, and turns also make the paper practically useful. The main significance hinges on the physical validity of ignoring Meissner screening; the paper itself acknowledges that this is not established for high-temperature superconducting tapes, which directly affects the credibility of the reactor-scale feasibility claims.","major_comments":[{"comment":"The central feasibility claim rests on neglecting Meissner effects, but the appendix's justification covers only a Type-I superconductor with a circular cross-section in a background field that is uniform on the scale R of the coil. For the reactor-scale designs in Secs. III.A–C, the field gradient length is comparable to the PSC radii (e.g., a minor radius of 1.7 m versus PSC radii of 1.14–3.26 m), so the assumption L0 ≫ R used in the appendix fails. For the proposed REBCO tapes, which are Type-II with penetration depth comparable to the layer thickness, the paper itself states that the accuracy of approximating the screening currents as surface currents is unclear. Since Eq. (2) determines PSC currents solely from flux-freezing, any additional Meissner or vortex-magnetization currents would alter the currents and the resulting magnetic field, so the presented B·n errors and Poincaré surfaces would not necessarily be reproduced by a physical array. This point must be addressed quantitatively, or the claims must be restricted to the flux-freezing-only model.","section":"Appendix A / Sec. I.A"},{"comment":"The reported two-term quasisymmetry error is inconsistent: the text says the error is 'at approximately 7 × 10−5 is reduced by an order of magnitude from the original value of 7 × 10−6', but 7 × 10−5 is one order of magnitude larger than 7 × 10−6, not smaller. Please correct the direction or the numbers; this affects the claim about how much quasisymmetry was preserved or degraded.","section":"Sec. III.A, Eq. (6)"},{"comment":"The manuscript correctly notes that long-timescale coupling between passive-coil currents and plasma currents cannot be removed by a controller, unlike the case for powered coils. Because the PSC currents in Eq. (2) are determined by the vacuum flux Ψ from the TF coils, any modification of the field by equilibrium plasma currents will induce additional currents in the PSCs (to maintain zero enclosed flux), and this is not accounted for in the optimized solutions. The paper should either demonstrate via a self-consistent equilibrium calculation that the magnetic surfaces remain good when plasma currents are included, or discuss why this effect is expected to be small. As written, the feasibility claim extends beyond the vacuum-field regime that is actually optimized.","section":"Sec. I.A, Eq. (2)"}],"minor_comments":[{"comment":"The statement that the Jacobian calculations 'were extensively verified against finite differences' would be more informative if the maximum relative error and the perturbation sizes were reported.","section":"Sec. II.B"},{"comment":"There is a typographical error: the text reads 'diagmagnetic' where 'diamagnetic' is intended.","section":"Appendix A"},{"comment":"The sentence about the two-term quasi-axisymmetry error is difficult to parse even after correcting the numerical inconsistency; please rewrite it to state clearly the original value, the final value, and whether the change is an increase or a decrease.","section":"Sec. III.A"},{"comment":"The phrase 'the average normalized two-term quasi-axisymmetry error of approximately 3.9 × 10−3 is fairly close to the original value of 1.4 × 10−3' is vague; please clarify the magnitude of the change and whether this is considered acceptable in the context of the larger B·n error.","section":"Sec. III.C"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically sound within the flux-freezing model, and the optimization results appear reproducible. The main issue is the physical justification for neglecting Meissner effects in the reactor-scale HTS regime; if the authors can provide a quantitative estimate or a clear scoping of their claims, the paper would be much stronger. I also recommend the authors re-examine the quasisymmetry error statement in Sec. III.A, which currently contains an order-of-magnitude inconsistency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is the first real optimization of passive superconducting coil arrays for stellarators, and it is a solid piece of work. The core idea is simple: currents in superconducting loops are set by flux-freezing, LI+Ψ=0, and the authors derive the full derivative chain through the inductance-coupled system so that coil positions, orientations, and shapes can be optimized jointly with the background TF coils. They then demonstrate the method on four stellarators, including reactor-scale designs, with normalized B·n errors in the 1e-3 range and Poincaré plots showing good surfaces. That is a genuine step beyond permanent-magnet and active dipole arrays, and the paper is honest about constraints: no net toroidal flux, currents clamped at a few MA, plasma coupling ignored, and more superconducting material than modular coils.\n\nThe math checks out. Equations (4) and (5) are the right chain rule, and the finite-difference verification is convincing. The use of SIMSOPT is a plus. I also credit the authors for including Appendix A at all—too many papers would have just said 'Meissner effects are negligible.' That appendix, however, is the main soft spot. It proves that a uniform background field produces no net Meissner current on a circular Type-I coil. That is a nice result, but it does not cover the actual designs: the reactor-scale PSCs have radii 1.1–3.3 m while the field varies on scales comparable to the minor radius (1.7 m), so the R/L0≪1 assumption is not satisfied; and REBCO tapes are Type-II with λ≈0.2 μm comparable to the 1–2 μm layer, a regime the paper explicitly leaves unresolved. If Meissner or vortex-magnetization currents produce fields at the plasma surface comparable to the claimed ~1e-3 B·n errors, the optimized currents and surfaces would shift. This is not a fatal flaw—the framework stands as a design tool—but the feasibility claim for HTS reactor arrays should be softened until someone quantifies these fields. Minor issue: the quasisymmetry numbers in Sec. III.A are inconsistent (7×10−5 vs 7×10−6, 'reduced by an order of magnitude'), which needs a fix. No runnable artifacts, but the SIMSOPT implementation should make reproduction possible.\n\nBottom line: this paper deserves serious refereeing. The core contribution is real, the derivations are clean, and the limitations are mostly stated. I would send it to peer review with a request for a more quantitative Meissner estimate and a fix to the quasisymmetry numbers.","headline":"First real PSC array optimization for stellarators; solid math, but Meissner effects are not adequately quantified for the HTS reactor-scale regime.","tokens_in":12805,"tokens_out":2141,"would_cite":true,"duration_ms":20333,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.55.Hc"],"model":"deepseek-v4-flash","headline":"The paper shows that passive superconducting coil arrays, with currents induced by the background field rather than by power supplies, can be optimized to shape stellarator magnetic fields to high precision across four stellarator designs.","keywords":["stellarator","passive superconducting coils","coil optimization","inverse magnetostatics","flux-freezing","quasi-symmetry","magnetic field shaping"],"falsifier":"A benchmark experiment: place one passive HTS coil (REBCO tape) in a uniform, slowly ramped background field, measure the external field with a Hall probe, and compare to the Biot-Savart prediction from the flux-freezing current alone; a discrepancy at the level of the surface-current approximation would falsify the assumption.","tokens_in":11829,"feed_emoji":"🧲","tokens_out":4683,"duration_ms":41904,"temperature":0.7,"pith_summary":"This paper tries to show that arrays of passive superconducting coils (PSCs) can replace powered shaping coils in stellarators. Because a superconductor's zero-resistance loops keep their enclosed magnetic flux at zero, a changing background field induces currents in the PSCs with no power supplies at all. The authors perform the first large-scale joint optimization of the positions, orientations, and radii of these coils together with the background field coils, and report normalized B·n errors as low as 5.9e-4 with good magnetic surfaces for four stellarator designs. If correct, this opens a path toward simpler, cheaper, and more accessible stellarator coil systems.","feed_headline":"Passive coils shape stellarator fields with no power supplies","feed_subtitle":"Joint optimization of coil positions and background fields reaches 0.06 percent field error in four stellarators.","key_machinery":"The load-bearing mechanism is the flux-freezing constraint LI + Ψ = 0, which converts the infinitely many possible passive coil currents into a unique vector determined by the coil geometry and the background field. Because the currents depend on all coil positions through the inductance matrix and the flux integrals, the gradients of the objective are computed by an autodifferentiation chain that includes derivatives of the PSC currents with respect to every coil degree of freedom. The practical object that carries the optimization is the surface-averaged normal-field error K, together with engineered constraints on coil distance, curvature, and force limits.","core_discovery":"The central discovery is that passive superconductor coil arrays can be optimized to shape stellarator magnetic fields to precision comparable to active coil designs. The currents in the PSCs are not free parameters; they are locked by the flux-freezing condition LI + Ψ = 0, where L is the inductance matrix and Ψ is the flux from the powered coils. The optimization therefore varies the geometry of the PSCs and the active coils, with PSC currents computed by solving that linear system, and minimizes the surface-averaged normal-field error K = (1/2)∫|B·n|^2 dS. Applied to four stellarators (two Landreman-Paul designs, a Schuett-Henneberg design, and the CSX university experiment), the method achieves ⟨B·n⟩/⟨B⟩ errors from 5.9e-4 to 2.4e-3, preserves quasi-symmetry at levels better than most existing devices, and produces good Poincaré surfaces.","pith_inferences":["If the Meissner effect is not negligible for high-temperature superconducting tapes, the optimized currents and field errors would change; a direct measurement of the field around a passive HTS coil in a known background field would settle this.","The sparse-solution technique used to prune low-current coils could be applied more aggressively, potentially reducing the number of PSCs needed and lowering superconducting material cost.","The same flux-freezing formulation could be transferred to other coil-design problems outside fusion, such as compact MRI magnet arrays or accelerator magnets, wherever a background field is already present.","The paper's assumption that a spatially uniform background field on the coil scale makes Meissner contributions vanish suggests a testable extension: computing the Meissner correction for realistic nonuniform fields and Type-II conductors and re-running the optimizations."],"forward_implications":["Stellarator shaping fields can be produced without power supplies for the shaping coils, since currents are induced by the background toroidal field coils.","The optimization method applies to any inverse magnetostatic problem with passive conductors, not just circular coils; the paper demonstrates nonplanar and shape-varying passive coils for CSX.","PSC arrays keep coil forces and torques within material tolerances without explicit force objectives in most cases, because the induced currents are naturally limited.","University-scale stellarators such as CSX could be built with passive coils, reducing the number of powered window-pane coils needed.","For quasi-helically symmetric designs, field ripple from small passive coils can degrade quasi-symmetry, so PSC arrays may not benefit every stellarator type."],"supporting_citations":[{"why":"Kaptanoglu et al. dipole-array optimization; the baseline active-coil approach whose force and torque minimization the PSC work extends.","marker":"[23]"},{"why":"Helander et al., establishes that permanent-magnet or passive arrays cannot produce net toroidal flux, motivating the need for TF coils.","marker":"[24]"},{"why":"SIMSOPT code, the open-source framework in which the optimizations are implemented.","marker":"[25]"},{"why":"Bromberg et al., earlier proposal of passive superconducting tiles/monoliths for flux-surface prescription, which this work numerically optimizes.","marker":"[26]"},{"why":"Landreman-Paul precise quasisymmetry equilibria, the target configurations for the QA and QH reactor-scale designs.","marker":"[31]"},{"why":"Schuett-Henneberg compact QA design, the third reactor-scale test case.","marker":"[32]"},{"why":"Baillod et al., CSX design study whose window-pane coil solution the PSC solution for CSX is compared against.","marker":"[33]"},{"why":"ARIES-CS parameters used to scale the reactor-scale configurations to B0 = 5.7 T and r0 = 1.7 m.","marker":"[34]"},{"why":"Wechsung et al., four-coil modular solutions whose coil length and curvature the PSC TF coils are compared with.","marker":"[41]"}],"fun_headline_variants":["Passive coils sculpt stellarator fields without power","Induced currents let passive coils shape stellarator fields","No-power superconducting arrays optimize stellarator fields","First joint optimization of passive coils and background fields","Passive superconductor coils bend stellarator magnetic fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The designs assume the Meissner screening field of the superconducting material is negligible, so only the flux-freezing surface currents matter; if that fails for realistic high-temperature tapes, the optimized coil currents and resulting field shapes would differ.","fun_headline_variants_meta":{"raw":{"variants":["Passive coils sculpt stellarator fields without power","Induced currents let passive coils shape stellarator fields","No-power superconducting arrays optimize stellarator fields","First joint optimization of passive coils and background fields","Passive superconductor coils bend stellarator magnetic fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000738,"raw_usage":{"total_tokens":3233,"prompt_tokens":815,"completion_tokens":2418,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":431,"completion_tokens_details":{"reasoning_tokens":2345}},"tokens_in":431,"tokens_out":2418,"duration_ms":19083,"temperature":1.0,"reasoning_tokens":2345,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:09:28.116980+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A benchmark experiment: place one passive HTS coil (REBCO tape) in a uniform, slowly ramped background field, measure the external field with a Hall probe, and compare to the Biot-Savart prediction from the flux-freezing current alone; a discrepancy at the level of the surface-current approximation would falsify the assumption.","supporting_citations":[{"cited_title":"Helander, M","cited_arxiv_id":null,"evidence_quote":"Helander et al., establishes that permanent-magnet or passive arrays cannot produce net toroidal flux, motivating the need for TF coils."},{"cited_title":"Landreman, B","cited_arxiv_id":null,"evidence_quote":"SIMSOPT code, the open-source framework in which the optimizations are implemented."},{"cited_title":"Bromberg, M","cited_arxiv_id":null,"evidence_quote":"Bromberg et al., earlier proposal of passive superconducting tiles/monoliths for flux-surface prescription, which this work numerically optimizes."},{"cited_title":"Landreman and E","cited_arxiv_id":null,"evidence_quote":"Landreman-Paul precise quasisymmetry equilibria, the target configurations for the QA and QH reactor-scale designs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Schuett-Henneberg compact QA design, the third reactor-scale test case."},{"cited_title":"Integrating Novel Stellarator Single-Stage Optimization Algorithms to Design the Columbia Stellarator Experiment","cited_arxiv_id":"2409.05261","evidence_quote":"Baillod et al., CSX design study whose window-pane coil solution the PSC solution for CSX is compared against."},{"cited_title":"Najmabadi, A","cited_arxiv_id":null,"evidence_quote":"ARIES-CS parameters used to scale the reactor-scale configurations to B0 = 5.7 T and r0 = 1.7 m."},{"cited_title":"Wechsung, M","cited_arxiv_id":null,"evidence_quote":"Wechsung et al., four-coil modular solutions whose coil length and curvature the PSC TF coils are compared with."}],"review_version":1}