## Featuring the topology with the unsupervised machine learning

Kenji Fukushima, 1, 2, ∗ Shotaro Shiba Funai, 3, † and Hideaki Iida 1, ‡

1 Department of Physics, The University of Tokyo,

7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan

2 Institute for Physics of Intelligence (IPI), The University of Tokyo,

7-3-1

Hongo, Bunkyo-ku, Tokyo 113-0033, Japan

3

Physics and Biology Unit, Okinawa Institute of Science and Technology (OIST), 1919-1 Tancha Onna-son, Kunigami-gun, Okinawa 904-0495, Japan

Images of line drawings are generally composed of primitive elements. One of the most fundamental elements to characterize images is the topology; line segments belong to a category different from closed circles, and closed circles with different winding degrees are nonequivalent. We investigate images with nontrivial winding using the unsupervised machine learning. We build an autoencoder model with a combination of convolutional and fully connected neural networks. We confirm that compressed data filtered from the trained model retain more than 90% of correct information on the topology, evidencing that image clustering from the unsupervised learning features the topology.

∗ fuku@nt.phys.s.u-tokyo.ac.jp

† shotaro.funai@oist.jp

‡ iida@nt.phys.s.u-tokyo.ac.jp

## I. INTRODUCTION

Brains ingeniously function with networks of neurons. For understanding of intrinsic brain dynamics, physicists would favorably decompose such an integral system into irreducible elements, so that we can analyze relatively simpler function of each building element that takes rather primitive actions. Then, numerical simulations on computer are handy devices to test if a postulated mechanism of brains should go as expected. Such modeling embodies non-equilibrium processes of brains, which is an approach acknowledged commonly as computational neuroscience. Besides, a hypothesis called quantum brain dynamics implements quantum fluctuations and Nambu-Goldstone bosons for brain sciences [1, 2] (for discussions for/against quantum phenomena in brain dynamics, see Ref. [3]), which bridges a devide between computational neuroscience and modern physics.

In contrast to such 'off-equilibrium' problems, in the language of physics, perception and recognition are 'static' problems. For the latter problems, model-independent research tools are available for computer simulations. That is, the machine learning enables us to emulate the neural structure of brains on computer. One of intriguing attributes of the machine learning, particularly with deep neural networks (NNs), is that any nonlinear mapping can be represented by data transmission through multiple hidden layers [4].

These days we have witnessed tremendous progresses in the field of image recognition and classification by means of the machine learning. In particular the progress has been driven by Convolutional Neural Network (CNN) [5], which was originally proposed as a multi-layer neural network imitating animal's visual cortex [6]. The CNN has become the most common approach for high-level image recognition since an overwhelming victory of AlexNet [7] at 'ImageNet Large Scale Visual Recognition Challenge 2012.' Challenges of minimizing interclass variability, reducing error rate, achieving large-scale image recognition, etc are ongoing improvements and they are all crucial steps for practical usages.

Physicswise, the image handling with the deep learning has proved its strength in identifying phase transitions. Some successful attempts are found in Ref. [8] for two-dimensional systems analyzed by the supervised learning, Ref. [9] for its extension to three-dimensional systems, and Ref. [10] for statistical systems studied by the unsupervised learning. Here, we point out an essential difference between the supervised and unsupervised learning; the former is useful for regression and grouping problems, while the latter efficiently makes fea-

ture extraction and clustering of data. Interestingly, similarity between the unsupervised learning and the renormalization group in physics has also been investigated, see Ref. [11].

In the context of image recognition, at the same time, a distinct direction toward more fundamental research would be as important to demystify blackboxed artificial intelligence, which may be somehow beneficial for so-called explainable artificial intelligence [12, 13]. The fundamental question of our interest in the present work is what would be the simplest element of images that categorizes those images into representative clusters. For the sake of image clustering, a useful mathematical notion, which underlies modern physics, has been developed known as the 'topology' theorized into the form of homotopy. The most wellknown example is that a mug with one handle and a torus-shaped donut belong to the same grouping class; the shape can be smoothly deformed from one to the other, and they are of the same homotopy type. In this work we report leading-edge results from our simulations with the CNN supporting an idea that the topology is critical information for image feature extraction and clustering.

## II. TOPOLOGY AND THE WINDING NUMBER

The topology is classified by the homotopy group in mathematics. The simplest example is what is called the fundamental homotopy group denoted as π 1 ( S 1 ) = Z associated with a mapping from S 1 (i.e., one dimensional unit sphere) to another S 1 and an integer n W ∈ Z corresponds to the winding number. To demonstrate the idea concretely, let us consider the following function on S 1 of U(1),

$$\phi ( x ) = \mathrm e ^ { \mathrm i \theta ( x ) } = \cos \theta ( x ) + \mathrm i \sin \theta ( x ) \,.$$

If x is a coordinate on a circle with period L , the above function represents a mapping from S 1 in coordinate space to S 1 on Gauss' plane with Euler's angle θ (which is also called the 'lift' in homotopy theory). While x travels around from 0 to L under a condition, φ (0) = φ L ( ), Euler's angle θ should return to the original position modulo 2 π . The winding number associated with the above function (or the 'degree' of this function) reads,

$$n _ { W } = \frac { \theta ( L ) - \theta ( 0 ) } { 2 \pi } = \frac { \ln [ \phi ( L ) / \phi ( 0 ) ] } { 2 \pi i } = \frac { 1 } { 2 \pi i } \int _ { 0 } ^ { L } \mathrm d x \, \phi ^ { - 1 } ( x ) \frac { \mathrm d \phi ( x ) } { \mathrm d x } \,.$$

Figure 1 schematically illustrates one winding configuration of φ x ( ) having n W = 1.

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

j

F

o

/

7

Q

c

8

N

I

0

B

W

J

H

k

w

3

&gt;

A

X

V

S

E

O

K

9

r

G

Y

2

M

d

L

u

R

y

C

+

z

v

Z

q

g

p

m

f

P

n

U

T

5

D

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

j

F

o

/

7

Q

c

8

N

I

0

B

W

J

H

k

w

3

&gt;

A

X

V

S

E

O

K

9

r

G

Y

2

M

d

L

u

R

y

C

+

z

v

Z

q

g

p

m

f

P

n

U

T

5

D

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

j

F

o

/

7

Q

c

8

N

I

0

B

W

J

H

k

w

3

&gt;

A

X

V

S

E

O

K

9

r

G

Y

2

M

d

L

u

R

y

C

+

z

v

Z

q

g

p

m

f

P

n

U

T

5

D

FIG. 1. Schematic illustration of π 1 ( S 1 ) realized by φ x ( ) = e i θ x ( ) . In the left representation the behavior of θ x ( ) is manifest and n W = 1 is easily concluded, but in our simulation, only (Re φ, Im φ ) as shown in the right is given to the NN model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-1c791d48-2ac4-4cb4-b85e-c7e8661883b1.jpg)

For us, human-beings, it would be an elementary-class exercise to discover a counting rule of n W . If we were given many different images with correct answers of n W , it would be just a matter of time for us to eventually find a right counting rule out. This description is nothing but the machinery of the supervised learning . Interestingly, it has been reported that the deep neural network trained by the supervised learning correctly reproduced n W of π 1 ( S 1 ) [14]. There, the machine discovered a nice fit of the formula (2) only from the information of given (Re φ, Im φ ) and n W , but not referring to θ x ( ) directly. In the present work, we are taking one step forward; we would like to see the NN model not only optimizing a fit from the supervised learning but featuring the topology from the unsupervised learning .

More specifically, we would like to think of classification of many images without giving the answers of n W . It would be very intriguing to ask such a question of how the CNN makes clustering of differently winding data, which would be a prototype model of how our brains categorize images based on the topological characterization. Thus, the unsupervised learning as adopted in this work should tell us surpassing information than the supervised learning for the purpose to dissect the topological contents.

## III. RESULTS AND DISCUSSIONS

In the numerical procedure we represent φ x ( ) by a sequence of numbers on discretized x with L = 128 grids, i.e., we generate 2 × 128 = 256 sized data of (Re φ , i Im φ i ) with

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

k

W

u

Y

c

o

w

n

B

v

z

D

+

F

R

5

g

M

A

J

U

&gt;

C

Z

H

V

N

L

G

q

K

I

E

3

S

T

Q

0

O

X

d

8

/

f

j

p

P

7

r

m

9

2

y

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

E

u

w

P

/

S

B

q

k

G

J

A

o

5

H

W

T

M

Q

&gt;

C

X

c

V

7

g

N

F

D

Z

3

O

0

K

9

m

n

8

v

2

I

+

Y

R

z

d

p

r

L

f

y

U

j

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

p

N

j

+

m

G

g

3

Y

f

Q

&gt;

A

C

X

c

V

F

L

w

B

H

v

q

o

W

y

E

S

O

D

k

R

n

K

I

0

Z

u

9

d

5

z

P

/

7

J

8

M

r

U

T

2

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

Z

D

T

u

Y

y

N

v

c

J

F

E

U

0

A

&gt;

C

X

V

H

L

g

R

B

3

M

w

2

m

r

j

I

Q

S

d

P

9

k

5

z

p

n

o

8

f

q

O

W

K

+

/

G

7

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

G

9

0

p

r

D

3

f

Y

V

c

E

+

U

C

o

k

&gt;

A

X

7

S

g

N

B

F

Z

/

R

q

I

w

j

m

L

K

y

u

z

2

O

H

M

8

Q

J

T

d

P

W

5

n

v

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

u

W

d

m

G

+

y

V

f

7

2

j

v

Z

N

D

&gt;

A

C

X

c

H

L

S

F

/

p

o

U

3

R

g

B

q

r

k

J

8

0

I

Y

O

E

n

Q

z

w

5

M

P

K

T

9

FIG. 2. One example of generated data for (+ + + , , , -, +) with n W = +3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-faa2223d-3a50-49e3-929e-436b088465d7.jpg)

i = 1 , . . . , L under the periodic boundary condition; (Re φ , L Im φ L ) = (Re φ , 1 Im φ 1 ). These 256 numbers consist of the input data onto the CNN side.

We prepare the training and test data randomly. We will give more detailed explanations on the numerical procedure in Method. Each data consists of distinct N s segments along x with either positive or negative winding, where the segment lengths, glyph[lscript] m , are chosen randomly, where N s ∑ m =1 glyph[lscript] m = L should be kept. For the data used in this work, we take account of N s = 0 1 2 3 4, and 5. , , , , We then assign positive and negative winding randomly to each segment, which is symbolically labeled as ( p , p 1 2 , . . . , p N s ) with p m = ± 1, where + (and -) stands for positive (and negative, respectively) winding. For example, (+ + + , , , -, +) for N s = 5 is a configuration with net winding, n W = N s ∑ m =1 p m = +3. In m -th segment we first postulate that θ x ( ) linearly changes its value by 2 πp m . Then, zigzag lines are distorted with random noises to enhance the learning quality and thus the adaptivity. Figure 2 depicts an example of generated data with a choice of (+ + + , , , -, +). With N s up to 5 2 0 + 2 1 + 2 2 + 2 3 + 2 4 + 2 5 = 63 winding patterns are possible, and n W can take a value from -5 to +5. This setup is for the moment sufficiently general for our goal to check performance of the topology detection.

The unsupervised learning utilizes the autoencoder [15]; it first encodes the data compressed into smaller number of neurons in the CNN (in our case, 16 sites × 4 filters from original 256 sites) and then decodes the compressed data with fully connected NN into the original size. We repeat such encoding and decoding processes to minimize the loss function

FIG. 3. Examples of feature maps on the deepest CNN layer. (Left) N s = 5 and all possible (i.e., 2 5 = 32) winding patterns averaged over 1,000 test data. (Right) N s = 5 with n W = +3 fixed and five winding patterns averaged over 1,000 test data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-6c3bdebe-a432-4b50-9370-90e5656b5a1e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-41e61831-3ab7-466c-8be4-b236405190f4.jpg)

measured by the squared difference between the original input data and the coarse-grained output data. The learning process optimizes the filters in the CNN encoder and simultaneously the weights in the NN decoder.

We shall see the results from our NN model that has been optimized by the unsupervised learning with the training dataset of 63 winding patterns times 1,000 randomly generated data (i.e., 63,000 data in total). We input the test data into the optimized NN model and observe feature maps of 16 neurons convoluted with 4 filters on the deepest CNN layer. Figure 3 summarizes sampled feature maps from the 1st filter. The left plot is for N s = 5 with all possible winding patterns averaged over 1,000 test data. The right plot particularly picks up the averaged feature maps for five different windings with N s = 5 and n W = +3. At a glance one may think that the behavior looks like coarse-grained Re φ x ( ) or Im φ x ( ) with segment lengths normalized.

Now, the most fascinated question is whether the compressed data on 16 sites convoluted with 4 filters could retain information on the topology or not, and if yes, how we can retrieve it. From the left of Fig. 3 it is obvious that the peak heights reflect different winding patterns. In fact, the averaged feature maps exhibit a clear hierarchy of four separated heights with one-to-one correspondance to the winding sequence; that is, the peak heights increase with sequential windings as

$$( +, + ) \ < \ ( +, - ) \ < \ ( -, + ) \ < \ ( -, - ),$$

and the height at the far right end is determined solely by ± , which is due to zero-padding

FIG. 4. Feature maps from the 1st filter (top-left), the 2nd filter (top-right), the 3rd filter (bottomleft), and the 4th filter (bottom-right) for (+ , + + , , -, +) data without taking the average. Four randomly selected data are shown by four different colors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-baa77ac5-cb7b-49c9-a1d1-073f8dec33f9.jpg)

in the convolution to keep the data size.

For example, if we see the far left peak around the 2nd neuron in the right of Fig. 3, ( -, + + + +) has the highest peak, , , , (+ , -, + + +) the second highest, , , and others are degenerated in accord with Eq. (3). Also, we notice in the right of Fig. 3 that, for (+ + + + , , , , -), three consecutive short peaks appear from the left, one middle peak follows, and then one tall peak sits at the right end. The short peaks correspond to (+ +) of , ((+ +) + + , , , , -), (+ (+ +) + , , , , -), and (+ + (+ +) , , , , -). The middle peak corresponds to (+ + + (+ , , , , -)) and the tall peak is sensitive only to -of (+ + + + , , , , -).

One might have thought that such a clear hierarchy as in Fig. 3 is visible only after taking the average. This is indeed the case, as exemplified in Fig. 4; the feature maps for one test data do not always show prominent peaks with well separated heights. Nevertheless, surprisingly, we can see that such fluctuating coarse-grained data in Fig. 4 still retain information on the topology!

It is impossible for our eyes to recognize any topological contents from Fig. 4, so we will

FIG. 5. Correct answer rate for n W guessed from the feature maps. For each of 63,000 test data we have checked where the correct winding number ranks (1st, 2nd, ... in the legend) in the probability output from the supervised NN model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-e1b5c059-a6ad-4ce2-bfe5-3eae5f5ac44a.jpg)

ask for a help of another machine learning device. For each of 63,000 training data, we have such feature maps like Fig. 4 and also the corresponding n W . We can then perform the supervised learning to train a fully connected NN (with one hidden layer) such that the output gives the probability distribution of guessed n W (out of -5 , . . . , 5) in response to the input of feature maps. Figure 5 is the correct answer rate. If we input the feature maps from only 1 filter, the most-likely n W hits the correct value at the rate of 50%, and the second-likely one at the rate of 31%. If we use the feature maps from 2 filters, the available information is doubled, and the correct answer rate of the most-likely n W becomes 78%. Amazingly, for the feature maps from 3 filters, it increases up to 90%! There is almost no difference between the 3-filter and the 4-filter results, and it seems that the correct answer rate is saturated at 90%.

To summarize, from these results and analyses, we can conclude that the coarse-grained data of the feature maps through the unsupervised learning do retain information on the topology. In other words, the unsupervised learning with autoencoder can provide us with a nontrivial machinery to compress data without losing the topological contents.

## IV. DISCUSSION

We demonstrated that the unsupervised machine learning of images makes feature extraction without losing information on the topology. This evidences that the winding number corresponding to the fundamental group, π 1 ( S 1 ), should be one of the essential indices that characterize image clustering. We trained an autoencoder with the CNN and the fully connected NN for the unsupervised learning using randomly generated data of functions from x on S 1 to a U(1) element. We found that the averaged feature maps on the deepest CNN layer show a clear hierarchy pattern of the configurations with one-to-one correspondence to the winding sequences. With help of the supervised learning technique we also revealed that feature maps for each image data look coarse-grained images and such compressed data retain information on the topology correctly.

The extension of the present work, i.e., the unsupervised learning of higher-dimensional images with nontrivial winding would be quite interesting. We note that the supervised machine learning has been utilized for π 2 ( S 2 ) [14], but as we illustrated in this work, the unsupervised machine learning would be more interesting. Implications from this extension include intriguing applications in quantum field theories. In fact, some classical solutions of the equations of motion in quantum field theory are topologically stabilized. In such cases the field configurations are classified according to the winding number. For representative examples, π 2 (SU(2) / U(1)) = Z for monopoles, π 3 (SU(2)) = Z for Skyrmions (with no timedependence), and π 3 (SU(2)) = Z for instantons (with the fields at infinity identified) in pure Yang-Mills theory [16, 17]. Actually, it is a long-standing problem how to visualize the topological contents in quantum field theories. The numerical lattice simulation is a powerful tool to solve quantum field theories, and field configurations should in principle contain information on the topology. Some algorithms to extract topologically nontrivial configurations such as monopoles, Skyrmions, and instantons have been proposed [18-20]. The most well-known approach, i.e., the cooling method has a serious flaw, however. If the cooling is applied too many times, the topology is lost and the field configuration would become trivially flat. Therefore, the cooling procedures should be stopped at some point, and this artificial termination of the procedures causes uncertainties. Alternatively, we would emphasize that the compression of field configuration images by means of the unsupervised machine learning is a promising candidate for the superior smearing algorithm not to lose

the topological contents. We are testing this idea in simple two-dimensional lattice model, namely, C P N -1 model that has π 2 ( S 2 ) instantons.

Mathematically, it would be also a very interesting question to consider not only the homotopy groups but also the homology H n ( X ) of a topological space X defined by a coset of cycles in n dimensions over boundaries of ( n + 1)-dimensional elements. For example, H X 0 ( ) counts the number of connected drawings of images, and H X 1 ( ) counts the number of loops of images, etc. This direction of research is now ongoing.

## ACKNOWLEDGMENTS

We thank Yuya Abe and Yuki Fujimoto for discussions. K. F. was supported by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant No. 18H01211.

## METHOD

## Input data of winding patterns

Input data φ i = e i θ i ( i = 1 , . . . , L ), including training data and test data, are generated in the following way. Note that since the data φ i are complex numbers, we actually input the combinations of their real and imaginary parts (Re φ , i Im φ i ).

First we impose the periodic boundary condition:

$$\text{Re} \, \phi _ { 1 } = \text{Re} \, \phi _ { L }, \quad \text{Im} \, \phi _ { 1 } = \text{Im} \, \phi _ { L }.$$

Then we divide L sites into N s segments, and length of each segment glyph[lscript] m ( m = 1 , . . . , N s ) is randomly chosen as

$$\lambda _ { 0 } = \frac { L - 1 } { N _ { s } } \left [ 1 + 0. 4 ( \xi _ { m } - \xi _ { m - 1 } ) \right ] _ { 0 } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where ξ m is a random number from a uniform distribution in the open interval ( -1 1). , This means the lengthes in all the segments satisfy 0 2 . L -1 N s &lt; glyph[lscript] m &lt; 1 8 . L -1 N s . We set ξ 0 = ξ N s = 0 so that the total length L is kept. To be exact, the length glyph[lscript] m should be an integer, so the right hand side is rounded to an integer.

In each segment m , the angle θ i is composed of a linear part from 0 to ± 2 π and a random noise:

$$\frac { \theta _ { i } } { 2 \pi } = p _ { m } \frac { i - i _ { 0 } } { \ell _ { m } } + 0. 1 \zeta _ { i }$$

FIG. 6. Schematic figure of the autoencoder used for our unsupervised learning.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-259737ce-f5ff-4d67-9c7f-1dcab9670da2.jpg)

where p m = ± 1 is the winding direction and i 0 = 1 + ∑ m -1 m ′ =1 glyph[lscript] m ′ . The random noise ζ i is from a Gaussian distribution with mean 0 and variance 1.

In our experiment, we set L = 128 and N s = 0 1 , , . . . , 5. Then all the combinations of winding directions p m have ∑ 5 N s =0 2 N s = 63 patterns. For each winding pattern, we generate 1,000 + 1,000 input data with the parameters ( ξ m , ζ i ) chosen randomly. The first 1,000 data (in total 63,000 data) is the training data, which is used for training our autoencoder and supervised NN. The other 1,000 data is the test data for analyzing the feature maps in the autoencoder (see Figs. 3 and 4) and the output from the supervised NN (see Fig. 5).

## Autoencoder

The autoencoder for our unsupervised learning consists of the CNN (encoder) and the fully-connected NN (decoder). A schematic figure of our autoencoder is shown in Fig. 6. We made a training code using TensorFlow [21].

The CNN encoder has two layers. In the both layers, we use the convolution with 8 × 1( × 2 channel) sized filters and stride 1. We also use the zero padding to keep the size of data, and the ReLU as an activation function. In the convolution part, difference between the first and second layers is only the number of filters. After the convolution, our encoder has the pooling part in each layer. We use the max pooling with 2 × 1 sized window and stride

FIG. 7. Loss function during training of our autoencoder.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00281v1-figure-5af075d2-4117-41ea-a178-89d2bdce8a1a.jpg)

2 for each channel in the first layer. This pooling compresses the data size into 1 / stride of the original size. The second layer has 4 × 1 sized window and stride 4, then as a result, the CNN encoder outputs the feature map with L/ / 2 4 = 16 sites (for each filter out of 4 filters), as we saw in Fig. 3.

The fully-connected NN decoder has two layers, too. In the first layer, we use the dropout method with probability 0.5 to avoid overlearning, then use the ReLU again. The final layer has no dropout and no activation function, then outputs coarse-grained data ( ϕ i, 1 , ϕ i, 2 ) with the same size as the input data (Re φ , i Im φ i ).

As the loss function for the training, we choose the squared difference between input data φ i and output data ϕ i,j , that is,

$$\frac { 1 } { 2 L } \sum _ { i = 1 } ^ { L } \left [ \left ( \text{Re} \, \phi _ { i } - \varphi _ { i, 1 } \right ) ^ { 2 } + \left ( \text{Im} \, \phi _ { i } - \varphi _ { i, 2 } \right ) ^ { 2 } \right ].$$

We prepared the training and test data, both of which contain 63 winding patterns times 1000 randomly generated data. Then we found the unsupervised learning with the learning rate 10 -7 and the mini-batch size 10 decreases the loss function of the test data to its minimum around 6000 epochs, as shown in Fig. 7.

- [2] M. Jibu and K. Yasue, Quantum Brain Dynamics and Consciousness: An Introduction , Advances in consciousness research (J. Benjamins Publishing Company, 1995).
- [3] P. Jedlicka, Frontiers in Molecular Neuroscience 10 , 366 (2017).
- [4] Y. LeCun, Y. Bengio, and G. Hinton, Nature 521 , 436 (2015).
- [5] Y. LeCun, B. Boser, J. S. Denker, D. Henderson, R. E. Howard, W. Hubbard, and L. D. Jackel, Neural Computation 1 , 541 (1989).
- [6] K. Fukushima, Biological Cybernetics 36 , 193 (1980).
- [7] A. Krizhevsky, I. Sutskever, and G. E. Hinton, in Advances in neural information processing systems (2012) pp. 1097-1105.
- [8] T. Ohtsuki and T. Ohtsuki, Journal of the Physical Society of Japan 85 , 123706 (2016).
- [9] T. Ohtsuki and T. Ohtsuki, Journal of the Physical Society of Japan 86 , 044708 (2017).
- [10] S. S. Funai and D. Giataganas, (2018), arXiv:1810.08179 [cond-mat.stat-mech].
- [11] S. Iso, S. Shiba, and S. Yokoo, Phys. Rev. E97 , 053304 (2018), arXiv:1801.07172 [hep-th].
- [12] W. Samek, T. Wiegand, and K. M¨ uller, CoRR abs/1708.08296 (2017), arXiv:1708.08296.
- [13] F. K. Doˇ silovi´, M. Brˇi´, c c c and N. Hlupi´, in c 2018 41st International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO) (2018) pp. 0210-0215.
- [14] P. Zhang, H. Shen, and H. Zhai, Phys. Rev. Lett. 120 , 066401 (2018).
- [15] G. E. Hinton and R. R. Salakhutdinov, Science 313 , 504 (2006).
- [16] A. A. Belavin, A. M. Polyakov, A. S. Schwartz, and Yu. S. Tyupkin, Phys. Lett. B59 , 85 (1975).
- [17] G. 't Hooft, Phys. Rev. D14 , 3432 (1976).
- [18] B. Berg, Phys. Lett. 104B , 475 (1981).
- [19] M. Teper, Phys. Lett. B171 , 86 (1986).
- [20] E.-M. Ilgenfritz, M. M¨ uller-Preu β ker, G. Schierholz, H. Schiller, et al. , Proceedings, 23RD International Conference on High Energy Physics, JULY 16-23, 1986, Berkeley, CA , Nucl. Phys. B268 , 693 (1986).
- [21] M. Abadi, P. Barham, J. Chen, Z. Chen, A. Davis, J. Dean, M. Devin, S. Ghemawat, G. Irving, M. Isard, M. Kudlur, J. Levenberg, R. Monga, S. Moore, D. G. Murray, B. Steiner, P. Tucker, V. Vasudevan, P. Warden, M. Wicke, Y. Yu, and X. Zheng, in 12th USENIX Symposium on Operating Systems Design and Implementation (OSDI 16) (USENIX Association, Savannah,

GA, 2016) pp. 265-283.
Kenji Fukushima Shotaro Shiba Funai Hideaki Iida
Featuring the topology with the unsupervised machine learning
Images of line drawings are generally composed of primitive elements. One of
the most fundamental elements to characterize images is the topology; line
segments belong to a category different from closed circles, and closed circles
with different winding degrees are nonequivalent. We investigate images with
nontrivial winding using the unsupervised machine learning. We build an
autoencoder model with a combination of convolutional and fully connected
neural networks. We confirm that compressed data filtered from the trained
model retain more than 90% of correct information on the topology, evidencing
that image clustering from the unsupervised learning features the topology.