pith. sign in

arxiv: 2605.18002 · v1 · pith:4XRAONZInew · submitted 2026-05-18 · ❄️ cond-mat.str-el

Exact Organization of Density Matrices and Entanglement Structure in the Kitaev Spin Liquid

Chen-Chih Wang , Yi-Ping Huang , Sungkit Yip This is my paper

Pith reviewed 2026-05-20 00:58 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords Kitaev spin liquidentanglement structuredensity matrixgauge symmetryentanglement spectrumZ2 gauge theoryhoneycomb latticestring operators
0
0 comments X

The pith

The density matrix of the Kitaev spin liquid organizes into blocks according to the symmetries of its emergent gauge theory, producing extensive degeneracy in the entanglement spectrum.

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

The paper derives an exact expression for the density matrix of the Kitaev honeycomb model using only the original spin operators. This matrix is structured by equivalence classes of string operators that reflect the underlying gauge structure. The authors combine this form with the model's exact Gauss law and one-form Wilson symmetry to prove that the reduced density matrix for a subsystem is block diagonal in symmetry sectors. This block structure accounts for the large number of degenerate levels observed in the entanglement spectrum and allows the total entanglement entropy to be split into separate gauge and matter contributions. The same framework also covers subsystems with an odd number of sites and connects the entanglement spectrum to the parity of the fermions.

Core claim

In the Kitaev spin liquid, the density matrix takes an explicit form in terms of spin operators and is organized by equivalence classes of string operators tied to the gauge structure. Together with the exact Gauss law of the emergent Z2 gauge theory and the exact 1-form Wilson symmetry, this organization establishes a symmetry-resolved block-diagonal structure in the reduced density matrix. The block-diagonal form produces extensive degeneracy throughout the entanglement spectrum and is responsible for the separability of the entanglement entropy into gauge and matter sectors. The formalism further extends to odd-sized subsystems, where it relates the entanglement spectrum to fermion parity

What carries the argument

The symmetry-resolved block-diagonal structure of the reduced density matrix, generated by equivalence classes of string operators and protected by the exact Gauss law and 1-form Wilson symmetry of the emergent gauge theory.

Load-bearing premise

The exact Gauss law and 1-form Wilson symmetry of the emergent gauge theory continue to hold without mixing when the density matrix is written solely in the original spin operators.

What would settle it

A direct computation of the reduced density matrix for a small Kitaev cluster that finds matrix elements connecting different symmetry sectors, or an entanglement spectrum lacking the predicted extensive degeneracy.

Figures

Figures reproduced from arXiv: 2605.18002 by Chen-Chih Wang, Sungkit Yip, Yi-Ping Huang.

Figure 1
Figure 1. Figure 1: FIG. 1. The left panel shows two equivalent edge configura [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Two [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. An illustration of the block-diagonal structure of the [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. For cases with odd [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. A ’t Hooft loop operator [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. The left panel illustrates the plaqeutte configurations [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. A triangular lattice and a honeycomb lattice are [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
read the original abstract

We give an exact form of the density matrix of the spin-1/2 Kitaev spin liquid represented in terms of spin operators and study the entanglement structures of the Kitaev honeycomb model within the spin framework. We show that the density matrix is naturally organized by equivalence classes of string operators associated with the underlying gauge structure of the model. With the explicit form of the density matrix, plus the exact Gauss law of the emergent gauge theory and the exact 1-form Wilson symmetry in the Kitaev model, we demonstrate the existence of the underlying symmetry-resolved block-diagonal structure of the reduced density matrix, which gives rise to the extensive degeneracy in the entanglement spectrum. The block-diagonal structure is then proven to be responsible for the separability of the entanglement entropy into the gauge and matter parts. Furthermore, we extend the formalism to subsystems with an odd number of lattice sites, revealing a relation between the entanglement spectrum and the fermion parity that is seldom mentioned in the literature.

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.

Referee Report

1 major / 1 minor

Summary. The paper provides an explicit expression for the density matrix of the Kitaev honeycomb model written in the original spin operators. It invokes the exact Gauss law and 1-form Wilson symmetry of the emergent Z2 gauge theory to establish a symmetry-resolved block-diagonal structure in the reduced density matrix on a spatial subsystem. This block structure is shown to produce extensive degeneracy in the entanglement spectrum and to permit exact separation of the entanglement entropy into gauge and matter contributions. The formalism is extended to subsystems containing an odd number of sites, where a relation between the entanglement spectrum and fermion parity is identified.

Significance. If the central derivation holds, the work supplies a direct, symmetry-based explanation for the entanglement properties of the Kitaev spin liquid entirely within the spin-operator language. The explicit density-matrix construction and the resulting separability of the entanglement entropy constitute a useful technical advance for studies of topological order and entanglement in exactly solvable models.

major comments (1)
  1. [section establishing the symmetry-resolved block-diagonal structure of the reduced density matrix] The demonstration that the Gauss law and 1-form Wilson symmetry continue to label exact blocks of the reduced density matrix after the partial trace (the step that produces the claimed block-diagonal structure, extensive degeneracy, and gauge/matter separability) must explicitly rule out mixing generated by boundary string operators. Please add the calculation or argument showing that any such boundary contributions cancel or vanish identically.
minor comments (1)
  1. [presentation of the explicit density matrix] Clarify the precise definition of the string operators used to organize the density matrix and ensure their notation is uniform across all sections.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive assessment of its significance. We address the major comment below and will revise the manuscript to incorporate the requested clarification.

read point-by-point responses

  1. Referee: [section establishing the symmetry-resolved block-diagonal structure of the reduced density matrix] The demonstration that the Gauss law and 1-form Wilson symmetry continue to label exact blocks of the reduced density matrix after the partial trace (the step that produces the claimed block-diagonal structure, extensive degeneracy, and gauge/matter separability) must explicitly rule out mixing generated by boundary string operators. Please add the calculation or argument showing that any such boundary contributions cancel or vanish identically.

    Authors: We thank the referee for this constructive suggestion. While the manuscript establishes that the Gauss law and 1-form Wilson symmetry label the blocks of the reduced density matrix, we agree that an explicit demonstration ruling out mixing from boundary string operators strengthens the argument. In the revised manuscript we will add a dedicated paragraph (or short subsection) immediately following the definition of the reduced density matrix. There we show that any string operator crossing the entanglement cut can be continuously deformed, using the exact 1-form Wilson symmetry, into a combination of a closed loop entirely in the complement and a path segment lying wholly inside the subsystem. The Gauss-law constraint then forces the flux through the deformed loop to match the sector label, so that the operator factors into a product of an intra-subsystem operator and a traced-out operator whose expectation value is identical within each equivalence class. Consequently, off-diagonal matrix elements between distinct symmetry sectors vanish identically after the partial trace. This calculation confirms that the block-diagonal structure survives without additional mixing terms. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation uses explicit density matrix plus standard symmetries to exhibit block structure

full rationale

The paper begins with an explicit construction of the density matrix in original spin operators, then invokes the known exact Gauss law and 1-form Wilson symmetry (standard in Kitaev literature and not redefined here) to demonstrate the symmetry-resolved block-diagonal form of the reduced density matrix. This step is presented as a demonstration rather than a redefinition or tautological relabeling; the symmetries label sectors independently of the partial-trace reduction once the explicit form is given. No equations reduce by construction to fitted inputs, self-citations, or ansatzes smuggled from prior work by the same authors. The central claims about degeneracy and gauge/matter separability follow from this organization without the result being presupposed in the inputs. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the exact solvability of the Kitaev model and the persistence of its emergent Z2 gauge symmetries when the density matrix is written in spin operators. No free parameters or new invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption The Kitaev honeycomb model possesses an exact Gauss law and an exact 1-form Wilson symmetry associated with its emergent Z2 gauge structure.
    Invoked in the abstract to establish the symmetry-resolved block-diagonal structure of the reduced density matrix.

pith-pipeline@v0.9.0 · 5700 in / 1483 out tokens · 34129 ms · 2026-05-20T00:58:29.751635+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

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

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

56 extracted references · 56 canonical work pages · 1 internal anchor

  1. [1]

    We add this here to emphasize that we restrict ourselves to physical space

    In general, the above relation implies the direct corre- spondence between Σ D andM D PG   Y ⟨ij⟩∈D ˆuij   MD = ΣDPG.(C4) The Gutzwiller projection of the last term is actually redundant since Σ D is already physical. We add this here to emphasize that we restrict ourselves to physical space. With this equation, the gauge part, composed of gauge-field...

    work page

  2. [2]

    Kitaev,Anyons in an exactly solved model and beyond, Annals of Physics321, 2 (2006), january Special Issue

    A. Kitaev,Anyons in an exactly solved model and beyond, Annals of Physics321, 2 (2006), january Special Issue

    work page 2006

  3. [3]

    Savary and L

    L. Savary and L. Balents,Quantum spin liquids: a review, Reports on Progress in Physics80, 016502 (2016)

    work page 2016

  4. [4]

    Y. Zhou, K. Kanoda, and T.-K. Ng,Quantum spin liquid states, Rev. Mod. Phys.89, 025003 (2017)

    work page 2017

  5. [5]

    Matsuda, T

    Y. Matsuda, T. Shibauchi, and H.-Y. Kee,Kitaev quan- tum spin liquids, Rev. Mod. Phys.97, 045003 (2025)

    work page 2025

  6. [6]

    Kasahara, T

    Y. Kasahara, T. Ohnishi, Y. Mizukami, O. Tanaka, S. Ma, K. Sugii, N. Kurita, H. Tanaka, J. Nasu, Y. Mo- tome, T. Shibauchi, and Y. Matsuda,Majorana quanti- zation and half-integer thermal quantum Hall effect in a Kitaev spin liquid, Nature559, 227 (2018)

    work page 2018

  7. [7]

    L. J. Sandilands, Y. Tian, K. W. Plumb, Y.-J. Kim, and K. S. Burch,Scattering Continuum and Possible Frac- tionalized Excitations inα−RuCl 3, Phys. Rev. Lett.114, 147201 (2015)

    work page 2015

  8. [8]

    J. Nasu, J. Knolle, D. L. Kovrizhin, Y. Motome, and R. Moessner,Fermionic response from fractionalization in an insulating two-dimensional magnet, Nature Physics 12, 912 (2016)

    work page 2016

  9. [9]

    B. Zeng, X. Chen, D.-L. Zhou, and X.-G. Wen,Quan- tum Information Meets Quantum Matter: From Quan- tum Entanglement to Topological Phases of Many-Body Systems, Quantum Science and Technology (Springer, New York, NY, 2019)

    work page 2019

  10. [10]

    Wen,Colloquium: Zoo of quantum-topological phases of matter, Rev

    X.-G. Wen,Colloquium: Zoo of quantum-topological phases of matter, Rev. Mod. Phys.89, 041004 (2017)

    work page 2017

  11. [11]

    Laflorencie,Quantum entanglement in condensed mat- ter systems, Physics Reports646, 1 (2016), quantum en- tanglement in condensed matter systems

    N. Laflorencie,Quantum entanglement in condensed mat- ter systems, Physics Reports646, 1 (2016), quantum en- tanglement in condensed matter systems

    work page 2016

  12. [12]

    Amico, R

    L. Amico, R. Fazio, A. Osterloh, and V. Vedral,Entan- glement in many-body systems, Rev. Mod. Phys.80, 517 (2008)

    work page 2008

  13. [13]

    Eisert, M

    J. Eisert, M. Cramer, and M. B. Plenio,Colloquium: Area laws for the entanglement entropy, Rev. Mod. Phys. 82, 277 (2010)

    work page 2010

  14. [14]

    Li and F

    H. Li and F. D. M. Haldane,Entanglement Spectrum as a Generalization of Entanglement Entropy: Identification of Topological Order in Non-Abelian Fractional Quantum Hall Effect States, Phys. Rev. Lett.101, 010504 (2008)

    work page 2008

  15. [15]

    Thomale, A

    R. Thomale, A. Sterdyniak, N. Regnault, and B. A. Bernevig,Entanglement Gap and a New Principle of Adi- abatic Continuity, Phys. Rev. Lett.104, 180502 (2010)

    work page 2010

  16. [16]

    Chandran, M

    A. Chandran, M. Hermanns, N. Regnault, and B. A. Bernevig,Bulk-edge correspondence in entanglement spectra, Phys. Rev. B84, 205136 (2011)

    work page 2011

  17. [17]

    F. J. Burnell and C. Nayak,SU(2) slave fermion solution of the Kitaev honeycomb lattice model, Phys. Rev. B84, 125125 (2011)

    work page 2011

  18. [18]

    Feng, G.-M

    X.-Y. Feng, G.-M. Zhang, and T. Xiang,Topological Characterization of Quantum Phase Transitions in a Spin-1/2Model, Phys. Rev. Lett.98, 087204 (2007)

    work page 2007

  19. [19]

    Chen and J

    H.-D. Chen and J. Hu,Exact mapping between classical and topological orders in two-dimensional spin systems, Phys. Rev. B76, 193101 (2007)

    work page 2007

  20. [20]

    Chen and Z

    H.-D. Chen and Z. Nussinov,Exact results of the Ki- taev model on a hexagonal lattice: spin states, string and brane correlators, and anyonic excitations, Journal of Physics A: Mathematical and Theoretical41, 075001 (2008)

    work page 2008

  21. [21]

    J. Fu, J. Knolle, and N. B. Perkins,Three types of rep- resentation of spin in terms of Majorana fermions and an alternative solution of the Kitaev honeycomb model, Phys. Rev. B97, 115142 (2018)

    work page 2018

  22. [22]

    Yao and X.-L

    H. Yao and X.-L. Qi,Entanglement Entropy and Entan- glement Spectrum of the Kitaev Model, Phys. Rev. Lett. 105, 080501 (2010)

    work page 2010

  23. [23]

    Jin and V

    B.-Q. Jin and V. E. Korepin,Quantum Spin Chain, Toeplitz Determinants and the Fisher—Hartwig Conjec- ture, J. Stat. Phys.116, 79 (2004)

    work page 2004

  24. [24]

    Vidal, J

    G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev,Entangle- 19 ment in Quantum Critical Phenomena, Phys. Rev. Lett. 90, 227902 (2003)

    work page 2003

  25. [25]

    J. I. Latorre, E. Rico, and G. Vidal, Ground state en- tanglement in quantum spin chains (2004), arXiv:quant- ph/0304098 [quant-ph]

    work page arXiv 2004

  26. [26]

    Peschel,Calculation of reduced density matrices from correlation functions, Journal of Physics A: Mathemati- cal and General36, L205 (2003)

    I. Peschel,Calculation of reduced density matrices from correlation functions, Journal of Physics A: Mathemati- cal and General36, L205 (2003)

    work page 2003

  27. [27]

    Peschel and V

    I. Peschel and V. Eisler,Reduced density matrices and entanglement entropy in free lattice models, Journal of Physics A: Mathematical and Theoretical42, 504003 (2009)

    work page 2009

  28. [28]

    Chung and I

    M.-C. Chung and I. Peschel,Density-matrix spectra of solvable fermionic systems, Phys. Rev. B64, 064412 (2001)

    work page 2001

  29. [29]

    Cheong and C

    S.-A. Cheong and C. L. Henley,Many-body density ma- trices for free fermions, Phys. Rev. B69, 075111 (2004)

    work page 2004

  30. [30]

    Mandal, M

    S. Mandal, M. Maiti, and V. K. Varma,Entanglement and Majorana edge states in the Kitaev model, Phys. Rev. B94, 045421 (2016)

    work page 2016

  31. [31]

    Meichanetzidis, M

    K. Meichanetzidis, M. Cirio, J. K. Pachos, and V. Lahti- nen,Anatomy of fermionic entanglement and criticality in Kitaev spin liquids, Phys. Rev. B94, 115158 (2016)

    work page 2016

  32. [32]

    W.-T. Xu, T. Rakovszky, M. Knap, and F. Pollmann, Entanglement Properties of Gauge Theories from Higher- Form Symmetries, Phys. Rev. X15, 011001 (2025)

    work page 2025

  33. [33]

    F. L. Pedrocchi, S. Chesi, and D. Loss,Physical solutions of the Kitaev honeycomb model, Phys. Rev. B84, 165414 (2011)

    work page 2011

  34. [34]

    Zschocke and M

    F. Zschocke and M. Vojta,Physical states and finite- size effects in Kitaev’s honeycomb model: Bond disorder, spin excitations, and NMR line shape, Phys. Rev. B92, 014403 (2015)

    work page 2015

  35. [35]

    E. H. Lieb,Flux Phase of the Half-Filled Band, Phys. Rev. Lett.73, 2158 (1994)

    work page 1994

  36. [36]

    Fluxes, Laplacians and Kasteleyn's Theorem

    E. Lieb and M. Loss, Fluxes, Laplacians and Kasteleyn’s Theorem (1992), arXiv:cond-mat/9209031 [cond-mat]

    work page internal anchor Pith review Pith/arXiv arXiv 1992

  37. [37]

    E. H. Lieb,The flux-phase problem on planar lattices, Helvetica Physica Acta65, 247 (1992)

    work page 1992

  38. [38]

    In our discussion, the gauge is chosen asA 0 = 0, which leads to the relation [A µ(x), Eν(x′)] = −iδµν δ(x−x ′) mentioned in the main text

    The quantization relation of the gauge theory depends on the gauge choice. In our discussion, the gauge is chosen asA 0 = 0, which leads to the relation [A µ(x), Eν(x′)] = −iδµν δ(x−x ′) mentioned in the main text

    work page

  39. [39]

    After that, the exponent becomes the familiar form∇ ·E=ρ

    In this Gauss law formula, we can see thec-Majorana fermion operator is exactly the matter fieldc i = e iπρi if theb-Majorana operators are represented asb αij i = eiπEij . After that, the exponent becomes the familiar form∇ ·E=ρ

    work page

  40. [40]

    Baskaran, S

    G. Baskaran, S. Mandal, and R. Shankar,Exact Results for Spin Dynamics and Fractionalization in the Kitaev Model, Phys. Rev. Lett.98, 247201 (2007)

    work page 2007

  41. [41]

    These four cases are equivalent for the analysis of the entanglement feature, but only in different topological sectors

    For a torus, there are four choices in total for non- contractible Wilson loops corresponding to⟨ ˜Wv⟩=±1 with⟨ ˜Wh⟩=±1. These four cases are equivalent for the analysis of the entanglement feature, but only in different topological sectors. In this work, we only consider one of the topological sectors

    work page

  42. [42]

    Levin and X.-G

    M. Levin and X.-G. Wen,Fermions, strings, and gauge fields in lattice spin models, Phys. Rev. B67, 245316 (2003)

    work page 2003

  43. [43]

    Note that the projectorP (p) becomes P (p) → Y p 1+ Q ⟨ij⟩∈p ˆuij 2 ! in the Majorana representation, where “⟨ij⟩ ∈p” stands forthose edges of the plaquettep

    work page

  44. [44]

    M. A. Levin and X.-G. Wen,String-net condensation: A physical mechanism for topological phases, Phys. Rev. B 71, 045110 (2005)

    work page 2005

  45. [45]

    Levin and X.-G

    M. Levin and X.-G. Wen,Colloquium: Photons and elec- trons as emergent phenomena, Rev. Mod. Phys.77, 871 (2005)

    work page 2005

  46. [46]

    Wen,Quantum order from string-net condensa- tions and the origin of light and massless fermions, Phys

    X.-G. Wen,Quantum order from string-net condensa- tions and the origin of light and massless fermions, Phys. Rev. D68, 065003 (2003)

    work page 2003

  47. [47]

    Lee and F

    J. Lee and F. Wilczek,Algebra of Majorana Doubling, Phys. Rev. Lett.111, 226402 (2013)

    work page 2013

  48. [48]

    One can choose an arbitrary dimer covering with a remaining site lat the boundary that covers the subregion

    The dimer covering is not uniquely defined. One can choose an arbitrary dimer covering with a remaining site lat the boundary that covers the subregion

    work page

  49. [49]

    H.-Y. Lee, R. Kaneko, T. Okubo, and N. Kawashima, Gapless Kitaev Spin Liquid to Classical String Gas through Tensor Networks, Phys. Rev. Lett.123, 087203 (2019)

    work page 2019

  50. [50]

    Lahtinen and J

    V. Lahtinen and J. K. Pachos,Non-Abelian statistics as a Berry phase in exactly solvable models, New Journal of Physics11, 093027 (2009)

    work page 2009

  51. [51]

    Lahtinen and J

    V. Lahtinen and J. K. Pachos,Topological phase transi- tions driven by gauge fields in an exactly solvable model, Phys. Rev. B81, 245132 (2010)

    work page 2010

  52. [52]

    Lahtinen, G

    V. Lahtinen, G. Kells, A. Carollo, T. Stitt, J. Vala, and J. K. Pachos,Spectrum of the non-abelian phase in Ki- taev’s honeycomb lattice model, Annals of Physics323, 2286 (2008)

    work page 2008

  53. [53]

    Lee, G.-M

    D.-H. Lee, G.-M. Zhang, and T. Xiang,Edge Solitons of Topological Insulators and Fractionalized Quasiparticles in Two Dimensions, Phys. Rev. Lett.99, 196805 (2007)

    work page 2007

  54. [54]

    Yu,Gauge symmetry in Kitaev-type spin models and index theorems on odd manifolds, Nuclear Physics B799, 345 (2008)

    Y. Yu,Gauge symmetry in Kitaev-type spin models and index theorems on odd manifolds, Nuclear Physics B799, 345 (2008)

    work page 2008

  55. [55]

    Yao and S

    H. Yao and S. A. Kivelson,Exact Chiral Spin Liquid with Non-Abelian Anyons, Phys. Rev. Lett.99, 247203 (2007)

    work page 2007

  56. [56]

    In this way, the fermion parity operator represented using the Majorana fermions becomes ˆπ≡(−1) f †f =ic 1c2

    Two distinct Majorana fermions, sayc 1 andc 2, can com- bine to form a normal fermion degree of freedomfin a way thatf= 1 2 (c1 +ic 2) andf † = 1 2 (c1 −ic 2). In this way, the fermion parity operator represented using the Majorana fermions becomes ˆπ≡(−1) f †f =ic 1c2. In a more general sense, the fermion parity operator is equal to the multiplication of...

    work page