Pith. sign in

REVIEW 3 major objections 5 minor 3 references

Numerical approach to the London Equation of superconductivity

T0 review · 3 major / 5 minor · reviewed 2026-07-30 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In elongated ellipsoidal superconductors, the major axis is not the first vortex path: off-axis U-shaped competitors win the isoflux ratio.

desk verdict Solid FEM–BEM London solver with clean ball validation; the elongated-ellipsoid isoflux claim is honest numerical evidence inside a restricted family, not a secured non-uniqueness theorem. read the letter →

arxiv 2607.27103 v1 pith:4COHGTZV submitted 2026-07-29 math.NA cs.NAmath-phmath.APmath.MP

Nicolás Barnafi , Ignacio Labarca-Figueroa , Carlos Román This is my paper · ORCID
classification math.NAcs.NAmath-phmath.APmath.MP MSC 35Q5665N2282D5565N3078M15
keywords LondonequationisofluxproblemfirstcriticalfieldFEM-BEMcouplingH(curl)finiteelementstype-IIsuperconductivityvortexnucleation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper builds a numerical method that solves the London equation for a superconducting sample in all of three-dimensional space, without cutting the exterior off at an artificial boundary. From that solution it recovers the vector field B0 that enters the isoflux problem: the variational problem that predicts where magnetic vortices first nucleate in the Ginzburg–Landau model. On the unit ball the method recovers the known fact that the diameter aligned with the applied field is optimal. On sufficiently elongated ellipsoids aligned with the same field, a two-parameter family of off-axis piecewise-linear curves already beats the major axis in isoflux ratio. Rotational symmetry then turns any such off-axis maximizer into a continuous family, so uniqueness fails and a degenerate rotational direction appears. The result matters because existing sharp asymptotic expansions for the first critical field assume a unique non-degenerate maximizer; elongated samples appear to live outside that regime and may nucleate U-shaped vortices instead.

What carries the argument

A nonstandard FEM–BEM transmission coupling for the magnetic potential U0 (interior Nédélec elements plus second-kind boundary integral operators on the exterior traces), followed by a Kikuchi mixed H(curl)–H1 discretization that recovers the divergence-free field B0 from the computed H0=curl U0.

What would settle it

Recompute B0 on a finer mesh (or with higher-order elements) for the same a=0.2 ellipsoid and re-evaluate the isoflux ratios of the reported best competitor Γ(x0,λ) versus the major axis; if the inequality reverses or disappears under controlled error, the non-optimality claim fails.

Watch

Extended reading notes

Core claim

For prolate ellipsoids under a constant applied field along the major axis, numerical evaluation of B0 shows that certain off-axis competitors attain a strictly larger isoflux ratio than the major axis once the minor semi-axis is small enough (demonstrated for a=0.2). The major axis is therefore not a maximizer; rotational symmetry produces a continuous family of equivalent configurations and a degenerate rotational direction in the isoflux problem.

Load-bearing premise

That beating the major axis inside a restricted two-parameter family of piecewise-linear curves is already enough to conclude the major axis is not optimal, without solving the full isoflux problem or proving that discretization error cannot reverse the inequality.

Editorial extensions

If this is right

  • Sharp expansions of the first critical field that assume a unique non-degenerate isoflux maximizer cannot be applied directly to sufficiently elongated ellipsoids.
  • Vortex nucleation in cigar-shaped samples is expected to select from a continuous rotational family rather than a single distinguished filament.
  • The geometry of the winning competitors points toward smooth U-shaped vortex lines analogous to those seen in rotating Bose–Einstein condensates.
  • The same FEM–BEM pipeline can be used to map the critical elongation at which the major axis loses optimality for other axisymmetric samples.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural next computation is a systematic scan in the aspect-ratio parameter a to locate the critical elongation where the major-axis isoflux ratio is first overtaken.
  • If U-shaped maximizers persist under mesh refinement, the next-order vortex interaction energy in the Ginzburg–Landau expansion will need a genuinely multi-curve or continuum formulation.
  • The same transmission formulation could test non-ellipsoidal elongated domains (e.g., rounded cylinders) to see whether the loss of uniqueness is geometry-generic.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a FEM–BEM discretization of the London equation on R^3 by reformulating the magnetic potential as a transmission problem for an auxiliary field U0, coupling an interior H(curl) variational equation to a second-kind exterior BIE with Buffa–Christiansen duality pairing and a reduced Steklov–Poincaré operator. From the resulting H0 it recovers B0 via Kikuchi’s mixed curl–curl formulation with a weakly enforced divergence constraint, proves discrete well-posedness and Céa/Strang-type estimates, and validates O(h) rates on manufactured solutions and the explicit ball solution. As an application it evaluates isoflux ratios on a two-parameter family of piecewise-linear competitors in prolate ellipsoids and reports numerical evidence that, for a sufficiently elongated ellipsoid (a=0.2), some off-axis curves beat the major axis, implying non-uniqueness and a degenerate rotational direction under axial symmetry.

Significance. If the discretization analysis holds, the work supplies a truncation-free, H(curl)-conforming 3D solver for the bulk London problem that is of genuine use for first-critical-field and vortex-nucleation studies beyond the ball. The sequential H0→B0 pipeline, second-kind BIE conditioning, and Strang estimate for approximate sources are concrete algorithmic contributions. The isoflux finding, if robust, would be scientifically important: it would show that the unique nondegenerate maximizer assumed in the refined Hc1 expansions of Román–Sandier–Serfaty fails in elongated geometries and would point toward U-shaped nucleation scenarios analogous to rotating BECs. Strengths include explicit benchmarking against the independent London ball formula, standard FE exterior-calculus stability for the Kikuchi system, and openly framed “numerical evidence” rather than a claimed theorem for the ellipsoid conclusion.

major comments (3)
  1. [Section 7, Figures 6–7] §7 and Figs. 6–7: the central applied claim—that for a=0.2 the major axis is not an isoflux maximizer—rests on R(Γ(x0,λ))>R(Γ(0,0)) computed from a single discrete field B0,h at h≈0.0221. Lemma 5.4 only gives ||B0−B0,h||curl=O(h); there is no a-posteriori or mesh-refinement control on the line-integral (or Stokes flux) difference that defines the isoflux gap. Visible oscillations of curl U0,h near ∂Ω (Fig. 5) sit exactly where off-axis competitors meet the boundary. A refinement study (or quantified quadrature error) showing that the reported strict inequality persists and stabilizes under h↓0 is needed before the abstract’s “therefore not a maximizer / non-uniqueness / degenerate rotational direction” language is justified.
  2. [Section 7] §7, construction of Γ(x0,λ): superiority inside a restricted two-parameter piecewise-linear family (symmetric, orthogonal boundary meeting, convex SΓ) is used to conclude that the major axis is not optimal among all admissible curves in X. The geometric reduction (Stokes + nonnegativity of curl B0·ŷ, length comparison) correctly rules out longer curves and motivates the family, but does not by itself make the family dense enough in the isoflux landscape. Either enlarge the competitor class (smooth U-shaped curves, free endpoint angles) or explicitly limit the claim to “the major axis is outperformed by admissible competitors in this family, hence is not a global maximizer,” which is logically sufficient once the inequality is mesh-robust, and avoid suggesting that the true maximizers have been identified.
  3. [Section 5.1, Proposition 5.2] §5.1, Proposition 5.2: the inf-sup for the reduced coupling q is only sketched—compactness of K0 reduces the problem to a principal part with 2W0, then a contraction argument is cited to [EEK21, Ste11] without checking that those references cover the static Maxwell / London transmission setting and the precise trace spaces used here (including the div-free subspace for λ). A self-contained statement of the needed ellipticity/contraction constants, or a precise theorem citation matching operators (27) and spaces (29), is load-bearing for the well-posedness claim of the FEM–BEM scheme.
minor comments (5)
  1. [Section 6] Table 1–3: report the polynomial degree k used in NEDk/Pk and whether the same k is used for volume and boundary spaces; rates alone do not identify the scheme order.
  2. [Sections 6–7] Figure 2 and 5: the oscillations in H0,h=curl U0,h near the boundary are noted in the text; a brief remark on whether post-processed (e.g. projected) magnetic fields are used for isoflux integrals would help reproducibility.
  3. [Section 5.1] Eq. (22) vs. discrete system: the continuous right-hand side uses λ0,ex while the discrete line writes λex; unify notation.
  4. [Appendix A] Appendix A: the ball formulas are the right benchmark; stating the radius R used in the numerics (unit ball) next to the general-R expressions would avoid a trivial mismatch when comparing constants C and M.
  5. Typos/notation: “sufficiently” and similar fi-ligature artifacts appear throughout; “div∂Ω 0” spacing in H−1/2(div∂Ω 0,∂Ω) is hard to parse—prefer H−1/2(div0,∂Ω) as in (29).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: numerical London/isoflux pipeline is independently discretized, benchmarked, and not forced by self-citation or fitted inputs.

full rationale

The paper’s load-bearing chain is (i) reformulation of the London equation as a transmission problem, (ii) FEM–BEM discretization of H0/U0 and Kikuchi mixed FE for B0, (iii) convergence/Strang estimates, (iv) validation on manufactured solutions and the classical ball formula (Lon50/ABM06), and (v) evaluation of isoflux ratios of an explicit competitor family on a computed B0,h for elongated ellipsoids. None of these steps defines the claimed output in terms of itself, fits a parameter to the isoflux conclusion and renames it a prediction, or imports a uniqueness theorem that forces the ellipsoid non-optimality result. Self-citations (Rom19, RSS23, RSS25) supply analytic motivation for R0 and Hc1 and the known ball maximizer; they are not used as numerical ground truth for the ellipsoid inequality. The competitor family Γ(x0,λ) is an ansatz used only to exhibit larger ratios than the major axis, which is a one-sided numerical comparison rather than a circular derivation. Correctness risks (mesh error possibly reversing a small gap; restricted family) are separate from circularity. The derivation is self-contained against external benchmarks.

Assumptions & free parameters 3 free parameters · 6 assumptions · 1 invented entities

Load-bearing content is standard Maxwell/FEM–BEM analysis plus the physical identification of B0 with the isoflux integrand from prior GL asymptotics. No new physical entities; free choices are discretization and the restricted competitor family used to probe optimality.

free parameters (3)
  • Ellipsoid aspect a (demo value 0.2) = a=0.2 in main isoflux plots
    Elongation at which major-axis failure is reported is chosen for the numerical study; ‘sufficiently elongated’ is not given a sharp threshold.
  • Competitor family parameters (x0, λ) and 500×500 grid = best reported x0≈0.1529, λ≈0.8862
    Hand-designed two-parameter piecewise-linear family and search grid determine which off-axis curves are compared to the major axis.
  • Mesh size h and FE degree k = e.g. h≈0.0221 for a=0.2 ellipsoid; k≥1 Nédélec/Lagrange
    Discretization parameters control approximation quality of H0,h and B0,h and thus of computed isoflux ratios.
assumptions (6)
  • domain assumption London equation (1) and B0 system (3) correctly encode the Meissner/isoflux data of 3D magnetic GL as ε→0
    Taken from Rom19/RSS23/RSS25 and §2 derivation; numerics assume this continuum identification.
  • domain assumption Ω bounded, simply connected, C² boundary; H0,ex normalized divergence-free applied field
    Stated in §1.1; needed for traces, Helmholtz–Hodge, and exterior representation.
  • standard math de Rham / Babuška–Brezzi stability of Kikuchi mixed curl–curl with Nédélec–Lagrange pair
    Invoked via BBF13 and FE exterior calculus for well-posedness of (18)/(24).
  • standard math Second-kind Maxwell boundary integral operators (K0 compact, stable RWG–BC pairing) yield a valid exterior Steklov–Poincaré map
    Appendix B and BC07/BHvPS03; underpins FEM–BEM coupling (22).
  • domain assumption curl B0 · ŷ ≥ 0 in meridional sections so longer-than-axis curves cannot beat the axis by the Stokes argument
    Cited from RSS23 §3 and used in §7 to restrict competitor geometry.
  • ad hoc to paper Restricted family Γ(x0,λ) is rich enough that beating the major axis inside the family implies the major axis is not a global isoflux maximizer
    Methodological leap in §7: only necessity of non-optimality of the axis is truly shown if some admissible curve wins; they argue Γ(x0,λ)⊂X, which is fair for non-optimality, but geometry of true maximizers remains conjectural.
invented entities (1)
  • Auxiliary potential U0 (interior A0−∇φ0, exterior A0−A0,ex−∇ψ) independent evidence
    purpose: Cast London as a transmission problem with homogeneous exterior curl-curl and convenient jumps for FEM–BEM
    Definition (14) in §4; standard gauge/reformulation device, not a new physical field.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Numerical approach to the London Equation of superconductivity." pith.science (2026). https://pith.science/paper/4COHGTZV

@misc{pith2026260727103,
  author       = {Pith},
  title        = {Pith review of: Numerical approach to the London Equation of superconductivity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4COHGTZV}},
  note         = {Machine review of arXiv:2607.27103}
}
abstract

In this work, we propose a general discretization strategy for solving the London equation for type-II superconductors in the whole space $\mathbb{R}^3$. To compute the magnetic field $H_0$, we reformulate the problem for the magnetic potential as a transmission problem and discretize it through a nonstandard FEM-BEM coupling. This formulation accounts for both the bounded interior domain and the unbounded exterior domain without introducing an artificial truncation. We then compute the vector field $B_0$, which arises from the Helmholtz-Hodge decomposition of the magnetic potential in the superconducting sample. This field enters the isoflux problem, which identifies the curves along which vortex nucleation first becomes energetically favorable in the Ginzburg--Landau model of superconductivity. We recast the equations for $B_0$ using the mixed formulation of Kikuchi, in which the divergence-free constraint is imposed weakly, and discretize the resulting problem using a classical $H(\operatorname{curl})$-conforming finite element discretization. We validate our discretization strategy through convergence tests and conclude with an application to the isoflux problem. For a ball under a constant applied magnetic field, the unique maximizer is the diameter aligned with the field. For ellipsoids under a constant applied magnetic field aligned with their major axis, our computations provide numerical evidence of a different behavior in sufficiently elongated, cigar-shaped geometries: off-axis competitors reminiscent of U-shaped vortex configurations attain a larger isoflux ratio than the major axis. Since the major axis is therefore not a maximizer, any off-axis maximizer generates, by rotational symmetry, a continuous family of equivalent configurations, implying non-uniqueness and the presence of a degenerate rotational direction in the isoflux problem.

Figures

Figures reproduced from arXiv: 2607.27103 by the authors.

Figure 1
Figure 1. Off-axis competitor attaining the largest isoflux ratio within the family considered, with an isoflux ratio larger than that of the major axis of the ellipsoid and a geometry resembling a U-shaped vortex configuration. See Section 7 for details [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Theoretical validation test on the unit ball. 7. Application to the isoflux problem We consider the prolate ellipsoid Ω =  (x, y, z) ∈ R 3 : x 2 + y 2 a 2 + z 2 < 1  , 0 < a < 1, and assume that the normalized applied magnetic field is H0,ex = ˆz. Thus, the major axis of the ellipsoid is aligned with the z-direction, while the two minor semi￾axes, in the x- and y-directions, have length a. Since both the domain an… view at source ↗
Figure 3
Figure 3. An oriented curve Γ lying on one side of the major axis and the region SΓ enclosed by Γ and the corresponding boundary curve Γbd. The boundary condition B0 × ν = 0 on ∂Ω is crucial here, since it implies that B0 has no tangential component along Γbd and hence Z Γbd B0 · dℓ = 0. We may therefore close Γ along the boundary without changing its line integral: Z Γ B0 · dℓ = Z ∂SΓ B0 · dℓ. Applying Stokes’ theorem in the… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Geometric construction of the family Γ(x0,λ) in the (x, z)-section of the ellipsoid. The figure shows the limiting configuration Γ(x0,0), an intermediate competitor Γ(x0,λ) , and the configuration Γ(x0,1), which reaches and follows the major axis. Fix x0 ∈ (0, a) and d…
Figure 5
Figure 5. Figure 5: Discrete fields H0,h and B0,h for the ellipsoid, used to evaluate the isoflux ratio. providing numerical evidence that the major axis is not optimal when the minor semi-axis is sufficiently small. The corresponding competitor is shown in [PITH_FULL_IMAGE:figures/full_…
Figure 6
Figure 6. Figure 6: Isoflux ratios of the competitors Γ(x0,λ) for a sample of x0-values. For the unit ball (left), none of the competitors considered outperforms the diameter aligned with the applied field. For the elongated ellipsoid (right), some off-axis competitors have a larger isofl…
Figure 7
Figure 7. Figure 7: Best off-axis competitor within the family considered for the ellipsoid with a = 0.2, with an isoflux ratio larger than that of the major axis. Appendix A. Explicit solution in a ball under a uniform normalized applied field We recall the explicit solution of the Londo…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references · 1 linked inside Pith

  1. [339]

    4 [Hip02b] R. Hiptmair, Symmetric coupling for eddy current problems , SIAM J. Numer. Anal. 40 (2002), no. 1, 41–65. MR 1921909

    MR 2009375 "4 [Hip02b] R. Hiptmair, Symmetric coupling for eddy current problems , SIAM J. Numer. Anal. 40 (2002), no. 1, 41–65. MR 1921909 "4, 24 [Kik87] F. Kikuchi, Mixed and penalty formulations for finite element analysis of an eigenvalue problem in electromagnetism, Proceedings of the first world congress on computational mechanics (Austin, Tex., 198...

    arXiv 2002

  2. [485]

    4 [Bra94] E. H. Brandt, Thin superconductors in a perpendicular magnetic ac field: General formulation and strip geometry , Physical Review B 49 (1994), no. 13, 9024–9040

    MR 2012928 "4, 8, 24, 25 [BJOS13] S. Baldo, R. L. Jerrard, G. Orlandi, and H. M. Soner, Vortex density models for superconductivity and superfluidity , Comm. Math. Phys. 318 (2013), no. 1, 131–171. MR 3017066 "4 [Bra94] E. H. Brandt, Thin superconductors in a perpendicular magnetic ac field: General formulation and strip geometry , Physical Review B 49 (1...

    arXiv 2013

  3. [2018]

    10, 11, 14 [BC07] A. Buffa and S. H. Christiansen, A dual finite element complex on the barycentric refinement , Math. Comp. 76 (2007), no. 260, 1743–1769. MR 2336266

    MR 3908678 "5 [BBF13] D. Boffi, F. Brezzi, and M. Fortin, Mixed finite element methods and applications , Springer Series in Computational Mathematics, vol. 44, Springer, Heidelberg, 2013. MR 3097958 "10, 11, 14 [BC07] A. Buffa and S. H. Christiansen, A dual finite element complex on the barycentric refinement , Math. Comp. 76 (2007), no. 260, 1743–1769. M...

    2013

Pith tools

Reviewed July 30, 2026 · model on record in the stance chip above.