Machine learning reveals features of spinon Fermi surface

With rapid progress in simulation of strongly interacting quantum Hamiltonians, the challenge in characterizing unknown phases becomes a bottleneck for scientific progress. We demonstrate that a Quantum-Classical hybrid approach (QuCl) of mining sampled projective snapshots with interpretable classical machine learning can unveil signatures of seemingly featureless quantum states. The Kitaev-Heisenberg model on a honeycomb lattice under external magnetic field presents an ideal system to test QuCl, where simulations have found an intermediate gapless phase (IGP) sandwiched between known phases, launching a debate over its elusive nature. We use the correlator convolutional neural network, trained on labeled projective snapshots, in conjunction with regularization path analysis to identify signatures of phases. We show that QuCl reproduces known features of established phases. Significantly, we also identify a signature of the IGP in the spin channel perpendicular to the field direction, which we interpret as a signature of Friedel oscillations of gapless spinons forming a Fermi surface. Our predictions can guide future experimental searches for spin liquids.


Introduction.
As our ability to simulate quantum systems increases, there is a corresponding need for determining how to characterize unknown phases realized in simulators.Going from measurements to the nature of the underlying state is a challenging inverse problem.Full quantum state tomography [1] of the density matrix is impractical.Although the classical shadow [2] scales better than full tomography, the approach does not prescribe to researchers the proper observables to evaluate.Viewing the inverse problem as a data problem invites adopting machine learning methods: a quantum-classical hybrid approach.Machine learning has been widely applied for characterizing quantum states [3].Such methods have been most fruitful with symmetry-broken states, with a diverse set of approaches increasingly bringing more interpretability and reducing bias [4][5][6].The characteristic features of ordered phases are ultimately local and classical, hence ML models tuned for image processing have readily learned such features.By contrast, past learning of quantum states defined without order parameters has relied on theoretically guided feature preparation [7,8].However, such reliance on prior knowledge blocks the researchers' access to new insights into unknown states: the ultimate goal of simulating quantum states.
To push the limits of the nascent quantum-classical hybrid approach, we need a setting known to host a non-trivial quantum phase of unknown nature.Recent investigations into extended Kitaev models [9][10][11][12] have led to the observation of a mysterious intermediate gapless * Corresponding author: kz345@cornell.eduphase (IGP) sandwiched between the Kitaev spin liquid and the trivial polarized state under a non-perturbative [111] magnetic field [13][14][15][16], whose identification presents an interesting and important puzzle away from the perturbative limit.However, the nature of this field-induced IGP has raised debate in the community.
Several theories have shown evidence that supports a gapless quantum spin liquid phase with an emergent U (1) spinon Fermi surface [15,[17][18][19], while there are also mean field theories indicating that the low energy effective theory of the intermediate phase is gapped with a non-zero Chern number [20,21].This tension between theories arises due to the challenge in determining the nature of the IGP that forms in a non-perturbative region under a magnetic field.Unlike the gapped topological phase adiabatically connected to the exactly solvable limit with known loop operators [7,22,23], the absence of measurable positive features for the possible candidate IGP states [18,19] also makes this problem a worthy challenge for machine learning.
We present a quantum-classical hybrid approach, QuCl, to reveal characteristic motifs associated with states without known signature features.We treat variational wavefunctions obtained from density matrix renormalization group (DMRG) [24][25][26] as output of a quantum simulator.Namely, we sample snapshots from the ground state and train an interpretable neural network architecture, i.e. the correlator convolutional neural network (CCNN) [4] (Figure 1(c)).Based on the trained network, we use regularization path analysis [27] to determine the distinct correlation functions learned by the CCNN as characteristic features of the state captured by snapshots.We benchmark the performance of this hybrid approach on the known phases and confirm that the CCNN learned features are consistent with the known characteristic features.Importantly, we reveal the signature feature of the IGP to imply the existence of a spinon Fermi surface, as proposed in [18,19].
Results and Discussion.Model.The Kitaev-Heisenberg model under an external field is defined by where γ = x, y, z enumerates the three colorings of bonds on the honeycomb lattice (Figure 1(a)), and S γ is the γ projection of spin-1/2 degrees of freedom on each site.We also add a uniform Zeeman field h along the [111] direction, i.e., the out-of-plane ê3 direction in the lab frame (see also Supplementary Note III), as well as a ferromagnetic Heisenberg term of strength J.Here we investigate the antiferromagnetic Kitaev intereaction (K > 0), for which a field in the [111] direction gives rise to an intermediate phase over a significant field regime.For the ferromagnetic Kitaev interaction, on the other hand, the intermediate phase is either absent or exists over a very small field regime [18].Starting from the exactly solvable point at [28], we consider three axes of the phase diagram that are controlled by the parameters h, J, and For the J axis of the phase diagram in Figure 1(b), the system undergoes a sequence of transitions through magnetically ordered states [11].For small values of J the system preserves time-reversal symmetry and the system remains a gapless Z 2 spin liquid.As J is increased, the system acquires a zigzag magnetic order (also experimentally observed in α-RuCl 3 [11,29,30]).At even larger values of J, the system eventually becomes a Heisenberg ferromagnet.On the other hand, a small magnetic field h ∥ [111] breaks the time reversal symmetry of the Kitaev model and opens a gap in the spectrum of free majorana fermions, resulting in a CSL [28].However, upon leaving the perturbative regime, numerical evidence through DMRG [19] and exact diagonalization [16,18] have shown that the system goes through an IGP before entering a partially polarized (PP) magnetic phase.Although the precise nature of the IGP is unknown, a U (1) spinon Fermi surface has been proposed recently [18,19], which, as we are to show, is in agreement with our CCNN results.Finally, the Kitaev model has an exact solution along the K z axis when J = h = 0 [28], where the system undergoes a transition to a gapped Z 2 spin liquid upon increasing K z .We use this axis for benchmarking the QuCl outcome to known exact results.
In order to generate a single snapshot from a wavefunction, we perform the following procedure sequentially on each site i of the lattice: 1. Find the reduced density matrix for site i, and exactly evaluate the expectation value of the spin operator projected along the chosen axis α i .
3. Collapse the wavefunction onto the associated eigenstate of site i using the projector |±α⟩ i ⟨±α| i .
5. Organize the snapshots into channels, one for each unique axis α i ; see Figure 1(c).
The choice of axis α i is random for the J and K z axes but tailored to the target phases for the h axis.The wavefunctions at phase space points of interest are obtained using DMRG on finite size systems composed of 6 × 5unit-cell (60 sites).In Supplementary Note III we also show results from an extended 12 × 3-unit-cell (120 sites) cylinder geometry.In both cases we used a maximum of 1200 states, giving converged results with a truncation error ∼ 10 −7 or less in all phases.Within a phase, we generated 10,000 snapshots for each wavefunction in question.Each resulting snapshot forms a three-dimensional array of bit-strings, with two spatial dimensions and a "channel" dimension (see Figure 1(c)).Such a collection of snapshots is a classical shadow of the quantum state [2].Since our goal is to characterize a quantum state without prior knowledge of the best operator to measure, we treat the snapshot collection as data rather than using them to estimate an operator expectation value as in Refs.[2,31].
For each axis of the phase diagram, we set up a binary classification problem between a pair of phase space points, |Ψ 0 ⟩ and |Ψ 1 ⟩, each deep within a phase.The machine learning architecture of choice, CCNN, was introduced in Ref. [4] as an adaptation of a convolutional neural network where a controlled polynomial non-linearity splits into different orders of correlators for the neural network to use (Figure 1(c)).Compared to the more standard CNN architecture, the CCNN has reduced expressibility due to using a low-order polynomial as the nonlinearity.However, at the expense of this reduction, we gain access to interpreting the network's learning that can be analytically connected to the traditional notion of correlation functions.Combined with regularization path analysis (RPA; see Methods and Supplementary Note I for details) [32], the CCNN can reveal spatial correlations or motifs that are characteristic of a given phase.
For a given channel α of filter k to be learned, f k,α , the CCNN samples correlators for each snapshot bit-string B α (x) through an estimate for the n-th order spatially averaged correlator associated with filter k c where the inner sum is over all n unique pairs of filter positions a and filter channels α.These correlator estimates are then coupled to coefficients β , where 0 ≤ ŷ ≤ 1 is the CCNN output for the given input snapshot.We reserved 1,000 samples from each wavefunction as a validation set, and used the remaining 9,000 for training.The orders of correlators were restricted to be between 2 and 6, inclusive.We allowed the neural network to learn up to 4 different filters, corresponding to 0 ≤ k ≤ 3.
The training optimizes the model parameters, namely the filters and the weights, by comparing the output ŷ to the training label (see Methods).
Once the CCNN is successfully trained for a given phase, we uncover the characteristic motif that is most informative for the contrast using RPA [32].For this, we fix the filters and relearn the weights of each learned correlation β

The β (n) k
that turns on at the lowest value of 1/λ points to the specific filter k and the correlation order (n) of that filter which is most informative for the contrast task.The sign of the onsetting β (n) k reveals whether the associated correlation is a feature of phase 0 (− sign) or of phase 1 (+ sign); see Figure 1(c).
Gapless Z 2 versus Heisenberg phases.As a benchmark, we first focus on the phases along the J axis (Figure 2(a)).At intermediate J, the system has zigzag order, while it is a Heisenberg ferromagnet for large J.We trained the CCNN to distinguish wavefunctions from the two points marked by stars in Figure 2(a), at J/K = 0 and K/J = 0, corresponding to the Kitaev spin liquid and Heisenberg ferromagnetic states, respectively.The snapshots were generated by choosing a random axis from x, y, or z for each site.The RPA shown in Figure 2(b) reveals that the most informative correlation functions are the twopoint functions of filters 2 and 3 presented in Figure 2(a).The negative sign of the onsetting β's means these features are positive indicators of the ordered phase (see Supplementary Note I).Given that the correlation length vanishes at the exactly solvable point at the origin (phase 0), the network's choice to focus on features of phase 1 is sensible.Moreover, the learned motif of phase 1 is clearly consistent with a ferromagnetic correlation.Hence this benchmarking confirms that the CCNN's learning is consistent with our theoretical understanding when both phases 0 and 1 are known.
Chiral spin liquid.Next, we contrast the CSL phase (phase 1) and the IGP (phase 0) along the h axis (Figure 3(a)).Neither of these phases is characterized by a local order parameter.However, the chiral phase is known to be a Z 2 quantum spin liquid characterized by non-local Wilson loop expectation values [28].To confirm that such non-local information can be learned with our architecture, we first use snapshots with a fixed basis shown in Figure 3(c) so that the architecture can access the necessary information.The RPA with the positive onset of β (6) 0 (Figure 3(d)) implies that a sixth-order correlator of the filters shown in Figure 3(a) is learned to be the key indicator of phase (1), the CSL phase.Remarkably, the relevant correlator ⟨σ z 1 σ x 2 σ y 3 σ z 4 σ x 5 σ y 6 ⟩ is exactly the expectation value of the Wilson loop associated with the plaquette p consisting of the six sites ⟨W p ⟩, shown in Figure 3(b).Theoretically, ⟨W p ⟩ ≈ 1 implies the state is well-described by the Z 2 gauge theory of the zero-field gapless phase [28].The fact that none other than ⟨W p ⟩ was learned to contrast the CSL phase from the intermediate gapless phase reveals that the latter is a distinct state.However, discovering the indicator of the intermediate phase requires a different approach, as we discuss below.
Intermediate Gapless Phase.We next discuss how we discover the physically meaningful features of the IGP (phase 0).Previous work has focused on mapping out the low energy excitations S(q, ω ≈ 0) in momentum space Also, in real space, the spin-spin correlations averaged over all directions shows power-law decay, indicating gapless spin excitations for intermediate fields.However, it has not been clear how to translate these correlations to positive signatures of a particular state that can be experimentally detected.While the QuCl approach has the potential to reveal such signatures, we have to first overcome a ubiquitous challenge accompanying using ML for scientific discoveries: the need to guide the machine away from trivial features.While the unbiased pursuit of representative feature in data is the benefit of using ML, a non-trivial cost is that the neural network's learning can be dominated by features that are trivial from the physicist's perspective.The neural network's propensity to make decisions based on what appears most visible to the network means it is essential that we guide the CCNN away from the trivial yet dominant difference between phase 0 and phase 1: the field-driven magnetization along the e 3 axis (see Supplementary Note II).This basic requirement for extracting meaningful information using ML led us to supply CCNN with snapshots in the basis orthogonal to the field direction, such as e 1 basis (see Figure 4(a)).This decision to guide the CCNN away from trivial features led to a sought-after discovery.
The RPA shown in Figure 4(b) unambiguously points to two-point correlators of the filter shown in Figure 4(c) as a signature feature of the IGP.As is clear from its Fourier transform shown in Figure 4(d), the filter implies the emergence of a length scale in the e 1 component of the magnetization.Given that the e 1 direction is perpendicular to the direction of h-field, the repeating arrangement of the motif the filter is detecting must be anti-ferromagnetic.One such ansatz we conjecture shown in Figure 4(e) will single out specific momenta points marked in Figure 4(d) from the Fourier intensity of the filter (see Supplementary Note IV for more details).To confirm this conjecture, we explicitly measure the per-site e 1 -magnetization, ⟨S e1 (r)⟩, of the two states.The measurement outcome (Figure 4(f)) and its Fourier transform (Figure 4(g)) confirms indeed the IGP state has a modulating e 1 -magnetization that we inferred from the CCNN learned filter and the ansatz tiling the filter.Furthermore, the contrast between Fourier transforms from the IGP Figure 4(g) and from the CSL Figure 4(h) establishes that the pattern and the associated length scale are unique features of the IGP.Remarkably, we find such modulation to be consistent with a conjecture [15-19, 22, 33] that the IGP is a U (1) spin liquid with a spinon Fermi surface.Note here that in a translationally invariant system the corresponding quantity is the twopoint spin-spin correlation function ⟨S e1 (0)S e1 (r)⟩ and its Fourier transform.
As detailed in Supplementary Note III, ⟨S e1 (r)⟩ can be mapped to fermionic spinon density in the Kitaev model.If spinons are gapless and deconfined to form a spinon Fermi surface, the Friedel oscillation of the spinon density due to the open boundary [34][35][36][37] will be reflected in the modulation of ⟨S e1 (r)⟩: where C and θ are constants, and r the distance measured from the boundary.We confirm the spinon Friedel oscillation origin of the observed modulations by fitting the ⟨S e1 (r)⟩ measured at different field strengths h to Eq. ( 3).The resulting excellent fit in Figure 4(i) shows that the modulation period increases with the increase in the perpendicular field h.This is consistent with a mean field picture in which the magnetic field plays the role of the chemical potential; the spinon bands successively get depleted upon increasing field until the system enters a trivial phase through a Liftshitz transition [38].Evaluation of S e1 on a longer 20 unit cell system of a 3-leg ladder shows a modulation pattern that agrees with the 6-leg ladder results in Figure 4(i) (see Supplementary Note III).
Gapless Z 2 vs gapped Z 2 .Finally, we contrast the gapless and gapped Z 2 phases along the K z axis as a sanity check in distinguishing two spin liquid phases.As one tunes K z , the model Eq. ( 1) is known to go through a phase transition between a gapless Z 2 and a gapped Z 2 spin liquid phases [28].However, since both phases have only short-range correlations in the ground state the distinction cannot be learned from the correlation lengths, unlike usual transitions between a gapless phase and a gapped phase.Hence it is a non-trivial benchmarking test for QuCl-based state characterization.Contrasting the two points marked by stars in Figure 5(a), again using the random basis snapshots, we find signature motifs consistent with exact solutions.Specifically, the RPA (Figure 5(b)) shows that nearest-neighbor correlation functions of x and y axes are a feature of the gapped Z 2 phase while the z axis nearest-neighbor correlation function is the feature of the gapped Z 2 phase.These results are consistent with the exact solution of the zero-field Kitaev model [39,40].
Conclusion.The significance of our findings is threefold.Firstly, we gained insight into the intermediate field spin liquid phase in the Kitaev-Heisenberg model.Confronted by two complementary predictions: a gapless spin liquid based on exact diagonalization and DMRG versus a gapped spin liquid in the same region from mean field theory, an identification of a positive signature for either possibilities was critical.The need for guiding the CCNN away from a trivially changing feature led to the discovery that it is critical to focus on snapshots taken along a direction e 1 perpendicular to the magnetic-field e 3 axis.Remarkably, the network then learned a geometric pattern characteristic of Friedel oscillations of spinons in the IGP.This observation strongly supports earlier theoretical proposals of a spinon Fermi surface in the IGP, thus advancing our understanding of this phase.
Secondly, our discovery translates to a prediction for experiments by providing a direct evidence of spinon FS in the modulated magnetization and the spin-spin correlations perpendicular to the field direction along e 1 .Such a feature in the computational data has been previously missed since the focus has been on isotropic spin correlation ⟨ ⃗ S i • ⃗ S j ⟩ which is dominated by the e 3 component.Our results can guide future experimental searches for spin liquids with spinon Fermi surfaces.
Finally, on a broader level, we have demonstrated that hidden features of a quantum many-body state can be discovered using QuCl: a data-centric approach to snapshots of the quantum states, employing an interpretable classical machine learning approach.Conventionally, quantum states have been studied through explicit and costly eval- uation of correlation functions.However, when the descriptive correlation function is unknown in a new phase, the conventional approach gets lost in the overwhelming space of expensive calculations.Although our method does not explicitly evaluate the correlation functions that it extracts, snapshots that can be readily treated with QuCl will enable computationally efficient identification of new phases associated with a quantum state, including topological states or states with hidden orders.Finally, our method is also broadly applicable to searches for physical indicators of states prepared on quantum simulators which are naturally accessed through projective measurements.

Methods.
In this section, we describe the architecture of the neural network and the training procedure.The CCNN, as first proposed in Ref. [4], consists of two layers: the correlator convolutional layer and the fully connected linear layer.We fed as input to the CCNN three-spinchannel (two-spin-channel for rotated basis measurements) snapshot data.Since the CCNN was originally applied to square lattice data at its conception, we reinterpreted our hexagonal lattice geometry as a rectangular grid with a 1 × 2 unit cell forming its two-site basis.We modified the convolutional layer to consist of 4 different learnable filters of dimension 2 × 2 unit cells, for a total receptive field of 8 sites each.To accommodate the 1 × 2 unit cell, we also introduced a horizontal stride of 2 in the convolution operation between filters and snapshots.
The filter weights are learnable nonnegative numbers indicated by f α,k (a), where 1 ≤ α ≤ 4 indexes the filter, 1 ≤ k ≤ 3 indexes the channel of the weight, and a is a spatial coordinate.The weights are convolved with the input snapshots using the recursive algorithm described in Ref. [4] to produce per-snapshot correlators as (n) k (x) is the position-dependent n-th order correlator of filter k, and B αj (x + a j ) indicates the snapshot value at location x + a j in channel α j .The correlator estimates are then defined as the spatially-averaged correlators, c of the linear layer, and summed to produce the logistic regression classification output ŷ = 1 so that they are constrained to the range 0 ≤ ŷ ≤ 1.For a visual overview of the architecture, see Figure 6.
During training, the weights of the network are updated with stochastic gradient descent to optimize the loss function where y ∈ {0, 1} is the ground truth label of the snapshot, and γ 1 and γ 2 are L1 and L2 regularization strengths, respectively.We took γ 1 = 0.005 and γ 2 = 0.002.The training was performed for 20 epochs consisting of 9,000 snapshots each with a learning rate of 0.006, using Adam stochastic gradient descent.For the regularization path analysis, the weights f are kept fixed, and the model is retrained in the same way with loss function where γ is the regularization strength to be swept over.Data Availability.Data is available upon request to the authors.4), then spatially averaged to obtain correlator estimates c Correlator estimates are then used in a logistic regression to predict the class of the input snapshot.
indicates the specific phase that the correlation functions are associated with: in this case, the signs being negative indicates that the filters point towards the ferromagnetic phase, which is consistent with the identification of the filters with the q = 0 spin-spin correlation function.

SUPPLEMENTARY NOTE II. TRAINING IN ROTATED BASIS Parallel axis
As a test case, we also trained the CCNN to distinguish the chiral spin liquid from the intermediate gapless phase using snapshots taken in the e 3 axis, which is parallel to the magnetic field.Since the magnetization of the gapless phase is larger, the CCNN should be able to use this information to trivially distinguish the two phases.Indeed, the two-point correlators of Filter 3 in Figure S2(a) were found to trivially characterize the gapless phase, as stated in the main text.Since the CCNN estimates unconnected correlation functions, the larger ⟨S e3 ⟩ expectation value would cause all two-point functions of S e3 variables to be larger in the IGP, which is consistent with the learned correlation function.

Perpendicular axis
Here, we present details of CCNN training using the e 1 basis snapshots Figure S3.Since the snapshots were only taken in the e 1 basis, the e 2 and e 3 channels of the filters are left unused.Regularization path analysis (Figure S3(b)) shows the most significant correlation function characterizing the IGP being the second-order correlators of Filter 0, as analyzed in the main text.
We note that the magnetization of the IGP (Figure S3(c)) is highly contrasted with the magnetization of the CSL phase, since the CSL magnetization is largely featureless except for edge effects.On the other hand, the IGP magnetization is rich, and we provide an explanation of how Filter 0 predicts the observed pattern both in the main text, and also in SM VI. which gives the new spin axis according to: where ⃗ R ∈ SO(3) ⊂ SU (2).Hence we can rewrite this in terms of S e1 , S e2 , S e3 .After a spin rotation, they are related to the original operators by It is readily seen that these new operators satisfies SU(2) Lie algebra: Spin operators are decomposed into four majorana particles in Kitaev's exact solution [1].In our rotated basis we label these four majoranas as c, b e1 , b e2 , b e3 .Now, let us write the majorana decomposition in the following way.Define two fermion modes ψ 1 , ψ 2 for each site, where such that each spin is now where by the onsite gauge constraint This fermionization is equivalent to both Kitaev's four-majorana representation as well as slave fermion representation [2].To appreciate this point, recall that the Abrikosov decomposition fermionizes a spin-1 2 into spinons according to with the onsite constraint Then if we identify we would find that Eq. 13 -15 are close to, but actually not identical to Eq. 9, Eq. 10 and Eq.11.Indeed, it is straightforward to find their relation: However, after a rotation and the application of constraint in Eq. 16 i.e. single particle per site, we immediately see their equivalence S ej i ≃ Sej i .Furthermore, the slave fermion constraint in Eq. 16 and the gauge constraint in Eq. 12, are also equivalent, which can be made explicit by noticing Therefore, the gauge transformed four-majorana representation used in the main text is equivalent to the Abrikosov spinons for spin-1 2 , which is responsible for the Friedel oscillation in the emergent gapless phase.In this picture, the external magnetic field can be interpreted as a chemical potential for the spinons, which are depleted with increasing field as sketched in Figure S5.This mapping is always faithful as long as the on-site constraint is met.In this representation, the Kitaev spin liquid can also be solved in a mean field picture whereby spinons form a weakpairing p-wave superconductor [2]; the gapless intermediate phase is then effectively a neutral metal via a spinon superconductor to spinon metal transition at h c1 , which is naturally accompanied by a U(1) gauge field and Fermi surface; and the trivial polarized phase is a trivial gapped insulator leading to confinement of spinons.Among all three aforementioned phases, only in the U(1) gapless phase does the fermion density (which is equivalent to on-site magnetization amenable to DMRG) oscillate with a length scale π/k F due to the emergence of a gapless Fermi surface [3,4], whereas the CSL and polarized phase do not due to the absence of it.The real and momentum space modulation of the local ⟨S e1 ⟩ is shown in Figure S4 for different lattice geometry and magnetic fields.In the intermediate gapless phase, we note that, the larger lattice with 3×20 unit cells in the lower panel of Figure S4 gives similar pattern to that of the smaller one (top panel and in the main text).For the h c1 < h = 0.28 < h c2 wavefunction, the location of the peak is approximately the same as that measured in the main text, i.e. q = (−1.5,1/ √ 3)π/a.For the wavefunction at larger h c1 < h = 0.32 < h c2 , the peak is shifted to larger k due to the increasing magnetic field (and thus smaller k when folded into the first Brillouin zone), which is consistent with the spinon picture described in the main text.

SUPPLEMENTARY NOTE IV. ANALYSIS OF MAGNETIZATION BRAGG PEAKS IN U (1) PHASE
Here we show how the learned filter characterizing the gapless phase predicts the amplitudes of the Bragg peaks in the ⟨S e1 (x)⟩ magnetization of the gapless phase.First, we assume that the magnetization pattern can be formed by a repeated, tiled, filter motif.Since the applied field is in a perpendicular direction, the tiling pattern must be antiferromagnetic; further, since we performed DMRG with open boundary conditions, the antiferromagnetic tiling must be striped along the open boundary.The simplest and densest such tiling pattern is shown in Figure S6( is the superlattice pattern.Then, by the convolution theorem, we have M (q) = f (q)S(q), which implies that the Bragg peaks in the magnetization should be located at the reciprocal superlattice vectors, and the peak height should be proportional to f (q * ) where q * is the peak location in momentum space.The size of our DMRG geometry is 6 × 5 unit cells, which does not fit an integer number of filter motifs in the vertical direction.To address this issue, we double the lattice vertically (attaching an identical copy of the same wavefunction) to form a 6 × 10 lattice (since we applied periodic boundary conditions in the vertical direction).The

FIG. 1 .
FIG. 1.The Kitaev-Heisenberg model and schematic of the Quantum-Classical (QuCl) approach.a) Honeycomb lattice of Kitaev model with bond-dependent interactions indicated by the three different colored bonds.b) Phase diagram of Kitaev-Heisenberg model in an external magnetic field (Eq (1)) long three axes of Kz, h, and J. c) A schematic description of QuCl: (i) From a pair of variational wavefunctions|Ψ0⟩ and |Ψ1⟩, labeled projective measurements ("snapshots") B0 and B1 are generated.(ii) The collection of labeled snapshots are used to train the correlator convolutional neural network (CCNN).(iii) The CCNN is configured with four filters (k = 0, • • • , 3) each with three channels, for the binary classification problem minimizing the distance between the prediction ŷ and the label.(iv) Once the training is completed, we fix the filters and use regularization path analysis to reveal signature motifs of the two phases, 0 and 1, under consideration.The correlator weight β onsetting upon reduction of the regularization strength λ to a negative (positive) value signals a feature of the phase 0 (1).

FIG. 2 .
FIG. 2. Gapless Z2 v.s.Heisenberg ordered phase.a) Ground state wavefunctions from the gapless Z2 phase and ferromagnetically ordered phase are obtained at the points on the J-axis marked by stars.The highlighted box shows the two most informative filters for the classification task.The pink and blue dots in filters correspond to projections in z and y basis.b) The regularization path analysis results pointing to the filters in panel a as signature motifs of the ordered phase.
layer (Figure 1(c); green arrows) according to ŷ (n) k with regularization that penalizes the magnitude of the β (n) k 's with strength λ (see Methods).

FIG. 3 .
FIG. 3. The chiral spin liquid v.s. the intermediate gapless phase, benchmarking the indicator of the chiral spin liquid phase.a) Ground state wavefunctions from the chiral spin liquid (CSL) phase and the intermediate gapless phase (IGP) are obtained at the points on the h-axis marked by stars.The highlighted box shows the most informative filter that signifies the CSL phase.Inset shows the bond anisotropy for the three colorings of bonds, e.g.x refers to a S x i S x j coupling.b) The plaquette operator Wp is a operator defined on the six sites around a hexagonal plaquette.c) A sample snapshot from the CSL phase showing the measurement basis that makes the plaquette operator Wp accessible.d) The regularization path analysis pointing to the six-point correlator of filter in panel a as indicator of the CSL.

FIG. 4 .
FIG. 4. The chiral spin liquid v.s. the intermediate gapless phase, discovering the indicator of the intermediate gapless phase.a) The rotated basis vectors in relation to the cardinal axes; the external field is along e3.The signature of the intermediate gapless phase (IGP) is targetted using e1 basis snapshots from the same pair of wavefunctions as in Figure 3(a).b) The regularization path analysis results pointing to the filter in panel c as an indicator of the IGP.c) Most significant filter learned by the correlator convolutional neural network to be associated with the IGP.d) Fourier transform of filter in panel c, with black lines indicating first and extended Brillouin zones.Green circles mark six Bragg peaks associated with the filter tiling pattern in panel e. e) Simplest possible tiling of filter shown in panel c, resulting in a superlattice of antiferromagnetic stripes.The Bragg peaks of this tiling pattern are marked by green circles in panels d and g. f) On-site magnetization ⟨S e 1 i ⟩ of the wavefunction |Ψ0⟩ in the IGP.g) Real part of Fourier transform of panel f, again with Bragg peaks of antiferromagnetic tiling marked (imaginary part is negligible).h) Real part of Fourier transform of ⟨S e 1 i ⟩ of the wavefunction |Ψ1⟩ in the CSL phase showing no discernable features.i) The perpendicular magnetization ⟨S e 1 (r)⟩ as a function of distance from boundary for various values of field strength h, showing decreasing modulation period with increasing field strength.Solid lines show fitted curve based on Equation 3.

FIG. 5 .
FIG. 5. Gapless v.s.gapped Z2 spin liquid, benchmarking distinguishing two spin liquids.a) Ground state wavefunctions from the gapless Z2 phase and gapped Z2 phase are obtained at the points on the Kz-axis marked by stars.The highlighted boxes shows the most informative filters that signifies the each phases.b) The regularization path analysis results associate two-point correlators of filters 0 and 2 to the gapped phase and that of filter 1 to the gapless phase.

FIG. 6 .
FIG. 6. Visualization of correlator convolutional neural network architecture.Input snapshots B α (x) (with α enumerating the spin channels) are shaped into rectangular data arrays with a 1 × 2 unit cell.The inputs are convolved with filters to obtain orders n = 2...6 of correlations C (n) k (x) (Equation4), then spatially averaged to obtain correlator estimates c FIG. S1.Interpretation of regularization path analysis in model distinguishing Z2 gapless spin liquid from Heisenberg ferromagnet.a) All filters learned by neural network.Each column represents a separate filter, while the rows represent the different channels in the filters.b) Regularization path analysis showing relevant correlation functions used by neural network.
FIG. S3.Interpretation of regularization path analysis in model distinguishing intermediate gapless phase from chiral spin liquid, using e1 basis snapshots.a) All filters learned by neural network.Each column represents a separate filter, while the rows represent the different channels in the filters.b) Regularization path analysis showing relevant correlation functions used by neural network.c,d) Real-space magnetization in e1 axis for gapless and CSL phases, respectively.

28 U( 1 )
FIG. S4.Top block: Results for CSL and gapless U(1) phase obtained by 60-site DMRG.Bottom block: Results for gapless U(1) phase at different fields.a) Momentum space resolution of the 60-site cluster.b,c) The magnetization or spinon density modulation of the CSL phase in real and momentum space;and d,e) showing that of the gapless U(1) phase in real and momentum space.f) Momentum space resolution of the 120-site cluster.g,h) Magnetization in the U(1) gapless phase at field h = 0.28; and i,j) at h = 0.32, respectively, as well as their momentum space pattern.
FIG. S6.Visualization of Bragg peaks, amplitude predictions, and observed data.a) Superlattice S(x) showing antiferromagnetic tiling pattern of filter.b) Fourier transform of panel a, with the six Bragg peaks within the extended Brillouin zone marked.c) Filter f (x) approximated by discrete 0 and 1 weights.d) Fourier transform of panel c. e) Numerically obtained magnetization from DMRG lattice, where the lattice has been doubled in the vertical direction as explained in the text.f) Real part of Fourier transform of panel e, showing Bragg peaks at the same locations as shown in panel b with varying amplitudes and phases.Imaginary part is negligible and thus not shown.