REVIEW 3 major objections 4 minor 28 references
Asymptotics for a class of planar orthogonal polynomials and truncated unitary matrices
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read This paper proves explicit asymptotics for planar orthogonal polynomials and derives the full moment expansion and central limit theorem for |det(B_n-x)| of truncated unitary matrices.
desk verdict A strong, likely correct Riemann–Hilbert paper with one real but fillable gap in Lemma 2.1 and a definite typo in the CLT proof; worth a serious referee. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The key step is Lemma 2.1: planar orthogonality with weight (1-|z|^2)^{$\alpha$-1}|z-x|^gamma is equivalent to contour orthogonality with weight w(z)=(1-$x^{2}$ z)^{$\alpha$+gamma/2}((z-1)/z)^{gamma/2} on a loop enclosing [0,1]. This contour orthogonality is encoded in a 2x2 Riemann-Hilbert problem. A g-function built from the potential V(z)=log z - c log(1-$x^{2}$ z), with level curves Re phi = constant, produces normalized jumps; after lens opening and a global parametrix, the remaining Riemann-Hilbert problem is small-norm, so the orthogonal polynomials are expressed as contour integrals. These integrals are evaluated through an incomplete gamma function identity, yielding explicit three-region asymptotics without constructing a local parametrix. A differential identity for d/dx log R_gamma(x) in terms of the Riemann-Hilbert matrix entries, together with the explicit radially symmetric value at x=0, completes the moment formula.
What would settle it
Simulate truncated unitary matrices for moderate n, for example n=50 to 500 with mu=1/2, x=0.5, and gamma=1, and compare the empirical E(|det(B_n-x)|^gamma) with the theorem's leading term, checking that the relative error decays like O(1/n). Separately, for gamma in (-2,0), numerically compare the planar integral defining the orthogonality in (1.2) with the contour integral in (2.2) for low-degree polynomials to verify Lemma 2.1 directly.
Extended reading notes
Core claim
The central claim is that in the strong non-unitarity regime, with fixed gamma satisfying Re gamma > -2 and fixed x in [0,$\sqrt$(mu)), the moments satisfy E(|det(B_n-x)|^gamma) = $n^{{gamma^2/8}}$ $mu^{{gamma n/2}}$ ((1-mu)/(1-$x^{2}$))^{$\alpha$ gamma/2} C_{gamma,mu}(x)(1+O(1/n)) with C_{gamma,mu}(x) explicitly given in terms of the Barnes G-function. This is proved through uniform asymptotics of the monic planar orthogonal polynomials in the exterior, interior, and near z=1; near z=1 the behaviour is governed by an incomplete gamma function. A corollary is the convergence in distribution of (log|det(B_n-x)| - kappa_1/2)/((1/2)$\sqrt$(log n)) to a standard normal, where kappa_1 is the explicit centering term. In a separate double-scaling regime $x^{2}$=1-v/n with $\alpha$ fixed, the same quantity is expressed through a $\sigma$-Painleve V function following the theory of Toeplitz determinants with Fisher-Hartwig singularities.
Load-bearing premise
The whole construction rests on Lemma 2.1, the reduction of planar orthogonality to contour orthogonality; for the singular range gamma in (-2,0) the Green's theorem step near z=x is stated but the details are omitted, so if that reduction fails, the later theorems lose their foundation.
Editorial extensions
If this is right
- The moment generating function of log|det(B_n-x)| now has a full leading-order expansion, so cumulants of the logarithm of the characteristic polynomial can be extracted order by order in the strong non-unitarity limit.
- A central limit theorem holds for log|det(B_n-x)| for every fixed x in the bulk of the limiting spectrum, generalizing the previously known determinant case x=0.
- The zeros of the planar orthogonal polynomials accumulate on the explicit level curve Gamma_1, and the polynomials have different algebraic, rational, and incomplete-gamma behaviours in the exterior, interior, and near-boundary regions.
- In the double-scaling weak regime, the moments are expressed in terms of a sigma-Painleve V function, giving a concrete connection between truncated unitary matrices and the Fisher-Hartwig/Toeplitz theory.
- The integer moment case gamma=2k reproduces results previously obtained by duality arguments, providing a consistency check for the new Riemann-Hilbert approach.
Reading between the lines
- One could test whether the same parametrix-free steepest descent scheme applies to planar weights with several point charges or with other radially symmetric base weights, which would give complete asymptotics beyond the Gaussian and unit-disc cases.
- The variance (1/4)log n seen in the central limit theorem is the signature of a log-correlated field; comparing higher cumulants with Gaussian multiplicative chaos predictions would be a natural numerical check.
- The differential identity used here might be iterated to compute subleading coefficients in the O(1/n) correction, and it may connect the strong-regime formula to the sigma-Painleve V expression across the critical curve x=sqrt(mu).
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops the large-n asymptotics of planar orthogonal polynomials for the measure dµ(z)=(1-|z|^2)^{α-1}|z-x|^γ 1_{|z|<1} d^2z in the strong non-unitarity regime n/N → µ̃∈(0,1). The authors convert the planar orthogonality into a contour orthogonality, set up a Riemann-Hilbert problem, and then apply Deift-Zhou steepest descent. The main results are: all-region asymptotics for the (rescaled) monic orthogonal polynomials (Theorem 1.1), an asymptotic expansion for the moments E(|det(B_n-x)|^γ) with explicit constants and uniform error O(1/n) (Theorem 1.4), a central limit theorem for log|det(B_n-x)| (Corollary 1.6), and a double-scaling result near x=1 expressed through a σ-Painlevé V solution (Theorem 1.8).
Significance. If the technical gaps are completed, this is a substantial contribution. The paper gives the first all-region asymptotic description for this non-Gaussian planar ensemble with a point charge, and it does so with explicit, parameter-free constants and no fitted parameters. The resulting moment formula and CLT for truncated unitary matrices are concrete and independently checkable. The method also avoids local parametrices in the strong regime, which is a useful structural improvement over earlier Gaussian-weight treatments, and the weak-regime Painlevé V connection extends the Toeplitz-determinant results of Claeys–Its–Krasovsky to the truncated-unitary partition function. The authors are explicit about error terms and about the range of parameters claimed, and they correctly identify the singular case γ∈(-2,0) as requiring separate treatment.
major comments (3)
- [Appendix A, Lemma 2.1] The proof of Lemma 2.1 explicitly omits the singular case γ∈(-2,0). After equation (A.4) the text states: 'for γ∈(-2,0) there is singular behavior near z=x, but one can still verify the statement. We omit the details.' This is a load-bearing gap, because Theorems 1.1 and 1.4 both claim uniform validity for Reγ>-2, which includes this interval. The Green's-theorem reduction requires a justification that the integration-by-parts step applies to the unbounded integrand p_j(z)(z-x)^{γ/2}h_k(z, \bar z, x) near z=x, and a proof that the small-circle boundary contribution around x vanishes. Since the contour orthogonality (2.2), the Riemann-Hilbert problem (2.10)-(2.13), and all subsequent steepest-descent analysis depend on Lemma 2.1, this omission affects the foundation of the paper. Please supply the missing limiting argument for -2<γ<0, or explicitly restrict the main theorems to Reγ>0 if the singular range is not intended to be covered.
- [§3.5, equation (3.23)] The expansion (3.23) for R(z) is only sketched. The Neumann series (3.33)-(3.34) is written down, but the step 'expanding up to order 1 and bounding the remaining terms' is not supported by explicit estimates for the iterated Cauchy operators C_{Σ_{r1}}, C_{Σ_{r2}}. The case r=r2=z0 is also delicate: there φ_r(r2)=0 and the RHP is not directly in the small-norm setting; the rescaling R_ε is mentioned but its analysis is not provided. The later estimates in §3.6-3.7, Proposition 3.4, and ultimately Theorem 1.1 rely on (3.36) and the claimed exponential smallness of the error. A complete proof of (3.23) with explicit norm bounds is therefore needed.
- [§4, proof of Theorem 1.4] The proof of uniformity in Theorem 1.4 is not fully justified. After treating γ=0 in (4.5)-(4.7), the text says 'The general case Reγ>-2 is essentially the same' and asserts that the added factors depend continuously on x and γ. This skips the behavior of the error as γ approaches -2 and as x approaches √µ̃-δ, exactly where the exponent difference φ(r1)-φ(z0) degenerates and the factors of h_γ grow. Since the theorem claims uniformity over compact subsets of Reγ>-2 and over [0,√µ̃-δ], the omitted continuity argument is load-bearing. Please give the explicit dependence of the error term in (4.7) on γ and x, or state uniformity only over fixed compact subsets of (γ,x) that avoid the degenerate limits.
minor comments (4)
- [§3.6, equation (3.42)] The error term O(e^{k(2φ(t)-φ(z0))}) appears to contain an undefined constant k; presumably this should be n, matching the rest of the section. Please clarify.
- [§4] There is a typo: 'preceeding' should be 'preceding'.
- [Theorem 1.4] The theorem states uniformity for 0≤x<√µ̃-δ, while the measure (1.1) and the polynomial construction are stated for x>0 in the introduction. The case x=0 is computed separately in the proof and in Remark 1.3, but the theorem statement should make explicit whether x=0 is included in the uniform statement.
- [Lemma 3.2] The domain of the coefficient functions c_m(z) is written as C\{γ_t}; this should presumably be C\γ_t, i.e. the complement of the curve, not a set minus a set of points.
Circularity Check
No circular reduction found; the only flagged weakness is an omitted singular-case verification in Lemma 2.1, which is a proof gap rather than circularity.
full rationale
The derivation chain is: planar orthogonality (1.2) → Lemma 2.1 contour orthogonality (2.2) with weight (2.1) → Riemann-Hilbert problem (2.10)–(2.13) → steepest descent analysis yielding Theorem 1.1 → differential identity (2.15) → integration from x=0 using the exact radial value R_γ(0) to obtain Theorem 1.4. Each link is explicit and parameter-free. The weight (2.1) is obtained by algebraically rewriting the planar integral in Appendix A, not by assuming the asymptotics it feeds. The RHP solution (2.13) is written in terms of the same orthogonal polynomials, and the small-norm/steepest-descent estimates depend only on the analytic properties of g and φ. The constant R_γ(0) is computed exactly from the rotationally symmetric case x=0, and no fitted constant is introduced. External inputs such as [14], [16], [25], and [26] are independent results whose assumptions do not include the strong-regime theorem proved here. Self-citations to [15] are used for a similar contour reduction and for product asymptotics, but Lemma 2.1's proof is reproduced in Appendix A and the product step is standard Barnes-G asymptotics, so these citations are not load-bearing circularity. One flagged passage is Appendix A, after Eq. (A.4): 'for γ ∈ (−2, 0) there is singular behavior near z = x, but one can still verify the statement. We omit the details.' This is a genuine completeness gap in a load-bearing lemma and a correctness risk for part of Re γ > −2, but it is not circularity: the claimed reduction is Green's theorem plus algebraic manipulation, and it does not presuppose Theorem 1.1 or Theorem 1.4. Accordingly, no significant circularity is found.
Assumptions & free parameters
assumptions (4)
- domain assumption Contour orthogonality reduction (Lemma 2.1) is valid for all Re gamma > -2, including gamma in (-2,0), where the Green's theorem step near z=x is asserted without proof.
- standard math The Deift-Zhou steepest descent method applies to the RHP with the g-function and level curves Gamma_r of Lemma 3.1.
- standard math The RHP (2.10)-(2.12) has a unique solution with the prescribed behavior at z=1.
- standard math Toeplitz determinant asymptotics of Claeys-Its-Krasovsky [14] apply to the symbol (5.3) with a=(gamma+alpha)/2, b=alpha/2, for the stated parameter ranges.
Cite this review
Pith. "Pith review of Asymptotics for a class of planar orthogonal polynomials and truncated unitary matrices." pith.science (2026). https://pith.science/paper/SD2O7EZL
@misc{pith2026250512633,
author = {Pith},
title = {Pith review of: Asymptotics for a class of planar orthogonal polynomials and truncated unitary matrices},
year = {2026},
howpublished = {\url{https://pith.science/paper/SD2O7EZL}},
note = {Machine review of arXiv:2505.12633}
}
abstract
We carry out the asymptotic analysis as $n \to \infty$ of a class of orthogonal polynomials $p_{n}(z)$ of degree $n$, defined with respect to the planar measure \begin{equation*} d\mu(z) = (1-|z|^{2})^{\alpha-1}|z-x|^{\gamma}\mathbf{1}_{|z| < 1}d^{2}z, \end{equation*} where $d^{2}z$ is the two dimensional area measure, $\alpha$ is a parameter that can grow with $n$, while $\gamma>-2$ and $x>0$ are fixed. This measure arises naturally in the study of characteristic polynomials of non-Hermitian ensembles and generalises the example of a Gaussian weight that was recently studied by several authors. We obtain asymptotics in all regions of the complex plane and via an appropriate differential identity, we obtain the asymptotic expansion of the partition function. The main approach is to convert the planar orthogonality to one defined on suitable contours in the complex plane. Then the asymptotic analysis is performed using the Deift-Zhou steepest descent method for the associated Riemann-Hilbert problem.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
- [22]
-
[26]
C. Webb and M. D. Wong. On the moments of the characteristic polynomial of a Ginibre random matrix. Proc. Lond. Math. Soc., 118(5):1017–1056, 2019
work page 2019
- [14]
-
[15]
A. Dea˜ no and N. Simm. Characteristic Polynomials of Complex Random Matrices and Painlev´ e Transcendents.Int. Math. Res. Not. , 2022(1):210–264, 05 2020
work page 2022
-
[1]
Y. Ameur, N.-G. Kang, and S.-M. Seo. The random normal matrix model: insertion of a point charge. Potential Anal., 58:331–372, 2021
work page 2021
- [2]
- [3]
-
[4]
S. Berezin, A. B. J. Kuijlaars, and I. Parra. Planar orthogonal polynomials as type I multiple orthogonal polynomials. SIGMA, 19:020, 2023
work page 2023
Show all 28 references
-
[5]
Bertola, J
M. Bertola, J. G. Elias Rebelo, and T. Grava. Painlev´ e IV critical asymptotics for orthogonal polynomials in the complex plane. SIGMA, 091:34 pages, 2018
2018
-
[6]
Byun, S.-Y
S. Byun, S.-Y. Lee, and M. Yang. Lemniscate ensembles with spectral singularity,
-
[7]
S.-S. Byun. Planar equilibrium measure problem in the quadratic fields with a point charge. Comput. Methods Funct. Theory , 24(2):303–332, 2024
2024
-
[8]
S.-S. Byun, P. J. Forrester, A. B. J. Kuijlaars, and S. Lahiry. Orthogonal polynomials in the spherical ensemble with two insertions, 2025. arXiv:2503.15732 [math.CA]
2025 arXiv
-
[9]
S.-S. Byun, P. J. Forrester, and S. Lahiry. Properties of the one-component coulomb gas on a sphere with two macroscopic external charges, 2025. arXiv:2501.05061 [math- ph]
2025 arXiv
-
[10]
Byun, N.-G
S.-S. Byun, N.-G. Kang, S.-M. Seo, and M. Yang. Free energy of spherical Coulomb gases with point charges, 2025. arXiv:2501.07284 [math.PR]
2025 arXiv
-
[11]
Byun, S.-M
S.-S. Byun, S.-M. Seo, and M. Yang. Free energy expansions of a conditional GinUE and large deviations of the smallest eigenvalue of the LUE, 2024. arXiv:2402.18983 [math.PR]
2024 arXiv
-
[12]
Byun and M
S.-S. Byun and M. Yang. Determinantal Coulomb gas ensembles with a class of dis- crete rotational symmetric potentials. SIAM J. Math. Anal. , 55(6):6867–6897, 2023
2023
-
[13]
Byun and E
S.-S. Byun and E. Yoo. Three topological phases of the elliptic Ginibre ensembles with a point charge, 2025. arXiv:2502.02948 [math-ph]
2025 arXiv
-
[16]
Deift, A
P. Deift, A. Its, and I. Krasovsky. Asymptotics of Toeplitz, Hankel, and Toeplitz+Hankel determinants with Fisher-Hartwig singularities. Ann. of Math. (2) , 174(2):1243–1299, 2011
2011
-
[17]
http://dlmf.nist.gov/, Release 1.1.6 of 2022-06-30
NIST Digital Library of Mathematical Functions. http://dlmf.nist.gov/, Release 1.1.6 of 2022-06-30. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds
2022
-
[18]
Gradshteyn and I.M
I.S. Gradshteyn and I.M. Ryzhik. Table of Integrals, Series, and Products . Else- vier/Academic Press, Amsterdam, 7th edition, 2007
2007
-
[19]
Hedenmalm
H. Hedenmalm. Soft Riemann-Hilbert problems and planar orthogonal polynomials. Comm. Pure Appl. Math. , 77(4):2413–2451, 2024
2024
-
[20]
Hedenmalm and A
H. Hedenmalm and A. Wennman. Planar orthogonal polynomials and boundary uni- versality in the random normal matrix model. Acta Math., 227(2):309–406, 2021
2021
-
[21]
J. P. Keating and N. C. Snaith. Random matrix theory and ζ(1/2+it). Comm. Math. Phys., 214(1):57–89, 2000
2000
-
[23]
Lee and M
S.-Y. Lee and M. Yang. Strong asymptotics of planar orthogonal polynomials: Gauss- ian weight perturbed by finite number of point charges. Commun. Pure Appl. Math. , 76(10):2888–2956, 2023
2023
-
[24]
L. D. Molag. Edge behavior of higher complex-dimensional determinantal point pro- cesses. Ann. Henri Poincar´ e, 24(12):4405–4437, 2023
2023
-
[25]
Serebryakov, N
A. Serebryakov, N. Simm, and G. Dubach. Characteristic polynomials of random truncations: Moments, duality and asymptotics. Random Matrices Theory Appl. , 12(01):2250049, 2023
2023
-
[27]
˙Zyczkowski and H-J
K. ˙Zyczkowski and H-J. Sommers. Truncations of random unitary matrices. J. Phys. A, 33(10):2045–2057, 2000. alfredo.deanho@uc3m.es. Departamento de Matem´aticas, Universidad Car- los III de Madrid, Spain. kmclaughlin@tulane.edu. School of Science & Engineering, Tulane Uni- ve...
2000
-
[2021]
arXiv:2107.07221 [math.PR]
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.