Singular Vector Finite Element Basis Functions for Tetrahedra in Complex Electromagnetic Geometries
Pith reviewed 2026-06-26 21:40 UTC · model grok-4.3
The pith
Singular vector basis functions for tetrahedra substantially improve accuracy in modeling electromagnetic fields near singularities while cutting computational costs.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors develop singular vector finite element basis functions for tetrahedra that are additive to the standard vector basis functions. These functions are constructed to be singular with respect to the nodes and edges of the tetrahedron, allowing combinations that accommodate multiple unique singular features in one element. Higher-order interpolatory versions are also provided. When used in electromagnetic FEM simulations, these basis functions substantially improve accuracy in representing singular fields and enable the performance of otherwise expensive simulations at far lower computational costs, as shown by reducing a superconducting qubit simulation from 21.27 hours and 6.23 TB m
What carries the argument
additive singular vector basis functions singular with respect to each node and edge, allowing per-feature combinations inside a single tetrahedron while remaining compatible with standard bases
If this is right
- Accuracy improves substantially relative to standard basis functions near singular features.
- Simulations that are otherwise prohibitively expensive become feasible at far lower cost.
- The same element can handle multiple distinct singular features without separate meshing treatment.
- Higher-order versions of the singular functions further increase modeling accuracy.
- Design optimization becomes practical for electromagnetic geometries dominated by singular behavior such as superconducting qubits.
Where Pith is reading between the lines
- The memory and processor reductions could make iterative design loops feasible for structures where singular fields control performance.
- The additive construction suggests the functions can be introduced incrementally into existing mixed-basis codes without requiring a full rewrite.
Load-bearing premise
The new singular functions remain linearly independent from the standard bases and their combinations for multiple singularities inside one tetrahedron do not introduce instabilities in the assembled system.
What would settle it
A side-by-side comparison on the superconducting-qubit geometry in which the singular-basis results for critical design quantities deviate from a converged reference solution obtained with an extremely refined standard mesh, or in which the global matrix becomes singular.
Figures
read the original abstract
Electromagnetic finite element method (FEM) implementations using traditional basis functions struggle to accurately represent field behavior near singular features such as conducting wedges. To combat this, specialized singular basis functions have been introduced to directly model the singular fields in these regions, leading to substantially improved performance. While these efforts have been pursued extensively in 2D, few functions have been developed for 3D elements. In this work, we develop basis functions for this in tetrahedra. Unlike prior functions, these basis functions are additive, meaning they are included alongside the standard vector basis functions to achieve more robust performance. Further, these functions are designed to be adaptable to tetrahedra touching several unique singular features by using combinations of basis functions singular with respect to each node and edge in the element, making them applicable to highly complex geometries. Higher-order interpolatory versions of the basis functions for modeling singular behavior with greater accuracy are also provided. These basis functions lead to substantial improvements in accuracy relative to the standard basis functions, and allow otherwise expensive simulations to be performed at far lower costs. As an application example, we perform simulations to extract critical quantities for designing superconducting qubits that significantly depend on the behavior of singular fields. In Ansys HFSS, this took 21.27 hours and a peak memory usage of 6.23 TB with 800 processors available, while using our singular basis functions achieved comparable results in 196 seconds while using 27.24 GB of memory and only 16 processors. Due to these benefits, our singular basis functions could be applied to enable design optimization of electromagnetic geometries with dominantly singular behavior, such as superconducting qubits.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops additive singular vector finite element basis functions for tetrahedral elements to accurately represent singular electromagnetic fields near conducting wedges and other features. Unlike prior 3D work, the functions are designed to be included alongside standard vector bases and to handle multiple distinct singular nodes/edges within one tetrahedron via per-feature combinations; higher-order interpolatory versions are also derived. The central demonstration is a superconducting qubit simulation in which the new basis functions achieve comparable accuracy to Ansys HFSS while reducing runtime from 21.27 hours to 196 seconds and memory from 6.23 TB to 27.24 GB.
Significance. If the linear-independence and stability claims hold, the work would be a meaningful advance for computational electromagnetics in geometries dominated by singularities, enabling routine high-fidelity modeling of devices such as superconducting qubits at far lower cost. The additive construction and explicit multi-singularity adaptability for tetrahedra address a documented gap relative to existing 2D and limited 3D singular bases.
major comments (1)
- [Basis function construction and assembly (sections describing the additive property and multi-feature combinations)] The central claim that the singular functions remain linearly independent from the standard vector bases and that combinations for multiple singularities per tetrahedron produce stable local matrices rests on the construction alone. No explicit verification—such as condition-number tables for the combined element matrix, null-space dimension checks, or eigenvalue spectra when two or more singular features coexist—is supplied in the sections describing the basis assembly or the global system. This property is load-bearing for the reported performance gains and for the assertion of applicability to complex geometries.
minor comments (2)
- [Numerical results] The abstract states numerical gains but the manuscript should include at least one table or figure quantifying the L2 or energy-norm error reduction versus polynomial degree and mesh density for a canonical singular test case (e.g., a conducting wedge) before the qubit example.
- [Mathematical formulation] Notation for the per-node and per-edge singular functions should be introduced with an explicit table listing the singularity exponents and the associated vector forms to improve readability.
Simulated Author's Rebuttal
We thank the referee for their careful review and constructive comments on our manuscript. We address the single major comment below.
read point-by-point responses
-
Referee: [Basis function construction and assembly (sections describing the additive property and multi-feature combinations)] The central claim that the singular functions remain linearly independent from the standard vector bases and that combinations for multiple singularities per tetrahedron produce stable local matrices rests on the construction alone. No explicit verification—such as condition-number tables for the combined element matrix, null-space dimension checks, or eigenvalue spectra when two or more singular features coexist—is supplied in the sections describing the basis assembly or the global system. This property is load-bearing for the reported performance gains and for the assertion of applicability to complex geometries.
Authors: We agree that the current manuscript presents the linear independence and stability properties as following from the additive construction without supplying explicit numerical verification (condition numbers, null-space checks, or eigenvalue spectra) for multi-singularity tetrahedra. The construction is explicitly designed to preserve independence by superposing per-feature singular modes that are orthogonal to the standard vector bases, but we acknowledge the absence of supporting numerical data. In the revised version we will add, in the basis-assembly sections, condition-number tables for the local matrices of elements containing one and multiple singular features, together with eigenvalue spectra of the combined element matrices to confirm stability. revision: yes
Circularity Check
New additive basis construction presented without reduction to inputs or self-citations
full rationale
The paper introduces singular vector basis functions for tetrahedra as an explicit new construction that is additive to standard vector bases and adaptable via per-node/edge combinations. No equations, derivations, or fitted parameters are shown that reduce any claimed prediction or uniqueness result to prior inputs by construction. The abstract and description frame the work as a methodological development for complex 3D geometries, with performance gains demonstrated via external comparison to Ansys HFSS rather than internal self-referential fitting. This is a self-contained construction; no load-bearing self-citation chains or self-definitional steps are present.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
The behavior of electromagnetic fields at edges,
J. Meixner, “The behavior of electromagnetic fields at edges,”IEEE Transactions on Antennas and Propagation, vol. 20, no. 4, pp. 442– 446, 1972
1972
-
[2]
The edge condition in electromagnetics,
R. Hurd, “The edge condition in electromagnetics,”IEEE Transactions on Antennas and Propagation, vol. 24, no. 1, pp. 70–73, 1976
1976
-
[3]
Field singularities at metal-dielectric wedges,
J. Van Bladel, “Field singularities at metal-dielectric wedges,”IEEE Transactions on Antennas and Propagation, vol. 33, no. 4, pp. 450– 455, 1985
1985
-
[4]
Wiley-IEEE Press, 1996
——,Singular Electromagnetic Fields and Sources. Wiley-IEEE Press, 1996
1996
-
[5]
The behavior of the electromagnetic field near the edge of a resistive half-plane,
I. Braver, P. Fridberg, K. Garb, and I. Yakover, “The behavior of the electromagnetic field near the edge of a resistive half-plane,”IEEE Transactions on Antennas and Propagation, vol. 36, no. 12, pp. 1760– 1768, 1988
1988
-
[6]
Quasi-TEM analysis of microwave trans- mission lines by the finite-element method,
Z. Pantic and R. Mittra, “Quasi-TEM analysis of microwave trans- mission lines by the finite-element method,”IEEE Transactions on Microwave Theory and Techniques, vol. 34, no. 11, pp. 1096–1103, 1986. IEEE TRANSACTIONS ON MICROW A VE THEORY AND TECHNIQUES 15
1986
-
[7]
Incorporation of edge conditions in moment method solutions,
D. Wilton and S. Govind, “Incorporation of edge conditions in moment method solutions,”IEEE Transactions on Antennas and Propagation, vol. 25, no. 6, pp. 845–850, 1977
1977
-
[8]
The calculation of stress intensity factors using the finite el- ement method with cracked elements,
E. Byskov, “The calculation of stress intensity factors using the finite el- ement method with cracked elements,”International Journal of Fracture Mechanics, vol. 6, no. 2, pp. 159–167, 1970
1970
-
[9]
Finite elements for determination of crack tip elastic stress intensity factors,
D. M. Tracey, “Finite elements for determination of crack tip elastic stress intensity factors,”Engineering Fracture Mechanics, vol. 3, no. 3, pp. 255–265, 1971
1971
-
[10]
Analysis of power type singularities using finite elements,
D. M. Tracey and T. S. Cook, “Analysis of power type singularities using finite elements,”International Journal for Numerical Methods in Engineering, vol. 11, no. 8, pp. 1225–1233, 1977
1977
-
[11]
Families of consistent conforming elements with singular derivative fields,
M. Stern, “Families of consistent conforming elements with singular derivative fields,”International Journal for Numerical Methods in En- gineering, vol. 14, no. 3, pp. 409–421, 1979
1979
-
[12]
Efficient singular element for finite element analysis of quasi-TEM transmission lines and waveguides with sharp metal edges,
J. Gil and J. Zapata, “Efficient singular element for finite element analysis of quasi-TEM transmission lines and waveguides with sharp metal edges,”IEEE Transactions on Microwave Theory and Techniques, vol. 42, no. 1, pp. 92–98, 1994
1994
-
[13]
Mixed finite elements inR 3,
J. C. Nedelec, “Mixed finite elements inR 3,”Numerische Mathematik, vol. 35, no. 3, pp. 315–341, 1980
1980
-
[14]
Two- dimensional singular vector elements for finite-element analysis,
Z. Pantic-Tanner, J. Scott Savage, D. Tanner, and A. Peterson, “Two- dimensional singular vector elements for finite-element analysis,”IEEE Transactions on Microwave Theory and Techniques, vol. 46, no. 2, pp. 178–184, 1998
1998
-
[15]
Two-dimensional curl- conforming singular elements for FEM solutions of dielectric waveg- uiding structures,
D.-K. Sun, L. Vardapetyan, and Z. Cendes, “Two-dimensional curl- conforming singular elements for FEM solutions of dielectric waveg- uiding structures,”IEEE Transactions on Microwave Theory and Tech- niques, vol. 53, no. 3, pp. 984–992, 2005
2005
-
[16]
Substitutive divergent bases for fem modelling of field singularities near a wedge,
J. Masoni, G. Pelosi, and S. Selleri, “Substitutive divergent bases for fem modelling of field singularities near a wedge,”Microwave and Optical Technology Letters, vol. 44, no. 4, pp. 327–328, 2005
2005
-
[17]
A new edge element for the modeling of field singularities in transmission lines and waveguides,
J. Gil and J. Webb, “A new edge element for the modeling of field singularities in transmission lines and waveguides,”IEEE Transactions on Microwave Theory and Techniques, vol. 45, no. 12, pp. 2125–2130, 1997
1997
-
[18]
Singular higher order complete vector bases for finite methods,
R. Graglia and G. Lombardi, “Singular higher order complete vector bases for finite methods,”IEEE Transactions on Antennas and Propa- gation, vol. 52, no. 7, pp. 1672–1685, 2004
2004
-
[19]
Singular hierarchical curl-conforming vector bases for triangular cells,
R. D. Graglia, A. F. Peterson, L. Matekovits, and P. Petrini, “Singular hierarchical curl-conforming vector bases for triangular cells,”IEEE Transactions on Antennas and Propagation, vol. 62, no. 7, pp. 3632– 3644, 2014
2014
-
[20]
Hierarchical additive basis functions for the finite-element treat- ment of corner singularities,
——, “Hierarchical additive basis functions for the finite-element treat- ment of corner singularities,”Electromagnetics, vol. 34, no. 3-4, pp. 171–198, 2014
2014
-
[21]
Singular basis functions and curvilinear triangles in the solution of the electric field integral equation,
W. Brown and D. Wilton, “Singular basis functions and curvilinear triangles in the solution of the electric field integral equation,”IEEE Transactions on Antennas and Propagation, vol. 47, no. 2, pp. 347– 353, 1999
1999
-
[22]
First-order singular basis functions for corner diffraction analysis using the method of moments,
A. K. Ozturk, R. Paknys, and C. W. Trueman, “First-order singular basis functions for corner diffraction analysis using the method of moments,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 10, pp. 3160–3168, 2009
2009
-
[23]
Singular higher order divergence- conforming bases of additive kind and moments method applications to 3D sharp-wedge structures,
R. D. Graglia and G. Lombardi, “Singular higher order divergence- conforming bases of additive kind and moments method applications to 3D sharp-wedge structures,”IEEE Transactions on Antennas and Propagation, vol. 56, no. 12, pp. 3768–3788, 2008
2008
-
[24]
Hierarchical divergence conforming bases for edge singularities in quadrilateral cells,
R. D. Graglia, A. F. Peterson, and P. Petrini, “Hierarchical divergence conforming bases for edge singularities in quadrilateral cells,”IEEE Transactions on Antennas and Propagation, vol. 66, no. 11, pp. 6191– 6201, 2018
2018
-
[25]
Hierarchical divergence conforming bases for tip singularities in quadrilateral cells,
——, “Hierarchical divergence conforming bases for tip singularities in quadrilateral cells,”IEEE Transactions on Antennas and Propagation, vol. 68, no. 12, pp. 7986–7994, 2020
2020
-
[26]
Singular tetrahedral finite elements for vector electro- magnetics,
J. P. Webb, “Singular tetrahedral finite elements for vector electro- magnetics,”IEEE Transactions on Microwave Theory and Techniques, vol. 55, no. 3, pp. 533–540, 2007
2007
-
[27]
Gradient-singular, hierarchical finite elements for vector electro- magnetics,
——, “Gradient-singular, hierarchical finite elements for vector electro- magnetics,”IEEE Transactions on Antennas and Propagation, vol. 60, no. 6, pp. 2814–2820, 2012
2012
-
[28]
Basis functions for vertex, edge, and corner singularities: A review,
A. F. Peterson and R. D. Graglia, “Basis functions for vertex, edge, and corner singularities: A review,”IEEE Journal on Multiscale and Multiphysics Computational Techniques, vol. 1, pp. 161–175, 2016
2016
-
[29]
S. T. Elkin, G. Khan, E. Forati, B. W. Langley, D. Timucin, R. Molavi, S. Sussman, and T. E. Roth, “Opportunities and challenges of compu- tational electromagnetics methods for superconducting circuit quantum device modeling: A practical review,”arXiv preprint arXiv:2511.20774, 2025
arXiv 2025
-
[30]
A quantum engineer’s guide to superconducting qubits,
P. Krantz, M. Kjaergaard, F. Yan, T. P. Orlando, S. Gustavsson, and W. D. Oliver, “A quantum engineer’s guide to superconducting qubits,” Applied Physics Reviews, vol. 6, no. 2, p. 021318, 2019
2019
-
[31]
Precision measurement of the microwave dielectric loss of sapphire in the quantum regime with parts-per-billion sensitivity,
A. P. Read, B. J. Chapman, C. U. Lei, J. C. Curtis, S. Ganjamet al., “Precision measurement of the microwave dielectric loss of sapphire in the quantum regime with parts-per-billion sensitivity,”Physical Review Applied, vol. 19, p. 034064, Mar 2023
2023
-
[32]
Surpassing millisecond coherence in on chip superconducting quantum memories by optimizing materials and circuit design,
S. Ganjam, Y . Wang, Y . Lu, A. Banerjee, C. U. Leiet al., “Surpassing millisecond coherence in on chip superconducting quantum memories by optimizing materials and circuit design,”Nature Communications, vol. 15, no. 1, p. 3687, 2024
2024
-
[33]
Characterization and reduction of microfabrication-induced decoher- ence in superconducting quantum circuits,
C. M. Quintana, A. Megrant, Z. Chen, A. Dunsworth, B. Chiaroet al., “Characterization and reduction of microfabrication-induced decoher- ence in superconducting quantum circuits,”Applied Physics Letters, vol. 105, no. 6, p. 062601, Aug 2014
2014
-
[34]
Investigating surface loss effects in superconducting transmon qubits,
J. M. Gambetta, C. E. Murray, Y .-K.-K. Fung, D. T. McClure, O. Dial, W. Shanks, J. W. Sleight, and M. Steffen, “Investigating surface loss effects in superconducting transmon qubits,”IEEE Transactions on Applied Superconductivity, vol. 27, no. 1, pp. 1–5, 2017
2017
-
[35]
Analysis and mitigation of interface losses in trenched superconducting coplanar waveguide resonators,
G. Calusine, A. Melville, W. Woods, R. Das, C. Stullet al., “Analysis and mitigation of interface losses in trenched superconducting coplanar waveguide resonators,”Applied Physics Letters, vol. 112, no. 6, p. 062601, Feb. 2018
2018
-
[36]
Determining interface dielec- tric losses in superconducting coplanar-waveguide resonators,
W. Woods, G. Calusine, A. Melville, A. Sevi, E. Golden, D. Kim, D. Rosenberg, J. Yoder, and W. Oliver, “Determining interface dielec- tric losses in superconducting coplanar-waveguide resonators,”Physical Review Applied, vol. 12, p. 014012, Jul 2019
2019
-
[37]
Surface loss simulations of superconducting coplanar waveguide resonators,
J. Wenner, R. Barends, R. C. Bialczak, Y . Chen, J. Kellyet al., “Surface loss simulations of superconducting coplanar waveguide resonators,” Applied Physics Letters, vol. 99, no. 11, p. 113513, Sep 2011
2011
-
[38]
Surface participation and dielectric loss in superconducting qubits,
C. Wang, C. Axline, Y . Y . Gao, T. Brecht, Y . Chuet al., “Surface participation and dielectric loss in superconducting qubits,”Applied Physics Letters, vol. 107, no. 16, p. 162601, Oct 2015
2015
-
[39]
Wiring surface loss of a superconducting transmon qubit,
N. S. Smirnov, E. A. Krivko, A. A. Solovyova, A. I. Ivanov, and I. A. Rodionov, “Wiring surface loss of a superconducting transmon qubit,” Scientific Reports, vol. 14, no. 1, p. 7326, 2024
2024
-
[40]
Z. Wang, S. Guan, F. Wu, X. Zhang, Q. Li, J. Chen, X. Wan, T. Xia, and H.-H. Zhao, “SesQ: A surface electrostatic simulator for precise energy participation ratio simulation in superconducting qubits,”arXiv preprint arXiv:2603.28524, 2026
Pith/arXiv arXiv 2026
-
[41]
J. D. Jackson,Classical Electrodynamics. John Wiley & Sons, 1998
1998
-
[42]
Jin,The Finite Element Method in Electromagnetics, 3rd ed
J.-M. Jin,The Finite Element Method in Electromagnetics, 3rd ed. John Wiley & Sons, 2015
2015
-
[43]
Higher order interpolatory vector bases for computational electromagnetics,
R. Graglia, D. Wilton, and A. Peterson, “Higher order interpolatory vector bases for computational electromagnetics,”IEEE Transactions on Antennas and Propagation, vol. 45, no. 3, pp. 329–342, 1997
1997
-
[44]
Hierarchal vector basis functions of arbitrary order for trian- gular and tetrahedral finite elements,
J. Webb, “Hierarchal vector basis functions of arbitrary order for trian- gular and tetrahedral finite elements,”IEEE Transactions on Antennas and Propagation, vol. 47, no. 8, pp. 1244–1253, 1999
1999
-
[45]
Moderate-degree tetrahedral quadrature formulas,
P. Keast, “Moderate-degree tetrahedral quadrature formulas,”Computer Methods in Applied Mechanics and Engineering, vol. 55, no. 3, pp. 339–348, 1986
1986
-
[46]
Coplanar waveguides for MMIC applications: Effect of upper shielding, conductor backing, finite-extent ground planes, and line-to-line coupling,
G. Ghione and C. Naldi, “Coplanar waveguides for MMIC applications: Effect of upper shielding, conductor backing, finite-extent ground planes, and line-to-line coupling,”IEEE Transactions on Microwave Theory and Techniques, vol. 35, no. 3, pp. 260–267, 1987
1987
-
[47]
Disentangling losses in tantalum superconducting circuits,
K. D. Crowley, R. A. McLellan, A. Dutta, N. Shumiya, A. P. M. Placeet al., “Disentangling losses in tantalum superconducting circuits,” Physical Review X, vol. 13, p. 041005, Oct. 2023
2023
-
[48]
Probing site- resolved current in strongly interacting superconducting circuit lattices,
B. Du, R. Suresh, S. L ´opez, J. Cadiente, and R. Ma, “Probing site- resolved current in strongly interacting superconducting circuit lattices,” Physical Review Letters, vol. 133, p. 060601, Aug. 2024
2024
-
[49]
Millisecond lifetimes and coherence times in 2D transmon qubits,
M. P. Bland, F. Bahrami, J. G. C. Martinez, P. H. Prestegaard, B. M. Smithamet al., “Millisecond lifetimes and coherence times in 2D transmon qubits,”Nature, vol. 647, no. 8089, pp. 343–348, Nov. 2025
2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.