Pith. sign in

REVIEW 3 major objections 6 minor 1 cited by

Exact critical exponents of the Motzkin and Fredkin Chains

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read The paper derives exact critical exponents for the q-deformed Motzkin and Fredkin spin chains: η=3/2 and ν=2/3.

desk verdict A credible TM-based derivation of η=3/2 with a plausible but not fully derived ν=2/3; the title's 'exact' is too strong but the paper deserves a serious referee. read the letter →

arxiv 2507.14656 v4 pith:FSK5YYO5 submitted 2025-07-19 cond-mat.str-el cond-mat.stat-mechhep-thmath-phmath.MPquant-ph

Olai B. Mykland , Zhao Zhang This is my paper · ORCID
classification cond-mat.str-elcond-mat.stat-mechhep-thmath-phmath.MPquant-ph
keywords MotzkinchainFredkinq-deformationcriticalexponentstransfermatrixproductstaterenormalizationgroupquantumphasetransition
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

The paper analyzes the q-deformed Motzkin and Fredkin spin chains, two frustration-free models whose ground states are superpositions of weighted random walks, and derives that their phase transition at $q=1$ is continuous with exactly determined critical exponents. Spin correlations at criticality decay with exponent $\eta=3/2$, equivalently the one-point function satisfies $\langle S^z_r\rangle\sim 1/\sqrt{r}$ for $1\ll r\ll L$. The correlation length on both sides of the transition diverges as $|\log q|^{-2/3}$, giving $\nu=2/3$, and the same exponent applies for $q<1$ and $q>1$. These values are obtained analytically from a transfer matrix built from the exact matrix product state, together with a renormalization-group argument, and are confirmed by MPS simulations and direct transfer-matrix diagonalization.

What carries the argument

The central object is the transfer matrix $T_{\max}$ constructed by contracting the physical legs of the exact matrix product state; it is a tridiagonal matrix acting in the height space of the random-walk configurations. At $q=1$ it is Toeplitz with sine-wave eigenvectors, which turns the correlation function into a saddle-point integral, while for $q\neq 1$ its eigenvalue gap determines the correlation length. A second ingredient is the dual zero-dimensional Hamiltonian $H=-\log T$, whose dispersion near the bottom of the spectrum directly gives the exponent $\eta$. The correlation-length exponent comes from an RG analysis of the effective free energy $F_L(\sigma,\tau)=\sum_r[\sigma(S^z_r)^2+\tau(L-r)S^z_r]$, together with the scaling dimension $S'^z=b^{1/2}S^z$ and the duality relating $T(q)$ and $T(1/q)$.

What would settle it

Measure $\xi(\tau)$ from the transfer-matrix gap at $q=e^{\pm\tau}$ over $\tau$ from $10^{-7}$ to $10^{-2}$; if the log-log slope of $\xi$ versus $|\tau|$ is not $-2/3$, the assumed RG closure is wrong. Separately, compute $\langle S^z_r\rangle$ at $q=1$ in the exact MPS for $L\gg r$ and fit $\log\langle S^z_r\rangle$ versus $\log r$; a slope different from $-1/2$ would falsify $\eta=3/2$.

Watch

Extended reading notes

Core claim

At the critical deformation $q=1$, the maximum transfer-matrix block for both chains becomes a tridiagonal Toeplitz matrix with known eigenvalues $2\cos(k\pi/L)$ (Fredkin) and $1+2\cos(k\pi/L)$ (Motzkin) and sine eigenvectors. Using the spectral decomposition of this matrix, the spin one-point function is evaluated by saddle-point integrals and found to be $\langle S^z_r\rangle\approx 1/\sqrt{2\pi(r-1)}$ plus a subleading oscillating term for the Fredkin chain and $\langle S^z_r\rangle\approx 4/\sqrt{3\pi(r-1)}$ for the Motzkin chain, so the correlation exponent is $\eta=3/2$. Away from criticality the $q$-deformation opens a gap in the transfer-matrix spectrum; the paper derives the correlation-length exponent from the scaling dimension of the spin operator together with an RG step, yielding $\nu=2/3$, and shows $\nu$ is the same in the ordered and disordered phases through the duality $q^{2L}T(1/q)=ST(q)S^{-1}$. The ordered phase is interpreted as bulk order induced by correlations with boundary spins under an effective edge field, with a domain-wall profile $\langle S^z_r\rangle=S\tanh(\tau(r-L))$ whose width scales as $|\log q|^{-1}$.

Load-bearing premise

The calculation of the correlation-length exponent depends on an assumed coarse-graining rule for how the deformation parameter and the block spin rescale together, rather than on a derivation from the microscopic transfer matrix; if that rule were replaced by a different closure, the exponent would change.

Editorial extensions

If this is right

  • At $q=1$, both chains are gapless with power-law spin correlations, so the ground state has logarithmic entanglement entropy and no finite correlation length, while the new exponents fix the scaling of correlation functions.
  • For every $q\neq 1$, correlations decay exponentially on both sides of the transition; the disorder-to-order transition is therefore continuous, with correlation length $\xi\sim|\log q|^{-2/3}$.
  • In the ordered phase $q>1$, the magnetization profile is a hyperbolic-tangent domain wall of width $w\sim|\log q|^{-1}$, so $w\sim\xi^{3/2}$, a scaling relation that differs from the Ginzburg-Landau expectation $w\sim\xi$.
  • The same transfer-matrix method applies to other translationally invariant gapless matrix product states: the exponent $\eta$ is read off from the low-energy dispersion of the dual zero-dimensional Hamiltonian $H=-\log T$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the RG step for $\nu$ is a scaling closure rather than a microscopic derivation, the cleanest independent test is to compute the transfer-matrix correlation length over a wide range of $|\tau|$; the analytic status of the closure could be settled by deriving the RG rule directly from coarse-graining the pentagonal tensor network, which the paper leaves as a future direction.
  • The unusual $w\sim\xi^{3/2}$ relation suggests the ordered phase is not described by the standard Ginzburg-Landau domain-wall theory; a modified coarse-grained free energy might reproduce both $\xi$ and $w$ from a single fixed point.
  • The transfer-matrix/dual-Hamiltonian spectral method is generic enough to be tested on colored or multi-species Motzkin and Fredkin chains, where the dispersion of the dual Hamiltonian would predict $\eta$ without solving the full model.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. The paper studies the q-deformed Fredkin and Motzkin spin chains, whose exact ground states are weighted superpositions of Dyck or Motzkin walks. It constructs a translationally invariant MPS representation and the associated transfer matrix (TM), then uses the TM spectrum at q=1 to derive the one-point function <S^z_r> ~ 1/sqrt(r) for 1 << r << L, yielding the critical exponent eta = 3/2 for spin correlations. It further introduces a dual zero-dimensional Hamiltonian whose dispersion gives eta generically, and combines the spin scaling dimension with an RG analysis of an effective free energy to obtain the correlation-length exponent nu = 2/3 for both sides of the transition. The analytic results are compared with MPS simulations and direct TM diagonalization, and an additional prediction for the domain-wall thickness w ~ |tau|^{-1} is derived and numerically confirmed.

Significance. If the claims hold, this is a valuable analytical treatment of an exotic frustration-free quantum phase transition. The construction of the TM from the exact MPS and the saddle-point evaluation leading to Eq. (27) is transparent and is validated by MPS data over a wide range of r. The observed duality between the ordered and disordered phases and the derivation of the domain-wall scaling w ~ |tau|^{-1} are interesting and go beyond previous work. The proposed dual-Hamiltonian perspective on critical exponents is a useful methodological idea. However, the analytic derivation of nu = 2/3 is significantly less controlled than the derivation of eta: it rests on an imported free-energy form and an assumed RG closure, while the numerical confirmation is limited to a relatively narrow range of tau without error bars. The paper would be strengthened by clearly separating rigorously derived statements from scaling hypotheses.

major comments (3)
  1. [Section 6, Eq. (34)-(38)] The derivation of nu = 2/3 depends on the free energy Eq. (34), which is imported from Ref. [29], and the paper itself states immediately after Eq. (34) that this free energy is only a good approximation for q != 1. Since nu is defined by xi ~ |tau|^{-nu} in the limit tau -> 0, the approximation is least controlled precisely in the regime used to extract the exponent. Moreover, Eq. (34) is a sum of single-site terms with no spatial gradient, so the partition function factorizes and this free energy contains no intrinsic correlation length; the xi appearing in Eq. (37) comes from the TM spectrum, not from Eq. (34). The authors should either derive the RG eigenvalue y_tau from the microscopic TM or explicitly state that nu = 2/3 is a scaling hypothesis supported by numerics rather than an exact analytic result.
  2. [Section 6, Eq. (36)] The RG closure tau' S'^z = b^2 tau S^z and sigma'(S'^z)^2 = b sigma (S^z)^2 is assumed, not derived from the coarse-graining of the TM. Once S'^z = b^{1/2} S^z is imposed, this closure yields y_tau = 3/2, but a different transformation of tau would give a different nu; the paper provides no microscopic justification for why this particular closure is the correct one. The numerical TM diagonalization in Fig. 6(b) supports the final value 2/3, but it does not test the analytic derivation. This is a load-bearing gap for the claim that nu is obtained exactly.
  3. [Section 4, Eq. (23)-(27)] The conversion from the discrete sum in Eq. (23) to the integral in Eq. (25) sets the fast oscillating factor (1-(-1)^{k1+k2}) to its average value 1, and the saddle-point evaluation in Appendix C relies on 'cancellations' of divergences in the cotangent factors. These steps are plausible but not rigorously justified, so the derivation of eta = 3/2 from the TM is controlled rather than exact. Given the title and abstract claim exact critical exponents, the authors should either supply a tighter error estimate for these approximations or qualify the wording.
minor comments (6)
  1. [Section 2, Eq. (6)] The sentence 'The Hamiltonians is written as' should read 'The Hamiltonians are written as'; this grammatical error appears near Eq. (6).
  2. [Figure 3 caption] The caption says 'the two Fredkin moves F^1_i and F^1_i'; the second move should be F^2_i.
  3. [Section 6, Eq. (37)-(38)] The notation tau = -log q makes tau negative for q > 1; Eqs. (37)-(38) and the fits in Figs. 6(b), 7(b), and 8(b) should consistently use |tau| or state the sign convention explicitly.
  4. [Figures 6-8] The statement that 'identical results are obtained for the Motzkin chain' is made for three figures without showing the Motzkin data; the figures should include the Motzkin data or this claim should be removed.
  5. [Eq. (34)] In Eq. (34), sigma is called the surface tension and S^z_r is used as a classical height variable, but earlier S^z_r denotes the spin operator; this double use of notation should be clarified.
  6. [Section 5 and footnote 1] The mapping T = e^{-H} is formal when T has non-positive eigenvalues, as acknowledged in the footnote; the main text should state more prominently that the relation (28) is used only to extract the asymptotic decay behavior, not as an exact spectral equivalence.

Circularity Check

1 steps flagged · score 4.0 of 10

RG derivation of ν=2/3 in Sec. 6 is built on a self-cited single-site free energy and an asserted closure (Eqs. 34, 36); independent TM/MPS numerics keep the claim from being fully circular.

  1. self citation load bearing [Section 6, Eqs. (34)-(37)]
    "An expression for the free energy of the ensemble of configurations in the GS has been derived for the scaling limit in Ref.[29]. ... In order for the free energy to stay invariant, its density must satisfy f ′(σ′,τ′) = b f(σ,τ), which demands the RG equations τ′S′z = b2τSz, σ′(S′z)2 = bσ(Sz)2. As Eq. (27) established the scaling dimension of the spin operator S′z = b1/2Sz, we must have τ′ = b3/2τ..."

    Eq. (34) is imported from Ref. [29], whose author list overlaps with the present paper (Z. Zhang), and is explicitly called only an approximation away from q=1. The analytic derivation of ν=2/3 then depends entirely on the asserted coarse-graining closure Eq. (36) (τ′S′z = b2τSz, σ′(S′z)2 = bσ(Sz)2) combined with the spin rescaling S′z = b1/2Sz obtained from Eq. (27). This fixes yτ = 3/2 and hence ξ ∼ τ^{−2/3}. Thus the claimed exact value of ν is not derived from the microscopic transfer matrix; it is the output of a scaling ansatz whose form is supplied by the authors' own prior work. The circularity is only partial: the numerical TM diagonalization in Fig. 6(b) and the MPS correlation lengths in Fig.

full rationale

The paper's main derivation of η=3/2 (Secs. 3-5) is self-contained: the TM is built from the exact MPS amplitudes (Eq. 12), the q=1 eigensystem is explicit (Eq. 21), and the saddle-point evaluation (Appendix C) yields Eq. (27) without fitting. The identification of η from the one-point boundary-to-bulk correlation is a direct use of the definition Eq. (10) with d=1. The duality ν+ = ν− follows from the TM relation Eq. (33), also derived in the paper. The only load-bearing item imported from the authors' prior work is the free energy Eq. (34) from Ref. [29], used in Sec. 6 to set the RG eigenvalue of τ; the RG closure Eq. (36) is asserted rather than derived from the TM. This makes the analytic path to ν=2/3 an RG scaling argument rather than a closed transfer-matrix computation. However, the numerical TM diagonalization (Fig. 6b) and MPS data (Fig. 7b) extract ν ≈ 0.67 without using Eq. (34), providing independent content. Accordingly, the circularity score is moderate (4), not higher.

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

The central claim rests on the exact GS/MPS forms (domain assumptions from earlier literature) and on the free-energy/RG input for nu. No parameters are fitted in the analytic derivations; sigma is a physical parameter that cancels in the RG. No new physical entities are postulated; the dual 0D Hamiltonian is a mathematical construct.

assumptions (5)
  • domain assumption The q-deformed ground state is exactly the area-weighted superposition of Dyck or Motzkin walks, Eq. (4).
    Standard result from Refs. [21,24,25] used throughout; not derived in this paper.
  • domain assumption The uniform MPS representation Eq. (12) is exact for the ground states at all q.
    Constructed in Section 3 from the height representation; used to build the transfer matrix.
  • domain assumption The effective free energy Eq. (34) is the correct scaling-limit description for q != 1.
    Taken from Ref. [29] by the same authors; this is the load-bearing input for the nu RG analysis.
  • ad hoc to paper RG closure tau' S'^z = b^2 tau S^z and sigma'(S'^z)^2 = b sigma(S^z)^2, Eq. (36).
    Asserted in Section 6 to keep the free-energy density invariant; not derived from the microscopic TM. This is the key step producing nu = 2/3.
  • standard math Saddle-point and continuum approximations in Section C are valid for 1 << r << L.
    Standard asymptotic analysis, verified numerically against MPS in Fig. 5.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exact critical exponents of the Motzkin and Fredkin Chains." pith.science (2026). https://pith.science/paper/FSK5YYO5

@misc{pith2026250714656,
  author       = {Pith},
  title        = {Pith review of: Exact critical exponents of the Motzkin and Fredkin Chains},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FSK5YYO5}},
  note         = {Machine review of arXiv:2507.14656}
}
abstract

The Motzkin and Fredkin chains are frustration-free spin models with exactly solvable ground states whose $q$-deformations realize an exotic quantum phase transition from a disordered phase to an ordered one under domain-wall boundary conditions. Previous work has mainly focused on their entanglement scaling and spectral gaps, particularly in color-enriched variants. Here we systematically characterize the critical behavior of this transition with three main advances. First, we interpret bulk magnetic order as arising from correlations with boundary spins subject to an effective edge field, directly relating the ordered phase to the domain-wall boundary condition. Second, using the transfer matrix (TM) constructed from the exact matrix product state (MPS) representation, together with a continuum and RG analysis, we derive the algebraic decay of spin correlations and obtain the critical exponent $\eta = \frac{3}{2}$. We further generalize this TM/dual-Hamiltonian method to translationally invariant gapless states with generic power-law correlations, via a zero-dimensional Hamiltonian dual to the one-dimensional TM. Third, the TM spectrum reveals a duality between the ordered and disordered phases which, combined with scale invariance at criticality and the scaling dimension of the spin operator, yields $\nu_\pm = \frac{2}{3}$ from an RG analysis of the $q$-deformed ground states. Both exponents are confirmed numerically by MPS simulations and by direct diagonalization of the TM.

Figures

Figures reproduced from arXiv: 2507.14656 by the authors.

Figure 1
Figure 1. Spin configurations and corresponding Motzkin walks for [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Phase diagram of the deformed Fredkin and Motzkin spin chain, showing [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Upper panel: the two Fredkin moves F 1 i and F 1 i related to the states in Eq. (8). Moving from left to right reduces the area under the walk segment by 2 units. Lower panel: the three Motzkin moves M1 i , M2 i and M2 i related to the states in Eq. (9). Moving from left to right reduces the area under the walk segment by 1 unit. The solid gray lines indicate the height. The fact that H F/M tot give the GS in Eq. (4… view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: (a) Diagrammatic depiction of the Fredkin TM (left) and its nonzero el [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Comparison of 〈S z r 〉 obtained from the MPS (L = 500) and from the analytic result in Eq. (27) obtained from the TM for the Fredkin and Motzkin chains. Panels (a), (b) correspond to the Fredkin chain, and panels (c), (d) to the Motzkin chain. Note that panels (a) and …
Figure 6
Figure 6. Figure 6: (a) The eigenvalue gap ∆(q, L) = (λ1 − λ2 )/λ1 as a function of the system size L for various qs. (b) Log-log plot of the correlation length ξ(τ) = −1/ ln| λ2 λ1 |, where τ = log q, for various q < 1 (marked by blue circles) and q > 1 (marked by red triangles), with L …
Figure 7
Figure 7. Figure 7: (a) Plot of log 〈S z r 〉 in the GS for various q < 1. The faint gray lines indicate the linear fits determining the correlation length ξ(q). The plot is shown for the first 500 sites with L = 1000 and F (resp. M) denotes Fredkin (resp. Motzkin). (b) Log￾log plot of the…
Figure 8
Figure 8. Figure 8: (a) Plot of 〈S z r 〉 for various q > 1, for a system of L = 1000. The faint gray lines indicate the fitted Fermi–Dirac distribution (almost indistinguishable from the data) giving the domain wall width w(q). Note that F (resp. M) denotes Fredkin (resp. Motzkin). (b) Lo…
Figure 9
Figure 9. Figure 9: (a) Contracting three tensors of the original Motzkin TN gives a uniform [PITH_FULL_IMAGE:figures/full_fig_p019_9.png]
Figure 10
Figure 10. Figure 10: (a) Two different ways to decompose a rank-4 tensors. (b) Five different [PITH_FULL_IMAGE:figures/full_fig_p020_10.png]
Figure 11
Figure 11. Figure 11: (a) Decomposition of a rank-5 tensor into the product of three rank-3 [PITH_FULL_IMAGE:figures/full_fig_p021_11.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Highly Entangled Quantum Spin Chains on Fermat's Spiral

    quant-ph 2025-06 conditional novelty 6.0 of 10

    A spiral-embedded Motzkin chain on a square lattice realizes a 2D ground state with volume-law entanglement and simple two-body interactions.

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry add.period pagebackref " Appearances: " * 'skip if write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = output.state before.all = or 'skip after.sentence 'output.state := if FUNCTION new.block.check...

    work page

  2. [2]

    , " * write output.state after.block =

    ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution isbn journal key month note number organization pages publisher school series title type url volume year label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 'af...

  3. [3]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.