Majorana Quantum Lego, a Route Towards Fracton Matter

Fracton topological phases host fractionalized topological quasiparticles with restricted mobility, with promising applications to fault-tolerant quantum computation. While a variety of exactly solvable fracton models have been proposed, there is a need for platforms to realize them experimentally. We show that a rich set of fracton phases emerges in interacting Majorana band models whose building blocks are within experimental reach. Specifically, our Majorana constructions overcome a principal obstacle, namely the implementation of the complicated spin cluster interactions underlying fracton stabilizer codes. The basic building blocks of the proposed constructions include Coulomb blockaded Majorana islands and weak inter-island Majorana hybridizations. This setting produces a wide variety of fracton states and promises numerous opportunities for probing and controlling fracton phases experimentally. Our approach also reveals the relation between fracton phases and Majorana fermion codes and further generates a hierarchy of fracton spin liquids.

Searching for and exploring exotic phases of matter is a principal goal of condensed matter physics. In the presence of strong interactions, quantum many-body systems composed of a limited number of elementary particles assume a remarkable variety of exotic phases whose low-energy degrees of freedom are much richer than suggested by their constituents. Prominent examples of such emergent quantum phases are topological phases, whose quasiparticle excitations carry fractional quantum numbers and obey anyonic statistics [1,2]. The low-energy properties of these topologically ordered states are characterized by topological quantum field theories (TQFT) [3,4].
Recently distinct long-range entangled states, transcending the TQFT paradigm and termed fracton phases, have been discovered and intensively studied in exactly solvable lattice models [5][6][7][8][9]. Fracton topological order shares many features with topological order, including nontrivial braiding statistics and symmetry fractionalization. At the same time, fracton phases have a subextensive ground-state degeneracy depending on system size in addition to lattice topology, and quasiparticles with restricted mobility, moving within lowerdimensional manifolds such as planes, lines, or fractals. The subdimensional nature of fracton excitations gives rise to unconventional features including glassiness and subdiffusive dynamics [9,10]. The restricted quasiparticle mobility makes fracton stabilizer codes interesting for quantum-memory and quantum-computation applications [5,11].
While theoretical aspects of fracton phases have been intensively explored via quantum stabilizer codes as well as higher rank gauge theories [12][13][14][15], their physical realization remains a key challenge of condensed matter physics and quantum information science [6,16,17]. The principal obstacle is that models exhibiting fracton physics tend to involve rather complicated spin cluster interactions. This raises the question whether and how such exotic fracton states emerge in models with more physical ingredients and interactions. These might then be amenable to experimental implementation, and assuming tunable tunable interaction parameters, allow for controlling and manipulating fracton phases and excitations, e.g., for quantum computing purposes.
Here, we show that many known fracton stabilizer codes can be obtained from Majorana band models with strong onsite interactions. There is currently a major push to develop the required technology for realizing such models in the context of Majorana-based quantum computing [18], and our work shows how to leverage this budding technology for investigating fracton phases of matter. Our constructions include fracton phases of both flavors, referred to as type-I and type-II fracton codes. While excitations of the former are created by line-or membrane operators, excitations of type-II fractons are generated by operators with fractal support [19].
The principal ingredients of our constructions are Majorana hybridization as well as local interactions which fix local fermion parities [20][21][22][23][24][25]. These interactions can be implemented using Majorana islands, also referred to as Majorana Cooper pair boxes, which underlie current designs for Majorana-based topological qubits [18,20,26,27]. Each island contains some number of Majoranas, e.g., at the ends of semiconductor wires proximity coupled to a superconductor [18,28,29]. The island's charging energy fixes its fermion parity, corresponding to a multi-Majorana interaction. In order to realize fracton phases, the fermion parity typically has to be fixed for overlapping sets of Majoranas. We show how such interactions can be realized by judiciously hybridizing Majoranas from multiple Majorana islands.
Our construction opens new avenues in the study of fracton phases from both conceptual and applied perspectives. On the theoretical side, our Majorana-based setups for fracton models imply that fractonic phases of matter can emerge from strongly interacting one-dimensional p-wave superconductors. Moreover, they expose the close relation between 3D Majorana fermion codes [30] and fracton codes, and further generate a hierarchy of 3D SO(N) invariant spin liquids with fractonic behavior. Analogous to Ref. [31], this correspondence also establishes a foliation-based equivalence relation between bosonic and fermionic fracton states. On the experimental side, our Majorana-based models not only provide a possible route towards realizing fracton phases of matter, but might also give access to interesting observables of fracton phases and their phase transitions. In particular, the Hamilto- nian would be amenable to tuning it away from the stabilizer limit, opening a path towards exploring confinement and disorder effects on fracton matter.

The planon-lineon code
We begin with a topological superconductor on a bodycentered cubic lattice. Each site contains eight Majoranas γ 1 , . . . , γ 8 which are each hybridized with a Majorana on a nearest-neighbor site as shown in Fig. 1. Thus, the topological superconductor has the Hamiltonian (1) i γ 3 i and the product of two Majoranas associated with any vertical edge as the Pauli-X operator σ x i , or vice versa. We make the first (second) choice for the red (green) sublattice in Fig. 1.
In the strong-U limit, we can treat the Majorana hybridizations as a perturbation. The effective low-energy Hamiltonian emerges from local terms which leave the fermion parities of the Majorana quartets unchanged, even though the latter are flipped by individual hybridization terms. The leadingorder Hamiltonian involves 16-Majorana terms for the octahedra shown in Fig. 1 (see Methods section below). Writing the Hamiltonian in the spin representation yields Here octa a and octa b refer to the two types of octahedra in Fig. 1 with red (green) sites at top and bottom and four green (red) sites in between. Thus, our construction exactly reproduces the planon-lineon model in Ref. [31] whose elementary quasiparticles are lineons and planons with mobility restricted to the z-direction and the xz (yz)-planes, respectively. Experiment implementation.-Our construction establishes a route towards building fracton codes from interacting Majoranas. To implement this construction, we must establish that one can realize the onsite interaction in Eq. (2).
To this end, we distribute the eight Majoranas of each site over two adjacent superconducting islands (SCI) as shown in Fig. 1. Each SCI could be made from two semiconductor quantum wires proximity coupled to the same superconductor. The proximity-coupled quantum wires effectively realize open Kitaev chains with two Majorana zero modes localized at their ends, so that there are a total of four Majoranas on each SCI. By virtue of their charging energy, each SCI can be tuned to have even fermion parity, effectively implementing the interaction terms U (γ 1 [26,27]. To generate the remaining four-Majorana interaction in Eq. (2), we turn on inter-island Majorana hybridization with amplitude t, Note that at the same time, there is no hybridization between Majoranas γ 1 i and γ 2 i as well as γ 5 i and γ 6 i . These inter-island hybridizations can in principle be implemented by direct tunnel coupling. Alternatively, and perhaps more flexibly, one can bridge between the two Majorana islands using a coherent link [26]. Such a link consists of an additional proximitycoupled quantum wire with its fermion parity fixed by charging energy. Its two Majorana end states would then be hybridized with the two Majorana end states of the Majorana islands which one wants to hybridize. Since the hybridization between the Majoranas on the coherent link and the islands can be realized through gate controlled tunnel junctions, the hybridization strength is tunable. We finally mention that the same hybridization via coherent link can be used to implement the Majorana hybridizations between different sites as given by Eq. (1).
In the limit of a large charging energy, which fixes the fermion parities of the SCI, a single Majorana tunneling between the islands is suppressed and the lowest order processes involve pairs of Majorana tunneling terms. At leading order, the resulting Hamiltonian is, where c is a number of order one. This produces an anisotropic version of the onsite interaction in Eq. (2). The anisotropy of the interaction coefficients does not modify the ground state manifold and thus the low-energy spin-1/2 degree of freedom. The additional weak hybridization t of Majoranas on nearest-neighbor sites as described by Eq. (1) generates the spin interactions on the octahedra. As a sufficient condition, the local spin-1/2 degrees of freedom remain intact in the limit U > t > t such that t 2 /U > t , and the resulting effective Hamiltonian realizes the planon-lineon code at leading order.

The fractal Chamon code
Following the spirit of this construction, one can realize many fracton stabilizer codes based on crossing Kitaev chains and strong onsite interactions. Specifically, we show now how a type-II fracton code -the fractal Chamon code -emerges in this manner. Unlike type-I fracton codes whose excitations are created by straight-line or planar membrane operators, the excitations of type-II fracton codes involve operators with fractal support [5,32]. In the case of the fractal Chamon code [33], excitations live at the corners of 2D Sierpinski triangles of side length 2 n [19]. These corner excitations are immobile in the xy plane unless one enlarges the Sierpinski triangle to side length 2 n+1 which is inhibited by a large energy barrier.
The configuration of crossing Kitaev chains illustrated in Fig. 2 generates a hexagonal close-packed lattice, with each site containing eight Majoranas. Six of these labeled γ 1 , . . . , γ 6 are paired to Majoranas on nearest-neighbor sites in the same xy plane, while the remaining two (labeled γ A and γ B ) are paired with Majoranas on neighboring sites along the ±z direction. Now consider the onsite interaction In the strong-coupling limit, the eight Majoranas retain a spin-1/2 degree of freedom. As for the planon-lineon model, we can define spin operators through σ With these definitions, the low-energy Hamiltonian takes the form of the fractal Chamon code with five-qubit stabilizers defined on prisms as illustrated in Fig. 2 [33]. The elementary excitations include lineons with restricted mobility along the z-direction, generated by a line of σ z operators, in addition to the fractal excitations within the xy planes which live at the corners of 2D Sierpinski triangles.
The onsite interaction in Eq. (6) can be implemented as for the planon-lineon code. The Majoranas involved in γ 1 γ 2 γ 3 γ 4 and γ A γ B γ 5 γ 6 are placed onto two floating SCI as shown in Fig. 2. The charging energy of each SC island fixes their fermion parity. Hybridizations it(γ 1 i γ 5 i + γ 2 i γ 6 i ) effectively generate the remaining four-Majorana interaction γ 1 i γ 5 i γ 2 i γ 6 i , thereby exactly reproducing the interaction in Eq. (6).

The octahedral Chamon code
We now consider a system of hybridized Majoranas on a face-centered cubic lattice with 12 Majoranas on each site. Each of these Majoranas, placed on the edges of a cube, pairs with their partner on one of the nearest-neighbor sites as illustrated in Fig. 3.
For each site, we place the 12 Majoranas onto three SCI as shown in Fig. 3. The charging energy of the SCI fixes the fermion parities η 1 η 2 η 9 η 12 = η 4 η 3 η 10 η 11 = η 7 η 6 η 5 η 8 = −1. Weak tunneling between the different SCI islands, generates the effective Majorana interactions η 1 η 4 η 12 η 11 and iη 9 η 5 η 7 η 10 η 12 η 11 . The full set of interactions projects each site into a single spin-1/2 degree of freedom and the corresponding Hamiltonian becomes the octahedral Chamon code [9], The six-spin terms associated with the octahedra originate from products of 12 Majorana pairs on the hinges (see SM for details). The Chamon code is a type-I fracton with quasiparticles of restricted mobility. The ground state degeneracy of this model depends on the greatest common divisor of the three spatial lengths.

Fractons from Majorana fermion codes
Our fracton construction from interacting hybridized Majoranas suggests that many bosonic fracton models are equivalent to 3D Majorana fermion codes [30]. By enlarging each site of the bosonic model into a cell, the stabilizers of the Hamiltonian can be expressed as products over the Majoranas located on the corners of certain cells or plaquettes, provided that similar to color codes [34], the plaquettes and cells share an even number of Majoranas. We make this argument explicit for the planon-lineon code in Fig. 1. We replace each site by a small cube, placing the eight Majoranas at the corners of the cube. The octahedra are then replaced by truncated octahedra as shown in Fig. 4 which share a face with cubes below or above and a hinge with a side cube. The spin Hamiltonian in Eq. (2) is replaced by the Majorana fermion code, The first term is a product over the 16 Majoranas at the corners of a truncated octahedron. The remaining terms involve products of four Majoranas on the side plaquettes (xz or yz) of the cube. This Majorana fermion code defines a commuting projector Hamiltonian which exactly reproduces the ground-state manifold of the planon-lineon code.
The same stratagem works for the other fracton codes which we discuss in this paper. While a wide variety of 3D Majorana fermion codes exhibiting Z 2 × Z 2 × Z 2 order have been proposed [35], our construction suggests that numerous fracton topological orders can also be represented by Majorana fermion codes akin to color codes. In particular, this reveals the equivalence between Majorana fracton codes and bosonic fracton codes.
Fracton spin liquid.-Now imagine we have four Majoranas per site, with each flavor forming a Majorana fermion code as given in Eq. (9). By imposing a strong onsite inter-flavor interaction U γ 1 γ 2 γ 3 γ 4 , each site is constrained to even fermion parity and the corresponding low-energy Hilbert space is reduced to a spin-1/2 degree of freedom. The resulting spin Hamiltonian is the SO(3) invariant fracton spin liquid One can also take 2N Majoranas per site and apply interflavor interactions to obtain an SO(2N-1) invariant Hamiltonian, yielding a hierarchy of spin liquids with fractonic behavior.

Discussion and experimental implementation
We have demonstrated the emergence of fracton phases of matter in systems of interacting Majorana band models. As discussed in Ref. [36], the onsite interactions gauge the subsystem fermion parity symmetry and the resulting degrees of freedom obey a higher-rank Z 2 gauge theory describing a fracton phase of matter. On the theory side, our construction exposes the relation between fracton models and Majorana fermion codes. By extension, such Majorana fermion codes enable us to generate a hierarchy of fracton spin liquids with SO(N) symmetry. In addition to the examples dis- cussed here explicitly, interacting Majorana wires can also realize additional models of fractons in particular and topological phases in general, including the checkerboard model proposed in Ref. [36], the 3D toric code, the X-cube model, the cuboctahedron code, and the 2D fractal spin glass (see SM for details).
On the experimental side, this provides a platform for exploring these novel strongly interacting phases. We conclude this paper by discussing opportunities for measuring and manipulating fracton phases of matter.
So far, we have focused attention on realizing exactlysolvable stabilizer Hamiltonians of fracton codes. Sufficiently strong transverse fields applied to the local spin degrees of freedom drive the fracton phase into a confined phase [37]. In our platform, such transverse fields are readily implemented by adding additional local Majorana hybridization terms as illustrated in Fig. 5. It would therefore be interesting to probe these phases and the intervening phase transition experimentally. While fracton phases lack a local order parameter, their phases and phase transitions are characterized by non-local string or membrane order parameters. For instance, the expectation values of certain planar Wilson loop operators may obey perimeter or area laws depending on the phase [37]. Such expectation values can be read out by repeated preparation of a ground state and subsequent projective measurements of, say, the σ z components of all involved sites. Averaging the results for the string operators over the repeated measurements then provides access to the desired expectation value. Such projective measurements of local spins correspond to the measurement of local two-or four-Majorana parities for which several readout schemes have been proposed [26,27]. A promising scheme to measure, say, a two-Majorana parity couples a single-level quantum dot to the two Majoranas. The resulting energy shift of the quantum dot level will then depend on the Majorana parity. This is illustrated in Fig. 5. The same scheme can be extended to measure four-Majorana parities [26].
When applied only to subsets of sites, transverse fields can also be used to implement synthetic twist defects [38]. For instance, adding intra-island Majorana hybridizations along a defect line, one can effectively create a twist defect of the fracton code which permutes different types of quasiparticles. This realizes quasiparticles of the 3D fracton phase which exhibit projective non-abelian statistics [38]. The experimental tunability of the Majorana hybridization enables control over their location, an important prerequisite for probing their nonabelian braiding properties. In addition to the static expectation values discussed above, one can also extract dynamic correlation functions of the fracton codes from linear-response measurements. Applying a time-dependent transverse field and reading out the spins at a later time provides access to spin-spin correlation functions of the fracton code. The dynamic properties of 3D fracton codes are particularly interesting due to their glassy dynamics which results from the restricted quasiparticle mobility. Moreover, disorder in the stabilizer flip energies or the transverse fields would would allow for studies of many body localization. Both kinds of disorder are easily tunable in this setting as inter-site hybridization and onsite stabilizer energies are controllable via gate-tuned coherent links.

Methods
Traditionally, any spin model can be converted into a Majorana model by expressing the Pauli operators associated with each spin in terms of Majorana bilinears, made up from four Majorana operators γ 1 , γ 2 , γ 3 , γ 4 , and fixing the Majorana parity γ 1 γ 2 γ 3 γ 4 = −1 on each site. In view of the Majorana anticommutation relations the Majorana bilinears satisfy the same algebra as the Pauli matrices. When applying this conventional fermionization to fracton codes, one again obtains complicated many-body interaction among fermion clusters. These are obviously difficult to implement, and the nontrivial challenge is thus to find a special Majorana representation which maps the spin cluster model to a Majorana model with local hybridization and merely onsite interaction implemented by charging energy. As we show above, this requires more involved Majorana representations of the spin degrees of freedom of the fracton codes, for instance involving also four-Majorana terms to represent Pauli matrices. We now sketch the procedure how the spin cluster interactions corresponding to the stabilizers of the fracton models can be obtained from local interactions involving only small numbers of Majoranas. In the strong coupling limit, these interactions define a nearly degenerate low-energy subspace, with nontrivial Hamiltonian terms appearing only in highorder perturbation theory in the inter-site Majorana hybridizations. We sketch the underlying Brillouin-Wigner perturbation theory using the example of how to obtain the planonlineon code from the strongly hybridized Majorana model in Eq. (1). Strong interactions in Eq. (2) enforce even fermion parity for each four-Majorana cluster. The effective Hamilto-nian in this even-parity Hilbert space is obtained in seventhorder perturbation theory and takes the form Defining spin operators through σ z i = γ 1 i γ 2 i γ 4 i γ 3 i and the product of two Majoranas associated with any vertical edge as σ x i , the Hamiltonian becomes the planon-lineon code in Eq. (3) with a six-spin cluster interaction on each octahedron. Any higher order perturbation terms are just products of the stabilizers in the Hamiltonian which do not change the ground state manifold. Analogous arguments apply for all other fracton codes which we obtain from interacting Majorana models.

Acknowledgments
This work was supported in part by a PCTS fellowship (YY) and by CRC 183 of the Deutsche Forschungsgemeinschaft (FvO), and was initiated at the Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611. FvO is grateful for sabbatical support from the IQIM, an NSF physics frontier center funded in part by the Moore Foundation. YY acknowledges the hospitality of IQIM-Caltech, where part of this work was completed. Competing interests: The authors declare that they have no competing interests. Author contributions: Both authors contributed to the theoretical results and experimental proposals presented here. Data availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data available from authors upon reasonable request.

SUPPLEMENTAL MATERIAL I. 3D TORIC CODE AND X-CUBE MODEL
Following the strategy in the main text, one can readily obtain the 3D toric code and the X-cube model from a similar Majorana construction. As pointed out in Ref.
[39], both models can be obtained via coupled layer constructions based on the 2D toric code. Moreover, the 2D toric code can be obtained via a Majorana network construction [20,36].
We briefly review the Majorana construction for the 2D toric code. Consider a Majorana model on a square lattice with four Majoranas per site. Each Majorana hybridizes with its closest neighbor, as shown in Fig. S1 and described by the Hamiltonian Here, the lattice sites are enumerated by j and connected by lattice vectors e r and e r . Each site can be view as a SCI whose charging energy fixes the site's fermion parity, γ 1 j γ 2 j γ 3 j γ 4 j = −1. The low-energy Hilbert space reduces to an effective spin-1/2 degree of freedom on each site for which one can define Pauli operators (see, e.g., [36]). Treating the Majorana hybridizations as perturbations, the leading-order Hamiltonian involves Majorana hopping terms around all elementary plaquettes. In terms of spin operators, the plaquette terms form a checkerboard pattern and become σ z i σ z j σ z k σ z l and σ x i σ x j σ x k σ x l for white and yellow plaquettes, respectively (see Fig. S1). This is just the toric (or surface) code.
We can now implement the coupled layer construction by placing two SCI on each link of a cubic lattice. Each SCI hosts four Majoranas forming a 2D Toric code on each ij layer. The model resembles decoupled 2D toric codes with each i-link holding qubits from the ij and ik layers as shown in Fig. S1.
We can now add additional tunneling terms between Majoranas on same site but different SCI, see Fig. S1. We consider two separate cases, described by the two Hamiltonians Here a, b labels the two SCI on the same site. These tunneling terms, together with the fixed fermion parities of each SCI, generate effective onsite XX or ZZ interactions between spins from perpendicular layers, but on the same link. This couples the perpendicular 2d toric code layers. In the strong coupling limit, the interaction H B leads just to the 3D toric code model, The Hamiltonian involves six-spin vertex interacts involving σ z and four-spin plaquette interactions via σ x . This model exhibits 3D Z 2 topological order with nontrivial particle loop braiding. Likewise, the interaction H A leads to the 3D X-cube model II. TRIANGLE ISING MODEL WITH FRACTAL SYMMETRY BREAKING Symmetries are indispensable for characterizing different phases of matter. Typically, one deals with global symmetries, whose operation acts extensively on an the volume of the system. Fractal subsystem symmetries, which act only on a subset of sites whose number scales with linear size L as L d with some fractal dimension d have attracted much attention with the recent developments on fracton topological order. Systems with such symmetries appear most notably in the context of glassiness. Examples are the triangle or tetrahedral Ising models, whose Hamiltonians have fractal Z 2 symmetry when flipping spins on arbitrary Sierpinski triangles. We now show how to utilize a crossing Majorana network to generate the 2D triangle Ising model whose low energy ground states exhibit fractal symmetry breaking [32]. To implement the onsite fermion interaction in Eq. (S5), we place four and two Majoranas on one SC island each as shown in Fig. S2. While the two-Majorana island is grounded, the left SC island with four Majoranas is floating so that its charging energy fixes the parity η 1 η 2 η 3 η 4 = −1. Weak tunneling it(η 3 η 5 + η 4 η 6 ) between the islands effectively creates the other four-Majorana interaction η 3 η 4 η 5 η 6 .

III. THE CUBOCTAHEDRON FRACTON CODE
We now consider the same structure of crossing Kitaev chains as in Eq. (7), yielding a face centered cubic lattice. At each site, there are twelve Majoranas, each of which pairs with a partner at a nearest neighbor site.
We place the 12 onsite Majorana onto three SCI as shown in Fig. S3. The charging energies of the three SC islands fixes the fermion parities η 5 η 6 η 7 η 8 = η 1 η 3 η 10 η 12 = η 2 η 4 η 11 η 9 = −1. Turning on tunneling between the SC islands, H t = it(η 5 η 1 + η 12 η 8 + η 3 η 4 + η 9 η 6 ) (S8) generates effective Majorana interactions η 1 η 5 η 8 η 12 and iη 3 η 4 η 8 η 12 η 9 η 6 . When combined with the parity fixing for each SCI, these interactions project each site into a single spin-1/2 degree of freedom and the Hamiltonian becomes the cuboctahedron code as shown in Fig. S3. Each stabilizer is a 12-spin interaction on the cuboctahedron. The cuboctahedron is made of six corner-sharing plaquettes from the xy,yz,xz planes. The 12-spin term on the cuboctahedron originates from the product of the 24 Majorana pairs on the hinges. This fracton code belongs to the family of type-I fracton codes whose quasiparticles have restricted mobility. Similar to the Chamon code, the number of ground-state degeneracy on a three-torus is dependent to the greatest common divisor of the system size. The construction of the cuboctahedron code shares many similarities with the checkerboard code which we discussed in Eq. 8. The geometry of this Majorana network corresponds to four intersecting plaquettes sharing a corner at a site. The four-Majorana projections γ 1 γ 3 γ 10 γ 12 , γ 5 γ 6 γ 7 γ 8 , γ 2 γ 9 γ 11 γ 4 produce three intersecting Wen-plaquette models on the xy, yz, and xz planes labelled by spin qubits (X 1 , Z 1 ), (X 2 , Z 2 ), (X 3 , Z 3 ). The remaining six Majoranas and four-Majorana interactions couple the three intersecting Wen-plaquette models via an anyon condensate. By imposing Y 1 Y 2 Y 3 = 1 and Z 1 X 2 = 1 on each site, the three spin qubits are reduced to one Pauli spin degree of freedom represented by σ x = X 1 Z 2 , σ y = X 2 Z 3 , σ z = X 3 Z 1 . The three corner-sharing intersecting Wen-plaquette models form an octahedron which exactly reproduces the Chamon code in Eq. 8. Likewise, if we impose Y 1 Y 2 Y 3 = 1 and X 1 Z 2 = 1 on each site, the three spin qubits are reduced to one Pauli spin degree of freedom represented byσ x = Z 1 X 2 ,σ y = Z 2 X 3 ,σ z = Z 3 X 1 . This exactly reproduces the cuboctahedron structure which ties the six Wen-plaquette models on the six plaquette-faces of the cuboctahedron. Our construction, en route, reveals the relation between 2D coupled toric code layer and 3D Chamon code from anyon condensate.

IV. SO(2N-1) FRACTON SPIN LIQUID
We place 2N Majoranas on each site, with each flavor forming a Majorana fermion code as given in Eq. (9). By imposing a strong onsite inter-flavor interaction U γ 1 ...γ 2N , each site is constrained to even fermion parity and the the corresponding low-energy Hilbert space is reduced to spin-(2N − 1)/2 degree of freedom. One can use the Clifford algebra for Sp(2N-1) representation, The resulting spin Hamiltonian is the SO(2N-1) invariant fracton spin liquid