## Variational wave-functions for correlated metals

Ankush Chaubey, 1 Harsh Nigam, 1 Subhro Bhattacharjee, 1 and K. Sengupta 2

1 International Center Of Theoretical Science, Tata Institute of Fundamental Research, Bangalore 560089, India. 2 School of Physical Sciences, Indian Association for the Cultivation of Science, 2A and 2B Raja S.C.Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.

We study a set of many-body wave-functions of Fermions that are naturally written using momentum space basis and allow for quantum superposition of Fermion occupancy, { n k } . This enables us to capture the fluctuations of the Fermi-surface (FS) - the singularly most important signature of a metal. We bench-mark our results in one spatial dimensions (1D) to show that these wave-functions allow for quantitative understanding of the Tomonaga-Luttinger liquid (TLL); computations of certain correlators using them can in fact be extended to larger systems sizes compared to conventional exact diagonalization (ED) allowing for a more systematic comparison with bosonization techniques. Finally we show that this basis may be useful for obtaining fixed-point wave-function for strongly correlated metals in dimensions greater that one. In particular, we study the case of coherent (equal) superposition of elliptical FS in continuum (2D) and on a square lattice. In case of the former, our variational wave-function systematically interpolates between the phenomenology of the Fermi liquid ground state, i.e. , finite single-Fermion residue at a sharp FS, to a non-Fermi liquid (NFL) with zero residue. In the NFL the jump in ⟨ n k ⟩ at the FS is replaced by a point of inflection (similar to a 1D TLL) whose contour is consistent with the Luttinger Theorem. In case of the square lattice, we find highly anisotropic distribution of the quasi-particle residue, which, at finite resolution has an uncanny resemblance to the Fermi-arcs, albeit at zero temperature, seen in the pseudo-gap state of the cuprates.

## I. Introduction

Our understanding of metallic phases squarely rest on the idea of Fermi-surface (FS). For non-interacting Fermions, it is sharply determined by the average single particle occupation in momentum space, ⟨ n k ⟩ , which in turn depends on the single-particle energy ε k -µ where µ is the chemical potential, and ε k is the single-particle dispersion. The shape and volume of the FS are determined by the underlying symmetry and filling respectively [1]. The corresponding ground state (GS) wave-function has a product structure in momentum space corresponding to non-trivial entanglement in real space [2].

In the presence of a moderate short-ranged four Fermion interaction, the above description survives and forms the backbone of Landau-Fermi liquid (FL) theory [3-5] in spatial dimensions d &gt; 1. Indeed the single particle residue, Z k , obtained from the jump in ⟨ n k ⟩ at the FS, and the associated Luttinger theorem [6, 7] provides a definition for the long-lived gapless electron-like Landau quasi-particles near the FS in such a FL which then can be described via an effectively weakly interacting framework. In a FL, 0 &lt; Z k F &lt; 1, such that the FS remains sharply defined via the position of the jump in occupation inspite of fluctuations in n k near the FS. One then expects that the corresponding GS is very well approximated via short-range entangled states in the momentum space across the non-interacting FS such that the real-space entanglement signatures are similar to that of the non-interacting case [8].

particle. We explicitly construct variational wavefunctions for such quantum fluctuating FS, which, we note, is different from the dynamical fluctuation of a classical membrane. Instead, in this case, the location of the membrane (the FS) has a uncertainty in the momentum space. This superposition is then closely related to the dimer coverings in Rokhsar-Kivelson (RK) models [9], but applied to momentum space for the FS. Such wave-functions can form a starting point for understanding the loss of the Landau quasi-particle leading to a socalled non-Fermi liquid (NFL) phase. Such phases generically arise in one dimensional (1D) interacting Fermion systems and are collectively known as Tomonaga Luttinger liquids (TLL) [10, 11]. The quasi-particle residue, for these TLLs, vanishes at the FS and is replaced by a point of inflection with ⟨ n k ⟩ ∼ | k -k F | p near the FS with p &gt; 0 being related to the Luttinger parameter [11] (also see below). Thus while the Landau quasi-particle is lost, the FS is still sharply defined in accordance with the Luttinger theorem [12]. In higher dimensions ( d &gt; 1) similar phenomena, i.e., loss of quasi-particles can occur with a well defined critical FS [13-15] which arises at Fermionic quantum critical point [13, 16-18] between the large and small FSs in heavy Fermion systems [18-22] or Mott transition in correlated metals [23, 24] where the entire FS disappears on approaching from the metallic side [13, 14].

In this paper, we show that quantum fluctuations arising from the coherent superposition of wave-functions with different momenta occupation { n k } can lead to a loss of FS along with the the loss of Landau quasi-

For d &gt; 1, effective field theory calculations [25] of Fermions (at finite density) interacting with gapless bosons near a metallic quantum critical point or in a gapless phase (such as in the U(1) quantum spin liquid with a spinon-FS [26, 27]) indicate that due to the presence of the tangential scattering with the bosons, even at ultra-low energies, the Landau quasi-particle is

killed possibly leading to a NFL [18, 28-31] particularly in 2D [25, 27, 32-38]. However, the fate of the infra-red (IR) fixed point of such theories is not quite settled due to lack of control over re-summation procedures; this is the case even in largeN f limit with N f being the flavour of Fermions. This calls for new ideas in understanding such a NFL state which is parent to a large number of diverse strongly correlated phenomena; this is the central motivation for the variational wave-function approach that we study here.

The variational wave-functions that we propose here explicitly capture the partial or complete loss of the single-Fermion residue signalling the loss of the Landau quasi-particle at low energies. In 1D, we show that we can re-interpret the existing numerical and fieldtheoretic understanding for the 1D TLL [11] for spinless Fermions within the above variational wave-functions (Figs. 1 and 2). This can be achieved by re-writing such wave-functions in the basis of fluctuating { n k } and comparing the average occupation, ⟨ n k ⟩ with predictions from bosonization. In addition, the entanglement signatures of such wave-functions can also be shown to match with those computed from TLL. Interestingly, following RK construction of dimer models, we can construct minimal projector Hamiltonians, albeit in momentum space, where the GS allows large fluctuations, both in 1D and 2D, in the occupation { n k } . These wavefunctions, at partial filling, do not break any symmetry and are devoid of sharp FS like signatures in momentum space. Therefore, in accordance with Lieb-SchultzMattis (LSM) [12, 39] theorem in 1D or their higher dimensional extensions [7, 40], they possibly correspond to GSs of gapless phases of correlated metals (in higher dimensions, in principle, gapped topological order is also allowed). We study the generic properties of such wavefunctions in 2D that correspond to smooth deformation of the FS as an example of fixed-point wave-functions for correlated metals without Landau quasi-particles and isotropically or anisotropically blurred FS. In particular, the anisotropic loss of the single-particle residue leads to Fermi-arc like phenomena, albeit at zero temperature, similar to underdoped cuprates and direction dependent Friedel oscillations.

In the rest of this paper, we provide the details of the construction of the variational wave-functions starting with the single particle distribution in momentum space and use it to provide initial insights into their nature via selected examples in 1D and 2D. In section II we provide the details of construction of the variational wavefunction which is a straightforward extension of superposing single-particle Slater determinants. Such superposition generically leads to fluctuations in momentum space occupation { n k } and this forms the right basis (Eq. 1) to capture the quantum fluctuation of the FS (Eq. 3) as is evident from the generic interacting Fermionic Hamiltonian, when written (Eq. 5) in this basis. In Sec. III we show that the new basis provides a useful momentumspace truncation scheme to understand the physics of

TLL in a system of spinless Fermions in 1D. In particular, we characterize the relation between the fluctuations of the FS and the loss of the Fermionic quasi-particle residue by quantitatively comparing with bosonization for larger system sizes (Figs. 2, 3 and 4). Having obtained a description of the TLL in d = 1, we turn to the general properties of the wave-functions of the type in Eq. 3. We mainly restrict ourselves to description of a class of effective projector models in Sec. IV whose GSs allow equal superposition of various { n k } distributions; we note that these wave-functions are qualitatively different from those describing a free Fermi gas GS. This naturally leads to possible fixed point Fermion wave-functions in Sec. V that are GSs similar to that of suitable RK-type projector Hamiltonians that describe strongly correlated NFLs; such NFL may be obtained near a metallic phase transition due to the Yukawa interactions with gapless bosons as discused above. Our study in Secs. IV and V may be relevant for some of these systems. We discuss various properties of a selected class of such wave-functions emphasising their NFL features. Various technical details and supporting plots are given in the different appendices.

## II. The variational wave-function

In order to define the wave-function, we consider a finite lattice of linear size N with periodic boundary conditions. In d -dimensions, the Brillouin zone (BZ) is d -tori with a grid of spacing 2 π/N along each direction. For concreteness we restrict ourselves to d = 1 2 while ex-, tension to higher dimensions are straightforward. With such a discretized momentum grid, we can now re-write the N p -particle occupation states (analogues to the Slater determinants) in the momentum-space occupation basis by defining an Ising variable θ k = ± 1 and using them to define

$$| \Psi [ \theta ] \rangle = ( c _ { k _ { B }$$

where the momenta are arranged from smallest to largest by assigning a particular ordering in which k 1 &lt; · · · &lt; k B over the entire BZ. Here c † k are spinless Fermion creation operators following usual Fermion algebra { c k , c † k ′ } = δ k k , ′ and { c k , c k ′ } = 0. In the above wave-function modes where θ k = +1 are occupied while θ k = -1 are empty. The above set of wave-functions form an orthonormal basis, i.e. ⟨ Ψ[ ] Ψ[ θ | θ ′ ] ⟩ = δ [ θ , θ ] [ ′ ] . The Slaterdeterminant form of the corresponding real-space representation and matrix elements of various many-body operator is easily represented in this basis and are shown in Appendix A for completeness. In particular, the total particle number and the net momenta are given by

$$\underset { \text{s of} } { \text{I we} } \sum _$$

Note that the second relation in Eq. 2 automatically implies that the momentum of a empty/filled band is zero. Generically, for a given filling ( ν = N /N p d ), the number of above basis states increases as N d C N p with N and is therefore of little value for numerical calculations. However, it provides an alternate way to construct N p -particle Fock space basis states in terms of their momentum space occupation. In terms of these new variables θ k , a single FS stands for a closed volume in BZ separating θ k = +1 (occupied) and θ k = -1 (unoccupied) regions such that we recover the familiar FS wave function.

Using the N p -particle basis states in Eq. 1, we can write down the generic state

$$| \psi \rangle = \sum _ { [ \theta ] | \sum _ {$$

where ψ θ [ ] is the amplitude for the state with the occupancy distribution [ θ ] and each basis state obeys Eq. 2. In this sense, Eq. 3 is a momentum space wave-functional which is equivalent to superposing Slater determinants. The above wave-functional is the central object in this work which captures the idea of quantum fluctuation of the FS. However, at this stage Eq. 3 is too generic; in the remaining sections, we therefore focus on its application to particular cases.

To this end, we consider an interacting four (spinless) Fermion Hamiltonian H = H 0 + H I , where

̸

$$H _ { 0 } \ = \ & \sum _ { k } E _ { k } c _$$

where E k = ( ε k -µ ) is the single-particle dispersion and V q is the interaction kernel for the density-density interaction. In the basis of Eq. 1, the above Hamiltonian can be written in the N p -particle sector as

$$\begin{array} {$$

where the first term in the expression of H denotes the non-interacting part which, as expected, is diagonal in the momentum space occupation basis. The details of the form of E θ [ ] is given by Eq. B1 in Appendix B. From Eq. 5, one can easily see that the [ θ ] configuration which yields minimum energy is given by (within the constraint of Eq. 2) choosing Min ( E θ [ ]). The corresponding distribution [ θ ] gives rise to the FS of the free electron gas which typically is a smooth but sharply defined surface in the thermodynamic limit ( N → ∞ ) dictated by the ultra-violet(UV) symmetry and energetics via the form of E k .

The second term in Eq. 5 corresponds to the fourFermion interaction which has off-diagonal contributions.

It's explicit form is given by Eq. B2. Thus interactions lead to generic superposition of distributions of [ θ ] having the same filling and total momentum (Eq. 2). It is easy to see that in the case of the four Fermion interactions, the matrix elements of H I are non-zero only between two distributions of [ θ ] that differ by occupation at four momenta (up to overall momentum conservation); this has been shown explicitly in Appendix B. The matrix corresponding to H I therefore generically leads to a sparse N d C N p × N d C N p matrix in this basis.

Other possible alternatives forms of the interaction stem from the electrons interacting with a gapless bosonic mode, ϕ q , (phonon/gauge field/collective modes) via typically a Yukawa type interaction ∼ ∑ k q , V ( q ) ϕ c q † k + q c k ; these interacting Hamiltonians also have a sparse offdiagonal structure once the bosons are integrated out. Such Yukawa interactions are renormalization group (RG) relevant and can lead to NFLs in d=2 [25, 27, 3238]. The central distinguishing feature of such Yukawa interactions in comparison to the short range four-Fermion interactions is the presence of tangential scatterings; their RG relevance indicates the collective effect of such scattering which would necessarily change the Fermion distribution about the non-interacting FS. This leads to the quantum fluctuation of the FS as suggested by the form of Eq. 5.

Below, we assume a generic form of ⟨ Ψ[ ] θ | H I | Ψ[ θ ′ ] ⟩ and understand the nature of the resultant GSs. We note that the simplest case where the four-Fermion interactions can already lead to a TLL [10, 11] in d = 1. So, before we proceed to study the generic properties of the wave-function in Eq. 3, we benchmark it to show that the TLL can indeed be understood in terms of fluctuating FS. Note however, in d = 1, fluctuations in the occupancy necessarily leads to disconnected segments of { n k } which is not the generic case in higher dimensions.

## III. Interacting Fermions in 1D: The TLL

̸

At finite but moderate V q the Hamiltonian in Eq. 4 supports a TLL. This is realized, in particular, in a microscopic model of spinless Fermions at half filling ( N /N p = ν = 1 2) consisting of nearest neighbour hopping such / that the single particle dispersion (in Eq. 4) is given by E k = -2 t cos k . In what follows, we shall use an interaction potential V q = V e 0 -q 2 / (2 q 2 0 ) , where q 0 denotes an upper cutoff for momentum around the FS. This, in real space corresponds to a density-density interactions, i.e. , ∑ i,j,i = j n n V i j | i - | j , where V ( x ) = V 0 exp[ -x / 2 (2 a 2 )], where a ∼ 1 /q 0 is the lattice spacing. We note that the nearest neighbor density-density interaction is obtained from this model in the limit a → 0.

We use exact diagonalization (ED) technique to study the GS properties of these interacting Fermions on a finite size chain (for details, see Appendix C) for nearest neighbour density-density interaction V 0 ∑ i n n i i +1 . The

FIG. 1. (a) The probabilities of the GS wave-function obtained from ED of the half filled Hubbard model ( N = 14 , t = 1 , V 0 = 0 3) for spinless Fermions (see main text and Appendix C for details) in the basis given by Eq. 1 (horizontal axis). . The relevant momenta in the first BZ are denoted by k = -π + απ/ 7 with α = 0 1 , , · · · , 13. In this labelling scheme, the six states with leading contributions are shown as [0] ≡ [00001111111000] being the non-interacting GS and the rest corresponding to the particle-hole excitations above it. These are [1] = [00010111110100], [2] = [00011011110010], [3] = [00100111101100], [4] = [00011101110001] and [5] = [01000111011100]. The states are arranged in terms of the increasing number of particle-hole excitations which can also be quantified via their root mean square momenta deviation from the FS as defined in Table I in Appendix C. Note that the vertical scale has been truncated for improved visibility of the contributing basis states. The normalised weights are given in Table I. (b) The exact diagonalization done using the momentum space cut-off q c ( n c = 3 ( q c = 2 π N n c ) , n 0 = 1 5 . ( q 0 = 2 π N n 0 ) (see main text for details). (c) The cumulative weights of the different basis states. (d) The IPR for the different basis states. The IPR has been normalised with respect to the maximum contribution (vertical axis truncated for better visibility) from [0].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-cd8d94b2-4ebe-45de-a188-280826ca01d0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-b22d737c-0943-4139-90e8-62d807797577.jpg)

GS for a finite chain can indeed be written as a superposition of the Fermionic wave-functions as envisaged in Eq. 3. This is shown in Fig. 1 for N = 14 and ν = 1 2 / where we plot the overlap probabilities, | ψ θ ( ) | 2 , of the exact GS wave-function, obtained using ED, with | ψ θ ( ) ⟩ . Fig. 1(a) shows the the six leading order contributions to | ψ θ ( ) | 2 obtained through real space ED. These states, as expected, consist of the non-interacting GS and particlehole excitations above it. Further details of these overlaps are shown in Table I in Appendix C which include up to 25 basis states ordered in terms of their increasing RMSmomenta and indicating number of particle-hole excitations about the FS. The cumulative weights of these states and their inverse participation ratios are shown in Fig. 1(c) and (d) respectively.

It is clear from Fig. 1 that only the basis states | Ψ[ ] θ ⟩ whose distribution in [ θ ] differs from that of the noninteracting FS (denoted by [0] in Fig. 1) only near its FS has appreciable overlap with the interacting GS wavefunction. In fact Figs. 1(c) and (d) show that the weight of the GS wave-function is saturated appreciably solely by basis states with particle-hole excitations near the non-interacting FS. Taking cue from this, we performed the ED directly in the occupation basis (Eq. 1) in momentum space (for details, see Appendix C). This is anal- ogous to momentum-space ED adapted to the FS by taking V q = V e 0 -q 2 / (2 q 2 0 ) with a soft ( q 0 ) as well as a hard cutoff on q = q c ˆ for x H . The resultant contributions of the basis states is shown in Fig. 1(b); these serve as a benchmark of the momentum space ED in terms of the weights of different configurations since they can be compared with their real-space counterparts. We can further benchmark the momentum space ED by comparing different correlators with those obtained from the real space ED. In this work we calculate the average momentum occupation, ⟨ n k ⟩ , the two-point connected correlator,

$$\langle \langle n _ { \mathbf r } n _ { 0 } \rangle \rangle \ = \ \langle n ( \mathbf r ) n ( 0 ) \rangle - \langle n ( \mathbf r ) \rangle \langle n ( 0 ) \rangle, \quad ( 6 )$$

and the bipartite von-Neumann entanglement entropy, S vN .

The plot of ⟨ n k ⟩ , shown in Fig. 2 (a), provides a comparison between the two methods. The results agree well and in fact the momentum space calculations with a cutoff about the FS can be used for much larger systems sizes as shown in Fig. 3 for N = 50 and q c = 2 π/ 10. This, therefore, can be compared with the standard bosonization results for ⟨ n k ⟩ which is devoid of a jump and has a form ⟨ n k ⟩ ∼ | k -k F | p with p = ( K +1 /K / ) 2 -1 with K being the Luttinger parameter [11]. An estimation of the point of inflection, obtained from such a comparison, is

FIG. 2. Comparison of (a) ⟨ n k ⟩ calculated from real space ED and k -space ED ( N = 14) and (b) the two-point connected correlator ⟨⟨ n n r 0 ⟩⟩ in real space ( N = 22). The other parameters are t = 1 0 . , V 0 = 0 3 . , n c = 3 ( q c = 2 π N n c ) , n 0 = 1 5 . ( q 0 = 2 π N n 0 ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-db211695-2b0d-4b54-899a-7b2f690d3798.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-29f55a76-45e2-4d4c-b3f9-04e0fb38a424.jpg)

FIG. 3. A plot of the Fermi distribution function ⟨ n k ⟩ as a function of k ; the inset shows the curve fitting of ⟨ n k ⟩ (near the Fermi point) with known bosonization result (near the Fermi point: ⟨ n k ⟩ = d sign( k f -k ) | k -k f | p +1 2 with / p = ( K + K -1 ) / 2 -1 [11] where K is the Luttinger parameter and d is a constant). For this plot, t = 1 0 . , V 0 = 0 6 . , n 0 = 1 5 . , n c = 5 ( q c = 2 π N n c ) and N = 50.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-03db1c42-e59b-4688-a35b-f91e4f273d04.jpg)

in accordance with the Luttinger theorem. Thus the basis in Eq. 1 effectively captures the TLL GS which has a sharp FS denoted by the point of inflection, but no Landau quasi-particle. This further confirms some advantage in using the basis in Eq. 1 and validates the fact that the low energy physics is dominated by FS fluctuations.

For the two-point correlator in real space (Fig. 2(b)) and the entanglement (Fig. 4), the advantage of our momentum space ED however is much less clear at this point. For example, while the two-point correlator can be calculated for larger system sizes, it suffers from spurious oscillations arising from the cut-off, q c , for parameters similar to Fig. 3. For the entanglement entropy, on the other hand, we are still restricted to smaller sizes due to calculational costs. In addition, also due to the cut-off, the estimation of the central charge (Fig. 4(b)) is somewhat crude. Therefore, while promising, our momentumspace ED scheme with cut-off around the FS requires more systematic investigation in future.

## IV. Exactly solvable toy model for GS with FS superposition

Having shown that the basis in Eq. 1 may be useful to capture the physics of TLL in 1D, in this section, we study another particularly useful question- are there parent Hamiltonians which allows massive FS fluctuation in the GS ? The positive answer is provided in a simple limit of the Hamiltonian in Eq. 5. This is obtained when, for a particular filling ν (= N /N P d ),

$$\frac { E [ \theta ] } { E _ { I } ( [ \theta ], [ \theta ^ { \prime } ] ) } = - 2 \left ( ^ { N ^ { d } } C _ { N _ { p } } - 1 \right ), \quad \quad ( 7 )$$

where E θ [ ] and E I ([ θ ′ ] , [ θ ]) are respectively the diagonal (non-interacting) and off-diagonal (interacting) matrix elements (Eq. 5). In this limit, the Hamiltonian then can be re-written as a sum of projectors h ([ θ , ] [ θ ′ ]) as

̸

$$H _ { \text{RK} } = \mathcal { C } \sum _ { [ \theta ], [ \theta ^ { \prime } ], [ \theta ] \neq [ \theta ^ { \prime } ] } h ( [ \theta ], [ \theta ^ { \prime } ] ), \quad \quad ( 8 )$$

where C &gt; 0 is a filling dependent positive constant and

$$\stackrel { \sim } { \stackrel { \sim } { 2 ( b ) ) } } \\ \text{r mo-} \\ ; \text{ this}$$

with h 2 ([ θ , ] [ θ ′ ]) = h ([ θ , ] [ θ ′ ]). We note that the condition in Eq. 7 may be obtained in systems with all-to-allinteractions.

FIG. 4. Comparsion of (a) Entanglement Entropy, (b) Central charge variation with 1 /N (we have done linear curve fit ( c = m (1 /N ) + c 1 ) for real space data) for the two ED methods. The system parameters are t = 1 0 . , V 0 = 0 3 . , n c = 3 ( q c = 2 π N n c ) , n 0 = 1 5 . ( q 0 = 2 π N n 0 ) and N = 22.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-97f539fb-0519-453b-a9b7-b6cc11bcb97f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-8e1362a3-ec3f-4d6b-af27-31e5580dca1c.jpg)

The Hamiltonian in Eq. 8 can be solved in the same spirit as the projector dimer model Hamiltonians by Rokhsar and Kivelson (RK) [9]. The difference between the two models comes from the fact that in the present case the [ θ ] basis is orthonormal unlike the hardcore dimer-covering basis of the quantum dimer models. The GS wave-function is then easily read off; it is given by the state which are annihilated by h ([ θ , ] [ θ ′ ]) for all pairs of ([ θ , ] [ θ ′ ]) and constitutes an equal and completely symmetric superposition

$$| \psi _ { \text{RK} } \rangle = \frac { 1 } { \sqrt { N _ { \text{RK} } } } \sum _ { [ \theta ] } | \Psi [ \theta ] \rangle, \quad \quad ( 1 0 )$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-bb79f118-cabc-4a91-8861-d5c62a4841c4.jpg)

with N RK being the number of states in the superposition. Notably ⟨ Ψ RK | ˆ n k | Ψ RK ⟩ = ν for all k ∈ BZ. Thus, an equal superposition of all | ψ θ ( ) ⟩ in momentum space leads to loss of any sharp structure and hence the FS. In fact, this is not very surprising because the condition in Eq. 7 mixes distributions of [ θ ] belonging to different momentum sectors and hence the resultant Hamiltonian does not have translation symmetry.

Hence, while it is interesting to study the properties of the generic projector Hamiltonian in Eq. 8 ( e.g. , we plot ⟨ n r ⟩ = ⟨ Ψ RK | ˆ n r | Ψ RK ⟩ in Fig. 17 in Appendix D along with few other details), we shall consider a sub-class of such superposed states where different total momentum sectors do not mix. This is obtained by considering distribution of [ θ ] k T with the same total momentum k T . We also put in an extra constraint that θ k = π = -1, i.e. , the k = π mode is always empty for convenience (see discussion in Appendix D). The GS is therefore an equal superposition of such basis states | Ψ[ ] θ k T ⟩ similar to Eq. 10 which we denote as | ψ k T RK ⟩ . It is fairly easy to show that ⟨ ψ k T RK | n r | ψ k T RK ⟩ = ν , i.e. , constant (Eq. D2); moreover, such a constrained RK GS is translation invariant unlike

FIG. 5. Plot of ⟨ n k ⟩ as a function of k for the constrained RK type projector GSs ( | ψ k T RK ⟩ ), here we have taken 20 sites for different k T for periodic boundary conditions (PBC). Note that n π (=0) is not shown as in the plot.

for the full RK case.

The average momentum-space occupation, ⟨ n k ⟩ , for such constrained RK GSs is plotted as a function of k in Fig. 5 for a 1D Fermion system with N = 20 and N p = 10, i.e. , ν = 1 2 for different / K T . The slope of ⟨ n k ⟩ is proportional to K T ; in particular it is flat for the zero momentum sector with a value tending to ν for large N . For finite systems ⟨ n k ⟩ shows oscillations (see Appendix D), but the amplitude of these oscillations vanishes with increasing system size. This indicates absence of sharp FS-like signature for | ψ k T RK ⟩ for large N even in presence of translation symmetry. Notably such symmetric wave-functions are not allowed to be unique at partial filling due to Lieb-Schultz-Mattis (LSM) theorem [39] and hence the above wave-function most likely corresponds to the GS of a gapless phase constituting a realization of a NFL in a toy model.

FIG. 6. Schematic Figure of an elliptical FS (Eq. E1) with major axis along k ∥ and eccentricity ϵ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-70d2fab8-e1e1-42b2-9292-2430d6a54de0.jpg)

## V. Fixed point wave-functions: Fermi surface superposition in two dimensions

The RK construction of the superposition of FS provides important clues to generating variational wavefunctions in higher dimensions. It shows that a generic superposition completely destroys all features of the metallic GS. This raises an interesting question - can we superpose FS but still partially retain properties of metals- both weakly and strongly correlated ones? In this section we investigate this question by constructing a series of variational wave-functions in 2D which interpolates between non-fluctuating FS in the free Fermigas limit and the systematically fluctuating FS (subset of Eq. 3); we expect it to capture the phenomenology of correlated metals. We note here that it is possible to write appropriate RK-like parent projector Hamiltonians for which some of the wave-functions (those with equal superposition) that we discuss below are exact GSs. Such wave-functions can be thought as a IR fixed point wavefunctions describing the GS of correlated metals.

One way to obtain such wave-functions is by choosing ψ θ [ ] such that various microscopic symmetries beyond translation are implemented manifestly via [ θ k ]. For example, under time reversal θ k → θ -k while for a generic spatial symmetry R , we get θ k → θ R [ k ] . To incorporate the symmetries we then need to choose ψ θ [ ] in Eq. 3 such that the distribution of [ θ ] and [ θ ′ ] that are connected by symmetry transformation are weighted equally, i.e. ψ θ [ ] = ψ ∗ [ θ ′ ] for TR and ψ θ [ ] = ψ θ [ ′ ] with spatial symmetries. We shall focus on the case of TR symmetric systems which is achieved by considering [ θ ] to be symmetric about k = 0 and ψ θ [ ] to be real. Other symmetries are implemented in a similar manner. We focus on two situations- (1) isotropic, and, (2) square lattice which have different spatial symmetries.

Finally, Eq. 1 contains arbitrary distributions of { n k } over the entire BZ subject to the constraints given by Eq.

2. We note that at low energies, various interaction terms can only cause scattering near the FS [4] similar to the case of TLL in the previous section. In contrast, in higher dimensional systems with Yukawa interactions that hosts NFLs, the RG relevant tangential scatterings [27] around the FS lead to a large number of Fermion-hole excitations. These excitations ultimately lead to the breakdown of the Landau quasi-particle. However, the proliferation of the such Fermion-hole excitation, particularly on integrating out the boson, would lead to change in the distribution of { n k } via low energy processes. For the rest of this paper, we assume that these low energy processes lead to smooth deformation of the FS and basis states of the form in Eq. 1 which represents such smooth deformation have non-zero amplitude in the GS wavefunction that we focus on.

Noting this, we consider the basic fluctuations of an elliptic distribution of [ θ ] which describes an elliptic FS at filling, ν , with eccentricity ϵ and major axis along ˆ in the k 2D momentum space. An example is shown in Fig. 6. We denote such a sub-set of these basis state wave-function (Eq. 1) as | ϵ, ν, ˆ k ⟩ (see Appendix E for details). The generic wave-function for superposition of such elliptical FS is given by

$$| \Psi _ { \nu } \rangle = \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon \, \int _ { 0 } ^ { \pi } d \hat { \mathbf k } \ \psi ( \epsilon, \hat { \mathbf k } ) | \epsilon, \nu, \hat { \mathbf k } \rangle, \quad ( 1 1 )$$

where the first integral is over the eccentricity between [0 , ϵ max ( ≤ 1)] for the elliptical distribution with major axis along ˆ and the second integral denotes the sum over k all possible directions of the major axis. To impose the filling constraint, all the participating ellipses for which the amplitude ψ ϵ, ( ˆ ) is non-zero, have same area and for k the rest of this work we set ν = 1 2. /

It is fairly easy to show that such wave-functions have uniform Fermion density in real space, i.e. , ⟨ Ψ e | n r | Ψ e ⟩ = ν , unlike the generic RK wave-functions (Eq. 10) (see Appendix E). However, not surprisingly, depending on the form of ψ ϵ, ( ˆ ), its single particle signatures can be rather k different from that of a FL. This is seen by calculating the single particle residue, Z ϵ, ˆ k , defined via the single particle correlator [41] for an elliptical FS, i.e. , a basis state | ϵ, ν = 1 2 / , ˆ k ⟩ ≡ | ϵ, ˆ k ⟩ (details in Appendix E 2),

$$Z _ { \epsilon, \hat { k } } = \lim _ { \delta \to + 0 } \frac { \delta ^ { 2 } } { \left ( \frac { k _ { \| } ^ { 2 } } { a ^ { 2 } } + \frac { k _ { \| } ^ { 2 } } { a ^ { 2 } ( 1 - \epsilon ^ { 2 } ) } - \mu \right ) ^ { 2 } + \delta ^ { 2 } }, \quad ( 1 2 )$$

where µ fixes the filling, ν , for the elliptic FS of eccentricity ϵ and area A = πµa 2 √ 1 -ϵ 2 and k ∥ ( k ⊥ ) are momentum components resolved along (perpendicular to) the major axis along ˆ k as shown in Fig. 6. For a single elliptic FS, the residue shows a jump at the perimeter of the ellipse as expected, signalling the presence of the FS which can then be used to define long-lived low energy quasi-particles. For the superposed elliptical FS (Eq. 11)

FIG. 7. (a) Plot of ⟨ n k ⟩ (in the vicinity of FS) as a function of k for Eq. 14 shown for ϵ 0 = 0 , ϵ max = 0 7 and different values . of Z ϵ 0 , ˆ k 0 = 0 8 (black), 0 5(red), and 0 2 (blue). . . . The Inset shows ⟨ n k ⟩ as a function of k ∈ (0 , π ). The jump in the occupancy of ⟨ n k ⟩ equals Z ϵ 0 , ˆ k 0 . (b) Plot of n k as a function of k for the isotropic NFL (Eq. 15) in the vicinity of FS ( k F = √ 2 π ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-214885fa-c7eb-4c80-8c15-6c62a32b47c3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-5e135482-d15b-4833-94d2-9cf76ce9196b.jpg)

then, the averaged residue is given by

$$Z = \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon \int _ { 0 } ^ { \pi } d \hat { \mathbf k } | \psi ( \epsilon, \hat { \mathbf k } ) | ^ { 2 } Z _ { \epsilon, \hat { \mathbf k } }. \quad \ ( 1 3 )$$

Importantly, Eq. 13 recovers the FL form for

$$\psi ( \epsilon, \hat { k } ) = & \sqrt { Z _ { \epsilon _ { 0 }, \hat { k } _ { 0 } } } \ \delta ( \epsilon - \epsilon _ { 0 } ) \ \delta ( \hat { k } - \hat { k } _ { 0 } ) + \sqrt { \frac { 1 - Z _ { \epsilon _ { 0 }, \hat { k } _ { 0 } } } { \pi \epsilon _ { \max } } }, \quad \widehat { \asymp }$$

for Z ϵ 0 , ˆ k 0 ∈ [0 , 1]. The resultant momentum space occupancy of the Fermions ⟨ n k ⟩ is shown in Fig. 7(a) where we chose ϵ 0 = 0 ( i.e. , circular FS) for convenience. The quasi-particle residue, by construction, is Z ϵ 0 , ˆ k 0 and represents the non-interacting limit for Z ϵ 0 , ˆ k 0 = 1. Further details are given in Appendix E 2.

The above form of the wave-function (Eq. 14) naturally leads us to investigate the case where the weight of the δ -function goes to zero such that

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-797df87a-9176-4f4a-8d37-f51b5349902c.jpg)

$$\psi ( \epsilon, \hat { \mathbf k } ) = \frac { 1 } { \sqrt { \epsilon _ { \max } } } \frac { 1 } { \sqrt { \pi } }. \quad \quad \ ( 1 5 ) \quad \ f u r$$

̸

The resultant plot of ⟨ n k ⟩ is shown in Fig. 7(b) for different extent of superposition of eccentricities obtained by changing ϵ max (further details in Appendix E 2 b). Notably for finite ϵ max , the jump in the occupancy is systematically smeared out as is evident by comparing the two panels of Fig. 7. This is further confirmed by calculating the residue using Eqs. 12 and 13 as shown in 21(b) for various ϵ max = 0 which shows that Z = 0. However, instead the residue is replaced by a point of inflection (see Fig. 20(d) in Appendix E 2 c) somewhat similar to the TLL in 1D. The position of the point of inflection as well as the finite resolution jump of ⟨ n k ⟩ is shown in Fig. 8. Notably the contour of the inflection point encloses an area in accordance with the Luttinger theorem in spite of the single particle residue being zero. Thus the presence of the inflection point sharply defines the position

FIG. 8. Plot of the single particle residue Z as a function of k x and k y for the superposed FS corresponding to the wavefunction in Eq. 15 for ϵ max = 0 2 and . δ = 0 12. Note that the . residue is uniform and is ≤ 1; it → 0 as δ → 0 for ϵ max = 0. The maximum value of Z (for finite δ ) is at the perimeter of circle of radius= √ 2 π (dashed contour) that also marks the position of the point of inflection of n k . See text for details.

̸

of the FS, while Z = 0 indicates the absence of Landau quasi-particles similar to the critical FS [13] or marginal FLs [42, 43]. Similar to the TLL, the average momentum mode occupancy is a power law, i.e., ⟨ n k ⟩ ∼ | k F -k | p with p ≥ 0 (see Fig. 22) and p increases monotonically with ϵ max . This power-law form is consistent with the scaling form suggested in Ref. [13].

Fig. 9 shows the two-point connected correlator W ( r , 0) = ⟨⟨ n n r 0 ⟩⟩ as a function of r = | r | (see Appendix E3 for detailed expressions) for the wave-functions given by Eq. 14 for ϵ max = 0 75 and different . Z ϵ 0 , ˆ k 0 (including

FIG. 9. Plot of the density-density correlation W ( r 0 , ) = ⟨⟨ n n r 0 ⟩⟩ as a function of r = | r | corresponding to state Eq. 14 for ϵ max = 0 75 and different values of . Z ϵ 0 , ˆ k 0 = 0 95 (black), . 0 80(red) and 0 00 (blue). . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-721f7e43-c9fd-4aca-80b0-4745037679b5.jpg)

FIG. 10. Plot of ⟨ n k ⟩ as a function of k , where k = ( k cos γ , k sin γ ) for the superposed state mentioned in Eq. 17 for different γ and ϵ max = 0 7. The inset shows . ⟨ n k ⟩ as a function of k in the vicinity of the FS.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-c01ef8f5-9125-40cb-bb14-7f56bf784c88.jpg)

Z ϵ 0 , ˆ k 0 = 0 whence we get Eq. 15). The only difference between the non-interacting limit ( Z ϵ 0 , ˆ k 0 → 1) and the completely superposed limit ( Z ϵ 0 , ˆ k 0 → 0) seems to be in the amplitude of oscillations while other features survive the loss of single-Fermion quasi-particle as signalled by the vanishing residue.

Finally, we move to the case of the square lattice where we only have a four-fold rotational symmetry. This should be reflected in the choice of the amplitudes in Eq. 11 via constraining them for different ellipses related by a four-fold rotation as discussed above. For illustrative purposes, in this case, we choose the major axis of the ellipses to lie along the Cartesian directions ˆ and ˆ such x y that we now have

$$\psi [ \epsilon, \hat { \mathbf k } ] = \sum _ { \alpha = \hat { \mathbf x }, \hat { \mathbf y } } \psi ( \epsilon ) \delta ( \hat { \mathbf k } - \hat { \alpha } ), \quad \quad ( 1 6 ) \quad \mathbf E _ { \text{var} }$$

where ψ ϵ ( ) is a smooth real (to ensure TR) function of ϵ such that Eq. 11 reduces to

$$| \Psi _ { s } \rangle = \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon \ \psi ( \epsilon ) \frac { ( | \epsilon, \hat { \mathbf x } \rangle + | \epsilon, \hat { \mathbf y } \rangle ) } { \sqrt { 2 } }. \quad ( 1 7 )$$

This is nothing but equally superposing ellipses of eccentricity up to ϵ max , but with major axis only restricted to horizontal and vertical directions with equal amplitude for any given ϵ . We note that dynamic lattice nematic fluctuations of the type discussed in context of cuprate superconductors [44] can possibly lead to such superpositions of the FS. Choosing different functions ψ ϵ ( ) result in superposing the above set of ellipses with varying amplitudes and leads to different wave-functions with manifest TR and lattice symmetries. Consider, for example, ψ ϵ ( ) = 1 √ ϵ max that correspond to an equal superposition of all such ellipses (see Appendix E 2 c for further details).

The momentum space occupation is shown in Fig. 10 for different directions γ (the angle with respect to k x ). As can be seen from the figure (see inset), in contrast to Fig. 7, both the position of the rapid change of ⟨ n k ⟩ and its form depends on γ (being sharpest for γ = ± π/ 4). This is reflected in the finite resolution quasi-particle residue (Eq. 13) plotted in Fig. 11 for different extent of superposition of eccentricities, ϵ max and finite resolution, δ (Eq. 12). Notably, on taking δ → 0 (see Fig. 21(c) in Appendix E 2) the residue goes to zero for all momentum points. However, unlike the isotropic case, it goes to zero anisotropically with the contrast of the anisotropy increasing with ϵ max ; this is revealed by comparing the finite resolution leftmost ( ϵ max = 0 1) and the rightmost . ( ϵ max = 0 75) plots of Fig. 11. . In other words, now we get ⟨ n k ⟩ ∼ | k F -k | p γ ( ) such that for a given ϵ max , the exponent p γ ( ) depends on the angle; this is in accordance with the scaling form expected in Ref. [13]. However, we note that the above simple power-law breaks down very close to the γ = π/ 4 (and equivalent points).

Concentrating on Fig. 11(c), we find that the largest residue at finite resolution, Z max ( δ ), always occur along the k x = ± k y lines with the maximum at ( √ π, √ π ) (Fig. 26) for all ϵ max (Figs. 11 and 24) and goes to zero parametrically more slowly compared to the other directions such that at finite resolution we get a smeared residue around these directions as shown in Fig. 11. One can ascribe a length to such segments of finite resolution residue as a function of both δ and ϵ max as mentioned in Fig. 27 of Appendix E 2 c. Finally turning to the density-density-correlator, ⟨⟨ n n r r ′ ⟩⟩ , the Friedel oscillations become highly direction dependent as shown in Fig. 12. This completes the phenomenology of the one and two particle static correlations of the many-Fermion wave-function obtained by superposition of ellipses.

## VI. Summary and outlook

Having described the construction of many-Fermion variational wave-functions via superposition of FS that

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-9ab0dbbf-5923-421a-9493-cf05544f202e.jpg)

FIG. 11. Plot of Z as a function of k x and k y for the superposed states (Eq. 17). The parameters used are δ = 0 01 and (a) . ϵ max = 0 10 , (b) . ϵ max = 0 50 , and (c) . ϵ max = 0 75. . The variations of Z with the resolution δ is shown in Fig. 23.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-d292eabb-30be-4e16-ba53-0144fb20dafa.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-e5a86716-220a-4873-9d3b-d3230c813f52.jpg)

FIG. 12. Plot of ⟨⟨ n n r r ′ = 0 ⟩⟩ (W( , r 0 )) as a function of r = | r | for the superposition of the elliptical FSs (Eq. 17) in different directions(tan γ = y x ) for (a) ϵ max =0.5 and (b) ϵ max =0.75.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-1a222ac5-3426-4817-b7cd-e16d217b3930.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-cf2fa4df-d684-4629-8d83-766c426d73bf.jpg)

can possibly describe correlated metals including NFLs, we now summarize our main results.

The variational wave-function of the form in Eq. 3 provide interesting insights into the two basic phenomenology of NFL metals(1) the lack of Landau quasiparticles, and (2) the well defined FS- that are obtained from quantum superposition of momentum space occupation. In view of the anisotropic loss of Landau quasiparticles on square lattice, it is interesting to re-plot Fig. 11 for the hole FS centered around k = ( π, π ). This is shown in Fig. 13 where the similarity with Fermi-arcs in underdoped cuprates [45] is apparent. However in cuprates such phenomena is observed at finite temperature while the present construction is strictly for the GS, i.e. , T = 0; in underdoped cuprates, a superconductor is obtained at lower temperatures. It has been suggested [13] that Fermi arcs can arise as a finite temperature effect of critical FS with angle-dependent scaling exponents. It would indeed be promising to investigate if similar low energy phenomenology can be obtained starting from the wave-function presented here.

a distribution [ θ ] representing a fully filled FS, | Ψ[ ] , at θ ⟩ filling ν . Now consider a related state, | Ψ[ θ ′ ] ⟩ , at the same filling obtained from [ θ ] by taking a particle at k 1 from inside the FS and putting it at an empty mode at k 2 . Then, superposing the above two results in an EPR state of a particle and a hole at the two momenta in the background of the rest of the Fermions. At low energies, both k 1 and k 2 lie within an annulus around the FS- hence local in momentum space. Thus, the GSs obtained by superposing the FS may be roughly thought as made up of local (in momentum space) particle-hole EPR pairs across the FS. It would be interesting to explore if this point provides advantage in understanding metallic phases using density-matrix-renormalization-group (DMRG) methods in momentum space [46].

Turning back to structure of the basis (Eq. 1), consider

Extending the present formulation to spinful Fermions can be obtained by defining a augmented momentumspace grid { k , σ } with σ = ↑↓ . Extension of the present formulation to the spinful case then provides for careful comparison with Gutzwiller wave-functions for correlated metals [47]. Finally, the present variational groundstates immediately lead to questions pertaining to the low

FIG. 13. Plot of Z as a function of k x and k y for superposed states (Eq. 17) for a hole FS centered around k = ( π, π ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-be17f686-448c-4008-84e3-12c3cf030994.jpg)

energy spectrum both within and beyond the RK type Hamiltonians. At this point we note that the variational wave-functions (Eq. 3 or 11), after implementing the symmetries, represent symmetric wave-functions at fractional filling, ν . Application of LSM theorem [12, 39] would suggest that this corresponds to the GS of a gapless phase as is indeed the case of the half filled Hubbard model as our calculations explicitly show. In higher dimensions, application of flux threading arguments of Ref. [7] indicate that it is possible to generate orthogonal states with momentum difference 2 πνN d -1 . For Hamiltonians which are short ranged (in real space) two dimensional generalization of LSM theorem [7, 40] would lead to a similar conclusion as in 1D about the fate of the present wavefunctions in regards to them being valid GSs of gapless phases. Addressing these issues would decide the feasibility of the present approach to understand the physics of correlated metals.

## Acknowledgments

The authors acknowledge useful discussions with F. Alet, K. Damle, R. Moessner, P. Majumdar, S. Mukerjee, M. Randeria, S. Roy, S. L. Sondhi, R. Sensarma and N. Trivedi. SB acknowledges funding by Swarna Jayanti fellowship of SERB-DST (India) Grant No. SB/SJF/202122/12; DST, Govt. of India (Nano mission), under Project No. DST/NM/TUE/QM-10/2019 (C)/7 and Max Planck Partner group Grant at ICTS. AC, HN and SB acknowledge the support of the Department of Atomic Energy, Government of India, under project no. RTI4001. KS thanks DST for support through SERB project JCB/2021/000030.

## A. The details of the [ θ ] basis

In this section, we chart out several properties of the | Ψ[ ] θ ⟩ (Eq. 1) basis. The real-space representation of the many-body state in Eq. 1, Ψ[ θ ] ( { } r ) = ⟨{ r }| Ψ[ ] , θ ⟩ can be written as

$$\left | \begin{array} { c c c c c } s _ { k _ { 1 } } \phi _ { k _ { 1 } } ( r _ { 1 } ) & s _ { k _ { 2 } } \phi _ { k _ { 2 } } ( r _ { 1 } ) & \cdots & \cdots & s _ { k _ { B } } \phi _ { k _ { B } } ( r _ { 1 } ) \\ s _ { k _ { 1 } } \phi _ { k _ { 1 } } ( r _ { 2 } ) & s _ { k _ { 2 } } \phi _ { k _ { 2 } } ( r _ { 2 } ) & \cdots & \cdots & s _ { k _ { B } } \phi _ { k _ { B } } ( r _ { 2 } ) \\ \cdots & \cdots & \cdots & \cdots & \cdots \\ \cdots & \cdots & \cdots & \cdots & \cdots \\ s _ { k _ { 1 } } \phi _ { k _ { 1 } } ( r _ { N _ { p } } ) & s _ { k _ { 2 } } \phi _ { k _ { 2 } } ( r _ { N _ { p } } ) & \cdots & \cdots & s _ { k _ { B } } \phi _ { k _ { B } } ( r _ { N _ { p } } ) \\ \end{array} \right |,$$

where s k = (1 + θ k ) / 2, { } r ≡ ( r 1 , r 2 , · · · , r N p ) is the position basis for particles corresponding to filling ν , and ϕ k ( r ) are the corresponding single-particle wavefunctions. Note that by changing of ν , the the number of rows and columns are appropriately changed by the constraint in Eq. 2.

Turning to the various many-body operators, consider the action of c † q on this state,

$$\text{ial} _ { m - } \quad c _ { q } ^ { \dagger } | \Psi [ \theta ] \rangle = c _ { q } ^ { \dagger } \prod _ { k } ( c _ { k } ^ { \dagger } ) ^ { \frac { 1 + \theta _ { k } } { 2 } } | 0 \rangle \quad \text{(A2)}$$

̸

$$\left | \text{al} \\ \text{se} \\ \text{as} \right | ^ { \text{gal} } \right | ^ { \text{el} } \left ( \prod _ { k ^ { \prime } \neq \mathbf q } \delta _ { \theta _ { k ^ { \prime } }, \theta ^ { \prime } _ { k ^ { \prime } } } \right ) \delta _ { \theta _ { \mathbf q }, - \theta ^ { \prime } _ { \mathbf q } } ( - 1 ) ^ { \sum _ { k > \mathbf q } \frac { 1 + \theta _ { k } } { 2 } } | \Psi [ \theta ^ { \prime } ] \rangle.$$

Similarly

̸

$$c _ { \mathbf q } | \Psi [ \theta ] \rangle \ = \ & \delta _ { \theta _ { \mathbf q }, 1 } \left ( \prod _ { k ^ { \prime } \neq \mathbf q } \delta _ { \theta _ { k ^ { \prime } }, \theta _ { k ^ { \prime } } ^ { \prime } } \right ) \delta _ { \theta _ { \mathbf q }, - \theta _ { \mathbf q } ^ { \prime } } \\ & \times ( - 1 ) ^ { \sum _ { k > \mathbf q } \frac { 1 + \theta _ { k } } { 2 } } | \Psi [ \theta ^ { \prime } ] \rangle. \quad \quad ( A 3 )$$

Therefore the matrix elements are given by

̸

$$\L y s i c s & \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \ \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \ Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambd \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \	\Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \LAMda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda\quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad\Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \ LAB \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \

 _ { h h f }$$

and

̸

$$\text{rjee}, \quad \text{````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````` ```````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````````
\text{d N. } \text{fel-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{\text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \\$$

These results are easily extended to matrix elements of multiple creation and annihilation operators. For two

particles, represented by corresponding Fermion creation operators,

̸

$$\text{operators}, \\ \left ( c _ { q _ { 1 } } ^ { \dagger } c _ { q _ { 2 } } ^ { \dagger } \right ) _ { [ \theta ], [ \theta ^ { \prime } ] } & = \left ( \prod _ { k ^ { \prime } \neq q _ { 1 }, q _ { 2 } } \delta _ { \theta _ { k ^ { \prime } }, \theta _ { k ^ { \prime } } ^ { \prime } } \right ) \delta _ { \theta _ { q _ { 1 } }, - \theta _ { q _ { 1 } } ^ { \prime } } \delta _ { \theta _ { q _ { 2 } }, - \theta _ { q _ { 2 } } ^ { \prime } } \\ & \quad \times \text{sign} ( q _ { 1 }, q _ { 2 } ) \left ( \prod _ { i = 1 } ^ { 2 } \frac { 1 - \theta _ { q _ { i } } } { 2 } \right ) \, \begin{array} { c } \delta _ { \theta _ { q _ { 1 } } } \\ \text{sign} \text{s} \\ \text{Fermi} \end{array} \right ) \, \text{$\quad (A6)$} \, \end{array}$$

reads

$$\text{sign} ( q _ { 1 }, q _ { 2 } ) = \begin{cases} + 1 & \forall \, q _ { 1 } > q _ { 2 } \\ 0 & \forall \, q _ { 1 } = q _ { 2 } \\ - 1 & \forall \, q _ { 1 } < q _ { 2 }. \end{cases} \quad ( A 7 )$$

This shows that the matrix elements encode the sign structure associated with the anticommutation of Fermion operators. Continuing for N -particle, we have

$$\left ( \prod _ { i = 1 } ^ { N } c _ { \mathbf q _ { i } } ^ { \dagger } \right ) = c _ { \mathbf q _ { 1 } } ^ { \dagger } c _ { \mathbf q _ { 2 } } ^ { \dagger } \dots c _ { \mathbf q _ { N - 1 } } ^ { \dagger } c _ { \mathbf q _ { N } } ^ { \dagger }, \quad ( A 8 )$$

where in the same ordering as used in Eq. 1, sign( q 1 , q 2 )

̸

$$\left ( \prod _ { i = 1 } ^ { N } c _ { \mathfrak { u } _ { \mathfrak { u } _ { 1 } } } ^ { i } \right ) = \left ( \prod _ { k ^ { \prime } \neq \mathfrak { u } _ { i } } \delta _ { \theta _ { k ^ { \prime } }, \theta _ { k ^ { \prime } } } \right ) \left ( \prod _ { i = 1 } ^ { N } \delta _ { \theta _ { \mathfrak { u } _ { 1 } } - \theta _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \leftrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak{ u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u  _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } ^ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { v } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u }_ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ {  \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \nu } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } | \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \comfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathf } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \mathfrak { u } _ { \math$$

We note that the above matrix elements are valid provided both [ θ ] and [ θ ′ ] obey the constraint relation of Eq. 2 with 0 ≤ ν ≤ 1. For particle-hole bilinears in momentum space the matrix elements are given by

̸

$$\left ( c _ { k + q } ^ { \dagger } c _ { k } \right ) _ { [ \theta ], [ \theta ^ { \prime } ] } = & \left ( \prod _ { k ^ { \prime } \neq k + q, k } \delta _ { \theta _ { k ^ { \prime } }, \theta _ { k ^ { \prime } } } \right ) \delta _ { \theta _ { k + q }, - \theta _ { k + q } } \delta _ { \theta _ { k }, - \theta _ { k } } \times \sin ( k + q, k ) \left ( \frac { 1 - \theta _ { k + q } } { 2 } \right ) \left ( \frac { 1 + \theta _ { k } } { 2 } \right ) \\ & \times \exp \left [ i \pi \sum _ { k + q > k ^ { \prime } > k } \frac { 1 + \theta _ { k ^ { \prime } } } { 2 } \right ] + \delta _ { q, 0 } \left ( \prod _ { k ^ { \prime } } \delta _ { \theta _ { k ^ { \prime } }, \theta _ { k ^ { \prime } } } \right ) \frac { 1 + \theta _ { k } } { 2 },$$

where the last term gives the diagonal contribution for the momentum space occupation n k = c † k c k .

## B. Hamiltonian in the [ θ ] basis

Having constructed the wave-function in the [ θ ] basis, we now consider the action of interacting Fermionic Hamiltonian H (Eq. 4) on these wave-functions. The matrix elements of the non-interacting part, H 0 , in this basis is given by

$$\langle \Psi [ \theta ^ { \prime } ] | H _ { 0 } | \Psi [ \theta ] \rangle \ = \ \left ( \sum _ { k } E _ { k } \left ( \frac { 1 + \theta _ { k } } { 2 } \right ) \right ) \langle \Psi [ \theta ^ { \prime } ] | \Psi [ \theta ] \rangle = E [ \theta ] \ \delta _ { [ \theta ^ { \prime } ], [ \theta ] }.$$

Thus H 0 is diagonal in the [ θ ] basis, as discussed in the main text.

Next, we consider the matrix elements for interacting part H I , given by the second term of Eq. 4. These matrix

elements can be written as

$$\mathbb { U } ^ { [ A ^ { \prime } ] } | H _ { \tau } | \mathbb { U } ^ { [ A ] } \mathbb { U } = F _ { \tau } ( [ A ^ { \prime } ] ) [ A ]$$

$$\text{ elements can be written as} \left ( \Psi [ \theta ^ { \prime } ] | H _ { 1 } | \Psi [ \theta ] \right ) & = E _ { 1 } ( [ \theta ^ { \prime } ], [ \theta ] ) \\ & = \sum _ { k _ { 2 } } \sum _ { q \neq 0 } V ( q ) \, \left ( \prod _ { k \neq k _ { 1 }, k _ { 2 }, k _ { 1 } + q, k _ { 2 } - q } \delta _ { k _ { 1 } \theta _ { k } } \right ) \left ( \delta _ { k _ { 1 } - 1 } - \theta _ { k _ { 1 } } \delta _ { k _ { 2 } - 1 } - \theta _ { k } \delta _ { k _ { 1 } + q } - \theta _ { k _ { 1 } + q } \delta _ { k _ { 2 } - q } - \theta _ { k _ { 2 } - q } } \right ) \\ & \quad \times \left ( \frac { 1 - \theta _ { k _ { 1 } } } { 2 } \right ) \left ( \frac { 1 - \theta _ { k _ { 2 } } } { 2 } \right ) \left ( \frac { 1 + \theta _ { k _ { 1 } + q } } { 2 } \right ) \left ( \frac { 1 + \theta _ { k _ { 2 } - q } } { 2 } \right ) \\ & \quad \times \sin ( k _ { 1 }, k _ { 2 } ) \, \sin ( k _ { 2 } - q, k _ { 1 } + q ) \times \exp [ i \pi \left ( \chi _ { k _ { 1 } } + \chi _ { k _ { 2 } } + \chi _ { k _ { 2 } - q } + \chi _ { k _ { 1 } + q } \right ) ] \quad ( B 2 ) \\ & = \sum _ { k _ { 3 } } \sum _ { k _ { 3 } } V ( q ) \, \left ( \prod _ { k \neq k _ { 1 }, k _ { 2 }, k _ { 1 } + q, k _ { 2 } - q } \delta _ { k _ { 1 } \theta _ { k } } \right ) \left ( \delta _ { k _ { 1 } - 1 } - \theta _ { k _ { 1 } } \delta _ { k _ { 2 } - 1 } - \theta _ { k } \delta _ { k _ { 1 } + q } - \theta _ { k _ { 1 } + q } \delta _ { k _ { 2 } - q } - \theta _ { k _ { 2 } - q } } \right ) \\ & \quad \times \left ( \frac { 1 - \theta _ { k _ { 1 } + q } } { 2 } \right ) \left ( \frac { 1 - \theta _ { k _ { 2 } } } { 2 } \right ) \left ( \frac { 1 + \theta _ { k _ { 1 } + q } } { 2 } \right ) \left ( \frac { 1 + \theta _ { k _ { 2 } - q } } { 2 } \right ) \\ & \quad \times \exp [ i \pi \left ( \chi _ { k _ { 1 } } + \tilde { \lambda } _ { k _ { 2 } } + \tilde { \lambda } _ { k _ { 2 } - q } + \tilde { \lambda } _ { k _ { 1 } + q } \right ) ], \quad ( B 3 ) \\ \text{where, } \chi _ { k _ { 1 } } = \sum _ { k > k _ { 1 } } \frac { 1 + \phi _ { k } } { 2 } \text{ for } k _ { 1 } \in \{ k _ { 1 }, k _ { 2 }, k _ { 2 } - q, k _ { 1 } + q \} \text{ and}$$

̸

̸

̸

̸

where, χ k i = ∑ k &gt; k i 1+ θ k 2 for k i ∈ { k 1 , k 2 , k 2 -q k , 1 + q } and

$$\tilde { \chi } _ { k _ { 1 } + q } = \sum _ { k _ { 1 } > k _ { 1 } + q } \frac { 1 + \theta _ { k _ { 1 } } } { 2 }, \ \tilde { \chi } _ { k _ { 2 } - q } = \sum _ { k _ { 1 } > k _ { 2 } - q } \frac { 1 + \theta _ { k _ { 1 } } } { 2 } - \frac { ( 1 - \text{sign} \left ( ( k _ { 2 } - q ) - ( k _ { 1 } + q ) \right ) ) } { 2 },$$

$$\tilde { \chi } _ { k _ { 1 } + q } & = \sum _ { k _ { 1 } > k _ { 1 } + q } \frac { 1 + \theta _ { k _ { 1 } } } { 2 }, \ \tilde { \chi } _ { k _ { 2 } - q } = \sum _ { k _ { 1 } > k _ { 2 } - q } \frac { 1 + \theta _ { k _ { 1 } } } { 2 } - \frac { ( 1 - \sin \left ( ( k _ { 2 } - q ) - ( k _ { 1 } + q ) \right ) ) } { 2 }, \\ \tilde { \chi } _ { k _ { 2 } } & = \sum _ { k _ { 1 } > k _ { 2 } } \frac { 1 + \theta _ { k _ { 1 } } } { 2 } - \frac { ( 1 - \sin \left ( k _ { 2 } - ( k _ { 2 } - q ) \right ) ) } { 2 } - \frac { ( 1 - \sin \left ( k _ { 2 } - ( k _ { 1 } + q ) \right ) ) } { 2 } - \frac { ( 1 - \sin \left ( k _ { 1 } - ( k _ { 2 } - q ) \right ) ) } { 2 } - \frac { ( 1 - \sin \left ( k _ { 1 } - ( k _ { 1 } + q ) \right ) ) } { 2 }.$$

## 1. Exact Diagonalization in Real space

The Hamiltonian given by Eq. 4 can be diagonalized by mapping it to the XXZ model using Jordan-Wigner transformation

$$\overset { 1 - } { N } \quad S _ { i } ^ { z } = \frac { 1 } { 2 } - c _ { i } ^ { \dagger } c _ { i }, \ S _ { i } ^ { + } = \frac { 1 } { 2 } \tilde { K } _ { i } c _ { i }, \ S _ { i } ^ { - } = \frac { 1 } { 2 } c _ { i } ^ { \dagger } \tilde { K } _ { i }, \quad ( C 1 )$$

(B4)

In Eq. B4, all the sign factors neutralize the effects of each other in the contribution. So we have,

$$\sum _ { i = 1 } ^ { 4 } \tilde { \chi } _ { k _ { i } } \, = \, \sum _ { i = 1 } ^ { 4 } \chi _ { k _ { i } } - \frac { ( 1 - \sin \left ( ( k _ { 2 } - \mathbf q ) - ( k _ { 1 } + \mathbf q ) \right ) ) } { 2 } - \frac { ( 1 - \sin \left ( k _ { 1 } - k _ { 2 } \right ) ) \right ) } { 2 },$$

which implies

$$( - 1 ) ^ { \sum _ { i = 1 } ^ { 4 } \tilde { x } _ { k _ { 1 } } } = & \sin \left ( ( k _ { 2 } - \mathbf q ) - ( k _ { 1 } + \mathbf q ) \right ) \quad \begin{pmatrix} \text{the} \, H \\ \text{tensor} \\ \end{pmatrix} \\ & \times \text{sign} \left ( k _ { 1 } - k _ { 2 } \right ) \right ) ( - 1 ) ^ { \sum _ { i = 1 } ^ { 4 } \tilde { x } _ { k _ { 1 } } }. \quad ( B 6 ) \quad \text{matri} \\ & \text{tonia} \quad \end{pmatrix}$$

Eq. B2 thus provides matrix elements for H I . The full Hamiltonian can now be written as in Eq. 5 which can be put into a more concise form as

$$H = \sum _ { [ \theta ], [ \theta ^ { \prime } ] } M \left ( [ \theta ], [ \theta ^ { \prime } ] \right ) | \Psi [ \theta ] \rangle \langle \Psi [ \theta ^ { \prime } ] |, \quad ( B 7 )$$

spins can be written as in Eq. C3. Now to construct the Hamiltonian matrix for the entire system, we use the tensor product (Kronecker product) to combine the Pauli matrices ( S α = σ α / 2) for different sites. Once the Hamiltonian matrix is constructed, we need to diagonalize it to find its eigenvalues and eigenvectors. We performed this process using numerical linear algebra libraries like NumPy in Python (we use sparse matrix technique to do it).

$$M ( [ \theta ], [ \theta ^ { \prime } ] ) = E [ \theta ] \ \delta _ { [ \theta ], [ \theta ^ { \prime } ] } + E _ { I } ( [ \theta ], [ \theta ^ { \prime } ] ) \left ( 1 - \delta _ { [ \theta ], [ \theta ^ { \prime } ] } \right ).$$

## C. Details of the Exact diagonalization in one spatial dimension.

To perform exact diagonalization (ED) on the XXZ spin-1/2 Hamiltonian, we begin by defining the Hamiltonian itself. The XXZ Hamiltonian for a chain of N

TABLE I. The probability density corresponding to different basis states representing particle-hole excitations over the FS as obtained through real space (XXZ) (Column 3) and the momentum-space with a fixed q c (Column 4) ED methods for the system discussed in Fig. 1 of the main text. The red numbers in column 1 indicate position of the FS in | 0 ⟩ and chart out the position of the particle-hole excitations near the FS for other states. Column 2 indicates the RMS momenta with respect to FS defined as k RMS = √ ⟨ k 2 ⟩ -√ ⟨ k 2 ⟩ FS , where ( √ ⟨ k 2 ⟩ = √∑ k k 2 (1 + θ k ) / 2), which measures deviation from | 0 ⟩ for the corresponding states.

| State                        |   k RMS |   | Ψ θ | 2 XXZ | (real space ED)   |   | Ψ θ | 2 q c ( k-space |
|------------------------------|---------|-----------------|-------------------|---------------------------|
| | 0 ⟩ = [0000 1 11111 1 000] |    0    |         0.93336 |                   |                   0.97935 |
| | 1 ⟩ = [000 10 11111 01 00] |    0.54 |         0.02284 |                   |                   0.01575 |
| | 2 ⟩ = [000 110 1111 001 0] |    0.99 |         0.00669 |                   |                   0.00124 |
| | 3 ⟩ = [00 100 1111 011 00] |    0.99 |         0.00669 |                   |                   0.00124 |
| | 4 ⟩ = [000 1110 111 0001 ] |    1.38 |         0.00296 |                   |                   8e-05   |
| | 5 ⟩ = [0 1000 111 0111 00] |    1.38 |         0.00296 |                   |                   8e-05   |

FIG. 14. (a) The probabilities of the GS wave-function obtained from ED of the half filled Hubbard model ( N = 10 , t = 1 , V 0 = 0 3) for spinless Fermions in the basis given by Eq. 1 (horizontal axis). . The relevant momenta in the first BZ are denoted by k = -π + απ/ 5 with α = 0 1 , , · · · , 9. In this labelling the four states with leading contributions are shown with [0] ≡ [0001111100] being the non-interacting GS and the rest corresponding to the excitations above the non-interacting GS. These are [1] = [0010111010], [2] = [0011011001], [3] = [0100110110], [4] = [0011010110] and [5] = [0100111001]. The states are arranged in terms of the increasing number of particle-hole excitations which can also be quantified via their root mean square momenta deviation from the FS. Note that the vertical scale has been truncated for improved visibility of the contributing basis states. (b) The exact diagonalization done using the momentum space cut-off q c ( n c = 2 ( q c = 2 π N n c ) , n 0 = 1 5 . ( q 0 = 2 π N n 0 ) , see main text for details). (c) The cumulative weights of the different basis states. (4) The IPR for the different basis states. The IPR has been normalised with respect to the maximum contribution (vertical axis truncated for better visibility) from [0].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-9db722cc-6c60-4d84-a6d3-964635c26782.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-57f84c9b-2e3a-4293-80d5-59cd1717ed7b.jpg)

where c i are spinless Fermions and

$$\tilde { K } _ { i } = \exp \left [ i \pi \sum _ { j = 1 } ^ { i - 1 } c _ { j } ^ { \dagger } c _ { j } \right ]. \quad \quad ( C 2 ) \quad \text{Sla} \\ \quad \text{Th} \\ \vdots. \.$$

The Hamiltonian (Eq. 4), using this transformation, can be written as

The GS of the XXZ Hamiltonian (Eq. C3) was obtained using ED in real space. The corresponding GS in momentum space is obtained by first computing its Slater determinant and a subsequent Fourier transform. The momentum-space GS so obtained can be understood in the [ θ ]-basis (Eq. 1) and can be used to compute the connected correlator and the entanglement entropy; these serve as a benchmark and are shown in Fig. 2.

$$H = - 2 t \sum _ { i } \left ( S _ { i } ^ { x } \ S _ { i + 1 } ^ { x } + S _ { i } ^ { y } \ S _ { i + 1 } ^ { y } \right ) + V _ { 0 } \sum _ { i } S _ { i } ^ { z } \ S _ { i + 1 } ^ { z }. \quad \text{$\text{$\text{inra}$ basis} \\ ( C 3 ) \quad \text{$\text{of the}$}$$

Since we are interested in understanding the role of interaction in the configuration space spanned by | Ψ[ ] θ ⟩ basis mentioned in Eq. 1. We looked at the amplitude of the each basis states in the computed GS as shown in

FIG. 15. (a) The probabilities of the GS wave-function obtained from ED of the half filled Hubbard model ( N = 14 , t = 1 , V 0 = 0 2) for spinless Fermions in the basis given by Eq. 1 (horizontal axis). . The relevant momenta in the first BZ are denoted by k = -π + απ/ 7 with α = 0 1 , , · · · , 13. In this labelling the four states with leading contributions are shown with [0] ≡ [00001111111000] being the non-interacting GS and the rest corresponding to the excitations above the noninteracting GS. These are [1] = [00010111110100], [2] = [00011011110010], [3] = [00100111101100], [4] = [00011101110001] and [5] = [01000111011100]. The states are arranged in terms of the increasing number of particle-hole excitations which can also be quantified via their root mean square momenta deviation from the FS. Note that the vertical scale has been truncated for improved visibility of the contributing basis states. (b) The exact diagonalization done using the momentum space cut-off q c ( n c = 3 ( q c = 2 π N n c ) , n 0 = 1 5 . ( q 0 = 2 π N n 0 ) , see main text for details). (c) The cumulative weights of the different basis states. (4) The IPR for the different basis states. The IPR has been normalised with respect to the maximum contribution (vertical axis truncated for better visibility) from [0].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-9b0041f5-d54d-415c-a6bf-81d5f527c223.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-13a98d2b-225e-4ec6-865e-ddd9ccaae2f9.jpg)

Fig. 1. The largest peak, as expected, corresponds to the state with non-interacting FS at a given ν . However, due to interaction between Fermions, there are other states with non-zero amplitudes. The details of the weights of different | Ψ[ ] θ ⟩ states with non-zero amplitude in the GS for the system parameters discussed in Fig. 1 is given in Table I. Results similar to Fig. 1 for N = 10 is shown in Fig. 14 and for N = 14 and V 0 = 0 2 in Fig. 15. .

## 2. Exact Diagonalization in momentum space

For the momentum space ED, we use a hard cutoff q c and also a soft cutoff q 0 (we choose q 0 = 2 πn /N 0 )for cutoff along k x and k y ) using ED in momentum space. This constitute a choice of the interaction potential (Eq. 4)

$$V ( { \mathbf q } ) = \begin{cases} \ V _ { 0 } e ^ { \frac { - q ^ { 2 } } { 2 q _ { 0 } ^ { 2 } } } \ \forall | { \mathbf q } | \leq | { \mathbf q } _ { c } |, \\ 0 & \text{otherwise}. \end{cases}$$

The Hamiltonian in Eq. 4 is then diagonalized by creating a matrix in the | Ψ[ ] θ ⟩ basis using the matrix elements mentioned in Eq. 5. For a given n c (we choose q c = 2 πn /N c for cutoff along k x and k y ), the number of states that are kept is given by N s = ( 2 n c C n c ) 2 in 2D; this yields a N s × N s matrix form of the Hamiltonian which is diagonalized using ED.

In 1D at ν = 1 2, there are two Fermi-points / k F1 = -π 2 ˆ x and k F2 = π 2 ˆ. x Then to calculate the different elements in Eq. B2, we note that the three momenta are restricted to the following domain with the possible values of k 1 , k 2 , q being:

$$k _ { 1 } \in \left ( k _ { F 2 } + \frac { 2 \pi } { N } \hat { x }, \ k _ { F 2 } + q _ { c } \right ), \quad \ \ ( C 5 )$$

$$k _ { 2 } \in \underset { \wedge } { \left ( k _ { F 1 } - q _ { c }, \ k _ { F 1 } - \frac { 2 \pi } { N } \hat { x } \right ) }, \quad \ \ ( C 6 )$$

$$\mathbf q \in \left ( \frac { 2 \pi } { N } \hat { x }, \ \mathbf q _ { \mathbf c } \right ),$$

$$k _ { 2 } \in \left ( k _ { F 2 } + \frac { 2 \pi } { N } \hat { x }, \ k _ { F 2 } + q _ { c } \right ), \quad \ \ ( C 8 )$$

$$k _ { 1 } \in \left ( k _ { F 1 } - q _ { c }, \ k _ { F 1 } - \frac { 2 \pi } { N } \hat { x } \right ), \quad \ \ ( C 9 )$$

$$\mathbf q \in \left ( \frac { 2 \pi } { N } \hat { x }, \ \mathbf q _ { \mathbf c } \right ).$$

The results for weights of different | Ψ[ ] θ ⟩ states in the momentum-space diagonalization are shown in Fig. 1(b), 14(b) and 15(b) for comparison.

and

FIG. 16. Difference of S vN calculated from real space ED and k-space ED (for N=22, V 0 =0.3, n 0 = 1.5). As we increase n c , it approaches zero, which means S q c vN starts to be in more agreement with S xxz vN .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-c4f1eff9-6df2-408e-a658-b8450cac061f.jpg)

## a. Details of calculation of entanglement entropy

For a 1D quantum system at a critical point, the Cardy-Calabrese [48] formula of the bipartite entanglement entropy for a system of size L with subsystem A of length l with periodic boundary conditions imposed on the whole system reads

$$S _ { v N } = \frac { c } { 3 } \log \left ( \frac { L } { \pi a } \sin \pi \frac { l } { L } \right ) + \text{constant}. \quad \text{(C11)}$$

We numerically calculate Von-Neumann Entropy

$$S _ { v N } = - T r [ \rho _ { A } \ln \rho _ { A } ], \quad \quad ( C 1 2 ) \quad \text{FIG}$$

where ρ A is the reduced density matrix of subsystem A obtained by tracing out the rest of the system ( ¯ ). A We employed standard singular value decomposition methods (SVD) to calculate the entanglement entropy, S vN , and have used Eq. C11 to obtain the central charge c as shown in Fig. 4.

For comparison, we show, in Fig. 16, the difference of the entanglement entropy calculated via real space ED and k-space ED for different n c . As n c is increased, the entanglement calculated by both methods agrees well showing the efficacy of the momentum-space method.

## D. Different Correlations for the RK state

For Eq. 10, we computed the average real space occupation ⟨ Ψ RK | ˆ n r | Ψ RK ⟩ in 1D for finite systems. This is plotted in Fig. 17. Useful insights are obtained by considering the single-particle (corresponding to ν = 1 16) / RK state (Eq. 10) ∼ ∑ k | k ⟩ ∼ | r = 0 . ⟩ The expectation value of n r to be sharply peaked at r = 0 for this state (see Fig. 17) exhibiting breaking of translation. However, this provides a way to generate Wannier

FIG. 17. Plot of ⟨ n r ⟩ as a function of r = r ˆ for x the RK ground states, where we have taken 16 sites for different fillings ν and used periodic boundary conditions (PBC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-0e8ebb81-73a7-4359-be87-674be95e4048.jpg)

FIG. 18. Plot of ⟨ n k ⟩ as a function of k for the constrained RK ground states ( | ψ k T RK ⟩ ), where we have taken 16 sites for different k T and used periodic boundary conditions (PBC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-6ee60704-8b56-46d5-ba8f-8eb735671573.jpg)

functions centered at other lattice sites by applying the lattice translation operator to the | r = 0 ⟩ state. The fully filled ( ν = 1) band, on the other hand has a flat average with the intermediate fillings interpolating between the two limits via a distribution of ⟨ Ψ RK | ˆ n r | Ψ RK ⟩ centered around r = 0. This distribution can be shifted to any lattice site via translation operator as mentioned above. The half-width of the distribution (not shown) increases monotonically with filling. It would be interesting to ask for the spectral properties of H RK which is an interesting future direction.

For the constrained RK state | ψ k T RK ⟩ , we always choose the momentum k = π to be empty for convenience in the main text. This is because on including this mode, the real space density shows oscillations with wave-vector π irrespective of the filling, ν . We think that this is an artifact of this wave-function. Turning back to | ψ k T RK ⟩ ,

FIG. 19. Plot of ⟨ n k ⟩ as a function of k for the constrained RK ground states ( | ψ k T RK ⟩ ), where we have taken 18 sites for different k T and used periodic boundary conditions (PBC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-3853bff9-4dc2-4a71-bbf9-ab39e7acfaf2.jpg)

we find

$$\langle \psi _ { R K } ^ { k _ { T } } | n _ { \mathbf r } | \psi _ { R K } ^ { k _ { T } } \rangle = \frac { 1 } { N } \sum _ { k _ { 1 }, k _ { 2 } } \langle \psi _ { R K } ^ { k _ { T } } | c _ { k _ { 1 } } ^ { \dagger } c _ { k _ { 2 } } | \psi _ { R K } ^ { k _ { T } } \rangle e ^ { i ( k _ { 1 } - k _ { 2 } ) \cdot \mathbf r }. \quad \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        } \text{                                                                                                                                                                                                    
   ( D 1 ) \quad \langle \psi _ { \nu } | n _ { \mathbf r } )$$

Since all states in the superposition carries same k T , this implies k 1 = k 2 such that the above equation gives

$$\langle \psi _ { R K } ^ { k _ { T } } | n _ { \mathbf r } | \psi _ { R K } ^ { k _ { T } } \rangle = \frac { 1 } { N } \sum _ { \mathbf k } \langle \psi _ { R K } ^ { k _ { T } } | c _ { \mathbf k } ^ { \dagger } c _ { \mathbf k } | \psi _ { R K } ^ { k _ { T } } \rangle = \nu. \quad ( D 2 )$$

The average momentum occupancy is shown in Fig. 5 in the main text. Similar plots for N = 16 and 18 are shown in Fig. 18 and 19. As mentioned in the main text, it is concluded that the oscillation in ⟨ n k ⟩ fades with increasing system size ( N ) and in the large N limit there is no sharp signature in the momentum space occupation.

## E. Superposition of Elliptical FS

## 1. The Elliptical FS

The area of BZ is 4 π 2 and for filling ν , the area of FS therefore is 4 π ν 2 . If the FS is elliptical, it's equation is given by (denoting states as | ϵ, ν, ˆ k ⟩ ),

$$\frac { k _ { \| } ^ { 2 } } { a ^ { 2 } } + \frac { k _ { \perp } ^ { 2 } } { a ^ { 2 } ( 1 - \epsilon ^ { 2 } ) } = \mu, \quad \quad \quad \ \text{(E1)}$$

where k ∥ ( k ⊥ ) are momentum components resolved along (perpendicular) to the major axis along ˆ . k For every ellipse (Fig. 6) the relation between the major axis (of length 2 a ), minor axis (of length 2 a √ 1 -ϵ 2 )for given eccentricity ϵ , such that the relation between the area of the ellipse and the filling is given by

$$\mu \, a ^ { 2 } \sqrt { 1 - \epsilon ^ { 2 } } \, = \, 4 \pi \nu, \quad \ \ \ ( \mathbf E 2 ) \quad \ s p \varepsilon$$

where a &lt; π for closed FS. For given eccentricity at half filling ( ν = 1 2), we have taken /

$$\mu = 1, \quad a = \left ( \frac { 4 \pi ^ { 2 } } { 1 - \epsilon ^ { 2 } } \right ) ^ { 1 / 4 }. \quad \quad ( E 3 )$$

Considering the states mentioned in Eq. 14, 15 and 17, we would like to emphasize that each individual states in the superposition carries same total momenta ( k T in Eq. 2). If we compute real space density ⟨ n r ⟩ , it will be uniform as shown below, and is equal to ν

$$\langle \psi _ { \nu } | n _ { \mathbf r } | & \psi _ { \nu } \rangle = \langle \psi _ { \nu } | c _ { \mathbf r } ^ { \dagger } c _ { \mathbf r } | \psi _ { \nu } \rangle \\ & = \frac { 1 } { N } \sum _ { k _ { 1 }, k _ { 2 } } \langle \psi _ { \nu } | c _ { k _ { 1 } } ^ { \dagger } c _ { k _ { 2 } } | \psi _ { \nu } \rangle e ^ { i ( k _ { 1 } - k _ { 2 } ) \cdot \mathbf r }.$$

Since all states carries same total momentum ( k T ), this implies k 1 = k 2 in eq. E4. Expanding | ψ ν ⟩ using Eq. 11, we finally obtain

$$& ( D 1 ) \quad \langle \psi _ { \nu } | n _ { \mathbf r } | \psi _ { \nu } \rangle \\ & = & \frac { 1 } { N } \sum _ { \mathbf k } \int _ { 0 } ^ { \pi } d \hat { \mathbf k } \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon | \psi ( \epsilon, \hat { \mathbf k } ) | ^ { 2 } \langle \epsilon, \nu, \hat { \mathbf k } | n _ { \mathbf k } | \epsilon, \nu, \hat { \mathbf k } \rangle \\ & = & \int _ { 0 } ^ { \pi } d \hat { \mathbf k } \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon | \psi ( \epsilon, \hat { \mathbf k } ) | ^ { 2 } \times \nu = \nu.$$

## 2. quasi-particle Residue

The quasi-particle residue has the following form [41]

$$Z _ { \mathbf K _ { \mathcal { F } } } ( \epsilon ) = n ( k _ { \mathbf F } - 0 ^ { + }, \epsilon ) - n ( k _ { \mathbf F } + 0 ^ { + }, \epsilon ), \quad ( \mathbf E 6 )$$

which captures the discontinuity at the FS. So for an ellipse that forms the basis state | ϵ, ˆ k ⟩ , the above formula can be cast into the form given by Eq. 12 using a Lorentzian regulator (which is also used in numerical calculations). For superposed ellipses, each basis state | ϵ, ˆ k ⟩ then produces its own jump with a strength of | ψ ϵ, ( ˆ ) k | 2 such that the net average residue is given by Eq. 13.

## a. The residue for the FL

The residue for the FL (Eq. 14) is obtained using Eq. 12 numerically for different resolution, δ and taking the limit δ → 0. This is shown in Fig. 20(a) while the variation of the residue at the FS with δ is shown in Fig. 21(a). Note that the residue is completely isotropic in momentum space.

FIG. 20. (a)Plot of the residue for FL state (Eq. 14) as a function of k , where ϵ 0 = 0 0, . ϵ max = 0 7, and . Z ϵ 0 , ˆ k = 0 8, (b) Plot . of Z k as a function of k for the isotropic NFL case, with different ϵ max and δ =0.001, (c) Plot of Z k as a function of k for the isotropic NFL case with the same ϵ max = 0 2 and different . δ , (d) Plot of d 2 ⟨ n k ⟩ /dk 2 corresponding to Fig. 7(b) as a function of k where ϵ max = 0 3. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-7c660ba1-46cc-4c84-903b-136ea2728bb7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-d39e1da8-1583-4fdb-a440-52ae63cd47a5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-0fdd9a6f-cebb-4602-9518-51fb3bbb5f2a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-9b2d49bb-0e01-4b6a-aa5a-3baebf05cbe4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-15e262df-2d12-4512-ad50-5b69b13b8dfa.jpg)

FIG. 21. Plot of Z k as a function of δ for (a) the FL state (Eq. 14) at k = √ 2 π , ϵ 0 = 0 0 and . ϵ max = 0 7, (b) for the Isotropic . NFL (c) for the anisotropic NFL (Eq. 17)) for different γ and ϵ max = 0 75. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-69a199d5-fac9-4d01-8c05-6ca01121e239.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-0aecf6a6-3f7c-48ae-ac17-3a14c43cb75e.jpg)

## b. Isotropic superposition of elliptic FS : Vanishing Residue

For the amplitudes in Eq. 15, the quantum state (Eq. 11) is given by

$$| \psi _ { \nu = 1 / 2 } \rangle = \frac { 1 } { \sqrt { \pi } \sqrt { \epsilon _ { \max } } } \int _ { 0 } ^ { \pi } d \hat { \mathbf k } \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon \ | \epsilon, \nu = 1 / 2, \hat { \mathbf k } \rangle, \quad \ \stackrel { \searrow } { \vdash }$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-db2a0588-3047-44f6-a027-f02d2a85c87a.jpg)

such that the expectation value of the occupancy of the momentum modes is given by Fig. 7(b). The corresponding residue calculated using Eq. 13 is then shown in Fig. 20(b) for a fixed resolution, δ = 0 001 for different . ϵ max and Fig. 20(c) for a fixed ϵ max and different resolutions.

The residue, shown in Fig 20(b) and 20(c), has its peak (at finite resolution ( δ )) at k = √ 2 π in all directions. This implies that the area of this circle is 2 π 2 which is consistent with Luttinger theorem. As shown in Fig. 21(b), Z k F → 0 as δ → 0. This is also concluded from the ⟨ n k ⟩ (in Fig 7(b)) as the discontinuity at the FS is replaced by the point of inflection. This is shown in Fig. 20(d), where the second derivative of ⟨ n k ⟩ is showing a jump at the FS( k = √ 2 π ) is plotted a s a function of k .

Finally in Fig. 22 we plot the fitting of the average momentum mode occupation to the power law form in agreement with the scaling theory of Ref. [13] as discussed in the main text.

FIG. 22. The scaling of ⟨ n k ⟩ ∼ | k F -k | p with p = 0 4828 . near the FS corresponding to Eq. 15 with ϵ max = 0 25. .

c. Superposition of ellipses on a Square Lattice : Anisotropic residue

For the state in Eq. 17, Fig. 11 shows that the superposition of FSs led to anisotropy in the quasi-particle residue (at finite resolution δ ). Here we present further details of the same as a function of the two parameters ϵ max and δ .

̸

From the plot of ⟨ n k ⟩ as a function of k , shown in Fig. 10, we note that for γ = π 4 , ⟨ n k ⟩ vanishes sharply, but smoothly. This is because in π/ 4 direction, radius of circle is max distance from origin. For ϵ = 0, it lies inside circle where it is smeared; such a smearing does not occur outside the circle.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-d39286c7-2c3a-4eab-985a-b0734aa80f53.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-73be14ec-a372-4cc9-897d-f742f3bbfe0a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-e16e5b2a-b4f0-4192-8bd1-659a056aec5c.jpg)

FIG. 23. Plot of Z as a function of k x and k y for the superposed state mentioned in Eq. 17 for ϵ max = 0 75 and (a) . δ = 0 005 . (b) δ = 0 01 . (c) δ = 0 02. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-b4d06684-1c96-4770-bed7-304c908ecf1e.jpg)

FIG. 24. Plot of Z k (for δ = 0 001) corresponding to state . Eq. 17 as a function of γ where k = ( √ 2 π cos( γ , ) √ 2 π sin( γ )) and γ ∈ { 0 , π 2 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-8531ad27-a32e-43b7-b6ae-0e70969f3af2.jpg)

FIG. 25. Plot of Z k F ( γ ) for superposed state mentioned in Eq. 17, where k F ( γ ) = ( √ 2 π cos( γ , ) √ 2 π sin( γ ))) as a function of ϵ max for different γ and δ = 10 -4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-f0167263-233d-4ad3-a8b4-0713d830a258.jpg)

The finite resolution residue (calculated using Eq. 13) is plotted in Fig. 11 for different ϵ max with a given resolution δ and as a function of δ in Fig. 23. Fig. 24 shows a plot of Z k ( δ = 0 001) as a function of . γ for different

FIG. 26. Plot of Z k ( k = ( √ π, √ π )) corresponding to Eq. 17 as a function of ϵ max for different δ .

FIG. 27. Plot of the arc length corresponding to Eq. 17 as a function of ϵ max with δ = 0 001 and for different . Z min .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00834v1-figure-9cad4e76-fc62-4e9d-86bc-7afa6e41fd05.jpg)

ϵ max indicating the direction dependence of Z k . As we move away from γ = π 4 direction, the value of the residue decreases rapidly. This is also seen by plotting Z k as a function of ϵ max for given δ (Fig. 25). At the momentum of highest smeared jump, we also plot the variation of the magnitude of jump as a function of ϵ max for different

resolution in Fig. 26. The plot shows that for any appreciable superposition, the residue indeed goes to zero. However, it goes to zero much more slowly compared to other angles as shown in Fig. 21(c).

The length of each of the segment of finite resolution residue has been computed as follows. We first normalize Z Norm = Z /Z k max . Then, while computing the length, we put a lower cut off, Z min on Z Norm and discard all Z Norm below that cutoff. Fig. 27 shows the width of segment, or the arc length, as a function of ϵ max for same δ but with different lower cutoffs Z min . We find that the maximum arc length is that of the perimeter of circle (2 π √ 2 π ).

## 3. Density-density Correlation in real space

The connected density-density correlator in real space, ⟨⟨ n r 1 n r 2 ⟩⟩ ≡ W ( r 1 , r 2 ) is given by

$$\langle \psi _ { \nu } | n _ { \mathbf r _ { 1 } } n _ { \mathbf r _ { 2 } } | \psi _ { \nu } \rangle - \langle \psi _ { \nu } | n _ { \mathbf r _ { 1 } } | \psi _ { \nu } \rangle \langle \psi _ { \nu } | n _ { \mathbf r _ { 2 } } | \psi _ { \nu } \rangle, \quad ( \mathbf E 8 ) \quad \text{$^{\pm$\dots$} } \text{or}$$

- [1] N. W. Ashcroft and N. D. Mermin, Solid state physics, (No Title) (1976).
- [2] D. Gioev and I. Klich, Entanglement entropy of fermions in any dimension and the widom conjecture, Phys. Rev. Lett. 96 , 100503 (2006).
- [3] L. D. Landau and E. M. Lifshit, Course of theoretical physics (Statistical Physics-2) (Elsevier, 2013).
- [4] R. Shankar, Renormalization-group approach to interacting fermions, Rev. Mod. Phys. 66 , 129 (1994).
- [5] J. Polchinski, Effective field theory and the fermi surface, arXiv preprint hep-th/9210046 (1992).
- [6] J. M. Luttinger and J. C. Ward, Ground-state energy of a many-fermion system. ii, Phys. Rev. 118 , 1417 (1960).
- [7] M. Oshikawa, Topological approach to luttinger's theorem and the fermi surface of a kondo lattice, Phys. Rev. Lett. 84 , 3370 (2000).
- [8] W. Ding, A. Seidel, and K. Yang, Entanglement entropy of fermi liquids via multidimensional bosonization, Phys. Rev. X 2 , 011012 (2012).
- [9] D. S. Rokhsar and S. A. Kivelson, Superconductivity and the quantum hard-core dimer gas, Phys. Rev. Lett. 61 , 2376 (1988).
- [10] F. Haldane, 'luttinger liquid theory'of one-dimensional quantum fluids. i. properties of the luttinger model and their extension to the general 1d interacting spinless fermi gas, Journal of Physics C: Solid State Physics 14 ,

for a many-Fermion state | ψ ν ⟩ at filling ν . For an elliptical FS given by | ϵ, ν, ˆ k ⟩ in the main text, we have

$$W ( { \mathbf r } _ { 1 }, { \mathbf r } _ { 2 } ) = \nu \delta ^ { 2 } ( { \mathbf r } _ { 1 } - { \mathbf r } _ { 2 } ) - | F ( \epsilon, { \mathbf r }, \nu, \hat { \mathbf k } ) | ^ { 2 }, \quad ( E 9 )$$

where,

$$F ( \epsilon, r, \nu, \hat { \mathbf k } ) = 2 \nu \left ( \frac { J _ { 1 } ( r t _ { \hat { \mathbf k } } \sqrt { 4 \pi \nu } ) } { ( r t _ { \hat { \mathbf k } } \sqrt { 4 \pi \nu } ) } \right ), \quad \quad ( \mathbb { E } 1 0 )$$

$$t _ { \hat { \mathfrak { k } } } = \frac { \sqrt { 1 - \epsilon ^ { 2 } \sin ^ { 2 } ( u - \gamma ) } } { ( 1 - \epsilon ^ { 2 } ) ^ { 1 / 4 } }, \quad \quad ( E 1 1 )$$

with

$$\tan u = \frac { \hat { k } _ { y } } { \hat { k } _ { x } }, \ r = | { \mathbf r } _ { 1 } - { \mathbf r } _ { 2 } |, \ \tan \gamma = \frac { y _ { 2 } - y _ { 1 } } { x _ { 2 } - x _ { 1 } }. \ \ ( E 1 2 )$$

Next, turning to superposed FS wave-functions (Eq. 11) the density-density correlator is given by

$$W = & \nu \delta ^ { 2 } ( { \mathbf r } _ { 1 } - { \mathbf r } _ { 2 } ) \\ + & \int _ { 0 } ^ { \epsilon _ { \max } } d \epsilon \int _ { 0 } ^ { \pi } d \hat { \mathbf k } \ | \psi ( \epsilon, \hat { \mathbf k } ) | ^ { 2 } | F ( \epsilon, r, \nu, \hat { \mathbf k } ) | ^ { 2 }.$$

This completes our discussion of computation of these correlators.

2585 (1981).

- [11] T. Giamarchi, Quantum Physics in One Dimension (Oxford University Press, 2003).
- [12] M. Yamanaka, M. Oshikawa, and I. Affleck, Nonperturbative approach to luttinger's theorem in one dimension, Phys. Rev. Lett. 79 , 1110 (1997).
- [13] T. Senthil, Critical fermi surfaces and non-fermi liquid metals, Phys. Rev. B 78 , 035103 (2008).
- [14] W. F. Brinkman and T. M. Rice, Application of gutzwiller's variational method to the metal-insulator transition, Phys. Rev. B 2 , 4302 (1970).
- [15] Q. Si, Quantum criticality and global phase diagram of magnetic heavy fermions, physica status solidi (b) 247 , 476 (2010).
- [16] F. Kr¨ger and J. Zaanen, Fermionic quantum criticality u and the fractal nodal surface, Phys. Rev. B 78 , 035104 (2008).
- [17] J. Zaanen, Quantum critical electron systems: the uncharted sign worlds, Science 319 , 1205 (2008).
- [18] H. v. L¨ ohneysen, A. Rosch, M. Vojta, and P. W¨ olfle, Fermi-liquid instabilities at magnetic quantum phase transitions, Rev. Mod. Phys. 79 , 1015 (2007).
- [19] H. v. L¨hneysen, T. Pietrus, G. Portisch, H. G. Schlager, o A. Schr¨der, o M. Sieck, and T. Trappmann, Non-fermiliquid behavior in a heavy-fermion alloy at a magnetic instability, Phys. Rev. Lett. 72 , 3262 (1994).

and

- [20] P. Gegenwart, Q. Si, and F. Steglich, Quantum criticality in heavy-fermion metals, nature physics 4 , 186 (2008).
- [21] J. Custers, P. Gegenwart, H. Wilhelm, K. Neumaier, Y. Tokiwa, O. Trovarelli, C. Geibel, F. Steglich, C. P´pin, e and P. Coleman, The break-up of heavy electrons at a quantum critical point, Nature 424 , 524 (2003).
- [22] P. Coleman and A. J. Schofield, Quantum criticality, Nature 433 , 226 (2005).
- [23] Y. Kurosaki, Y. Shimizu, K. Miyagawa, K. Kanoda, and G. Saito, Mott transition from a spin liquid to a fermi liquid in the spin-frustrated organic conductor κ -(ET) cu (CN) , Phys. Rev. Lett. 2 2 3 95 , 177001 (2005).
- [24] M. Plat´ e, J. D. F. Mottershead, I. S. Elfimov, D. C. Peets, R. Liang, D. A. Bonn, W. N. Hardy, S. Chiuzbaian, M. Falub, M. Shi, L. Patthey, and A. Damascelli, Fermi surface and quasiparticle excitations of overdoped tl 2 ba cuo 2 6+ δ , Phys. Rev. Lett. 95 , 077001 (2005).
- [25] S.-S. Lee, Recent developments in non-fermi liquid theory, Annual Review of Condensed Matter Physics 9 , 227 (2018).
- [26] J. Polchinski, Low-energy dynamics of the spinon-gauge system, Nuclear Physics B 422 , 617 (1994).
- [27] S.-S. Lee, Low-energy effective theory of fermi surface coupled with u(1) gauge field in 2 + 1 dimensions, Phys. Rev. B 80 , 165102 (2009).
- [28] J. A. Hertz, Quantum critical phenomena, in Basic Notions of Condensed Matter Physics (CRC Press, 2018) pp. 525-544.
- [29] A. J. Millis, Fluctuation-driven first-order behavior near the t = 0 two-dimensional stripe to fermi liquid transition, Phys. Rev. B 81 , 035117 (2010).
- [30] G. R. Stewart, Non-fermi-liquid behavior in d -and f -electron metals, Rev. Mod. Phys. 73 , 797 (2001).
- [31] S. Sachdev and J. Ye, Gapless spin-fluid ground state in a random quantum heisenberg magnet, Phys. Rev. Lett. 70 , 3339 (1993).
- [32] M. Y. Reizer, Relativistic effects in the electron density of states, specific heat, and the electron spectrum of normal metals, Phys. Rev. B 40 , 11571 (1989).
- [33] B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Low-energy properties of fermions with singular interactions, Phys. Rev. B 50 , 14048 (1994).
- [34] P. A. Lee, Gauge field, aharonov-bohm flux, and highT c
- superconductivity, Phys. Rev. Lett. 63 , 680 (1989).
- [35] M. A. Metlitski and S. Sachdev, Quantum phase transitions of metals in two spatial dimensions. ii. spin density wave order, Phys. Rev. B 82 , 075128 (2010).
- [36] D. F. Mross, J. McGreevy, H. Liu, and T. Senthil, Controlled expansion for certain non-fermi-liquid metals, Phys. Rev. B 82 , 045121 (2010).
- [37] D. Dalidovich and S.-S. Lee, Perturbative non-fermi liquids from dimensional regularization, Phys. Rev. B 88 , 245106 (2013).
- [38] S. Sur and S.-S. Lee, Chiral non-fermi liquids, Phys. Rev. B 90 , 045121 (2014).
- [39] E. Lieb, T. Schultz, and D. Mattis, Two soluble models of an antiferromagnetic chain, Annals of Physics 16 , 407 (1961).
- [40] M. B. Hastings, Lieb-schultz-mattis in higher dimensions, Phys. Rev. B 69 , 104431 (2004).
- [41] L. Landau, E. Lifshitz, and L. Pitaevskii, Course of theoretical physics: Statistical physics, part 2 : by e.m. lifshitz and l.p. pitaevskii, (No Title) (1980).
- [42] P. W. Anderson, Singular forward scattering in the 2d hubbard model and a renormalized bethe ansatz ground state, Phys. Rev. Lett. 65 , 2306 (1990).
- [43] C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abrahams, and A. E. Ruckenstein, Phenomenology of the normal state of cu-o high-temperature superconductors, Phys. Rev. Lett. 63 , 1996 (1989).
- [44] S. A. Kivelson, E. Fradkin, and V. J. Emery, Electronic liquid-crystal phases of a doped mott insulator, Nature 393 , 550 (1998).
- [45] M. R. Norman, H. Ding, M. Randeria, J. C. Campuzano, T. Yokoya, T. Takeuchi, T. Takahashi, T. Mochiku, K. Kadowaki, P. Guptasarma, et al. , Destruction of the fermi surface in underdoped high-t c superconductors, Nature 392 , 157 (1998).
- [46] T. Xiang, Density-matrix renormalization-group method in momentum space, Phys. Rev. B 53 , R10445 (1996).
- [47] M. C. Gutzwiller, Correlation of electrons in a narrow s band, Phys. Rev. 137 , A1726 (1965).
- [48] P. Calabrese and J. Cardy, Entanglement entropy and conformal field theory, Journal of Physics A: Mathematical and Theoretical 42 , 504005 (2009).
Ankush Chaubey Harsh Nigam Subhro Bhattacharjee K. Sengupta
Variational wave-functions for correlated metals
We study a set of many-body wave-functions of Fermions that are naturally
written using momentum space basis and allow for quantum superposition of
Fermion occupancy, $\{n_{\bf k}\}$. This {enables} us to capture the
fluctuations of the Fermi-surface {(FS)} -- the singularly most important
signature of a metal. We bench-mark our results in one spatial dimensions (1D)
to show that these wave-functions allow for quantitative understanding of the
Tomonaga-Luttinger liquid (TLL); computations of certain correlators using them
can in fact be extended to larger systems sizes compared to conventional exact
diagonalization (ED) allowing for a more systematic comparison with
bosonization techniques. Finally we show that this basis may be useful for
obtaining fixed-point wave-function for strongly correlated metals {in
dimensions greater that one}. In particular, we study the case of coherent
(equal) superposition of elliptical FS {in continuum (2D) and on a} square
lattice{. In case of the former, our variational wave-function systematically
interpolates between the phenomenology of the Fermi liquid ground state, i.e.,
finite single-Fermion residue at a sharp FS, to a non-Fermi liquid (NFL) with
zero residue. In the NFL the jump in $\langle n_{\bf k}\rangle$ at the FS is
replaced by a point of inflection (similar to a 1D TLL) whose contour is
consistent with the Luttinger Theorem. In case of the square lattice, we} find
highly anisotropic distribution of the quasi-particle residue, which, at finite
resolution has an uncanny resemblance to the Fermi-arcs{, albeit at zero
temperature,} seen in the pseudo-gap state of the cuprates.