Closed-Form Analytical Solution for Effective Resistance in Finite 2D Anisotropic Resistor Grids via Jacobi Theta Functions

Ruichao Liu

arxiv: 2604.14192 · v1 · submitted 2026-04-01 · 📡 eess.SP

Closed-Form Analytical Solution for Effective Resistance in Finite 2D Anisotropic Resistor Grids via Jacobi Theta Functions

Ruichao Liu This is my paper

Pith reviewed 2026-05-13 22:22 UTC · model grok-4.3

classification 📡 eess.SP

keywords resistor gridseffective resistanceJacobi theta functionsanisotropicclosed-form solutioninfinity mirrorcircuit analysisVLSI

0 comments

The pith

A closed-form expression using Jacobi theta functions computes effective resistance in finite 2D anisotropic resistor grids to machine precision.

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

The paper derives a closed-form analytical solution for effective resistance between nodes in finite 2D resistor grids where resistances differ along rows and columns. Starting from the infinity mirror technique for boundaries, it finds an exact primitive for the singular integral in the resistance formula and then recasts the infinite mirror sum as a compact expression involving the Jacobi theta function. This yields high accuracy with minimal terms for typical cases. The approach matters for applications in circuit analysis and network science because it replaces slow numerical simulations with direct evaluation.

Core claim

We transform the doubly infinite mirror series into a compact expression using the Jacobi theta function ϑ1 after deriving an exact analytical primitive for the singular integral term R2. This achieves machine precision with only a few terms, and a hybrid localized grid integration corrects residual errors under high anisotropy to maintain mean relative errors below 0.04% against SPICE simulations.

What carries the argument

The Jacobi theta function ϑ1 that compacts the doubly infinite mirror series of the resistance integral after the analytical primitive for R2 is obtained.

If this is right

The analytical expression avoids numerical truncation or polynomial fitting for the resistance calculation.
Under high anisotropy a hybrid method with dynamic numerical cache and localized grid integration is used to eliminate cross-shaped error artifacts.
Numerical experiments confirm mean relative errors below 0.04% compared to SPICE simulations.
The implementation provides an open-source calculator for 2D resistor grids.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This method could allow rapid evaluation of resistance in large grids without computational cost scaling with size.
Similar theta-function summations might apply to other lattice problems in electrostatics or diffusion on networks.
The hybrid remediation suggests a general strategy for combining analytical periodicity with local corrections in anisotropic media.
Applications in VLSI power networks could benefit from exact formulas for optimization.

Load-bearing premise

The infinity mirror technique supplies correct boundary conditions for finite grids and the singular integral R2 admits an exact analytical primitive that remains valid under the subsequent theta-function summation.

What would settle it

Computing the effective resistance for a small 2x2 or 3x3 anisotropic grid using the theta function formula and verifying it matches the exact solution obtained by solving the linear system of Kirchhoff's laws directly.

Figures

Figures reproduced from arXiv: 2604.14192 by Ruichao Liu.

read the original abstract

Computing the effective resistance between nodes in finite discrete resistor grids is a classical problem in circuit analysis with applications in VLSI power delivery network analysis, graph theory, and network science. Recent advances, particularly the infinity mirror technique, provide an elegant physical interpretation for boundary conditions in finite grids. Building upon this foundation, this paper presents a closed-form analytical expression that avoids numerical truncation or polynomial fitting. Our theoretical development proceeds in two steps. First, we derive an exact analytical primitive for the singular integral term $R_2$ within the integral operator $\Omega_\alpha$. Second, we transform the doubly infinite mirror series into a compact expression using the Jacobi theta function $\vartheta_1$. This transformation achieves machine precision with only a few terms. However, under high anisotropy, the pure analytical approximation exhibits a distinct "cross-shaped" residual error. To address this, we introduce a hybrid engineering remediation: a dynamic numerical cache that performs localized grid integration (LGI), combining $O(1)$ speed with exact near-field accuracy. Numerical experiments demonstrate mean relative errors below 0.04% compared to SPICE simulations, eliminating axis-localized error artifacts. To facilitate further research, the implementation of our proposed 2D resistor grid calculator is available at: https://github.com/SeaTheDestiny/2D-Resistor-Grid-Calculator.git.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper turns the mirror series for anisotropic grids into a theta-function form but needs a hybrid numerical patch for high anisotropy, so the closed-form claim is only partial.

read the letter

The new piece is the step that folds the doubly infinite mirror sum into a Jacobi theta_1 expression after finding an exact primitive for the R2 singular integral. That move is not in the earlier mirror-technique papers they cite, and it does cut the number of terms needed for machine precision in the isotropic or mildly anisotropic cases. The code release and the SPICE comparisons are also useful; mean relative error below 0.04 % is a practical result for VLSI-style grids. The hybrid LGI cache they add for the cross-shaped residuals shows they tested the formula carefully enough to notice where it breaks. The soft spot is exactly where the stress-test flagged: when anisotropy gets large the pure analytic expression leaves visible axis-aligned error, so they switch to localized numerical integration near the source. That tells me the claimed exact primitive does not commute with the mirror summation under strong directional imbalance, or the mirror boundary conditions themselves are not fully accurate in that regime. Either way, the headline result is not a truncation-free closed form for arbitrary anisotropy. The derivation steps are not shown in the abstract, but the full text appears to follow the standard integral-operator route without obvious circularity or invented parameters. A reader who needs fast resistance lookups in 2-D grids will get value from the formula plus the code. The work is coherent on its own terms and engages the right prior literature, so it is worth sending out for review even though the analytic claim needs qualification in the high-anisotropy limit.

Referee Report

3 major / 1 minor

Summary. The manuscript claims a closed-form analytical solution for effective resistance between arbitrary nodes in finite 2D anisotropic resistor grids. Building on the infinity-mirror technique, it first derives an exact primitive for the singular integral term R2 inside the operator Ω_α and then converts the resulting doubly infinite mirror series into a compact expression involving the Jacobi theta function ϑ1. The pure analytic form is asserted to reach machine precision with only a few terms; under high anisotropy a hybrid numerical cache (localized grid integration) is added to remove cross-shaped residuals, yielding mean relative errors below 0.04 % versus SPICE.

Significance. A rigorously derived, parameter-free closed-form expression for this classical problem would be a substantial advance for VLSI power-grid analysis and network science, replacing truncation or fitting with an exact, rapidly evaluable formula. The open-source implementation is a positive reproducibility feature. The reported need for a hybrid correction in the high-anisotropy regime, however, indicates that the central analytic claim does not hold universally as stated.

major comments (3)

[Abstract] Abstract: the headline claim of a 'closed-form analytical expression that avoids numerical truncation' is directly contradicted by the introduction of a post-hoc 'hybrid engineering remediation' (dynamic numerical cache + LGI) required to eliminate 'cross-shaped' residuals precisely when anisotropy is high. This regime-dependent breakdown shows that the ϑ1 transformation of the R2 primitive does not remain exact under strong directional imbalance.
[Theoretical development] Theoretical development (first and second steps): no derivation is supplied for the asserted exact antiderivative of R2 nor for the subsequent mapping of the mirror series onto ϑ1. Without these steps or a convergence/error analysis, the machine-precision claim cannot be verified, especially given the documented residuals.
[Abstract] Abstract and numerical experiments: the infinity-mirror boundary conditions are assumed to remain valid for finite grids after substitution of the R2 primitive, yet the cross-shaped artifacts appear exactly when anisotropy is large; this suggests either the mirror technique fails to enforce the correct finite-grid conditions or the analytic continuation does not commute with the summation in that limit.

minor comments (1)

[Numerical experiments] The GitHub link is provided, but the repository contents (test cases, anisotropy ranges, grid sizes) are not described in the text; explicit documentation of the validation suite would strengthen the numerical claims.

Simulated Author's Rebuttal

3 responses · 0 unresolved

Thank you for your thorough review. We address the major comments point by point, agreeing to revisions where the presentation can be improved to better reflect the scope of our analytic solution and the practical hybrid correction.

read point-by-point responses

Referee: [Abstract] Abstract: the headline claim of a 'closed-form analytical expression that avoids numerical truncation' is directly contradicted by the introduction of a post-hoc 'hybrid engineering remediation' (dynamic numerical cache + LGI) required to eliminate 'cross-shaped' residuals precisely when anisotropy is high. This regime-dependent breakdown shows that the ϑ1 transformation of the R2 primitive does not remain exact under strong directional imbalance.

Authors: We thank the referee for highlighting this potential inconsistency in the abstract. The core contribution is the closed-form expression using Jacobi theta functions for the effective resistance, which is exact in the sense that it represents the infinite mirror sum without truncation for the R2 term. However, for practical computation in high anisotropy, the series convergence slows, leading to residuals that we address with a hybrid method for accuracy. This does not contradict the analytic nature but supplements it for robustness. We will revise the abstract to state that the analytic form achieves high precision with few terms, and the hybrid ensures sub-0.04% error in all cases. revision: partial
Referee: [Theoretical development] Theoretical development (first and second steps): no derivation is supplied for the asserted exact antiderivative of R2 nor for the subsequent mapping of the mirror series onto ϑ1. Without these steps or a convergence/error analysis, the machine-precision claim cannot be verified, especially given the documented residuals.

Authors: The derivations for the exact antiderivative of R2 and the mapping to ϑ1 are provided in the theoretical development sections of the full manuscript. To address the concern, we will include more explicit step-by-step derivations, including the integral evaluation and the Poisson summation or equivalent transformation leading to the theta function. Additionally, we will add a convergence analysis demonstrating that the series converges to machine precision with a small number of terms for moderate anisotropy, with the hybrid handling the high-anisotropy case. revision: yes
Referee: [Abstract] Abstract and numerical experiments: the infinity-mirror boundary conditions are assumed to remain valid for finite grids after substitution of the R2 primitive, yet the cross-shaped artifacts appear exactly when anisotropy is large; this suggests either the mirror technique fails to enforce the correct finite-grid conditions or the analytic continuation does not commute with the summation in that limit.

Authors: The infinity-mirror technique correctly enforces the finite boundary conditions by construction, and the substitution of the R2 primitive maintains this validity. The cross-shaped artifacts are due to the directional dependence in the convergence rate of the theta function series under high anisotropy, not a breakdown of the boundary conditions or non-commutativity. The analytic expression remains valid, but numerical evaluation requires care in those regimes, which the hybrid LGI addresses by exact computation near the nodes. We will expand the numerical experiments section to include an analysis of these artifacts and their mitigation. revision: yes

Circularity Check

0 steps flagged

Theta-function compaction of mirror series is a direct identity; minor foundation citation on infinity mirror does not load-bear the final expression

full rationale

The derivation first obtains an exact antiderivative for the singular integral R2 inside Ω_α, then substitutes the known Jacobi ϑ1 summation identity into the doubly infinite mirror series. Both steps are presented as standard mathematical reductions rather than fitted parameters or self-referential definitions. The infinity-mirror boundary condition is invoked from prior literature as an external physical interpretation; even if that citation overlaps with the author, it supplies only the starting series and does not force the subsequent theta compactification. The reported cross-shaped residual under high anisotropy is treated by an auxiliary numerical cache, confirming that the pure analytical claim is not asserted to be exact in every regime and therefore does not collapse into its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the infinity-mirror boundary model and the existence of an exact analytic primitive for the singular integral R2; no free parameters or invented entities are mentioned.

axioms (1)

domain assumption Infinity mirror technique supplies exact boundary conditions for finite grids
Invoked as the foundation for the integral operator Ω_α

pith-pipeline@v0.9.0 · 5543 in / 1205 out tokens · 37980 ms · 2026-05-13T22:22:52.878336+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We derive an exact analytical primitive for the singular integral term R2 ... R2(α)=√α/4π ln(16/π²(α+1))
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

transform the doubly infinite mirror series into a compact expression using the Jacobi theta function ϑ1 via the Jacobi triple product identity

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

[1]

On oscillation hysteresis in a triode generator with two degrees of freedom,

B. van der Pol and N. Wiener, “On oscillation hysteresis in a triode generator with two degrees of freedom,”Philos. Mag., vol. 7, no. 43, pp. 227–258, 1929

work page 1929
[2]

Application of the lattice Green’s function for calculating the resistance of an infinite network of resistors,

J. Cserti, “Application of the lattice Green’s function for calculating the resistance of an infinite network of resistors,”Am. J. Phys., vol. 68, no. 10, pp. 896–906, 2000

work page 2000
[3]

Fast algorithms for effective resistance computation,

S. Köse and E. G. Friedman, “Fast algorithms for effective resistance computation,” inProc. IEEE Int. Symp. Circuits Syst., 2012, pp. 2533–2536

work page 2012
[4]

Efficient algorithms for fast IR drop analysis,

S. Köse and E. G. Friedman, “Efficient algorithms for fast IR drop analysis,” inProc. ACM/IEEE Int. Symp. Low Power Electron. Des., 2012, pp. 255–260

work page 2012
[5]

Effective resistance of finite two-dimensional grids based on infinity mirror technique,

E. Bairamkulov, K. Xu, M. S. Bakir, and E. G. Friedman, “Effective resistance of finite two-dimensional grids based on infinity mirror technique,”IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 67, no. 12, pp. 4581–4593, 2020

work page 2020
[6]

Atkinson and S

K. Atkinson and S. Han,Theoretical Numerical Analysis: A Functional Analysis Frame- work. Springer, 2012

work page 2012
[7]

AnaIR: Green’s function inspired neural network for IR drop prediction,

A. Pandeyet al., “AnaIR: Green’s function inspired neural network for IR drop prediction,” IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., 2025

work page 2025
[8]

Resistance distance,

D. J. Klein and M. Randić, “Resistance distance,”J. Math. Chem., vol. 12, no. 1, pp. 81–95, 1993. 14

work page 1993

[1] [1]

On oscillation hysteresis in a triode generator with two degrees of freedom,

B. van der Pol and N. Wiener, “On oscillation hysteresis in a triode generator with two degrees of freedom,”Philos. Mag., vol. 7, no. 43, pp. 227–258, 1929

work page 1929

[2] [2]

Application of the lattice Green’s function for calculating the resistance of an infinite network of resistors,

J. Cserti, “Application of the lattice Green’s function for calculating the resistance of an infinite network of resistors,”Am. J. Phys., vol. 68, no. 10, pp. 896–906, 2000

work page 2000

[3] [3]

Fast algorithms for effective resistance computation,

S. Köse and E. G. Friedman, “Fast algorithms for effective resistance computation,” inProc. IEEE Int. Symp. Circuits Syst., 2012, pp. 2533–2536

work page 2012

[4] [4]

Efficient algorithms for fast IR drop analysis,

S. Köse and E. G. Friedman, “Efficient algorithms for fast IR drop analysis,” inProc. ACM/IEEE Int. Symp. Low Power Electron. Des., 2012, pp. 255–260

work page 2012

[5] [5]

Effective resistance of finite two-dimensional grids based on infinity mirror technique,

E. Bairamkulov, K. Xu, M. S. Bakir, and E. G. Friedman, “Effective resistance of finite two-dimensional grids based on infinity mirror technique,”IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 67, no. 12, pp. 4581–4593, 2020

work page 2020

[6] [6]

Atkinson and S

K. Atkinson and S. Han,Theoretical Numerical Analysis: A Functional Analysis Frame- work. Springer, 2012

work page 2012

[7] [7]

AnaIR: Green’s function inspired neural network for IR drop prediction,

A. Pandeyet al., “AnaIR: Green’s function inspired neural network for IR drop prediction,” IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., 2025

work page 2025

[8] [8]

Resistance distance,

D. J. Klein and M. Randić, “Resistance distance,”J. Math. Chem., vol. 12, no. 1, pp. 81–95, 1993. 14

work page 1993