The ferromagnetic symmetric-group (SQS_Q) permutation model is the replica statistical model whose ordering transition encodes the entanglement transition of random tensor networks. A recent star-triangle (finite-basis) duality analysis by Ohzeki (arXiv:2607.14917; published as Prog. Theor. Exp. Phys. 2026, ptag162) estimates its critical bond dimension on the honeycomb lattice at Dchex=2.634929344884⋯D_c^{hex}=2.634929344884\cdots and, via a single-bond mutual-duality relation, the paired triangular estimate Dctri=1.475661534848⋯D_c^{tri}=1.475661534848\cdots. Independently, the direct triangular block condition reduces in the replica limit Q→0Q\to 0 to γ+S(D)=0\gamma+\mathcal{S}(D)=0 with S(D)=∑r≥1brD−2r\mathcal{S}(D)=\sum_{r\geq 1}b_r D^{-2r} a factorially divergent asymptotic series whose Farahat-Higman coefficients brb_r are known exactly through r=24r=24; Borel-Pade resummation gives only a diagnostic root 1.48463(10)1.48463(10), which the paper explicitly declines to adopt as a critical-point estimate. The open problem is the full resummation of this series: a closed form or exact representation of S(D)\mathcal{S}(D), or a rigorous summability result with certified error, yielding the definitive root of the direct condition and deciding whether the direct projection independently confirms the mutual-dual paired estimate or whether the roughly 0.6% offset of the current diagnostic reflects a genuine systematic discrepancy between the two finite-basis duality projections. No resolution or adjacent result addressing this question is known; the problem is open.

Background

Replica analyses of random tensor networks and of monitored random quantum circuits map disorder-averaged entanglement moments onto classical statistical models whose elementary degrees of freedom are permutations. In the resulting permutation model, each site ii of a graph G=(V,E)G=(V,E) carries a variable σi\sigma_i in the symmetric group SQS_Q (the group of all permutations of QQ objects, where QQ is the replica index), and the partition function is ZG(D,Q)=∑{σi}∏⟨ij⟩fD(σiσj−1)Z_G(D,Q)=\sum_{\{\sigma_i\}}\prod_{\langle ij\rangle}f_D(\sigma_i\sigma_j^{-1}) with fD(g)=DC(g)f_D(g)=D^{C(g)}, where C(g)C(g) is the number of cycles of the permutation gg. The parameter DD plays a double role: it is a ferromagnetic coupling (one writes D=eKD=e^K, so larger DD means lower effective temperature) and, in the quantum-information application, the bond dimension, i.e. the Hilbert-space dimension carried by each bond of the underlying random tensor network. The replica limit analytically continues QQ from positive integers, and the physical quantity is the coefficient linear in QQ at Q=0Q=0, exactly as in replica treatments of spin glasses. The model's ordering transition separates a disordered phase, corresponding to area-law entanglement of the network, from an ordered phase, corresponding to volume-law entanglement, which is further interpreted as a quantum-error-protecting phase. At Q=2Q=2 the model reduces to the ferromagnetic Ising model with D=e2KIsingD=e^{2K_{\rm Ising}}, a useful normalization check: the Q=2Q=2 critical bond dimensions reproduce the exact triangular and honeycomb Ising values 3\sqrt{3} and 2+32+\sqrt{3}.

Duality is the classical analytic tool for locating two-dimensional critical points: a Fourier-transform relation maps a planar lattice model to a model on the dual lattice, and in self-dual cases the equality of the original and dual principal Boltzmann factors (the weight of the configuration in which neighboring variables agree) identifies the critical point, the Kramers-Wannier argument for the Ising model. For random and replicated systems the principal-factor condition is only a scalar projection of the exact duality relation: predictive, but not exact. The established refinement, developed for random-bond Ising systems on the triangular-hexagonal pair by Nishimori and Ohzeki and extended to regular lattices by Ohzeki, is to impose the condition on a finite basis larger than a bare bond; for the triangular/honeycomb pair the appropriate block is the star-triangle block (a triangle of bonds on one lattice maps to a three-legged star on the dual lattice, the local move underlying the exact triangular-honeycomb Ising duality). For the permutation model, the single-bond replica condition log⁡D=ψ(D)+γ/2\log D=\psi(D)+\gamma/2, with ψ\psi the digamma function (the logarithmic derivative of the Gamma function) and γ\gamma the Euler-Mascheroni constant, gives the bare estimate D0=1.882008326950⋯D_0=1.882008326950\cdots for self-dual lattices; exact real-space renormalization on self-dual hierarchical lattices (self-similar fractal lattices on which the renormalization transformation closes exactly) with extrapolation toward the replica limit supports that estimate.

The source paper extends the finite-basis analysis to the non-self-dual triangular/honeycomb pair. On the honeycomb side, the star-triangle (Y-block) principal-factor condition collapses in the replica limit to the closed equation ψ(D3)=3ψ(D)+γ\psi(D^3)=3\psi(D)+\gamma, whose positive solution is Dchex=2.634929344884⋯D_c^{\rm hex}=2.634929344884\cdots; the single-bond mutual-duality relation between the two mutually dual lattices, log⁡Dtri+log⁡Dhex=ψ(Dtri)+ψ(Dhex)+γ\log D_{\rm tri}+\log D_{\rm hex}=\psi(D_{\rm tri})+\psi(D_{\rm hex})+\gamma, then yields the paired triangular estimate Dtripair=1.475661534848⋯D_{\rm tri}^{\rm pair}=1.475661534848\cdots. The direct triangular face condition is subtler. Its dual principal factor is the constrained three-permutation sum

YQ(Q,D)=1Q!∑g1g2g3=eDC(g1)+C(g2)+C(g3)=∑λ⊢Q(D)λ3hλ2,\mathcal{Y}_Q(Q,D)=\frac{1}{Q!}\sum_{g_1g_2g_3=e}D^{C(g_1)+C(g_2)+C(g_3)}=\sum_{\lambda\vdash Q}\frac{(D)_\lambda^3}{h_\lambda^2},

where the first sum runs over triples g1,g2,g3∈SQg_1,g_2,g_3\in S_Q with g1g2g3=eg_1g_2g_3=e, and the second is a sum over Young diagrams λ\lambda of QQ boxes, with (D)λ=∏□∈λ(D+c(□))(D)_\lambda=\prod_{\square\in\lambda}(D+c(\square)) the content product (c(□)c(\square) is the column index minus the row index of the box) and hλh_\lambda the product of hook lengths. Factoring out the leading behavior gives YQ=D3QΓ(Q+1)[1+∑r≥1cr(Q)D−2r]\mathcal{Y}_Q=\frac{D^{3Q}}{\Gamma(Q+1)}\left[1+\sum_{r\geq 1}c_r(Q)D^{-2r}\right], and the classical Farahat-Higman stability theorem (the structure constants of the center of the SQS_Q group algebra are polynomials in QQ) implies that each cr(Q)c_r(Q) is a polynomial in QQ. Writing br=cr′(0)b_r=c_r'(0), the replica limit of the direct triangular block condition D3Q=YQ(Q,D)D^{3Q}=\mathcal{Y}_Q(Q,D) becomes

γ+S(D)=0,S(D)=∑r≥1brD2r,\gamma+\mathcal{S}(D)=0,\qquad\mathcal{S}(D)=\sum_{r\geq 1}\frac{b_r}{D^{2r}},

with b1=−3/2b_1=-3/2, b2=5/4b_2=5/4, b3=−41/6b_3=-41/6, b4=2751/40b_4=2751/40, b5=−2057/2b_5=-2057/2, and the coefficients known exactly through r=24r=24. The series is factorially divergent at the physical value of DD. Its Borel and Borel-Leroy resummations (the Borel transform divides the coefficients by Γ(r+β+1)\Gamma(r+\beta+1), is extrapolated by Pade rational approximants, and is integrated against e−te^{-t}; approximants with poles on the positive Borel axis are discarded) cluster near D≈1.4846D\approx 1.4846, giving the diagnostic Dtridirect=1.48463(10)D_{\rm tri}^{\rm direct}=1.48463(10). The paper explicitly declines to use this diagnostic as the critical-point estimate, because a Pade diagnostic does not settle the value; the full resummation of the series is left open.

Problem Statement

Determine the full resummation of the direct triangular block condition for the Q→0Q\to 0 symmetric-group permutation model, and thereby the definitive direct estimate of the triangular-lattice critical bond dimension. Concretely, for the asymptotic series S(D)=∑r≥1brD−2r\mathcal{S}(D)=\sum_{r\geq 1}b_rD^{-2r}, where br=cr′(0)b_r=c_r'(0) are the replica-limit coefficients of the constrained three-permutation sum YQ(Q,D)=∑λ⊢Q(D)λ3/hλ2\mathcal{Y}_Q(Q,D)=\sum_{\lambda\vdash Q}(D)_\lambda^3/h_\lambda^2 (exactly known through r=24r=24): obtain a closed form or exact representation for S(D)\mathcal{S}(D), equivalently for ∂Qlog⁡YQ(Q,D)∣Q=0\left.\partial_Q\log\mathcal{Y}_Q(Q,D)\right|_{Q=0}, or failing that a resummation with a certified, rigorous error bound, sufficient to determine the root of γ+S(D)=0\gamma+\mathcal{S}(D)=0 to definitive accuracy (here a root pinned to within 10−310^{-3}, an order of magnitude below the roughly 0.6%≈9×10−30.6\%\approx9\times10^{-3} offset at issue); a determination by resummation must resolve the status of the positive-axis poles seen in Pade approximants of the Borel transform, establishing whether S(D)\mathcal{S}(D) is Borel (or Borel-Leroy) summable at the relevant D≈1.48D\approx 1.48. The answer must state whether this direct estimate confirms the mutual-dual paired estimate Dtripair=1.475661534848⋯D_{\rm tri}^{\rm pair}=1.475661534848\cdots, or whether the roughly 0.6%0.6\% offset already visible in the Borel-Pade diagnostic Dtridirect=1.48463(10)D_{\rm tri}^{\rm direct}=1.48463(10), about ninety times the quoted Pade scatter, reflects a genuine systematic discrepancy between the two finite-basis duality projections: the direct estimate confirms the paired estimate if and only if the definitive direct root agrees with DtripairD_{\rm tri}^{\rm pair} within 10−310^{-3}, and otherwise the offset is a genuine systematic discrepancy. Completion of this problem is the definitive direct determination together with the consistency verdict; the source's broader suggested extensions (an explicit duality analysis of the statistical-mechanical models of random quantum circuits on square, triangular, and honeycomb lattices; other non-self-dual lattices and higher replica numbers; the connection to quantum-error-correction thresholds) are open research directions without sharp completion boundaries, and are outside the scope of this problem.

References

  1. M. Ohzeki, 'Star-triangle duality estimates for triangular and honeycomb permutation models', arXiv:2607.14917 (2026); published as Progress of Theoretical and Experimental Physics 2026, ptag162, DOI 10.1093/ptep/ptag162. https://arxiv.org/abs/2607.14917
  2. M. Ohzeki, 'Duality Analysis in a Symmetric Group and its Application to Random Tensor Network Models', Progress of Theoretical and Experimental Physics 2024, 121A01, DOI 10.1093/ptep/ptae171; preprint arXiv:2310.14140 (2023). https://doi.org/10.1093/ptep/ptae171
  3. R. Ito, T. Matsuo, M. Ohzeki, 'Analysis of Critical Points in a Permutation Model on Hierarchical Lattices by Real-Space Renormalization Group', Journal of the Physical Society of Japan 95, 094003 (2026), DOI 10.7566/jpsj.95.094003; arXiv:2605.25683. https://arxiv.org/abs/2605.25683
  4. M. Ohzeki, 'Cavity method for permutation models on Cayley trees', arXiv:2606.17751 (2026). https://arxiv.org/abs/2606.17751
  5. M. Ohzeki, 'Finite-Basis Duality Estimate for the Surface-Code Threshold under Correlated Bit-Flip Errors', Journal of the Physical Society of Japan 95, 095002 (2026), DOI 10.7566/jpsj.95.095002; arXiv:2606.04723. https://arxiv.org/abs/2606.04723
  6. M. Ohzeki, 'A Minimal Duality Estimate for the Surface-Code Threshold under Nearest-Neighbor Correlated Errors', arXiv:2606.02052 (2026). https://arxiv.org/abs/2606.02052
  7. H. K. Farahat, G. Higman, 'The centres of symmetric group rings', Proceedings of the Royal Society of London Series A 250, 212-221 (1959), DOI 10.1098/rspa.1959.0060. https://doi.org/10.1098/rspa.1959.0060
  8. T. Zhou, A. Nahum, 'Emergent statistical mechanics of entanglement in random unitary circuits', Physical Review B 99, 174205 (2019), DOI 10.1103/PhysRevB.99.174205. https://doi.org/10.1103/PhysRevB.99.174205
  9. R. Vasseur, A. C. Potter, Y.-Z. You, A. W. W. Ludwig, 'Entanglement transitions from holographic random tensor networks', Physical Review B 100, 134203 (2019), DOI 10.1103/PhysRevB.100.134203. https://doi.org/10.1103/PhysRevB.100.134203
  10. Y. Bao, S. Choi, E. Altman, 'Theory of the phase transition in random unitary circuits with measurements', Physical Review B 101, 104301 (2020), DOI 10.1103/PhysRevB.101.104301. https://doi.org/10.1103/PhysRevB.101.104301
  11. S. Choi, Y. Bao, X.-L. Qi, E. Altman, 'Quantum Error Correction in Scrambling Dynamics and Measurement-Induced Phase Transition', Physical Review Letters 125, 030505 (2020), DOI 10.1103/PhysRevLett.125.030505. https://doi.org/10.1103/PhysRevLett.125.030505
  12. H. Nishimori, M. Ohzeki, 'Location of the Multicritical Point for the Ising Spin Glass on the Triangular and Hexagonal Lattices', Journal of the Physical Society of Japan 75, 034004 (2006), DOI 10.1143/JPSJ.75.034004. https://doi.org/10.1143/JPSJ.75.034004
  13. M. Ohzeki, 'Locations of multicritical points for spin glasses on regular lattices', Physical Review E 79, 021129 (2009), DOI 10.1103/PhysRevE.79.021129. https://doi.org/10.1103/PhysRevE.79.021129