Pith. sign in

REVIEW 2 major objections 5 minor 29 references

Surface Variables Description of Axion Topological Materials

T0 review · 2 major / 5 minor · reviewed 2026-07-14 · grok-4.5

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

Pith's one-line read The entire quasi-static magnetoelectric response of a topological insulator is carried by induced surface charge and current densities that obey linear integral equations on the material boundary.

desk verdict Clean transplant of soft-matter surface variables to quasi-static axion electrodynamics; solid analytics and a working BEM for non-trivial shapes, with the usual sharp-θ idealizations openly stated. read the letter →

arxiv 2607.10509 v1 pith:RX6P5FOA submitted 2026-07-12 cond-mat.mtrl-sci cond-mat.otherhep-thphysics.app-phquant-ph

Francisco J. Solis , Hugo García-Compeán This is my paper · ORCID
classification cond-mat.mtrl-scicond-mat.otherhep-thphysics.app-phquant-ph
keywords axionelectrodynamicstopologicalinsulatorssurfacevariablesmagnetoelectricresponseboundaryintegralequationsquasi-staticregimefinite-elementmethod
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 shows that, for a topological insulator whose axion field jumps discontinuously from a constant value inside to zero outside, the whole macroscopic quasi-static electromagnetic response under external static fields can be reduced to two surface densities: an induced surface charge and an induced surface current. Those densities satisfy a closed set of linear surface integral equations. Once the surface densities are known, the electric and magnetic fields everywhere in space are recovered from ordinary free-space Green-function integrals. Analytic solutions exist for spheres and other axial geometries; for arbitrary shapes the same equations are discretized by a boundary finite-element method and solved iteratively. The surface formulation is also recovered from a variational principle written solely in the surface variables. The reduction lowers the computational dimension from three to two and makes non-trivial geometries (tori, hollow spheres, multiple interfaces) routine.

What carries the argument

The surface integral operator X that maps the pair of surface densities (ω_s, s_s) into the self-fields they produce on the same surface; the governing equation is the linear relation w + Xw = w_D, where w_D is the direct contribution of the external free sources.

What would settle it

Solve the same sphere-plus-external-field problem by an independent Green-function method that retains bulk polarization; any systematic discrepancy larger than the reported one-percent numerical error would falsify the pure-surface reduction.

Watch

Extended reading notes

Core claim

A complete macroscopic description of the quasi-static magnetoelectric response of a topological insulator is given by the induced surface charge and current densities ω_s and s_s that satisfy the linear surface integral equations ω_s + Δ heta n · ( abla imes ∫ G s_s dS) = ω_D and s_s + Δ heta n imes ( abla ∫ G ω_s dS) = s_D. Once those surface densities are known, the bulk fields are recovered everywhere by free-space Coulomb and Biot-Savart integrals.

Load-bearing premise

The axion field is assumed to be strictly piecewise constant, so that every induced source collapses exactly onto a sharp surface with vacuum permittivity and permeability on both sides.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The manuscript reformulates the quasi-static magnetoelectric response of axion topological insulators in terms of induced surface charge and current densities ω_s and s_s. Starting from a Hamiltonian that includes the axion term θ E·B, the authors show that the modified Maxwell equations reduce to free-space Poisson problems sourced by free charges/currents plus induced surface densities localized where ∇θ is nonzero. These densities obey the linear surface integral equations (40)–(41) (abstractly w + Xw = w_D). Analytic solutions are obtained for spheres under axial external fields and recover the known Green-function results of Martín-Ruiz et al. A boundary-element discretization with two empirical refinements (near-field analytic Green integration and neighbor smoothing) is developed and applied to spheres, hollow spheres, and tori, with reported agreement to analytics better than 1 %. A surface-variable Hamiltonian is also shown to recover the same integral equations by extremization.

Significance. If the surface reduction is accepted, the work supplies a practical, dimensionally reduced computational route for the magnetoelectric response of topological insulators of arbitrary shape, together with an energy functional that may be useful for force calculations. The analytic recovery of published spherical solutions and the explicit numerical validation against those solutions are genuine strengths; the method is free of adjustable parameters once Δθ is fixed. The approach is a natural extension of established soft-matter surface-variable techniques to axion electrodynamics and is therefore of interest to both the topological-materials and computational-electromagnetism communities. The restriction to piecewise-constant θ and vacuum ε, μ is standard in the literature the paper cites, so the results are immediately comparable to existing Green-function calculations while opening the door to more complex geometries.

major comments (2)
  1. Section 6, paragraph after Eq. (68): the two empirical refinements (analytic near-field Green integration and neighbor smoothing) are essential for the reported <1 % accuracy, yet no systematic convergence study with respect to triangle number, smoothing stencil, or near-field cutoff is provided. A short table or plot quantifying residual error versus N_T for the sphere (where the exact solution is known) would make the numerical claims load-bearing rather than anecdotal.
  2. Section 3, Eqs. (27)–(30) and the subsequent energy evaluation (34): the surface-variable Hamiltonian contains several quadratic terms with negative signs that arise from the elimination of Lagrange multipliers. While the authors correctly note that the on-shell energy recovers the standard 1/2(E^{2}+B^{2}) form, it is not shown that the functional is bounded from below off-shell or that the iterative scheme of Sec. 6 is a descent method for that functional. Clarifying the variational status of the iteration would strengthen the claim that a variational principle underlies the surface equations.
minor comments (5)
  1. Abstract and Sec. 1: the phrase “axion topological materials” is used interchangeably with “topological insulators”; a single consistent terminology would avoid confusion with Weyl semimetals (where θ is not piecewise constant).
  2. Eq. (1) and the units paragraph of Sec. 2: the choice c = ε_{0} = μ_{0} = 1 is standard, but the numerical value θ_{0} ≈ 1/137 is introduced only in Sec. 7; stating the convention once in Sec. 2 would help the reader.
  3. Figs. 1–6: the color-scale bars and arrow lengths are not quantified; adding a single representative magnitude (e.g., max |ω_s| or |s_s|) would make the plots more informative.
  4. Typographical: “we show that that a variational principle” (abstract); “paramagnetic constants” should be “permeability” (Sec. 2); occasional missing spaces after periods.
  5. References: the soft-matter surface-variable papers [21,22] are cited, but a brief pointer to the classic boundary-element literature for magnetostatics would help readers outside that community.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: surface integral equations are derived from Hamiltonian extremization and independently recover published Green-function analytics; numerics are validated rather than fitted.

full rationale

The derivation chain is self-contained. Section 2 starts from a quasi-static Hamiltonian (1) whose extremization yields the modified Maxwell equations (8)–(9). Section 3 defines induced sources ω = - abla heta·B and s = abla heta imes E, rewrites the fields via free-space Green functions (20)–(24), and obtains the integral equations (25)–(26) for those sources; the same equations reappear as stationarity conditions of a surface-variable Hamiltonian (27)–(30). For piecewise-constant heta the equations collapse to the linear surface system (40)–(41). Analytic solutions for spheres (Sec. 5) are obtained by Legendre expansion and are shown to coincide with the independent Green-function results of Martín-Ruiz et al. (ref. [15]), not with any prior result of the present authors. The boundary-element discretization (Sec. 6) is validated by recovering those same analytics to <1 % error before being applied to tori and hollow spheres; no free parameters are adjusted to force agreement. Self-citations to the authors’ earlier soft-matter surface-variable papers supply only the methodological template and are not load-bearing for the axion-specific claims. No self-definitional loop, fitted-input-as-prediction, uniqueness theorem imported from the authors, or renaming of a known empirical pattern is present.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central claim rests on standard axion electrodynamics for topological insulators plus the modeling choice that θ is piecewise constant and material constants equal vacuum. No new particles or free fit parameters enter the analytic claim; the numerical examples simply adopt the conventional value θ_0 ≈ α.

free parameters (1)
  • θ_0 (numerical examples) = 0.0073 ≈ 1/137
    Set by hand to the fine-structure constant 0.0073 for all numerical demonstrations; not fitted to data but chosen as a conventional small-contrast value.
assumptions (5)
  • domain assumption Axion field θ is piecewise constant (θ = θ_0 inside the topological insulator, 0 outside), so ∇θ is a surface delta function.
    Stated in Introduction and Sec. 3–4; converts bulk axion terms into pure surface sources.
  • domain assumption Relative permittivity and permeability of the insulator equal vacuum values (ε_r = μ_r = 1).
    Explicitly adopted in Sec. 2 so that only axion-induced surface sources appear.
  • domain assumption Quasi-static regime: all fields, charges and currents are time-independent; radiation is absent.
    Stated throughout; allows Hamiltonian extremization without time derivatives.
  • standard math Free-space Green function G = 1/(4π|x−y|) and Coulomb/Biot-Savart reconstruction of fields from total sources.
    Used in Eqs. (20)–(24) and after the surface solve.
  • domain assumption Hamiltonian density (1) correctly encodes quasi-static axion electrodynamics with matter kinetic term.
    Starting point of Sec. 2; derived from a Lagrangian in the appendices.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Surface Variables Description of Axion Topological Materials." pith.science (2026). https://pith.science/paper/RX6P5FOA

@misc{pith2026260710509,
  author       = {Pith},
  title        = {Pith review of: Surface Variables Description of Axion Topological Materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RX6P5FOA}},
  note         = {Machine review of arXiv:2607.10509}
}
read the original abstract

We study the response of axion topological materials under the presence of external static electric and magnetic fields. We focus on the macroscopic quasi-static magnetoelectric response of topological insulators. We use techniques based on surface variables that have been previously employed in soft condensed matter problems for the description of heterogeneous systems composed of multiple homogeneous materials. A complete description of the whole system is written in terms of surface degrees of freedom, which in this case correspond to effective surface charge and surface current densities. We obtain a set of integral equations satisfied by these variables. We present exact analytic solutions for axial-symmetric cases. We develop a numerical method for the solution of the integral equations by means of a boundary finite element method. We apply the numerical method to topological insulators with different geometries such as spheres, hollow spheres and toroidal surfaces. We show that that a variational principle can be used to recover the surface variables equations.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references · 4 canonical work pages

  1. [1]

    Reviews of modern physics82(4), 3045–3067 (2010)

    Hasan, M.Z., Kane, C.L.: Colloquium: topological insulators. Reviews of modern physics82(4), 3045–3067 (2010)

    2010

  2. [2]

    Qi, X.L., Zhang, S.C.: Topological insulators and superconductors. Rev. Mod. Phys.83(4), 1057–1110 (2011) https://doi.org/10.1103/RevModPhys.83.1057 arXiv:1008.2026 [cond-mat.mes-hall]

    arXiv doi:10.1103/revmodphys.83.1057 2011

  3. [3]

    Yan, B., Felser, C.: Topological Materials: Weyl Semimetals. Ann. Rev. Condensed Matter Phys.8, 337–354 (2017) https://doi.org/10.1146/ annurev-conmatphys-031016-025458 arXiv:1611.04182 [cond-mat.mtrl-sci]

    arXiv 2017

  4. [4]

    Reviews of Modern Physics88(3), 035001 (2016)

    Witten, E.: Fermion path integrals and topological phases. Reviews of Modern Physics88(3), 035001 (2016)

    2016

  5. [5]

    La Rivista del Nuovo Cimento39(7), 313–370 (2016)

    Witten, E.: Three lectures on topological phases of matter. La Rivista del Nuovo Cimento39(7), 313–370 (2016)

    2016

  6. [6]

    arXiv preprint arXiv:2111.07290 (2021)

    Planelles, J.: Axion electrodynamics in topological insulators for beginners. arXiv preprint arXiv:2111.07290 (2021)

    arXiv 2021

  7. [7]

    Journal of Applied Physics129(14) (2021)

    Sekine, A., Nomura, K.: Axion electrodynamics in topological materials. Journal of Applied Physics129(14) (2021)

    2021

  8. [8]

    Physical review letters 58(18), 1799 (1987)

    Wilczek, F.: Two applications of axion electrodynamics. Physical review letters 58(18), 1799 (1987)

    1987

Show all 29 references
  1. [9]

    Physics Letters B86(3-4), 283–287 (1979)

    Witten, E.: Dyons of charge eθ/2π. Physics Letters B86(3-4), 283–287 (1979)

    1979

  2. [10]

    Physical Review D—Particles, Fields, Gravitation, and Cosmology84(10), 105008 (2011)

    Canfora, F., Rosa, L., Zanelli, J.:θterm in a bounded region. Physical Review D—Particles, Fields, Gravitation, and Cosmology84(10), 105008 (2011)

    2011

  3. [11]

    Physical Review D—Particles, Fields, Gravitation, and Cosmology85(8), 085024 (2012)

    Huerta, L., Zanelli, J.: Optical properties of aθvacuum. Physical Review D—Particles, Fields, Gravitation, and Cosmology85(8), 085024 (2012)

    2012

  4. [12]

    Physical Review D55(11), 6760 (1997)

    Colladay, D., Kosteleck` y, V.A.: CPT violation and the standard model. Physical Review D55(11), 6760 (1997)

    1997

  5. [13]

    Physical Review D58(11), 116002 (1998) 24

    Colladay, D., Kosteleck` y, V.A.: Lorentz-violating extension of the standard model. Physical Review D58(11), 116002 (1998) 24

    1998

  6. [14]

    Physical Review D92(12), 125015 (2015)

    Mart´ ın-Ruiz, A., Cambiaso, M., Urrutia, L.: Green’s function approach to Chern- Simons extended electrodynamics: An effective theory describing topological insulators. Physical Review D92(12), 125015 (2015)

    2015

  7. [15]

    Physical Review D93(4), 045022 (2016)

    Mart´ ın-Ruiz, A., Cambiaso, M., Urrutia, L.: Electro-and magnetostatics of topo- logical insulators as modeled by planar, spherical, and cylindricalθboundaries: Green’s function approach. Physical Review D93(4), 045022 (2016)

    2016

  8. [16]

    Physical Review D94(8), 085019 (2016)

    Mart´ ın-Ruiz, A., Cambiaso, M., Urrutia, L.F.: Electromagnetic description of three-dimensional time-reversal invariant ponderable topological insulators. Physical Review D94(8), 085019 (2016)

    2016

  9. [17]

    Physical Review B107(24), 245103 (2023)

    Bonilla, D.A., Casta˜ no-Yepes, J.D., Mart´ ın-Ruiz, A., Munoz, E.: Electromag- netic coupling and transport in a topological insulator–graphene heterostructure. Physical Review B107(24), 245103 (2023)

    2023

  10. [18]

    Physical Review B109(19), 195136 (2024)

    Barredo-Alamilla, E., Bobylev, D.A., Gorlach, M.A.: Axion electrodynamics without Witten effect in metamaterials. Physical Review B109(19), 195136 (2024)

    2024

  11. [19]

    Computer Physics Communications241, 80–91 (2019)

    Nguyen, T.D., Li, H., Bagchi, D., Solis, F.J., Cruz, M.O.: Incorporating surface polarization effects into large-scale coarse-grained molecular dynamics simulation. Computer Physics Communications241, 80–91 (2019)

    2019

  12. [20]

    Journal of electromagnetic waves and applications3(1), 1–15 (1989)

    Harrington, R.F.: Boundary integral formulations for homogeneous material bodies. Journal of electromagnetic waves and applications3(1), 1–15 (1989)

    1989

  13. [21]

    Physical review letters 109(22), 223905 (2012)

    Jadhao, V., Solis, F.J., De La Cruz, M.O.: Simulation of charged systems in het- erogeneous dielectric media via a true energy functional. Physical review letters 109(22), 223905 (2012)

    2012

  14. [22]

    Physical Review E—Statistical, Nonlinear, and Soft Matter Physics88(2), 022305 (2013)

    Jadhao, V., Solis, F.J., De La Cruz, M.O.: Free-energy functionals of the elec- trostatic potential for Poisson-Boltzmann theory. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics88(2), 022305 (2013)

    2013

  15. [23]

    Physical Review E—Statistical, Nonlinear, and Soft Matter Physics88(5), 053306 (2013)

    Solis, F.J., Jadhao, V., Cruz, M.: Generating true minima in constrained variational formulations via modified lagrange multipliers. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics88(5), 053306 (2013)

    2013

  16. [24]

    Mart´ ın-Ruiz, A., Cambiaso, M., Urrutia, L.F.: Electromagnetic description of three-dimensional time-reversal invariant ponderable topological insulators. Phys. Rev. D94(8), 085019 (2016) https://doi.org/10.1103/PhysRevD.94.085019 arXiv:1610.00576 [cond-mat.str-el]

    2016 doi arXiv

  17. [25]

    Mart´ ın-Ruiz, A., Cambiaso, M., Urrutia, L.F.: Electromagnetic fields induced by an electric charge near a Weyl semimetal. Phys. Rev. B99(15), 155142 (2019) https://doi.org/10.1103/PhysRevB.99.155142 arXiv:1809.04210 [hep-th] 25

    2019 doi arXiv

  18. [26]

    The European Physical Journal E47(9), 60 (2024)

    Solis, F.J., Jadhao, V.: Conduction in heterogeneous systems in the low-frequency regime: variational principles and boundary integral equations. The European Physical Journal E47(9), 60 (2024)

    2024

  19. [27]

    Ibarra-Meneses, M., Mart´ ın-Ruiz, A.: Electromagnetic radiation mediated by topological surface states. Phys. Rev. D113(7), 076002 (2026) https://doi.org/ 10.1103/2y93-2fmn arXiv:2603.26495 [hep-th]

    2026 doi

  20. [28]

    Ibarra-Meneses, M., Mart´ ın-Ruiz, A.: Radiation of a short linear antenna above a topologically insulating half-space. Eur. Phys. J. Plus139(2), 144 (2024) https:// doi.org/10.1140/epjp/s13360-024-04931-8 arXiv:2401.11285 [cond-mat.mes-hall]

    2024 doi arXiv

  21. [29]

    CRC Press, Boca Raton, USA (2019) 26

    Schwinger, J., DeRaad Jr, L.L., Milton, K., Tsai, W.-y.: Classical Electrodynam- ics. CRC Press, Boca Raton, USA (2019) 26

    2019

Pith tools

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