{"id":"e22a5c7f-8190-4c5d-8add-9e567c9b7dab","arxiv_id":"2607.10509","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"Axion topological insulators are fully described by surface charge and current densities satisfying integral equations solvable analytically for spheres and numerically for complex shapes.","lead":"This paper reformulates the static magnetoelectric response of topological insulators entirely in terms of induced surface charge and current densities that obey coupled integral equations. The approach yields analytic solutions for spheres and a boundary finite-element scheme for tori and hollow spheres, giving a practical route to fields in complex geometries.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the reader's already-flagged modeling idealizations.","rationale":"The paper’s strongest claim is a clean, inspectable reduction of quasi-static axion electrodynamics to surface integral equations plus free-space reconstruction. Analytic limits match prior Green-function results, the variational surface Hamiltonian recovers the same equations, and the numerical scheme reproduces those limits to <1 % on spheres before being applied to tori and hollow spheres. The only material limitation is the idealization already identified by the reader (sharp θ jump, vacuum material constants). Because that limitation is explicit and does not undermine the internal correctness of the surface formalism as written, no further adjustment to the CONDITIONAL verdict is warranted. The suggested concrete test simply reconfirms the analytic core without relying on the series machinery of Sec. 5.","tokens_in":15438,"tokens_out":478,"duration_ms":4621,"concrete_test":"Re-derive the n=1 amplitudes (63)–(64) for a sphere in uniform E_0, B_0 from the continuous integral equations (40)–(41) without using the Legendre expansion of Sec. 5; if the closed-form D = 1 + (2/9)Δθ^{2} is recovered, the surface reduction is confirmed independently of the series method.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the quasi-static magnetoelectric response is completely captured by surface densities ω_s, s_s obeying the linear integral equations (40)–(41), from which bulk fields are recovered via the free-space Green function—holds under the paper’s stated premises. The derivation from the Hamiltonian (Sec. 2–3), the reduction to surface sources via ∇θ = δ(x_n)Δθ n (Eq. 35), the recovery of known spherical analytic solutions (Sec. 5), and the boundary-element discretization (Sec. 6) are internally consistent. The reader’s weakest assumption (strictly piecewise-constant θ and vacuum ε, μ) is the genuine modeling boundary of the pure-surface reduction; it is openly stated and standard in the cited Green-function literature. No additional load-bearing inconsistency, hidden singularity, or unstated assumption that would invalidate the surface equations themselves was found.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","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.","tokens_in":15607,"tokens_out":1155,"duration_ms":10833,"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":[{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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).","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"Typographical: “we show that that a variational principle” (abstract); “paramagnetic constants” should be “permeability” (Sec. 2); occasional missing spaces after periods.","section":null},{"comment":"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.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, self-contained contribution that sits comfortably within the scope of a condensed-matter or computational-physics journal. The modeling idealizations (piecewise-constant θ, vacuum ε/μ) are the same ones used by the Green-function papers it recovers, so novelty lies in the surface-variable reformulation and the numerical capability for non-spherical shapes rather than in new physics. I see no citation or priority issues. Minor revision is appropriate; the two major points are fixable without new theory."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a practical methods paper that does what it claims. The authors take the surface-charge/current integral-equation machinery they (and others) already use for heterogeneous dielectrics and apply it to piecewise-constant axion electrodynamics. The new pieces are the coupled surface operator X that mixes ω_s and s_s, the closed-form multipole amplitudes for spheres under point-charge or uniform-field drive, and a boundary-element solver that they actually run on tori and hollow spheres.\n\nWhat works: the Hamiltonian reduction to surface sources is clean, the spherical analytics recover the known Green-function results of Martín-Ruiz et al. to <1 %, and the numerical examples (off-center monopole, torus under E or B) show the method can handle disconnected surfaces and non-trivial topology without drama. The variational principle that recovers the same integral equations is a nice consistency check. Citations to both the soft-matter literature and the axion-TI Green-function papers are appropriate; no circularity.\n\nSoft spots are real but limited. Everything rests on a sharp θ jump and vacuum ε, μ, so the pure-surface reduction is exact only under those idealizations (standard in the papers they cite, and they say so). The two empirical fixes (near-field analytic Green integration + neighbor smoothing) are ad-hoc; no code or systematic convergence study is shipped. Those are engineering details, not conceptual holes.\n\nWho it is for: people who need magnetoelectric fields for irregular TI geometries or who want a lower-dimensional alternative to volume Green functions. It is not a foundational advance, but it is usable progress inside the subfield.\n\nI would send it to referees. The math is inspectable, the central claim holds under the stated premises, and the numerical demonstrations are concrete enough to be useful.","headline":"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.","tokens_in":16229,"tokens_out":488,"would_cite":true,"duration_ms":7051,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"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.","keywords":["axion electrodynamics","topological insulators","surface variables","magnetoelectric response","boundary integral equations","quasi-static regime","finite-element method"],"falsifier":"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.","tokens_in":16325,"feed_emoji":"⚡","tokens_out":670,"duration_ms":6613,"temperature":0.7,"pith_summary":"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.","feed_headline":"Surface densities alone capture a topological insulator's magnetoelectric response","feed_subtitle":"Linear integral equations on the boundary recover all bulk fields for spheres, tori and hollow shells","key_machinery":"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.","core_discovery":"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 + Δ\theta n · (\nabla \times ∫ G s_s dS) = ω_D and s_s + Δ\theta n \times (\nabla ∫ 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.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Surface charge and currents fully encode TI magnetoelectric response","Boundary densities alone recover bulk fields in axion topological materials","Linear surface equations capture full magnetoelectric response of TIs","Induced surface densities determine all quasi-static fields in topological insulators","Surface variables alone describe axion materials under electric and magnetic fields"],"cache_read_input_tokens":128,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Surface charge and currents fully encode TI magnetoelectric response","Boundary densities alone recover bulk fields in axion topological materials","Linear surface equations capture full magnetoelectric response of TIs","Induced surface densities determine all quasi-static fields in topological insulators","Surface variables alone describe axion materials under electric and magnetic fields"]},"model":"grok-4.5","effort":"low","cost_usd":0.004602,"raw_usage":{"total_tokens":1293,"prompt_tokens":741,"num_sources_used":0,"completion_tokens":85,"cost_in_usd_ticks":46020000,"prompt_tokens_details":{"text_tokens":741,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":467,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":741,"tokens_out":85,"duration_ms":5060,"temperature":1.0,"reasoning_tokens":467,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T11:10:13.238797+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}