Mutual Information and Correlations across Topological Phase Transitions in Topologically Ordered Graphene Zigzag Nanoribbons

Graphene zigzag nanoribbons, initially in a topologically ordered state, undergo a topological phase transition into crossover phases distinguished by quasi-topological order. We computed mutual information for both the topologically ordered phase and its crossover phases, revealing the following results: (i) In the topologically ordered phase, A-chirality carbon lines strongly entangle with B-chirality carbon lines on the opposite side of the zigzag ribbon. This entanglement persists but weakens in crossover phases. (ii) The upper zigzag edge entangles with non-edge lines of different chirality on the opposite side of the ribbon. (iii) Entanglement increases as more carbon lines are grouped together, regardless of the lines’ chirality. No long-range entanglement was found in the symmetry-protected phase in the absence of disorder.


I. INTRODUCTION
The topological order of graphene [1][2][3] zigzag nanoribbons [4][5][6][7] (ZGNRs) manifests in the weak disorder regime through the presence of e/2 solitonic fractional charges [8][9][10][11][12] localized at opposite edges of the ribbon and spin-charge separation [13][14][15].In this regime, interacting disordered ZGNRs exhibit a universal value of topological entanglement entropy [16], providing evidence for the existence of long-range entanglement and, consequently, topological order in these systems [17,18].Additionally, in interacting disordered ZGNRs, the shape of the entanglement spectrum closely resembles the density of states of the edge states [19], aligning with the expectations for a topologically ordered system, as proposed by Li and Haldane [20].Furthermore, a phase diagram of the topological order of ZGNRs predicts crossover phases that exhibit non-universal values of topological entanglement entropy, each characterized by distinct properties in the regimes of intermediate on-site repulsion and intermediate disorder strength [21].In one of the crossover phases, crossover phase I, fractional charges and spin-charge separation are notably absent.However, charge transfer correlations persist between the zigzag edges.In the other crossover phase crossover phase II, fractional charges are present, but no correlations are observed between the opposite zigzag edges.The experimental verification of these theoretical predictions has an excellent chance of success since the fabrication of ZGNRs has achieved atomic precision [22][23][24].
In this study, we investigate how the entanglement pattern changes during the transition from either a crossover phase or a symmetry-protected phase to the topologically ordered phase.Our investigation is motivated by several factors.Firstly, near the usual quantum-critical point, electrons cease to behave independently and, instead, become entangled in large numbers.Entanglement patterns exhibit intriguing features near a quantum critical point, as previously described in works by Sachdev et al. [25,26].This quantum phase transition shares some similarities with the smooth topological phase transitions from crossover I/II phases to the topologically ordered phase.It is worthwhile to explore how the entanglement patterns evolve between these two phases.Secondly, the transition from disorder-free interacting ZGNRs [27] (symmetry-protected phase) to their interacting disordered counterparts (topologically ordered phase) represents an abrupt topological phase transition.In other words, disorder acts as a singular potential that significantly alters the system's physics.It remains unclear how the entanglement pattern changes during this transition.
To address these issues, we have computed the mutual information between different carbon lines in zigzag ribbons, and our findings are summarized in Figure 1.In the topologically ordered phase, A-chirality carbon lines are strongly entangled with B-chirality carbon lines on the opposite side of the zigzag ribbon, but not with other lines.In crossover phases, they are also entangled but to a lesser degree than in the topologically ordered phase.Furthermore, the upper zigzag edge entangles with carbon lines of opposite chirality in the lower part of the ribbon, as indicated by the solid lines in Figure 1.However, entanglement increases as more carbon lines are grouped together, regardless of the lines' chirality.Notably, we did not observe long-range entanglement in the symmetry-protected phase in the absence of disorder, but it increases suddenly with the addition of disorder.

Crossover phase
FIG. 1.A zigzag ribbon consists of horizontal carbon lines with A (red sites) and B (blue sites) chirality, as shown in the upper figure (disorder is not shown).The upper zigzag edge consists of A-type carbon atoms, while the lower zigzag edge consists of B-type carbon atoms (mirror carbon lines have opposite chirality.)In the ground state, the upper edge sites have net spin-up while the lower edge sites have net spin-down.Below is illustrated the entanglement of carbon lines represented by red and blue straight lines, in both topological order and crossover phases.The thickness of the solid curved lines represents the strength of the entanglement.Dashed lines present weak entanglement between lines not connected by mirror symmetry.
The rest of this paper is organized as follows.Section II introduces the Hamiltonian model of ZGNRs.Section III presents the method for calculating mutual information and provides relevant notations.Sections IV and V are devoted to presenting our numerical results of mutual information in topologically ordered and crossover phases, respectively.Comments and conclusions drawn from our results will be given in Section VI.

II. HAMILTONIAN MODEL
We employ the self-consistent Hartree-Fock approximation for our numerical calculations, as it has been proven to work effectively for graphene systems [28], and its results are consistent with those obtained through density functional theory and density matrix renormalization group methods [19,29,30].
Disordered ZGNRs can be described within the Mott-Anderson Hamiltonian [31][32][33][34] and may be approximated using the Hartree-Fock approximation: Here, c † m,σ , n m,σ , and s mx(y) represent creation operators, occupation numbers, and spin operators at site m with spin σ (σ =↑ (↓) for spin-up (down)).The summation with 'n.n.' specifies nearest neighbor hopping.The second term represents a short-range disorder potential [35], with 10% of the total ribbon's sites containing impurities.The magnitude of the potential at each site is drawn from a uniform distribution [−Γ, Γ].
There is nothing special about the value of 10%, i.e., one can also use different values.We examined ZGNRs with other values of disorder concentration n imp and disorder strength Γ.These results are presented in a previous paper [13].Note that, in the self-consistent Born approximation, the disorder strength is not solely characterized by a single parameter, either Γ or n imp , but by their product Γ √ n imp [36].Additionally, we have investigated ZGNRs with off-diagonal disorder in Ref. [15].The third and fourth terms arise from on-site repulsion U .The fourth term accounts for self-consistent "magnetic fields", where ⟨h mx ⟩ = −2U ⟨s mx ⟩ and ⟨h my ⟩ = −2U ⟨s my ⟩.Note that this term is zero in the undoped case because the ground state spins are collinear.The decoupling of the onsite term used in the Hartree-Fock Hamiltonian is derived from where the spin rasing and lowering operators are with site spin operators (A detailed derivation of this term using Hartree-Fock approximation can be found in chapter 6 of Ref. [5].)The relevant parameters in our study include on-site repulsion, disorder strength, and the number of doped electrons, denoted as (U, Γ, δN ).The total number of carbon lines, denoted as W , and the total number of sites on each line, represented by L, are assigned as the width and length of the ZGNR, respectively.Therefore, the total number of carbon sites is given by N s = L × W .A Hartree-Fock single-particle state, denoted as |k⟩ (k = 1, 2, . . ., 2N s ), can be expressed as a linear combination of site states |m, σ⟩, with the coefficients A k,m,σ representing the self-consistent Hartree-Fock eigenvectors.In the language of second quantization, this relation can be represented as The eigenstates of mean-field Hamiltonian in Equation (1) are given by With such representation of the eigenstates, the expectation values of site spins and site occupation numbers can be written as follows: and where n α = 1(0) for occupied (unoccupied) Hartree-Fock states.
For the Hartree-Fock calculations, the input consists of an antiferromagnetic initial state for the disorder-free case at half-filling (Γ = 0, δN = 0), and a paramagnetic initial state otherwise.The paramagnetic initial state exhibits a small imbalance between spin-up and spin-down occupations at every site, resulting in HF eigenvectors with spin splitting [13,31].The paramagnetic initial state produces a state with fractional charges.The dimension of the Hartree-Fock matrix can scale up to 16,000, which corresponds to the longest ribbon considered in this study.The Hartree-Fock eigenstates and eigenenergies are computed self-consistently through 20 iterations to achieve convergence.We employed GPUs to expedite the solution of the Hartree-Fock matrix.Additionally, due to the extensive nature of the calculations involving numerous disorder realizations, GPU computations in this paper were conducted on a supercomputer.Note that the ground state is doubly degenerate, with sites reversed.In undoped ribbons, site spins are collinear, while in doped ribbons, they are non-collinear.

III. MUTUAL INFORMATION
In topologically ordered ZGNRs, two peaks in the probability density of fractional charges are located on the zigzag edge with nearly the same horizontal coordinates, as previously reported [13].Additionally, occupation numbers are correlated between pairs of carbon lines.These properties have motivated us to examine mutual information between carbon lines in a ZGNR.Mutual information is defined as follows: where i and j can represent single or multiple carbon lines (domains).Here, S(i) and S(i, j) denote the entanglement entropy of spin-up electrons within their respective areas.The entanglement entropy of domain i is defined as the von-Neumann entropy [37] S where ρ i is the reduced density matrix.Within the Hartree-Fock approximation, it is comparatively convenient to adopt Peschel's method and compute the correlation matrix between sites m and n restricted within domain i [38] where |Ψ⟩ is one of the degenerate ground states.The reduced density matrix and spin correlation matrix C mn are related by the following identity: Thus, the entanglement entropy is equal to the sum of binary entropies [39] of each eigenvalue These measures can also be applied to spin-down electrons.Mutual information, M (i, j), quantifies the total amount of correlations between domains i and j, and it obeys boundary laws [40].The notation for S(i, j) is based on the following convention: For the purposes of our analysis, we maintain a fixed ribbon width of W = 32.We employ the following notation: (i, j) = (1, 32), which constitutes the first pair, comprising a line and its mirror counterpart in the opposite half of the ribbon.Continuing in this manner, (i, j) = (2, 31) corresponds to the second pair, signifying that it is the closest pair to the outermost one, while (i, j) = (3, 30) forms the third pair, and so forth.It is important to note that the two lines in these pairs are symmetrically positioned along the width of the ribbon.A visual representation of this notation is provided in Figure 2.

IV. RESULTS FOR TOPOLOGICALLY ORDERED PHASE
Let us delve into the entanglement patterns of the topologically ordered phase.To accurately quantify the longrange entanglement, it is essential that the two lines under consideration be separated by a distance greater than the correlation length.Figure 3 displays the dependence of the average mutual information between two symmetrically positioned lines on their distance, computed for ribbons with a width of W = 32.The mutual information values for the outermost pair differ significantly between the topologically ordered phase and symmetry-protected phase, implying the presence of long-range entanglement in the topologically ordered phase.Otherwise, the trends of mutual information in both the disorder-free and weakly disordered cases exhibit similarities, suggesting a comparable correlation length in both scenarios.As we approach the innermost pair, the mutual information values abruptly increase, indicating that the distance between two lines decreases to less than the correlation length.These results imply that the two lines, denoted as i and j, should be separated by at least 10 units of the lattice constant to ensure accurate mutual information calculations.Consequently, all mutual information calculations hereafter will be performed with a fixed ribbon width of W = 32, focusing primarily on the first few outermost pairs.There are several other noteworthy properties.Additional properties of mutual information for a weakly disordered ZGNR with (U, Γ, δN ) = (t, 0.1t, 0) are presented in Figure 4a,b.Firstly, mutual information follows the area law, showing a linear scaling with the length of the ribbon.This behavior arises from the fact that the number of mixed chiral states, including those carrying fractional charges, is proportional to the length of the ribbon [15].We also observe that, when both i and j consist of two lines, the mutual information follows the same linear scaling, as depicted in Figure 4c.The entanglement entropies S(i) in the double-line configuration are smaller than those in the singleline configuration, as shown in Table I.However, intriguingly, the mutual information in the double-line configuration is larger than that observed when i and j each consist of a single line.Specifically, the mutual information for (i, j) = (1 ∪ 2, 31 ∪ 32) is twice as large as (i, j) = (1, 32), while the mutual information for (i, j) = (2 ∪ 3, 30 ∪ 31) is one and a half times bigger than (i, j) = (3, 30) (the third pair).(a) Distances between symmetrically positioned lines are illustrated graphically.Numerical results for mutual information between symmetrically positioned lines (along the ribbon's width) as a function of their distance in the topologically ordered phase (b) with (U, Γ, δN, ND) = (t, 0.1t, 0, 20) and in the symmetry-protected phase (c) with (U, Γ, δN, ND) = (t, 0, 0, 1).
We find that mutual information vanishes for the even-numbered second, fourth, sixth, and so on, outermost pairs, while correlations are observed among the odd-numbered first, third, fifth, and so forth counterparts.As the pair number increases, the mutual information between lines in the pair decreases.Additionally, there is a non-zero correlation between one edge and the 30th carbon line, and its mutual information (M (1, 30)) is even higher than that of the third pair (M (3,30)).These properties can be explained by the mixed chiral states depicted in Figure 5.These states consist of two nonlocal parts residing on the edges of the A and B sublattices, with their probability densities peaking at these edges and rapidly decaying within the ribbon.Consequently, the density distributed on the first line is greater than that on the third line, so we can deduce that M (1, 30) > M (3,30).Note also that The ribbon is in the topologically ordered phase.(a) M (i, j) for (i, j) = (1, 32), (3,30), (1,30).(b) M (i, j) for (i, j) = (2, 30), (1,31), (2,31).Due to the fact that the corresponding values of mutual information are nearly zero, even numbered carbon lines in the upper half of ZGNR barely entangle with their symmetric lines in the lower half of the ribbon.(c) In the double-line configuration, we depict M (i, j) for (i, j) = (1 ∪ 2, 31 ∪ 32) and (2 ∪ 3, 30 ∪ 31) by using the values of M (i, j) for (1, 32) and (3, 30) for comparison.Note that the scales of y-axis differ among (a-c).
TABLE I. Numerical values of entanglement entropy, denoted as S(i) and S(j), were computed for (U, Γ) = (t, 0.1t) using a single disorder realization (ND = 1).In the case of double lines, values of entanglement entropy are smaller than those in the single-line case, whereas the mutual information M (i, j) in the double-line case is larger.odd-numbered pairs exhibit substantial probability density on both their A and B carbon lines, while even-numbered pairs exhibit diminishing probability density on their A and B carbon lines.

< l a t e x i t s h a 1 _ b a s e 6 4 = " +
h z X q x 3 q 0 P 6 3 P W u m L l M 0 c w B + v r F 8 x p n M g = < / l a t e x i t > 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " m s u g T r L 9 t X P f n

V. RESULTS FOR CROSSOVER PHASES
We discuss how the behavior of mutual information is related to the topological transitions into crossover I and crossover II phases of interacting disordered ZGNR.Before we proceed, we briefly mention that at the critical point Γ = 0, where the symmetry-protected phase transitions to the topological order phase, mutual information exhibits a singularity, as seen in Figure 6b.In the limit as Γ/U → 0 with a ribbon length of L = 300, the mutual information between two opposite edges is M (1, 32) ≈ 0.78 for (U, Γ) = (t, 0.1t) and M (1, 32) ≈ 0.20 for (U, Γ) = (t, 0).In the crossover I and crossover II phases, mutual information decreases upon increasing the strength of disorder and on-site repulsion, as seen in Figure 6.By increasing the strength of the disorder and on-site repulsion, the interacting disorder ZGNR undergoes a phase transition into a crossover phase and eventually a non-topologically ordered phase.Furthermore, the strength of the disorder affects mutual information in a manner similar to the impact of on-site repulsion, not only for the first pair but also for other lines: mutual information decreases as the disorder strength Γ increases.
It is noteworthy that the behavior of mutual information mirrors that of topological entanglement entropy.Topological entanglement entropy experiences an abrupt increase from zero to a universal value when transitioning from the symmetry-protected phase to the topologically ordered phase [16].However, as the phase of the interacting disordered ZGNR undergoes a crossover from topologically ordered phase to non-topologically ordered phase [21], topological entanglement entropy gradually decreases to zero and exhibits increasing variance.Additionally, in the topologically ordered phase, the mean value of topological entanglement entropy slightly decreases (with unchanged variance) when both on-site repulsion and disorder increase [21].A similar behavior is observed with mutual information, as shown in Figure 6: mutual information and topological entanglement entropy for (U, Γ) = (2t, 0.1t) and (U, Γ) = (t, 2t) are smaller than those for (U, Γ) = (t, 0.1t), and all these values correspond to the topologically ordered phase.
We provide a detailed description of how carbon lines become entangled during the transition in the crossover phases, as depicted in Figure 7.As the topologically ordered phase transitions into crossover with an increase in on-site repulsion (U ), the mutual information of the first pair M (1, 32), the third pair M (3, 30), and M (1, 30) monotonically decreases (see Figure 7a).The mutual information for double-line configurations also follows the same decreasing trend (phase transition into the crossover II phase shows similar results to those of crossover I).This decline in mutual information with increasing U is associated with the decrease in the number of well-localized fractional charges at the edge.Mutual information values associated with pairs with the next outmost lines, such as M (2, 31), M (2, 30), and M (1, 31), also decrease for sufficiently large U > 3, as shown in Figure 7b.These values in the crossover phase show larger fluctuations compared to those shown in Figure 7a (note that, when U ≳ 4, the ribbon is not in the crossover phase).In addition, the phase of interacting disordered ZGNRs can transit to the crossover II phase upon increasing the number of doped electrons δN .The insets in Figure 7c,d show that mutual information has a decreasing trend with higher doping concentration.

VI. DISCUSSION
Mutual information provides a quantitative measure of the entanglement between carbon lines.We calculated mutual information for both the topologically ordered phase and crossover phases.The entanglement correlations between carbon lines in an interacting disordered ZGNR are summarized in Figure 1 as follows: (i) Mutual information scales linearly with the ribbon's length in all cases.(ii) In both the topologically ordered and crossover phases, correlations exist among the first, third, fifth, and so on, pairs, while they are absent in the second, fourth, sixth, and so on, pairs.These properties can be well explained by mixed chiral states.Additionally, correlations exist between a zigzag edge and the even-numbered lines in the opposite half of the ribbon, which have different chirality.Our numerical results also show that all the non-vanishing correlations weaken as on-site repulsion or disorder strength increases.(iii) The mutual information among the outermost four lines (two from the upper and the other two from lower parts of the ribbon) appears to be approximately twice as large as that among outermost pairs (first pair).However, it is worth noting that mutual information alone is not sufficient to distinguish between the crossover I and crossover II phases.(iv) We also found that entanglement increases as more carbon lines are grouped together,

FIG. 2 .
FIG.2.Schematic figure for mutual information between carbon lines in ZGNR.The dashed horizontal lines specify the order of the carbon lines, and red (blue) dots represent A (B) type carbon atoms according to chiral symmetry.The ribbon's width is fixed at W = 32, as explained in the main text.(a) Mutual information between two zigzag edges (or between lines in the pair) with i = 1 and j = 32.(b) Mutual information between the top two lines and the bottom two lines with i = 1 ∪ 2 and j = 31 ∪ 32.

FIG. 5 .
FIG.5.Illustration of the probability distribution of a mixed chiral state in the absence of disorder.The horizontal carbon lines run perpendicular to the page.Solid (dashed) line represents probability density distributed on A (B) sites.Odd-numbered pairs of lines are indicated.

FIG. 6 .
FIG. 6.(a) Mutual information between zigzag edges M (1, 32) with respect to on-site repulsion U .The other parameters are (Γ, δN, ND) = (0.1t, 0, 50).The ribbon is in the undoped crossover I phase.(b) Mutual information between zigzag edges M (1, 32) with respect to disorder strength Γ.The other parameters are (U, δN, ND) = (t, 0, 50).The ribbon is in the undoped crossover II phase.The error bars represent the standard deviation of each data point.A reduction in the mean value is observed with the increase in on-site repulsion and disorder strength, leading to a decrease in the size of the error bars.The ribbon size used in calculations is (L, W ) = (300, 32).