## Machine learning topological energy braiding of non-Bloch bands

Shuwei Shi, 1,2 Shibing Chu, 1,2, ∗ Yuee Xie, 1,2 and Yuanping Chen 1,2, † 1 School of Physics and Electronic Engineering, Jiangsu University, Zhenjiang, Jiangsu 212013, People's Republic of China 2 Jiangsu Engineering Research Center on Quantum Perception and Intelligent Detection of Agricultural Information, Zhenjiang 212013, People's Republic of China

(Dated: August 5, 2024)

Machine learning has been used to identify phase transitions in a variety of physical systems. However, there is still a lack of relevant research on non-Bloch energy braiding in non-Hermitian systems. In this work, we study non-Bloch energy braiding in one-dimensional non-Hermitian systems using unsupervised and supervised methods. In unsupervised learning, we use diffusion maps to successfully identify non-Bloch energy braiding without any prior knowledge and combine it with k-means to cluster different topological elements into clusters, such as Unlink and Hopf link. In supervised learning, we train a Convolutional Neural Network (CNN) based on Bloch energy data to predict not only Bloch energy braiding but also non-Bloch energy braiding with an accuracy approaching 100% . By analysing the CNN, we can ascertain that the network has successfully acquired the ability to recognise the braiding topology of the energy bands. The present study demonstrates the considerable potential of machine learning in the identification of non-Hermitian topological phases and energy braiding.

## I. INTRODUCTION

In nature, braids and knots are prevalent structural forms across various fields[1]. These concepts are often employed to describe topological phase transitions, where systems with braiding structures exhibit rich and stable topological properties. Unlike Hermitian systems, non-Hermitian systems feature complex eigenenergies, leading to unique phenomena such as exceptional points[2-7] and non-Hermitian skin effect[8-14]. Studies have shown that the complex energy bands in one-dimensional non-Hermitian systems can form intricate and diverse braiding topologies, which are classified through the conjugacy classes of braid groups[11, 1522]. These complex braiding structures have been experimentally verified on platforms including photonic systems[19, 23], acoustic systems[22, 24, 25] and quantum circuits[26]. However, identifying complex energy braiding[27], especially non-Bloch braiding[28], remains challenging. An efficient identification method is urgently needed.

Machine learning methods have rapidly developed and found applications in fields such as image recognition[29, 30], natural language processing[23], recommender systems[31], and motion data analysis[32]. Their powerful ability to analyze large-scale data and uncover hidden patterns has led to advances in condensed matter physics, quantum physics, and astrophysics[33]. In physics, supervised learning has been used to identify phase transitions and topological invariants[34, 35]. Unsupervised methods, such as diffusion maps combined with k-means, can identify phase transition points and Bloch braiding quickly and in a purely data-driven manner without prior physical knowledge[27, 36-40]. However, most current research has focused on Bloch bands under periodic boundary conditions, yielding idealized results.

∗ c@ujs.edu.cn

† ypchen@ujs.edu.cn

Non-Hermitian systems under open boundary conditions are more representative of realistic systems, exhibiting phenomena such as the non-Hermitian skin effect, where eigenenergies collapse. This disrupts the topological braids of the Bloch bands and invalidates the topological invariants defined in the Brillouin zone, making non-Bloch braiding identification more challenging[28].

In this work, we utilize diffusion maps and Convolutional Neural Networks (CNN) to predict energy braiding in nonHermitian systems. We employ diffusion maps to identify non-Bloch braiding without prior knowledge and train a CNN to predict Bloch energy braiding, achieving high accuracy even for topological elements not included in the training set. Our study demonstrates the potential of machine learning for studying phase transitions and energy braiding in nonHermitian systems.

## II. ALGORITHMS

## A. unsupervised algorithm

When trying to identify topological phase transitions in non-Hermitian systems, linear downscaling methods (e.g., PCA) may not be effective in capturing the topology structure because they assume that the main direction of change in the data is linear. In contrast, nonlinear dimensionality reduction methods (e.g., diffusion maps) can capture the topology structure by constructing a neighbourhood map that preserves this nonlinear structure in the low-dimensional space[36].

For example, consider a two-band non-Hermitian system with the Hamiltonian expressed as:

$$H ( k ) = d ( k ) \hat { \sigma } = d _ { x } ( k ) \hat { \sigma } _ { x } + d _ { y } ( k ) \hat { \sigma } _ { y } + d _ { z } ( k ) \hat { \sigma } _ { z }, \quad ( 1 )$$

where ˆ σ x,y,z are the Pauli matrices and k is the Bloch wave vector. We can construct such a dataset X = { x 1 , x 2 , · · · x l } , where x i = d i ( k / E ) ( i + ( k ) -E i -( k )) rep-

resenting the i -th normalized Hamiltonian sample with E ± ( k ) = ± √ d 2 x ( k ) + d 2 y ( k ) + d 2 z ( k ) .

The distance between x i and x j is defined as: M i,j = ∥ x i -x j ∥ 2 p . In recent research, there has been discussion regarding the utilization of p = 1 , 2 , ∞ cases in unsupervised clustering topological phases[37] . In this paper, we set p = 1 , unless otherwise specified. The similarity matrix is derived from the Gaussian kernel

$$A _ { i, j } = \exp \left ( - \frac { M _ { i, j } } { 2 \epsilon N ^ { 2 } } \right ),$$

where ϵ is the Gaussian kernel coefficient( 0 &lt; ϵ ≪ 1 ). N is an adjustable parameter related to the specific model. A i,j ∈ [0 , 1] represents the similarity between x i and x j . The probability transition matrix is defined as P i,j = A i,j / ∑ N j =1 A i,j . After 2 t steps of random walk, we can use diffusion distance to describe the connectivity between x i and x j

$$D _ { t } ( i, j ) = \sum _ { k = 1 } ^ { l } \frac { \left ( P _ { i, k } ^ { t } - P _ { j, k } ^ { t } \right ) ^ { 2 } } { \sum _ { q } A _ { k, q } } = \sum _ { k = 1 } ^ { l - 1 } \lambda _ { k } ^ { 2 t } [ ( \psi _ { k } ) _ { i } - ( \psi _ { k } ) _ { j } ] ^ { 2 }, \quad \text{$\stackrel{i,\searrow$ breaks} \\ ( 3 ) } \quad \text{$\text{terizes} \\ \text{where $\lambda$, and $\sqrt{i,i}$, are the $i$th eisonvalue and right to onevector}$} \quad \text{$\text{the}$}$$

where λ k and ψ k are the k -th eigenvalue and right eigenvector of P , respectively.

After prolonged diffusion, only the few components ψ k with λ k ≈ 1 dominate the diffusion process. Therefore, the manifold diffusion distance information of the samples is encoded in these components, thus enabling the downscaling of high dimensional data[36]. By further applying clustering algorithms (e.g., k-means) to low-dimensional data, it is possible to cluster different topological phases without prior knowledge. Previous research has demonstrated that the number of λ k ≈ 1 corresponds to the the number of topological clusters[37].

## B. supervised algorithm

For supervised learning, we use a CNN as shown in Fig. 1.

FIG. 1. Schematic of the CNN. Here d r and d i represent the real and imaginary parts of the input data, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01141v1-figure-d98ce9b3-4917-4b2c-8724-4c13084837b6.jpg)

## III. MODEL AND RESULTS

## A. unsupervised learning results

Let us begin with a two-band non-Hermitian lattice as shown in Fig. 2(a). The Hamiltonian of this lattice model takes the form

$$H & = \sum _ { n } \sum _ { \sigma \in \{ A, B \} } \left ( \Delta ^ { - } c _ { n, \sigma } ^ { \dagger } c _ { n + 1, \sigma } + \Delta ^ { + } c _ { n + 1, \sigma } ^ { \dagger } c _ { n, \sigma } \right ) \\ & \quad + \frac { m } { 2 } \left ( c _ { n + 1, A } ^ { \dagger } c _ { n, B } - c _ { n + 1, B } ^ { \dagger } c _ { n, A } + H. c. \right ) \\ \vdots \\.$$

here c † n,σ ( c n,σ ) denotes the creation (annihilation) operator of the sublattice σ ∈ { A,B } on the n th site. The terms ∆ ± = (∆ ± δ / ) 2 represent the non-reciprocal nearest-neighbour hopping within the same sublattice. The nonzero parameter δ breaks the Hermicity of the model. The parameter m characterizes the reciprocal nearest-neighbor hopping between distinct sublattices, while γ denotes the intracell coupling. All parameters are real[28].

The form of the real-space Hamiltonian expressed by Eq.(4) can be written in the form of the Bloch Hamiltonian in momentum space as H k ( ) = ∆cos k -iδ sin k + ( γ + m sin k σ )ˆ y . We can obtain the non-Bloch Hamiltonian H β ( ) by extending the BZ into GBZ with β = re ik , r = √ | (∆ + δ + im / ) (∆ -δ -im ) | , k ∈ R . The energies of the non-Bloch bands take the form E ± ( β ) = ± γ +0 5(∆ . -δ ∓ im β ) +0 5(∆+ . δ ± im β ) -1 . The non-Bloch winding number is defined as

$$W = \frac { 1 } { 2 \pi i } \oint _ { C _ { \beta } } d \ln \det \left [ H ( \beta ) - \frac { 1 } { 2 } \text{Tr} [ H ( \beta ) ] \right ], \quad ( 5 )$$

which describes the braid degree of two non-Bloch bands[28].

In order to explore the non-Bloch braiding in this model as shown in Fig. 2(b), we consider a dataset X = { x 1 , x 2 , · · · x l } with m = 0 5 . , ∆ = 1 , δ = 0 7 . and variable γ ∈ [0 , 1 2] . , where x i = d i ( β / E ) ( i + ( β ) -E i -( β )) . The results of unsupervised learning are shown in Fig. 2(c-d). From Fig. 2(c), the similarity matrix is separated into two blocks, which corresponding to the two largest eigenvalues λ ≈ 1 of the diffusion matrix shown in Fig. 2(d). As a result, the diffusion maps identifies two distinct topological phases. In addition, we can get the phase transition point γ ≈ 0 58 . form Fig. 2(c), which is consistent with the theoretical value γ c = m r ( 2 +1) 2 / r ≈ 0 579 . .

We compare the results of unsupervised learning with those calculated in Eq. 5. For γ &lt; γ c , W = 2 , corresponding to Hopf link, and for γ &gt; γ c , W = 0 , corresponding to Unlink, as shown in Fig. 2(b). In Fig. 2(e-f), we display the braiding of non-Bloch bands in (Re [ E ] ,Im [ E k ] , ) space with γ = 0 3 . and γ = 0 8 . , corresponding to the blue and red topological clusters in the inset of Fig. 2(d), respectively.

FIG. 2. (a) Sketch of the two bands non-Hermitian lattice model. (b) Phase diagram of the non-Hermitian system under open boundary condition with ∆ = 1 0 . and δ = 0 7 . . (c) The visualisation of the Gaussian kernel. The breakpoint indicates the predicted phase transition point γ ≈ 0 58 . . (d) The first eight eigenvalues of the diffusion matrix P . The inset shows the result of k-means. The first two largest eigenvalues correspond to the two topological clusters in the inset. The samples are generated uniformly of γ ∈ [0 , 1 2] . with ∆ = 1 , δ = 0 7 . and m = 0 5 . . The Gaussian kernel coefficient ϵ = 0 01 . . (e-f) Complex energy braiding of non-Bloch bands in (Re[ E , ] Im[ E , k ] ) space. (e) and (f) are topologically distinct and correspond to Hopf link and Unlink, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01141v1-figure-6ba85dd3-6f57-4e34-8561-b943f96775c9.jpg)

## B. supervised and transfer learning results

## 1. learning Bloch braiding

In this section, we consider a two-band non-Hermitian lattice as shown in Fig. 3(a). The Hamiltonian of this lattice model takes the form

$$H _ { ( m, n ) } = \begin{pmatrix} 0 & t _ { 0 } \\ t _ { 0 } + t _ { m } e ^ { i m k } + t _ { n } e ^ { i n k } & 0 \end{pmatrix}, \quad ( 6 ) \quad \breve { \text{the C} }$$

1 kernel of size ( 40 1 1 , , ), followed by a fully connected layer before the output layer. Similar to but different from previous studies[41, 42], only the first convolutional layer is followed by the nonlinear activation function ReLU. We find that the prediction accuracy of the CNN depends on the first convolutional layer and its activation function. Shallow networks are more effective than deep networks in learning braid degree. Specifically, simply increasing the number of layers of the CNN can instead cause it to lose its generalisation ability.

where t 0 is reciprocal intracell coupling. t m and t n are nonreciprocal intercell hoppings that spans m and n lattices, respectively[22].

̸

The complex energy spectra reads E ± ( k ) = ± √ t 0 ( t 0 + t m imk e + t n e ik ) . For any two energy bands satisfy the separable band condition( E + ( k ) = E -( k ) for all k ∈ [0 , 2 π ] ), a topological invariant W can be defined to classify their braiding topology[22],

$$W = \frac { 1 } { 2 \pi i } \oint _ { 0 } ^ { 2 \pi } \left ( \frac { d \ln E _ { + } } { d k } + \frac { d \ln E _ { - } } { d k } \right ) d k. \quad ( 7 ) \quad \text{where}$$

We trained a generic CNN as shown in Fig. 1 to predict the Bloch energy braiding in this model. The CNN in our research has two convolution layers with 40 kernels of size ( 1 2 2 , , ) and

The input data are normalised spectral correlation configurations d ( n ) = { Re[ E 2 (2 πn/L )] , Im[ E 2 (2 πn/L )] } , E = E + ( k ) -E -( k ) . In the following, we set L = 201 , which is large enough for the CNN to capture the topology of the complex energy bands.

The labels corresponding to each configuration are calculated by the Eq. 7. The output of the CNN is a real number w representing the braid degree. The objective function to be optimized is defined by

$$J = \frac { 1 } { N } \sum _ { i = 1 } ^ { N } ( w _ { i } - W _ { i } ) ^ { 2 },$$

where w i and W i are, respectively, the predicted and true braid degree of i th sample, and N represents the size of the training data set. We uniformly sample 9 × 10 4 configurations from H (1 2) , and H (1 3) , , containing Unlink ( W = 0) ,

FIG. 3. (a) Sketch of the two-band non-Hermitian lattice model H ( m,n ) . (b) The loss with each training and validation. (c-g) The predicted phase diagrams, along with the corresponding braid degrees across different parameter spaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01141v1-figure-03e16e65-1fbf-4741-9b34-37b7a002d665.jpg)

Unknot ( W = 1) , Hopf link ( W = 2) and Trefoil ( W = 3) , respectively. There are five different types of test data, (i,ii) two test sets with 1 × 10 4 configurations generated from H (1 2) , and H (1 3) , , respectively, (iii,iv,v) 1 × 10 4 configurations sampled from H (2 3) , , H (3 4) , and H (3 5) , , respectively. The batch size for training is set to 50 and the epoch is set to 50.

After training, the losses on the training and validation sets are shown in Fig. 3(b). It can be observed that the neural network converges very quickly and the losses stay at a low level and no overfitting occurs. We test with the test data and round the output to plot the phase diagram of the predictions on the different test sets as shown in Fig. 3(c-g). In Fig. 3(c) and (d), the CNN successfully predicts the phase diagram using the test sets (i) and (ii), which seems natural since our train-

FIG. 4. The phase diagram predicted by the pre-trained CNN in the γ -m plane with ∆ = 1 0 . and δ = 0 9 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01141v1-figure-fb436eda-0d40-46de-adc4-64d565d0023a.jpg)

ing data are also sampled with the same parameters. However, we can verify the effectiveness of the CNN in Fig. 3(e-g). In Fig. 3(e), the CNN successfully predicts the braid degrees under the H (2 3) , case, which are contained in the training set. In Fig. 3(f) and (g), the CNN also predicts the braid phases with high accuracy, even predicting two braid degrees, W = 4 and W = 5 , that are not included in the training set. At this point, we have demonstrated the ability of the CNN to predict the braid degree under different parameter spaces of the same model.

## 2. transfer learning non-Bloch braiding

Further, we use the above pre-trained CNN to recognise non-Bloch braiding in the non-Hermitian system as shown in Fig. 2(a). The input data can be chosen as d ( n ) = { Re[ E 2 (2 πn/L )] , Im[ E 2 (2 πn/L )] } , E = E + ( β ) -E -( β ) with ∆ = 1 0 . and δ = 0 9 . , varying γ ∈ [0 , 1 2] . and m ∈ [0 , 1] . The predicted phase diagram with respect to the parameters γ and m is shown in Fig. 4. The dashed line is the theoretical value of the phase boundary. Without retraining, the CNN almost perfectly predicts the presence of Unlink and Hopf link elements accompanied by the exact phase boundary in the phase diagram.

To further understand the workflow of the CNN, we analyse the training details. The first layer of the CNN is a convolutional layer with 40 convolutional kernels to perform convolutional operations on the input data,

$$\text{al net-} _ { \boldsymbol x \text{ level} } \quad B ^ { i } ( n ) = f [ A _ { 1 1 } ^ { i } d _ { r } ( 2 \pi ( n - 1 ) / L ) + A _ { 1 2 } ^ { i } d _ { i } ( 2 \pi ( n - 1 ) / L ) \\ \text{round} _ { \text{on the} } \quad \quad + A _ { 2 1 } ^ { i } d _ { r } ( 2 \pi n / L ) + A _ { 2 2 } ^ { i } d _ { i } ( 2 \pi n / L ) + A _ { 0 } ^ { i } ],$$

where A is a kernel of size 2 × 2 , i = 1 , ..., 40 , n = 1 , ..., L and f ( x ) is the activation function. The second layer is a convolutional layer without an activation function that takes

the form,

$$D ^ { i } ( n ) = \sum _ { i = 1 } ^ { 4 0 } C ^ { i } B ^ { i } ( n ) + C _ { 0 } ^ { i }, \quad \ \ ( 1 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where C i is the 1 × 1 kernel. D n i ( ) depend on the winding angle of the two bands. Finally, all 40 neurons are mapped to a single output neuron to obtain the braid degree,

$$W = \sum _ { n = 1 } ^ { L } M _ { n } D ( n ) + M _ { 0 }. \quad \quad ( 1 1 ) \quad \text{ple}$$

During the training process, the neural network successfully fit all the parameters above. Through the above analysis, we find that the Eq. 9 used to extract the braid information is universal, and then the output of the Eq. 9 can be linearly summed to obtain the braid degree.

## IV. CONCLUSION

In summary, we applied unsupervised and supervised methods to identify the energy braiding in non-Hermitian sys-

- [1] C. C. Adams, The knot book (American Mathematical Soc., 1994).
- [2] K. Kawabata, T. Bessho, and M. Sato, Phys. Rev. Lett. 123 , 066405 (2019).
- [3] W. B. Rui, Y. X. Zhao, and A. P. Schnyder, Phys. Rev. B 99 , 241110 (2019).
- [4] X. Zhang and J. Gong, Phys. Rev. B 101 , 045415 (2020).
- [5] L. Xiao, T. Deng, K. Wang, Z. Wang, W. Yi, and P. Xue, Phys. Rev. Lett. 126 , 230402 (2021).
- [6] S. Sayyad, M. St˚ alhammar, L. Rødland, and F. Kunst, SciPost Physics 15 (2023), 10.21468/SciPostPhys.15.5.200.
- [7] B.-B. Wang, Y. Ge, S.-Q. Yuan, D. Jia, and H.-X. Sun, Progress In Electromagnetics Research 176 , 1 (2023).
- [8] S. Yao and Z. Wang, Physical Review Letters 121 , 086803 (2018).
- [9] F. Song, S. Yao, and Z. Wang, Phys. Rev. Lett. 123 , 170401 (2019).
- [10] L. Li, C. H. Lee, S. Mu, and J. Gong, Nature Communications 11 , 5491 (2020).
- [11] L. Zhang, Y. Yang, Y. Ge, Y.-J. Guan, Q. Chen, Q. Yan, F. Chen, R. Xi, Y. Li, D. Jia, S.-Q. Yuan, H.-X. Sun, H. Chen, and B. Zhang, Nature Communications 12 , 6297 (2021).
- [12] Q. Liang, D. Xie, Z. Dong, H. Li, H. Li, B. Gadway, W. Yi, and B. Yan, Phys. Rev. Lett. 129 , 070401 (2022).
- [13] K. Kawabata, T. Numasawa, and S. Ryu, Phys. Rev. X 13 , 021007 (2023).
- [14] Q. Zhou, J. Wu, Z. Pu, J. Lu, X. Huang, W. Deng, M. Ke, and Z. Liu, Nature Communications 14 , 4569 (2023).
- [15] H. Shen, B. Zhen, and L. Fu, Phys. Rev. Lett. 120 , 146402 (2018).
- [16] K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Phys. Rev. X 9 , 041015 (2019).

tems. The unsupervised learning algorithm, diffusion maps, can identify the phase transition points and separate different braiding elements into clusters. In supervised learning, we trained a CNN to capture the braiding topology of complex energy in non-Hermitian systems. It is worth noting that the CNN trained under Bloch conditions can predict not only Bloch braiding but also non-Bloch braiding. Our methods can be extended to identify other knots and links in non-Bloch bands. Our study demonstrates the potential of machine learning in identifying non-Hermitian topological phases and complex energy braiding.

## V. DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available upon reasonable request from the authors.

## VI. ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (No. 12074150, 12174157 and 11904137).

- [17] C. C. Wojcik, X.-Q. Sun, T. c. v. Bzduˇ sek, and S. Fan, Phys. Rev. B 101 , 205417 (2020).
- [18] Z. Li and R. S. K. Mong, Phys. Rev. B 103 , 155129 (2021).
- [19] K. Wang, A. Dutt, C. C. Wojcik, and S. Fan, Nature 598 , 59 (2021).
- [20] K. Wang, A. Dutt, K. Y. Yang, C. C. Wojcik, J. Vuˇ ckovi´ c, and S. Fan, Science 371 , 1240 (2021).
- [21] Y. S. S. Patil, J. H¨ oller, P. A. Henry, C. Guria, Y. Zhang, L. Jiang, N. Kralj, N. Read, and J. G. E. Harris, Nature 607 , 271 (2022).
- [22] Q. Zhang, Y. Li, H. Sun, X. Liu, L. Zhao, X. Feng, X. Fan, and C. Qiu, Phys. Rev. Lett. 130 , 017201 (2023).
- [23] J. Hirschberg and C. Manning, Science (New York, N.Y.) 349 , 261 (2015).
- [24] H. Qiu, Q. Zhang, T. Liu, X. Fan, F. Zhang, and C. Qiu, Nature Communications 14 , 1261 (2023).
- [25] Q. Zhang, L. Zhao, X. Liu, X. Feng, L. Xiong, W. Wu, and C. Qiu, Phys. Rev. Res. 5 , L022050 (2023).
- [26] Y. Yu, L.-W. Yu, W. Zhang, H. Zhang, X. Ouyang, Y. Liu, D.L. Deng, and L.-M. Duan, npj Quantum Information 8 , 116 (2022).
- [27] J. Chen, Z. Wang, Y.-T. Tan, C. Wang, and J. Ren, Communications Physics 7 , 209 (2024).
- [28] Y. Li, X. Ji, Y . Chen, X. Yan, and X. Yang, Phys. Rev. B 106 , 195425 (2022).
- [29] K. He, X. Zhang, S. Ren, and J. Sun, in 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR) (2016) pp. 770-778.
- [30] D. Zhang, Y. Liu, Y. Zhao, J. Liang, B. Sun, and S. Chu, Applied Sciences 13 (2023), 10.3390/app132312649.
- [31] G. Adomavicius and A. Tuzhilin, IEEE Transactions on Knowledge and Data Engineering 17 , 734 (2005).

- [32] J. Yin, M. Chen, Y. Ge, Q. Song, and H. Zheng, Journal of Sensors 2022 , 3907002 (2022).
- [33] G. Carleo, I. Cirac, K. Cranmer, L. Daudet, M. Schuld, N. Tishby, L. Vogt-Maranto, and L. Zdeborov´ a, Rev. Mod. Phys. 91 , 045002 (2019).
- [34] B. Narayan and A. Narayan, Phys. Rev. B 103 , 035413 (2021).
- [35] L.-F. Zhang, L.-Z. Tang, Z.-H. Huang, G.-Q. Zhang, W. Huang, and D.-W. Zhang, Phys. Rev. A 103 , 012419 (2021).
- [36] R. R. Coifman and S. Lafon, Applied and Computational Harmonic Analysis 21 , 5 (2006), special Issue: Diffusion Maps and Wavelets.
- [37] Y. Che, C. Gneiting, T. Liu, and F. Nori, Phys. Rev. B 102 , 134213 (2020).
- [38] L.-W. Yu and D.-L. Deng, Phys. Rev. Lett. 126 , 240402 (2021).
- [39] Q. Zhou, J. Wu, Z. Pu, J. Lu, X. . J. Huang, Y. Yu, L.-W. Yu, W. Zhang, H. Zhang, X. Ouyang, Y. Liu, D.-L. Deng, and L.M. Duan, npj Quantum Information 8 , 116 (2022).
- [40] Y. Long, H. Xue, and B. Zhang, Nature Machine Intelligence (2024).
- [41] P. Zhang, H. Shen, and H. Zhai, Phys. Rev. Lett. 120 , 066401 (2018).
- [42] L.-F. Zhang, L.-Z. Tang, Z.-H. Huang, G.-Q. Zhang, W. Huang, and D.-W. Zhang, Phys. Rev. A 103 , 012419 (2021).
Shuwei Shi Shibing Chu Yuee Xie Yuanping Chen
Machine learning topological energy braiding of non-Bloch bands
Machine learning has been used to identify phase transitions in a variety of
physical systems. However, there is still a lack of relevant research on
non-Bloch energy braiding in non-Hermitian systems. In this work, we study
non-Bloch energy braiding in one-dimensional non-Hermitian systems using
unsupervised and supervised methods. In unsupervised learning, we use diffusion
maps to successfully identify non-Bloch energy braiding without any prior
knowledge and combine it with k-means to cluster different topological elements
into clusters, such as Unlink and Hopf link. In supervised learning, we train a
Convolutional Neural Network (CNN) based on Bloch energy data to predict not
only Bloch energy braiding but also non-Bloch energy braiding with an accuracy
approaching 100%. By analysing the CNN, we can ascertain that the network has
successfully acquired the ability to recognise the braiding topology of the
energy bands. The present study demonstrates the considerable potential of
machine learning in the identification of non-Hermitian topological phases and
energy braiding.