Quantum reservoir computation utilising scale-free networks Akitada Sakurai,1, 2 Marta P. Estarellas,2, ∗ William J. Munro,3, 2 and Kae Nemoto2, 1, † 1School of Multidisciplinary Science, Department of Informatics, SOKENDAI (the Graduate University for Advanced Studies), 2-1-2 Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan 2National Institute of Informatics, 2-1-2 Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan 3NTT Basic Research Laboratories & Research Center for Theoretical Quantum Physics, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa, 243-0198, Japan Today’s quantum processors composed of fifty or more qubits have allowed us to enter a computa- tional era where the output results are not easily simulatable on the world’s biggest supercomputers. What we have not seen yet, however, is whether or not such quantum complexity can be ever useful for any practical applications. A fundamental question behind this lies in the non-trivial relation between the complexity and its computational power. If we find a clue for how and what quantum complexity could boost the computational power, we might be able to directly utilize the quantum complexity to design quantum computation even with the presence of noise and errors. In this work we introduce a new reservoir computational model for pattern recognition showing a quan- tum advantage utilizing scale-free networks. This new scheme allows us to utilize the complexity inherent in the scale-free networks, meaning we do not require programing nor optimization of the quantum layer even for other computational tasks. The simplicity in our approach illustrates the computational power in quantum complexity as well as provide new applications for such processors. I. INTRODUCTION some of which have been demonstrated on gate-based quantum computers [13, 14]. It has been suggested that The recent realization of quantum processors with fifty some quantum neural networks may offer an advantage plus qubits is undoubtedly a key milestone for the nascent over classical neural networks [11], however the optimiza- quantum technology field [1–3]. The quantum advantage tion of a parameterized quantum circuit with a classical has been estimated through its complexity generated by feedback loop would be difficult to scale [15]. Quantum these quantum processors, benchmarking it against the reservoir computation (QRC) can also be considered as necessary run time for a classical computer to simulate a type of QNNs, which are introduced to mainly anal- it. Such a complexity has been believed to be deeply yse temporal/sequential data processing with both qubits connected to its quantum computational power, yet it [16, 17] and continuous variables [18]. In those works still remains unclear in what way we could extract us- quantum neural networks are considered in the Hilbert able computational power for real applications from the space, whereas another approach uses a network of phys- quantum complexity itself. ical qubits [19], which has been used to construct a uni- Quantum neural networks (QNNs) have been consid- versal quantum reservoir computer [20]. ered as one of the potential directions to exploit quantum In this paper, we propose a new quantum reservoir complexity that the current quantum processors (includ- computation (QRC) model based on scale-free networks ing Noisey Intermidate Scale Quantum (NISQ) devices) that emerge during the melting of a discrete time crystal can generate. There has been a number of approaches [21] and apply it to a pattern recognition task. These and definitions for QNNs [4, 5]. In the early days of networks represent the effective Hamiltonian of a period- QNNs, the core ideas behind perceptron (a model of ar- ically driven system obtained through the application of tificial neural networks) had been extended to QNNs [6], tools in Floquet theory [22]. In our model, the classical where the perceptron activation function was replaced by reservoir, that is a large nonlinear system, is replaced by an operator. Despite the simple notion of the quantum a quantum hidden layer characterized by the non-trivial arXiv:2108.12131v2 [quant-ph] 31 Aug 2021 perceptron given by Altaisky, its implementation remains complexity of scale-free networks. Non-trivial complex- highly nontrivial. The main focus of the early works ity as well as randomness has played an important role on quantum perceptron is to integrate the nonlinearity in both quantum and classical neural networks [17, 23]. present in the classical model into quantum mechanics Our QRC model presents several advantages. First [7–10]. More recently QNNs have been extensively inves- and foremost the network grows exponentially with the tigated as a subclass of variational quantum algorithms number of qubits N; 50 qubits could generate a network (VQAs) [11, 12]. VQAs can be modeled with a feature as large as the neural network in a human brain. Second, map and a variational model with classical feedfoward, despite this rapid growth in complexity, the parameter setting for the quantum hidden layer remains simple. In fact, there is no programing nor optimization necessary ∗ Current Affiliation: Qilimanjaro Quantum Tech, Carrer dels for the quantum hidden layer once it has been set, even Comtes de Bell-Lloc, 161, 08014 Barcelona, Spain for different computational tasks. Third, by eliminating † Corresponding author: [email protected] the costly parameter optimization over the variational 2 (a) (b) (d) Dropout 0 1 0 1 2 3 4 5 6 7 M-layer 9 Output layer Scale-free network (c) MNIST Images Input layer Quantum hidden layer M-layer Output layer 0 1 PCA Map 9 Quantum dynamics Measurement ONN + data processsing FIG. 1. A schematic illustration of the structure of the quantum reservoir computer. (a) shows an example of scale-free networks represented in the Hilbert space. Each dot indicates a computational basis state, while each edge is weighted with the hopping strength between the two states. When the hopping strength is smaller than a threshold given by the percolation rule, there is no edge between the states. (b) indicates a quantum system used to implement the quantum hidden layer. Here the model of a discrete time crystal is used. (c) illustrates the total scheme of our QRC model. There are four layers: the input layer, the quantum hidden layer, the M-layer and the output layer. The input layer encodes the input data on to the quantum hidden layer, while the M-layer is to measure the quantum system in the computation basis, converting the quantum data (state) to the classical data. The one-layer neural network (ONN) is responsible for the training. model in VQAs, we can expect a high feasibility of our that the encoding can be done easily and efficiently. For QRC model and a significant speed up in the training this we chose the technique of Principle Component Anal- process. This advantage in the learning speed is similar ysis (PCA), as indicated in Fig. 1(c). After some time to what we expect with the extreme learning machine in in the dynamics of the quantum hidden layer the output the classical machine learning [24]. (k) quantum state |ψf i is measured in the computational (k) basis at the M-layer, converting the quantum data |ψf i to the classical data ~x (k). The classical data is then reno- II. THE MODEL malized for the one-layer neural network (ONN), that is a classical neural network only with the input and output Our quantum reservoir computation (QRC) model layers. Training is performed only at the output layer. consists of four layers as illustrated in Fig.1. The compu- Since the quantum hidden layer is the solo quantum (k) tation starts with encoding the classical data Iij for the component in our scheme, how it is designed is key to k-th image onto the initial state of the quantum hidden achieve a quantum advantage. The quantum hidden layer layer. The quantum hidden layer is the reservoir of this (Fig.1 (b)) and its network (Fig.1 (a)) are connected computational model. It may not be necessary, however through the visualization of the Hamiltonian [22]. The it is desirable to preprocess the classical input data so Hamiltonian of an N-qubit system in the computational 3 basis {|ii} where i is a binary digit from 0 to 2N − 1 may we set for the quantum hidden layer is shorter than half P P be represented as H = i Ei|iihi| + i,j(i6=j) Wij|iihj| of the typical coherence time of an ion trap experiment where Ei and Wij correspond to the energy of the com- with 10 qubits [25]. The discrete time crystal can also putational basis |ii and the transition energy between |ii be implemented on a gate-based quantum processor by a and |ji respectively. The network has an edge between gate decomposition of the Floquet operator or the entire the i-th node (state) and j-th node (state) with weight unitary map for the quantum dynamics [28]. Hence the Wij when the percolation condition |Ej − Ei| < |Wij| is model can be adapted to various quantum computational satisfied [22]. The percolation rule is to eliminate the systems. off-resonant transitions from the network. In the case of the periodically driven system, the network is generated from an effective Hamiltonian using the Floquet theory. III. PATTERN RECOGNITION WITH QRC There are now two approaches to setting the quantum hidden layer in the targeted regime. Using a network Here we illustrate how the QRC model performs pat- with the appropriate complexity (Fig.1 (a)) as a recipe, tern recognition. We use the MNIST handwritten digit we can generate a Hamiltonian with the scale-free nature data set publicly available [29] to benchmark the per- by assigning each computational basis state to each node formance.