Exact parent Hamiltonians of bosonic and fermionic Moore-Read states on lattices and local models

We introduce a family of strongly-correlated spin wave functions on arbitrary spin-1/2 and spin-1 lattices in one and two dimensions. These states are lattice analogues of Moore-Read states of particles at filling fraction 1/q, which are non-Abelian Fractional Quantum Hall states in 2D. One parameter enables us to perform an interpolation between the continuum limit, where the states become continuum Moore-Read states of bosons (odd q) and fermions (even q), and the lattice limit. We show numerical evidence that the topological entanglement entropy stays the same along the interpolation for some of the states we introduce in 2D, which suggests that the topological properties of the lattice states are the same as in the continuum, while the 1D states are critical states. We then derive exact parent Hamiltonians for these states on lattices of arbitrary size. By deforming these parent Hamiltonians, we construct local Hamiltonians that stabilize some of the states we introduce in 1D and in 2D.


I. INTRODUCTION
The Fractional Quantum Hall (FQH) effect is one of the most fascinating phenomena in strongly correlated electronic systems, in which the electrons of a twodimensional electron gas subject to a strong magnetic field form an incompressible quantum liquid supporting fractionally charged quasiparticle excitations. The understanding of this paradigm of topological order was in large part made possible by the discovery of analytical wave functions, such as the Laughlin's wave function [1], describing the electrons in a partially filled Landau level.
Since its discovery in 1987 [2], one FQH state has attracted a lot of attention : Unlike the states at filling factors with odd denominators, the ν = 5/2 FQH state with electrons occupying the second Landau level cannot be explained by a hierarchical construction based on the Laughlin's states [3]. This opens the door to the possibility of electron pairing and emergence of non-Abelian quasiparticle excitations. Indeed the leading candidate for the description of the ν = 5/2 FQH state is the Moore-Read "Pfaffian" state at filling 1/2 [4][5][6], describing the wave function of the electrons in the second Landau level. Moore-Read states have a wave function defined by [4] ψ(w 1 , . . . , where w i are the positions of the particles on the complex plane, 1/q is the filling factor and the magnetic length has been set to one. They support fractionally charged non-Abelian anyons possessing Majorana fermion states at zero energy [7,8]. These non-Abelian anyons have attracted a lot of attention due to their applications to topological quantum computation [9,10].
There has been a lot of recent interest in finding FQH physics in other systems. Indeed, it is a fundamental problem to understand the origin and properties of states supporting non-Abelian anyons. Finding different systems exhibiting the same physics is an important step towards an explanation of these phenomena, for which our understanding is still far from complete. Moreover, experimental realizations of FQH states and manipulation of their quasiparticles are a challenge. As such it is interesting to search for new possibilities for realising FQH physics and variations thereof experimentally. In particular, lattice models without Landau levels may pave the way towards experimental realization and investigation of FQH-like states, for example in optical lattices.
Two main approaches have been followed to recreate the FQH effect in lattices. In the first approach, one tries to mimic electrons in a magnetic field by replacing the fractionally filled Landau level by a nearly flat fractionally filled Chern band and adding local interactions [11][12][13][14][15]. Non-Abelian FQH states where found in such lattice models with topological flat bands [16][17][18][19].
In the second approach, instead of trying to reproduce the interactions between the electrons and a magnetic field, the focus is on reproducing the FQH wave functions on lattices. This approach first started with the introduction by Kalmeyer and Laughlin of the bosonic Laughlin state at filling fraction 1/2 on a square lattice [20], for which a parent Hamiltonian was later derived on the torus [21,22] and for more general lattices in the thermodynamic limit [23][24][25] (see also [26,27] for related results). Hamiltonians were also derived for non-Abelian chiral spin liquids with excitations with SU (2) k statistics [24], but only in the thermodynamic limit. For k = 2 this corresponds to a bosonic lattice Moore-Read state at filling fraction 1 on a spin-1 lattice [28]. Moore-Read states of bosons have also been considered on one dimensional lattices, were parent Hamiltonians have been obtained [29]. Other filling fractions and Moore-Read states on spin-1/2 lattices have not been introduced before.
In the continuum, a useful description of FQH states [4] uses wave functions expressed in terms of correlators of the related edge Conformal Field Theory (CFT). This description was extended to lattice systems in one and two dimensions [23,30] by writing the wave functions of spin systems as CFT correlators. The states obtained can be seen as an infinite dimensional version of Matrix Product States (MPS) which gives a unified treatment of 1D and 2D lattices systems. An alternative way of implementing MPS ideas to 2D FQH systems has been developed in [31]. For the Laughlin states, this construction of wave functions written as CFT correlators provides states that are close to, but not exactly the same as the Kalmeyer-Laughlin states on a lattice of finite size, but that become the same in the thermodynamic limit [23]. This modification of the lattice construction has made it possible to construct exact parent Hamiltonians for strongly interacting lattice spin systems of arbitrary sizes [23,[32][33][34][35], including topological FQH states such as Laughlin states of hardcore bosons and fermions [25]. It was shown in some cases that these states could be stabilized by a local Hamiltonian [26,30,35,36]. This description has been applied to find a Hamiltonian for the SU (2) 2 spin models but the focus was on the one dimensional spin chain [32] and the corresponding two dimensional non-Abelian FQH state has not been analyzed so far.
In this paper we fill this gap by extending the construction of lattice wave functions from CFT correlators to non-Abelian FQH states. Using this approach we construct a family of lattice versions of the Moore-Read state at filling fraction 1/q. This family of states allows us to interpolate between the continuum limit, where all states become continuum Moore-Read wave functions, and a lattice limit. The states are defined on arbitrary spin-1/2 and spin-1 lattices in one or two dimensions, where the value of each spin can be mapped to an occupation number of bosonic (odd q) or fermionic (even q) particles. We provide numerical evidence that the states are critical states in one dimension and that in two dimensions they are topologically ordered states. It is shown that states defined on a spin-1/2 square lattice (for q = 2) and states defined on a spin-1 square lattice (for q = 1) have a topological entanglement entropy which is constant along the interpolation to the continuum. This suggests that the lattice states have the same topological properties as Moore-Read states in the continuum.
Since we use an approach based on the wave functions, it is then relevant to ask whether these states are ground states of physical Hamiltonians. Using properties from CFT, we derive parent Hamiltonians for the wave functions in the lattice limit. These parent Hamiltonians have long-range interactions, are exact on any lattice of arbitrary size and we find numerically that they have a nondegenerate ground space for q ≤ 2. By deforming these parent Hamiltonians, we then show numerical evidence that the state at filling fraction 1 on a spin-1 square lat-tice can be stabilized by a local Hamiltonian in one and two dimensions, while the state at filling fraction 1/2 on a spin-1/2 square lattice can be stabilized by a local Hamiltonian in one dimension, which is a first step towards an experimental realization of these states.
The paper is organized as follows : In Sec. II, lattice Moore-Read states are defined from correlators of conformal fields. It is shown in Sec. III that these states reduce to Moore-Read states of particles in the continuum. Properties of the states are computed in Sec. IV, where evidence that the states are critical in one dimension and that the topological properties remain the same along the interpolation between the continuum and the lattice states in two dimensions is presented. Parent Hamiltonians are derived in Sec. V. Finally Sec. VI provides numerical evidence that some of these Hamiltonians can be deformed into local Hamiltonians stabilizing the lattice states.

II. LATTICE STATES FROM CORRELATORS OF CONFORMAL FIELDS
Let us consider a lattice with N sites at positions z j , j ∈ {1, 2, . . . , N } in the complex plane. We will refer to one dimensional models when the z j are restricted to the unit circle and two dimensional models otherwise. Let a be the average area per site (in 1D a is the average distance between two sites on the unit circle) and η ≡ a 2π a positive real number. The local basis at site j is |n j , where n j is an integer that will be interpreted as the number of particles at site j. We consider two classes of models, labelled by S ∈ { 1 2 , 1} : • Models with S = 1 2 are defined on a Hilbert space of size 2 N : the two states in the local basis are n j ∈ {0, 1} and correspond to the absence/presence of a particle at site j. These states can be expressed in terms of spin 1 2 variables at each site : s i = 2n i − 1.
• Models with S = 1 are defined on a Hilbert space of size 3 N : the local basis is n j ∈ {0, 1, 2}, which we interpret as the presence of 0/1/2 particles at site j. These states can be expressed in terms of spin 1 variables at each site : In general, a wave function defined on one of these two spaces will have the form |ψ S = n1,...,n N ψ(n 1 , . . . , n N )|n 1 , . . . , n N S . (2) In the following, we will consider wave functions for which the coefficients ψ(n 1 , . . . , n N ) can be expressed as the correlator of some conformal operators. Let us introduce the operators where φ(z) is a chiral bosonic field from the c = 1 Conformal Field Theory (CFT), χ is a Majorana fermion field, : . . . : denotes normal ordering, q is a positive integer and δ nj is 1 if n j = 1 and 0 otherwise. We also define a phase coefficient where the ξ j are phase factors to be specified. We propose to consider the wave functions defined by the correlator of the previous operators : where These wave functions are labelled by three parameters : q, η, S and will be referred to as the (q, η) S CFT states.
Re(zj ) Illustration of a square lattice on the complex plane in the continuum limit (η → 0, N → ∞) and in the lattice limit(η → 1). At each site there can be 0 (blue circle), 1 (blue disk) or 2 (red disks) particles. The interpolation is performed by fixing the number of particles M = η N q and by varying η = a 2π between 0 and 1, which changes the number of lattice sites per particle between infinity and q.
Let us now evaluate the correlator. Note that we do not need to add a background charge in this setting. The correlator is zero unless N i=1 n i = ηN/q. This condition fixes the total number of particles in the system to M ≡ N i=1 n i = ηN/q. η/q is therefore the lattice filling fraction and η is a parameter which interpolates between the continuum limit (η → 0, N → ∞, number of particles conserved) with infinitely many lattice sites per particle and the lattice limit (η = 1) at which the lattice filling fraction is equal to 1/q, which corresponds to the Landau level filling fraction in the FQH effect (Fig. 1). Let us write δ n = 1 if N i=1 n i = M and δ n = 0 otherwise. The evaluation of the correlator yields [37] ψ S (n 1 , . . . , n N ) ∝δ n where f N (z l ) ≡ ξ l j( =l) (z l − z j ) −η and the Pfaffian is evaluated at the coordinates where n i = 1. The Pfaffian is antisymmetric, so these are bosonic states when q is odd and fermionic states when q is even. Note that this formula defines different states for S = 1/2 than for S = 1, the difference being that in the second case the n i can take the value 2. The states (q, η) 1/2 are therefore projections of the states (q, η) 1 onto the Hilbert space allowing only for single occupancy at each site, while states with S = 1 and q odd (resp. even) can have sites with two bosons (resp. fermions of different types). The state (q = 1, η = 1) 1/2 is trivial since the number of particles is fixed to N and there are N sites, so in the following we restrict ourselves to states (q, η) 1/2 with q ≥ 2 and (q, η) 1 with q ≥ 1 (Table I).

III. THE CFT STATES BECOME MOORE-READ STATES IN THE CONTINUUM LIMIT
In this section we consider a two dimensional lattice defined on a disk D of radius R → ∞ and show that the CFT states we have introduced reduce to Moore-Read states of particles in the continuum, that is Eq. (1), when η → 0, N → ∞ and the number of particles M is fixed. We restrict ourselves to lattices where the area per site a i is constant equal to a, but the derivation remains true for any lattice if we make η position dependent [25].
Let us first compute l f N (z l ) n l . Notice that |f N | = exp(− j( =l) η ln(|z l − z j |)) and since η = a 2π , in the continuum limit this sum can be replaced by an integral D ln(|z l − z|))dz 2 /2π. This integral evaluates, in the thermodynamic limit, to |z l | 2 + constant [4], so that where g l = Im j( =l) η ln(z l − z j ) is a real number. It was found numerically in Ref. [25] and Ref. [23] that this formula was an accurate approximation even for moderately large N. We thus get that In the rest of this section we set the phase factors such that ξ l = e ig l to get rid of the overall gauge factor, which does however not change properties like the particleparticle correlation function and entanglement entropy of the state.
A. Continuum limit of the S = 1 2 states Let us now write the complete wave function in the continuum limit : where the gauge factor has been set to one. It is not straightforward to take the continuum limit in this basis, since one has to define the limit of the Hilbert space on which the wave functions are defined. However, since the number of particles is conserved, we can rewrite the wave function in the basis spanned by the positions w 1 , . . . , w M of the particles. For S = 1/2 there is at most one particle per site so the wave function can be simply expressed as where the w i are restricted to positions in the lattice. In the limit of infinitely many lattice sites per particle the lattice becomes a continuous plane and the positions w i become positions in the plane. This state then coincides with the Moore-Read state (1). The number of particles on the lattice is M = η N q = a 2π N q , so if the flux is Area 2π , then we can express 1 q = M flux . This explains that this quantity corresponds to the filling fraction in the continuum, defined as the number of particles per magnetic flux.
B. Continuum limit of the S = 1, q = 1 state For S = 1, q = 1, the state also has the form (10), however since the n i can take the value 2, it is not straightforward to take the continuum limit. We first have to write the wave function in the basis spanned by the position of the particles. For a basis element |n 1 , . . . , n N , let w r , r ∈ {1, . . . , M } be the positions of the particles.
Since we interpret the state |2 as the presence of two particles, positions z i where n i = 2 are listed with two different indices r in the set {w r }. We now write the wave function in the basis given by the sets {w r }.
As a starting point, observe that We will now prove that Suppose first that only the first site is doubly occupied. Then, Moreover, where we have used in the last line that all other positions have at most one particle. Combining (14) and (16) we see that the (z 1 − z 1 ) term cancels, so that (13) holds. The same derivation with a recursion on the number of sites with two particles proves that (13) holds for all basis elements.
Putting together (13) and (12), the wave function can therefore be written as where the w r can be repeated twice to allow for states with double occupation, in which case this expression does not vanish because of a cancellation between the Jastrow factor and the Pfaffian. If we now take the continuum limit, the positions of the particles can be anywhere on the plane and this becomes the bosonic Moore-Read state (1) at filling fraction 1, which also does not vanish when two particles are at the same site.
One may ask what happened to the doubly occupied site. In the continuum the ensemble of states with two particles at the same positions has measure zero compared to states with at most one particle at each position, therefore they are irrelevant and do not contribute to the wave function. Similarly, when we take the continuum limit of the lattice, states with two particles at the same site do not contribute to the wave function, as we will now explain in more details. Intuitively, this comes from the fact that there is a finite number i n i = M of particles, so configurations with at least one doubly occupied site are rare when the number of sites goes to infinity. More quantitatively, let us denote P m the number of basis elements |n 1 , . . . , n N satisfying n i = M with m doubly occupied sites. P 0 = N M is the number of basis elements with no doubly occupied sites and In the continuum limit, M and m are fixed and N goes to infinity, so P m /P 0 ∼ K/N m where K is a constant and This shows that the number of basis elements with at least one doubly occupied site is small compared to the number of basis elements with no doubly occupied sites. Let us now observe that the wave function has to be normalized by a factor Q = n1,...,n N |ψ(n 1 , . . . , n N )| 2 . We can decompose Q between the basis of states with zero doubly occupied sites and the basis of states with at least one doubly occupied site : where the first sum contains P 0 elements and the second sum contains M/2 m=1 P m P 0 elements. We now observe that elements appearing in these two sums are of the same order. Indeed, let us take a configuration with M particles at positions w r such that 2 particles are at the same site (w 1 = w 2 ), corresponding to an element ψ(n 1 , . . . , n N ) in the second sum. Then by slightly moving one of the particles (w 1 = w 1 + δ) we obtain a new configuration that is in the first sum, where all particles are at distinct positions. Since the number of particles is fixed, continuity of the coefficients of the wave function (17) implies that the new element obtained is close to ψ(n 1 , . . . , n N ). Therefore elements in the second sum are of the same order as elements in the first sum and since there are M/2 m=1 P m P 0 elements in the second sum, in the continuum limit this sum does not contribute to the normalization factor Q. This shows that the contribution of configurations with two particles at one or more sites can be neglected in the continuum limit.
Note that the previous derivation also shows that the (q = 1, η = 1) 1 CFT state in the thermodynamic limit is equivalent to the spin 1 non-Abelian chiral spin liquid introduced in Ref. [28], but the two states are different on finite lattices. C. Continuum limit of the S = 1, q ≥ 2 state When q ≥ 2 the state can be written as We have already derived the continuum limit of the state on the right at q = 1 and the remaining factor can be expressed as where In the continuum limit, the wave function can therefore be written as When no two particles are at the same site, this wave function is the same as the (q, η) 1/2 CFT state in the continuum limit. As explained in Sec. III B, configurations with two particles at the same site do not contribute to the wave function in the continuum limit, therefore the (q, η) 1 CFT states for q ≥ 2 have the same continuum limit as the (q, η) 1/2 CFT states, which is the Moore-Read state at filling fraction 1/q.

D. One dimensional continuum limit
So far we have focused on two dimensional states. In the one dimensional setting, when z j = e 2πij/N , the same results enable us to perform an interpolation between the lattice and the continuum. The only difference with the 2D case is that it is now possible to compute analytically so the wave functions in the continuum can be expressed as which is a one-dimensional version of the Moore-Read state.

IV. PROPERTIES OF THE CFT STATES
In two dimensions, the Moore-Read states in the continuum are topological states which support non-abelian quasi-particle excitations. It is of high interest to check whether the lattice CFT states we have introduced share these properties. In one dimension, we expect that the CFT states display critical behaviour related to the conformal operators used to construct the wave function. In this section we focus on the states with q = 1 and q = 2, and numerically compute some of the properties of the states we have introduced.

A. One dimensional critical states
We now look at one-dimensional chains such that z j = e 2πij/N . Since we will find local Hamiltonians for the states (q = 1, η) 1 and (q = 2, η) 1/2 in one dimension in the lattice limit in Sec. VI, we focus on these states. First we compute the Renyi entropy S (2) L = − ln Tr ρ 2 L , where ρ L is the density matrix of the CFT states restricted to a subsystem of size L of the chain. This computation can be performed by using a Metropolis-Hastings algorithm with two independent spin chains [30,[38][39][40].
The results are shown in Fig. 2. The entropy scales logarithmically with the size of the subchain for all values of η. Moreover, the scaling is approximately the same for different values of η and fits of the form yield a value of the central charge c approximately equal to 1.36 for the (q = 1, η) 1 CFT states (the main source of errors is here the finite size of the lattice considered). This value is in agreement with the value of 1.395 found for the state in the lattice limit in Ref. [32], where it was also shown that a value of 1.5 for the central charge, as expected for the SU (2) 2 WZW model, could not be excluded. For the (q = 2, η) 1/2 CFT states we find a value of 0.98 which is compatible with a central charge equal to 1. Another quantity that can be computed using Monte-Carlo techniques is the particle-particle correlation function C L = n 1 n L − n 1 n L . Results in Fig. 3 confirm that the states are critical since the correlation functions decay polynomially with the distance L. For the (q = 1, η = 1) 1 CFT state, the critical exponent is found to be 0.70, which is in agreement with the value of 0.69 found in Ref. [32], where it was observed that such a value can be influenced by a multiplicative logarithmic correction which could explain the difference with the expected value of 0.75 [41][42][43]. Moreover for the (q = 2, η) 1/2 states at different values of η, the correlations are very close once rescaled by a factor of 1/η 2 , which confirms that properties of the state do not change along the interpolation. For a Tomonaga-Luttinger liquid, the expected L of a subsystem of L consecutive sites for the 1D (q = 1, η)1 CFT states (a) and (q = 2, η) 1/2 CFT states (b) for different values of η. The number of particles M = ηN/q is fixed so the sizes of the chain are N = 40, 80, 160, 320 for η = 1, 1/2, 1/4, 1/8 respectively. The lines are linear fits of the points for η = 1/8 (blue) and η = 1 (red).
behaviour of the particle-particle correlation function is [44] where K is the Luttinger parameter, k F = ηπ/q is the Fermi momentum and A is a non-universal constant. For the (q = 2, η) 1/2 state we find a good agreement of this formula for K = 0.494, A = 0.123. This suggest that this state in one dimension is well described by a Tomonaga-Luttinger liquid with central charge c = 1 and Luttinger parameter K = 0.5, which corresponds to the properties of a free-boson CFT with radius √ 2, as was the case for the corresponding one-dimensional Laughlin state [25].
To compute the topological entanglement entropy (TEE) of the CFT states, we map a L x × L y square lattice on the cylinder to the complex plane by choosing the positions of the lattice sites to be z j = e 2π((xj −Lx−1/2)+iyj )/Ly , where x j ∈ {1, . . . , L x }, y j ∈ {1, . . . , L y } (Fig. 4). We then cut the cylinder in two halves and compute the Renyi entropy of the first half. The size along the cut is L y and we use the behaviour of to extract the topological entanglement entropy γ (Fig.  5). The results for the states (q = 1, η) 1 (resp. (q = 2, η) 1/2 ) are in agreement with the topological entanglement entropy of a Moore-Read state at filling fraction 1 (resp. 1/2), γ 0 (1) = 1 2 ln(4) ≈ 0.69 (resp. γ 0 (2) = 1 2 ln(8) ≈ 1.04). Moreover the value of the TEE does not change with η, so topological properties of the states remain the same along the interpolation between the continuum and the lattice limit.
We observe however that the state (q = 2, η) 1 has a TEE close to zero and different than γ 0 (2). The TEE does not stay constant when η is changed, which is compatible with the expectation that its value is γ 0 (2) in the continuum limit : there must be a phase transition along the interpolation between the continuum and the lattice. The states with S = 1, q ≥ 2 can therefore define distinct lattice states from the states with S = 1/2, q ≥ 2, while having the same continuum limit. The insets are enlarged views confirming that the topological entanglement entropy stays the same when η is varied and that its value corresponds to γ0(1) (resp. γ0(1/2)) in the first two cases, while the topological entanglement entropy of the (q = 2, η)1 CFT state is close to zero in the lattice limit and close to γ0(1/2) in the continuum limit.

V. PARENT HAMILTONIANS
So far we have considered wave functions for lattice states. It is also relevant to ask whether these states are ground states of some Hamiltonians and whether these Hamiltonians can be realized in nature or implemented in experiments. We turn now to the construction of parent Hamiltonians for which the CFT states are ground states. The CFT states are constructed from correlators of conformal fields. As has been shown in Ref. [32], this enables the construction of parent Hamiltonians from null fields of the considered CFT. Null fields are fields such that when inserted in a correlator of primary fields, the evaluation of the correlator gives zero. The procedure is as follows : 1. Find null fields χ a (z i ) labelled by a and acting at position z i of the considered CFT.
2. The vacuum expectation value of a product of primary chiral conformal fields is zero if one of the fields is a null vector, therefore we get equations 3. These equations are rewritten in the form Λ a i |ψ = 0, where Λ a i are operators acting on the degrees of freedom of the wave function.

a,i Λ a †
i Λ a i is then a positive semi-definite Hermitian operator annihilating the wave function.
In this section we apply this procedure to construct parent Hamiltonians for the CFT states in the lattice limit. In the rest of this work the phase factors are fixed to ξ l = 1.
As shown in Ref. [32], the null fields in this case can be parametrized by where repeated indices are summed over, ϕ n is a chiral spin 1 primary field, J b −1 are the −1 modes of the SU (2) 2 current operators and where abc is the Levi-Civita symbol and t a i are the spin 1 operators acting on site i. These operators can be written in the spin basis at site i as We then exploit that the correlator with a null field inserted is zero : The Ward identity enables to transfer the action of a current operator to fields at other positions : Using (36) it is possible to rewrite (35) as [32] Λ a i |ψ = 0, where where we have defined ω ij = zi+zj zi−zj . This can be used to construct a parent Hamiltonian H = a,i Λ a † i Λ a i , which gives This Hamiltonian is SU (2) invariant and numerical diagonalization on small systems confirm that it has the (q = 1, η = 1) 1 CFT state as a unique ground state. This Hamiltonian is similar to the one obtained in Ref. [24] for the spin 1 non-Abelian chiral spin liquid state introduced in Ref. [28]. However the Hamiltonian we just constructed is valid for any choice of lattice and not only in the thermodynamic limit as is the case in Ref. [24].
B. Parent Hamiltonians for the (q, η = 1)1, q ≥ 2 CFT states For q ≥ 2, the (q, η = 1) 1 CFT states do not display an SU (2) symmetry. However it is still possible to find null fields and construct operators annihilating the wave function. Let us define the operators G ± (z) = χ(z)e ±i √ qφ(z) , J(z) = i √ q ∂φ(z). We use the following q + 1 null fields : where Eq. (41) gives q−1 null fields χ p (w) with p running from 0 to q − 2. Let us define the operators d † , d to be creation and annihilation operators (bosonic for q odd, fermionic for q even) acting between states |0 and |1 , and d † , d to be creation and annihilation operators acting between states |1 and |2 . The number of particles at site i is thus n i = n Using the previous null fields, we derive the following operators annihilating the wave function (see Appendix A for the complete derivation).
This leads to a three-body Hamiltonian annihilating the wave function where the last term fixes the number of particles. This Hamiltonian annihilates the wave function (q, η = 1) 1 , however we find numerically that the ground space of this Hamiltonian is degenerate when q ≥ 3 and that the degeneracy does not depend on the number of sites.
Other simple null fields of the theory constructed with the same current operators do not lead to operators acting on the wave function that would reduce this degeneracy.
C. Parent Hamiltonians for the (q, η = 1) 1/2 , q ≥ 2 CFT states We can now use the previous results to construct parent Hamiltonians for the (q, η = 1) 1/2 , q ≥ 2 CFT states. They are projections of the (q, η = 1) 1 CFT states in the subspace allowing only for single occupation at each site, that we denote H 1 . Let us also define the Hilbert space H 2 spanned by basis elements containing at least one site with two particles. We will now project the operators annihilating the (q, η = 1) 1 CFT states onto H 1 in order to get operators annihilating the (q, η = 1) 1/2 CFT states (see Appendix B for the detailed derivation). We start by multiplying the operators Λ a i obtained previously on the left by d i : these operators continue to annihilate the (q, η = 1) 1 CFT states. Since that act on H 1 and are zero on H 2 , so they also annihilate the (q, η = 1) 1/2 CFT states. We can then use the fact i . The resulting operator d i Λ q i then acts separately on H 1 and H 2 . We finally get operators annihilating the (q, η = 1) 1/2 wave function : This leads to a five-body Hamiltonian As in the previous case, this parent Hamiltonian has a single ground state only when q = 2 and this ground state is the (q = 2, η = 1) 1/2 CFT state.

VI. LOCAL HAMILTONIANS
The parent Hamiltonians we have derived involve three or five-body interactions between all sites on the lattice. These Hamiltonians would therefore be very difficult to implement in experiments. However in some cases it has turned out that states constructed from correlators of conformal fields had very high overlaps with ground states of local Hamiltonians [26,30,35,36]. This has lead to a protocol to implement one of these states in experiments [26,49]. In this section we show that there is a local Hamiltonian for which the ground state is close to the (q = 1, η = 1) 1 CFT state in one and in two dimensions and that this result is also true for the (q = 2, η = 1) 1/2 CFT state in one dimension.
A. Local Hamiltonians for the (q = 1, η = 1)1 CFT state In one dimension the case of the (q = 1, η = 1) 1 CFT state was studied in Ref. [32]. It was shown that this state has a high overlap with the ground state of the bilinear-biquadratic spin 1 Hamiltonian with periodic boundary conditions, when β = 0.3213. Note that this Hamiltonian includes the 2-body terms present in the parent Hamiltonian (40).
We now study the two dimensional case and build a local Hamiltonian from the parent Hamiltonian (40). The operators Λ a i contain 2-body interactions between sites i and j. We cut these operators by keeping only terms for which the sites i and j are nearest-neighbours on the square lattice. This leads to a local Hamiltonian with three-body interactions. In addition to these terms, we include the two-body interactions between next-nearest neighbours present in the parent Hamiltonian. All six terms included in our trial Hamiltonian, that we denote H 2D , are shown in Table II. Note that the coefficients of these terms in the exact parent Hamiltonian are not position-independent. In our local Hamiltonian, however, we choose them to be position-independent and invariant under rotations.
By exact diagonalization and optimization on these coefficients, we find that there is a local Hamiltonian for which the overlap | ψ H |ψ CFT | between the ground state and the (q = 1, η = 1) 1 CFT state on a 4 × 4 square lattice on the plane is 97.36%. Considering the size of the Hilbert space 3 16 ≈ 4 × 10 7 , this overlap is very high. Note that with the same parameters, the overlap is also above 98% on a 4 × 3 or on a 4 × 2 lattice. On a cylinder geometry, i.e. periodic boundary conditions in one direction, the overlap on a 4 × 4 square lattice is 97.21%.
Compared to the local Hamiltonian found in Ref. [28], which is for a state that is equivalent to the (q = 1, η = 1) 1 CFT state in the thermodynamic limit, but different on finite lattices, the Hamiltonian we find has less free parameters (5 instead of 11) to fine-tune, which might make it easier to implement. Moreover, the very good scaling with lattice sizes lets us expect that a good agreement will persist on larger lattices. In Fig. 6, we show the low-energy spectrum of this local Hamiltonian. This figure is compatible with having a gap in the thermodynamic limit, but the limitations on the system sizes that we can consider prevent us from making a reliable extrapolation. The state (q = 2, η = 1) 1/2 CFT state has a five-body parent Hamiltonian. Let us consider the one dimensional case. If we cut the Λ a i operators by keeping only terms for which the sites i and j are nearest-neighbours, we get a local Hamiltonian with several terms. We find however numerically that even a smaller number of terms is already sufficient to get a good overlap. Specifically, we choose to keep only the simplest two-body and the simplest three-body terms, to obtain a local Hamiltonian with periodic boundary conditions where d † i , d i are fermionic creation and annihilation operators at site i and n i = d † i d i . For κ = 0.274 + 0.052i, we find that the overlap between the ground state of H (2) 1D and the (q = 2, η = 1) 1/2 CFT state is 97.71% for a chain with 20 spins (Fig. 7(a)). with κ fixed to 0.274 + 0.052i, as a function of the number N of lattice sites. The dotted line is a linear fit with equation y = 1.0075 − 0.00152N . If the overlap continues to follow this behavior at larger sizes, it will still be above 85% for a spin chain with 100 lattice sites. (b) Overlap per site | ψCFT|ψH | 1/N between the same two states.
In two dimensions, since the SU (2) symmetry is not present in this model, cutting the parent Hamiltonian leads to a local Hamiltonian with up to five-body interactions with many different coefficients. In addition, the fact that the five-body terms involve more sites means that each of them stretches over a larger part of the lattice. With the limited lattice sizes that we can consider with exact diagonalization, this is problematic because the local regions need to be small compared to the total size of the lattice (otherwise it would not be expected that the same local Hamiltonian would also work for other lattice sizes). This suggests that even if a local Hamiltonian that is related to the exact Hamiltonian exists, we may not be able to find it with exact diagonalization. Instead of cutting the exact parent Hamiltonian, we therefore asked whether, by chance, a local Hamiltonian can be obtained if we restrict the range of the interactions to all interactions preserving the number of particles on all possible configurations inside a plaquette of the lattice. Optimizing the coefficients in this Hamiltonian, we did not, however, find a set of coefficients for these interactions for which the ground state of this Hamiltonian is close to the (q = 2, η = 1) 1/2 CFT state. Whether there exists a more complicated, but still local, Hamiltonian stabilizing this state therefore remains an interesting open problem.

VII. CONCLUSION
We have introduced a three-parameter family (q, η) S of strongly-correlated spin states on arbitrary lattices in one and two dimensions. It was shown that these states reduce to continuum Moore-Read states of bosons (odd q) and fermions (even q) in the continuum limit η → 0. Numerical evidence that these states are critical states in one dimension and topological states in two dimensions was provided, and the topological entanglement entropy was shown to remain the same along the interpolation for the (q = 1, η) 1 and (q = 2, η) 1/2 states. Parent Hamiltonians of the states in the lattice limit were derived using analytical tools from CFT and in some cases it was shown that these states could be stabilized by local Hamiltonians in one and two dimensions.
There is currently a lot of interest in finding models possessing topological properties. Given the complexity of quantum many-body systems, the analysis of phenomena, like e.g. topology, can be greatly facilitated by having models in which at least the ground state can be found analytically. The results of the present paper show that CFT is a valuable tool to derive analytical models, also in the context of non-Abelian FQH states, and in addition can be used as a starting point to identify simpler models that are more realistic to realize physically. which shows that χ q (w) is a null field. Similarly, there are more simple null fields χ p (w) p=0,1,...,q−2 = w dz 2πi 1 (z − w) p G + (z)V 1 (w), (A10) Let us now use the fact that by replacing the field at site i by a null field, the correlator vanishes : 0 = V n1 (z 1 ) · · · χ a (z i ) · · · V n N (z N ) . (A12) We will transform this equation into an equation involving the wave function by deforming the contour integral and moving the operators in the null fields at different positions. Let us do it for the null field χ q (w). We will use the OPEs as well as the commutation relations : V nj (z j )G + (z) = (−1) (q+1)nj −1 G + (z)V nj (z j ), (A13) where we have used : e iαφ(z) :: e iβφ(w) : = (−1) αβ : e iβφ(w) :: e iαφ(z) : and χ(z)χ(w) = (−1)χ(w)χ(z), which adds a minus sign only when n j = 1.
We then have, starting with the first term involving χ q 2 (w) :  We also define the particle number operators corresponding to these operators as n (1) j = d † j d j and n (2) j = d † j d j , such that n j = n (1) j + 2n (2) j . We multiply Eq. (A14) by |n 1 , . . . , n j−1 , n j , n j+1 . . . , n i−1 , 2, n i+1 . . . , n N , (A19) and sum over all n k , k = i, to get an expression involving the wave function Let us now look at the second term involving χ q 3 (w): = q j( =i) zj dz 2πi