Quantum–Classical Hybrid Information Processing
via a Single Quantum System

Recent advances in machine learning (ML) and quantum computing have revolutionized the methodology of processing complex and large-scale data. While merging these fields, classical or quantum systems can generate a massive amount of time series data, such as sensing data or quantum states that flow through multiple quantum channels in a network of quantum devices [1, 2, 3]. This context leads to the requirement of a novel learning paradigm to process these data efficiently, such as the easy manipulation used in training and deployment, while maintaining rich representation capability. Currently, algorithms are being designed on specific homogeneous data, such as quantum-native or classical-native data. However, most quantum devices rely on classical controls [4, 5], such as temperature or signals from electronic controllers [6, 7]. The outputs of these devices are not simply derived from quantum channels and are also considered a function of classical controls and quantum input. A representative example is a quantum switch with classical control, which simulates the indefinite causal order between two operations [8, 9, 10, 11, 12] [Fig. 1(a)]. Therefore, the research on hybrid quantum and classical data processing can lead to broader and near-term applicability for quantum devices. For example, we can use the same resource to learn the tomography of devices receiving both classical and quantum data without doing it separately for each control setting.

Contrary to ML models such as artificial recurrent neural networks on a digital computer, a physical system with rich dynamics can be a good candidate for a learning system within the framework of physical reservoir computing (PRC) [13, 14]. In PRC, the input is fed into a dynamical system called a reservoir to create nonlinear dynamics of input data via sufficiently complex and high-dimensional trajectories [15, 16, 17, 18, 13]. A readout, which outputs a linear combination of the accessible observables in the reservoir, is the only part that needs to be trained without interfering with the reservoir’s internal parameters. Accordingly, the success and efficiency of PRC rely on good physical realizations of the reservoir, which has attracted considerable interest from diverse research fields [14]. The seminal work [19] uses a disordered ensemble quantum dynamics system as a quantum reservoir (QR) to process classical data, with the possibility of having a large number of degrees of freedom. QRs have been developed in various platforms, such as nuclear magnetic resonance (NMR) systems [19, 20, 21], superconducting quantum processors [22, 23], fermions and bosonic models [24, 25, 26, 27], quantum harmonic oscillators [28, 29], arrays of Rydberg atoms [30], and photonic quantum memristors [31]. Several studies have focused on the processing of data in the form of quantum states [24, 26, 25, 32, 33, 34], which provide certain advantages over classical ML methods. However, a QR is yet to be treated as a homogeneous data-driven model because it lacks the ability to deal with hybrid forms of quantum-classical data. Therefore, an unified architecture for hybrid quantum-classical processing is required from theoretical and applied perspectives.

In this study, we establish a framework that considers a QR as an analog processor to process hybrid quantum-classical data. Inspired by Refs. [24, 25], our QR is a network of quantum dots with random inter-site couplings. Classical inputs are encoded in classical controls, such as coherent pumps in the network, and quantum inputs are incident to the QR in the form of optical fields. For temporal processing, each quantum input interacts with the QR for a short duration before being replaced by another input. The time evolution of the interactions provides a high-dimensional nonlinear mapping of the input via the correlations in the QR, which can be extracted by classical or quantum readouts on accessible nodes. This enables us to learn the function of input sequence, leading to diverse applications in classical and quantum data processing.

Quantum–Classical Hybrid Information Processing via a Quantum Reservoir. When we describe a quantum device processing quantum data in a realistic scenario, we must incorporate classical control into the model. In this case, a quantum device is in fact a function of quantum input $\beta$ and classical control $u$ as $F(u,\beta )$, where we consider the scalar $u$ for ease of explanation. For a device processing the sequence of hybrid inputs $(u_1,{\beta }_1)$, $(u_2,{\beta }_2)$, …, we can describe it using the temporal map $F\left(\right\{(u_l,{\beta }_l)\left\}\right)$ of input history [33]. Our target is to develop a trainable framework to emulate $F$.

The proposed framework contains three main parts: an input part containing input modes to receive the data, a QR processor to interact with inputs in a quantum evolution, and a readout for further processing [Fig. 1(b)]. We consider the QR processor as a two-dimensional lattice of $N$ quantum dots, represented by the Hamiltonian

where ${\hat{c}}_i$, $E_i$, $h_{ij}$, $Q_i$, and $P\left(t\right)$ are the field operators, onsite energies, hopping amplitudes between the nearest neighbor sites, nonlinearity strengths, and uniform time-dependent coherent field strengths, respectively. $P\left(t\right)$ can be used to encode the classical input $u\left(t\right)$ as $P\left(t\right)=P+Wu\left(t\right)$, where $P$ and $W$ are the constant coefficient and input scaling, respectively.

The dynamics of the combined quantum state $\rho$ of the QR as well as the input modes can be described by the quantum master equation (we use the unit where Plank constant $\hslash =1$).

where $\Omega \left(t\right)=1$ for $t\ge 0$ and $0$ otherwise. Here, $\hat{A}\rho ={\sum }_k\frac{{\gamma }_k}{\gamma }L\left({\hat{a}}_k\right)\rho +{\sum }_{k,j}{W}_{jk}^{\text{in}}([{\hat{a}}_k\rho ,{\hat{c}}_j^{†}]+[{\hat{c}}_j,\rho {\hat{a}}_k^{†}])$ represents the cascade coupling between the input modes ${\hat{a}}_k$ and the QR [35]. The Lindblad superoperator $L\left(\hat{x}\right)$ is defined for any arbitrary operator $\hat{x}$ by $L\left(\hat{x}\right)\rho =\hat{x}\rho {\hat{x}}^{†}-\frac{1}{2}({\hat{x}}^{†}\hat{x}\rho +\rho {\hat{x}}^{†}\hat{x})$.

We explain quantum-classical hybrid processing using the proposed platform. First, the QR is excited only with the uniform $P$ for $0\le t<t_{\text{init}}$ and no incident quantum inputs. We choose $t_{init}$ such that the QR at time $t_{\text{init}}$ reaches a steady state. This setting ensures the echo state property [15] for the reproducible computation, where the response to the same input sequence is independent of the QR’s initial state. Then, the quantum input $\beta$ (described by the input modes ${\hat{a}}_k$) is incident to the reservoir, and the classical input $u\left(t\right)=u$ is activated at the same time. At time $t_1=t_{init}+\tau$ for time interval $\tau$, an appropriate and practical readout from the QR is performed for nontrivial transformations of input data (see Supplementary Information for detailed settings of $h_{ij}$, ${\gamma }_k$, $W_{jk}^{\text{in}}$, $\tau$, and $t_{init}$). We consider two readout schemes: a linear combination of measurement results on the accessible observables (classical readout) and the other with a linear combination of quantum modes (quantum readout). The former is associated with a measurement process, while the latter has been considered in a quantum neuromorphic platform for quantum state preparation [26].

For a non-temporal processing task, we repeat the above procedure for every hybrid data instance $(u,\beta )$. For a temporal processing task, at $t_l=t_{init}+(l-1)\tau$ ($l=1,2,\dots$), the classical input is switched to $u\left(t\right)=u_l$, and the quantum state ${\beta }_l$ replaces the partial state in the input modes. Since the input information is transferred into the QR during the interaction, this scheme enables the memory ability, which is required in temporal processing tasks.

In the classical readout, measuring the expectation values of the occupation numbers $n_j=⟨{\hat{c}}_j^{†}{\hat{c}}_j⟩$ can extract the information from the QR to reconstruct $F$. A representative application is quantum tomography, which reconstructs the density matrix output of $F$ via the linear regression model: $W^{out}n+b\approx Y_F$ [25, 33]. Here, $n=(n_1,\dots ,n_K{)}^{\top }$ is the $K$-dimensional reservoir state for readout, $Y_F$ is the real vector form to stack the real and imaginary elements of $F$, and $W^{out}$ and $b$ are the weight and bias parameters to be optimized via the training (see Methods). In the classical readout, multitasking is possible since the training cost is minimal for independent training with different $W^{out}$ for different tasks. If the measurement is performed after an interaction time $\tau$ for the current input and right before the next input, the dimensionality $K$ is equal to the number of quantum dots $N$. One can increase this dimensionality by performing measurements at different timings in the interval $\tau$, which is known as the temporal multiplexing technique. Between two inputs, we perform measurements at equal interval $\tau /V$, forming the dimensionality $K=NV$. Here, $V$ is called the measurement multiplexity. Another technique to increase the dimensionality $K$ is spatial multiplexing [20], where readout reservoir states in different QRs are combined to learn the target.

In the quantum readout, the standard toolbox of linear optical elements [36] enables us to generate $M$ quantum output modes ${\hat{C}}_m={\sum }_j{o}_{mj}{\hat{c}}_j$ with complex coefficients $o_{mj}$. The output modes must satisfy the commutation relations $[{\hat{C}}_m,{\hat{C}}_n^{†}]={\delta }_{mn}$, which impose the condition ${\sum }_j{o}_{mj}{o}_{nj}^{\ast }={\delta }_{mn}$. Since the target is the quantum state, the training process is not as simple as the one used for linear regression on the accessible observables in the classical readout. Consider the separation of non-adjustable and adjustable parameters in PRC, we assume that the parameters of Hamiltonian in Eq. (Results) are random and not trainable. Instead, we train interaction ($W_{jk}^{in}$) and readout ($\left\{o_{mj}\right\}$) coefficients such that the quantum state described via $\left\{{\hat{C}}_m\right\}$ becomes the same as the output of $F$ (see Methods).

Quantum Tomography and Channel Equalizer. We present an application of QR to hybrid tasks in which quantum tomography and noise-free reconstruction of classical data are performed simultaneously. Consider a temporal map $F\left\{\right(s_l,{\beta }_l\left)\right\}$ where $\left\{s_l\right\}$ and $\left\{{\beta }_l\right\}$ are the sequences of classical controls and quantum inputs, respectively. We assume that the output state $F_l=F\left\{\right(s_l,{\beta }_l\left)\right\}$ is accessible at $l=1,\dots ,L$ for training. The tomography task learns the relation between $F_l$ and $\left\{\right(s_l,{\beta }_l\left)\right\}$ for $l\le L$ and reconstructs $F_l$ with $l>L$. Obviously, the QR cannot learn this hybrid task without the information contained in $\left\{s_l\right\}$. Therefore, we further assume that the classical control data are also accessible, although only in a distorted form of a nonlinear transformation $s_l\to u_l$. Since multitasking is feasible in the classical readout, we can also reconstruct $\left\{s_l\right\}$ from $\left\{u_l\right\}$.

In the following example, we consider $F$ as a quantum switch with classical control [Fig. 1(a)]. Technically, a quantum switch includes two quantum channels $N_A$ and $N_B$ representing the operations by Alice and Bob, respectively, and an independent switch state ${\rho }_s$. Signal communication between Alice and Bob is only restricted to a partial order. However, the quantum switch can send the information under the indefinite causal order of quantum channels [8, 9, 10, 11, 12]. Given a state $\beta$ on which these channels act, the quantum switch produces an output $N_A\circ N_B\left(\beta \right)$ if ${\rho }_s=|0⟩⟨0|$ and $N_B\circ N_A\left(\beta \right)$ if ${\rho }_s=|1⟩⟨1|$. When the switch state is in a superposition of $|0⟩$ and $|1⟩$, such as ${\rho }_s=|{\psi }_s⟩⟨{\psi }_s|$ with ${\psi }_s=\sqrt{s}|0⟩+\sqrt{1-s}|1⟩$ ($0\le s\le 1$), the output becomes a quantum superposition of two alternative orders $N_A\circ N_B\left(\beta \right)$ and $N_B\circ N_A\left(\beta \right)$. Here, the quantum switch $S(N_A,N_B)$ can be considered a function of hybrid input $(s,\beta )$.

We use our QR to mimic the behavior of the quantum switch applied to the input sequence. Given a delay $d$, we demonstrate that the QR with current inputs ${\beta }_l$ and $u_l$ can utilize memory effects to reconstruct ${\sigma }_l=S(N_A,N_B)(s_{l-d},{\beta }_{l-d})$ and $s_{l-d}$. We consider $N_A,N_B$ as two depolarizing quantum channels and the reconstruction of $\left\{s_{l-d}\right\}$ from $\left\{u_l\right\}$ as the nonlinear channel equalization task (see Methods). Here, $\left\{{\beta }_l\right\}$ is an i.i.d. sequence of one-qubit density matrices, and $\left\{s_l\right\}$ is an i.i.d. discrete sequence of symbols, which are selected from $\{-3,-1,1,3\}$ with equal probability. The switch state at each $l$ is ${\rho }_{(3+s_l)/6}$, and the distorted input $\left\{u_l\right\}$ is transformed from $\left\{s_l\right\}$ via both linear and nonlinear channels [37] (see Methods). If $d\ge 1$, it requires a QR with the nonlinear effect and a memory of both quantum and classical inputs.

The QR’s output is divided into two parts: the tomography result $Y_l$ in the real vector form and the equalized result $y_l$. $Y_l$ is then transformed in the density matrix form ${\hat{\sigma }}_l$ with the consideration of a projection to obtain a positive semidefinite matrix (see Methods). $y_l$ is converted back into a nearest symbol ${\hat{s}}_l\in \{-3,-1,1,3\}$. The training is performed at $l=1,\dots ,L$ ($L=800$), and the tomography performance is evaluated via the root mean square of fidelities (RMSF)

where $T=200$ and $F(\rho ,\sigma )=Tr\left[\sqrt{\sqrt{\sigma }\rho \sqrt{\sigma }}\right]$. The equalization performance is evaluated via the symbol error rate (SER)

Figure 2(a) illustrates a sequence of one-qubit quantum input in the evaluation phase (upper panel) and a result for the channel equalizer (bottom panel) with delay $d=1$ (see Supplementary Information for results with other $d$ values). Here, the predict and target sequences for the reconstruction of classical symbols $\left\{s_l\right\}$ are overlapped at almost all time steps. The density matrix at each time step is represented as a real vector by stacking the real and imaginary parts. Figure 2(b) depicts that the quantum target sequence can be reconstructed well.

We systematically evaluate the performance of the tomography and channel equalizer tasks via the RMSF (left axis) and SER (right axis) in Fig. 3(a) for different $N$ and $W$. A large value of $W$ compared with $h_{ij}$ and $W_{jk}^{\text{in}}$ leads to non-ergodic behavior in the QR, i.e., a strong and qualitative dependence on the initial state at $t_{init}$ (Fig. S1 in Supplementary Information). In addition, in Supplementary Information, we further investigate the effects of the classical input in the reconstruction of the quantum input. With a large $W$, the input state is incident with weak coupling ($|W_{jk}^{in}|\ll |P\left(t\right)=P+Wu\left(t\right)|$) under a strong effect of the classical input to the QR’s dynamics, which means that not much information regarding quantum inputs can be retained in the QR. In contrast, a small $W$ reduces the memory effect in reconstructing the previous classical input. This explains the existence of a region of $W$ for an optimal performance ($W/\gamma \approx 1.0$).

The left panel of Fig. 3(b) displays the RMSF of the tomography task when we increase the number $N$ of reservoir sites. In the right panel of Fig. 3(b), we further compare the performance in the equalization task with the Echo State Network (ESN) in classical reservoir computing under the condition of the same number of computational nodes (see Methods). Here, we set the input scaling as $W/\gamma =1.0$ and use the QR with the measurement multiplexity $V=8$; therefore, the QR containing $16,24,32,40,48$ computational nodes corresponds to $N=2,3,4,5,6$ sites in the reservoir. We confirm that with appropriate setting of the constant coherent field $P$, we can obtain almost the same performance with the ESN.

Continuous Variable Tomography and Closed Loop. We modify the situation in the tomography task where, after the training phase, we were unable to access the information from the classical control $s_l$. Surprisingly, owing to the advantages of multitasking, our QR can autonomously generate $s_l$ in a closed-loop manner while performing the tomography task with the hybrid input. In the training phase, $s_l$ is learned in an open loop where we predict the next step $s_{l+1}$ given the input $u_l=s_l$. After training, the prediction is used as the classical input for the next step, forming a closed-loop control without any external interventions. This model-free prediction is well established in classical reservoir computing, for example, to predict low-dimensional chaotic systems [37] or large spatiotemporally chaotic systems [38]. However, to the best of our knowledge, our demonstration is the first to combine the closed-loop setting with the quantum tomography task, which is only effective in the QR setting.

We consider the quantum tomography of continuous variable states. The target is to reconstruct the output $F_l=F\left\{\right(s_l,{\beta }_l\left)\right\}$ in the Wigner function form $W(F_l;x_i,p_j)$ defined on a grid of continuous variables $x_i$ and $p_j$ (see Methods). We use 300 randomly generated one-mode thermal states ${\beta }_l$ and the periodic signals $s_l=0.5+0.5\sin \left(\frac{l\pi f}{510}\right)$ in the training phase. The target $F_l$ is created by applying one-mode squeezing operator to ${\beta }_l$ as

where ${\xi }_l=s_l{e}^{i\pi /4}$. In Supplementary Information, we consider another encoding: ${\xi }_l=0.3e^{i2\pi s_l}$. Here, we consider the cutoff Fock space dimension (the effective dimension) of these continuous variables states is $D_{eff}=9$.

Figure 4(a) shows examples of the control signals in the training and closed-loop phase for $f=60$. With $W/\gamma =0.8$ and $N=3$ sites, the control signal is almost reconstructed perfectly for all time steps in the closed-loop phase. This QR can efficiently reconstruct the Wigner function even without accessing the control signal [Fig. 4(b)]. We further investigate the stability of the closed-loop trajectories plotted in the $(s_l,s_{l+10})$ plane [Fig. 4(c)]. The QR presents a stable embedding of sinusoidal classical inputs if the trajectory can return to the target after adding a small perturbation (green line) into a predicted value, suggesting that our system successfully learned the target attractor. We observe an appropriate setting of input scaling $W$ to obtain stable closed loops ($W/\gamma \approx 0.8)$. Intriguingly, if we increase $W/\gamma$, for example to $W/\gamma =1.8$, the closed loop fails to reconstruct the trajectory of the sinusoidal input in the evaluation stage but can produce chaotic-like behavior in the embedding space. In this case, the generated trajectory is not elliptical as the trajectory of sinusoidal inputs but still robust with respect to a small perturbation. We also observe the dependency of the performance of closed-loop control and the production of chaotic-like behavior on time scales $f$ of the control signals, which are investigated in detailed in Supplementary Information.

Quantum Readout and Depolarizing Channel. Finally, we present an application using the quantum readout scheme to output quantum states. We use the QR to prepare a depolarizing quantum channel $F\left\{\right(s_l,{\beta }_l\left)\right\}=s_{l}I/D+(1-s_l){\beta }_l$, where $\left\{{\beta }_l\right\}$ are randomly generated in a $D$-dimensional Hilbert space and $\left\{s_l\right\}$ is a random sequence in $[0,1]$.

First, we consider a sequence of $200$ one-qubit quantum states for the training and $100$ states for the evaluation. The baseline is computed when we set the output as the same as the input. We use the Nelder–Mead simplex algorithm [39] (see Methods) to minimize the fidelity error

where ${\sigma }_l$ and ${\hat{\sigma }}_l$ are the target and preparing quantum states, respectively. In Fig. 5(a), the evaluated fidelity errors with different readout and training configurations are presented for the QR with $N=2$ sites, $P=1.0$, and $W=2.0$. The interquartile range is contained within the box, and the 5th and 95th percentiles are marked by whiskers. The median is the line across the box, and the outliers are located outside the whiskers of each box plot. Here, IN, RV, and ALL correspond to the setting where only input modes ${\hat{a}}_k$, only reservoir modes ${\hat{c}}_j$, or both of them are considered as the readout nodes, respectively. Wo and Wio correspond to the situation where only readout weights or both readout weights and interaction coefficients $W_{jk}^{in}$ are considered as the training parameters, respectively. The result implies that the consideration of both input and reservoir modes as $N_R$ readout modes and both interaction coefficients and readout weights for training leads to the best performance. Under this setting, we display the variation in fidelity errors EF with the input scaling $W/\gamma$ and $N_R$ in Fig. 5(b). Even with a small QR ($N_R=3,4$) we can prepare the target channel with a relatively low error ($<2\%$), which is significantly better compared with the baseline ($\approx 8\%$). Furthermore, increasing $W/\gamma$ basically leads to a better performance where more information regarding the classical input is integrated.

Finally, we prepare the depolarizing channel using the input quantum states as random squeezed thermal states in the continuous variable form. We minimize the cost function taken in the Wigner representation as follows:

Owing to the scale limitation, we only simulate the continuous variable states of the effective dimension $D_{eff}=3$, where $D=D_{eff}^2=9$. Figure 5(c) presents the errors in 50 training and 50 evaluating data varying with $W$. We can observe a similar trend in Fig. 5(b), that is, with sufficient classical information ($W/\gamma \ge 1.0$), the error EW ($\approx 0.16$) with $N_R=3$ readout nodes is significantly lower than the baseline’s error ($\approx 0.29$). This result is still below a considerably good preparation ($EW<0.1$), but it demonstrates that hybrid inputs can be effectively considered for training the quantum readout.

We proposed a framework for an analog QR processor with hybrid inputs and classical and quantum readouts for learning heterogeneous quantum-classical data. This aligns well with scenarios where one wishes to model a quantum device to process quantum input but must rely on classical control signals in physical experiments. Our framework, therefore, has the potential to be physically implemented in quantum network systems where classical control and quantum sources can interact with nonlinear quantum systems to form a quantum channel. It can help realize quantum adaptive systems capable of quantum information processing. These agents can be used to interpret and memorize both classical and quantum signals from their environment and to respond accordingly to the actions of their surroundings [40].

Processing hybrid quantum and classical data is a promising idea to facilitate future innovative use cases for quantum computers. This concept aims to leverage the advantages of quantum mechanics in ML with an unconventional computing framework and intriguing applications. It is not limited to the conventional discussion on practical quantum advantages, such as the “beating speedup” of quantum to classical ML methods [41]. For example, classical readouts lead to interesting applications of multitasking where quantum data can be processed in a closed loop of the classical control. Furthermore, adding this closed-loop mechanism allows us to utilize the unique coherence properties of quantum systems to generate unique classical dynamics. We consider the quantum readout to avoid the measurement process of preparing the quantum output. However, optimization can be challenging and requires improvement, since we need to simulate or drive the quantum system and evaluate the cost function for a wide range of parameters.

A further enticing discussion would be the case of the correlation between the processing of quantum data and classical data in a hybrid setting of the QR. We can consider a QR to simultaneously process quantum data and classical data as separate tasks. An intriguing research question arises: Can this multitasking mechanism induce positive or negative effects on information processing? For example, if we repeatedly modify the coherent field strengths of the QR via a classical input with a large magnitude, it can limit the short-term memory properties of quantum data processing (see Supplementary Information). However, one can also expect positive effects and not only negative ones. There may exist a situation where simultaneously processing different modals of data can actually bring an optimal regime rather than solely solving a single task. We can start by investigating relations between hybrid input protocols with the dynamics of the QR, such as the classical input may induce the dynamical phase transition in the QR [42]. We can also study how classical and quantum data are processed via the QR’s dynamics, such as by decomposing the readout reservoir states in terms of basis polynomials for input history [43, 44]. Along with this research line, one can refer to a recent study demonstrating that quantum noise in real quantum processors can induce the information processing capability when using classical data [45].

Reservoir Computing (RC). RC is based on the modeling for the transformation of the input sequence in high-dimensional trajectories [15, 16, 17, 18]. RC contains three main parts: the input part to store the input sequence, the reservoir, and the readout part. A reservoir is a dynamical system driven by an input to encode recurrent relations and nonlinear dynamics of the input history. The readout part with connections to the reservoir can extract reservoir states, which are useful features retained by the reservoir for emulating the target sequence. In general, we only need to train the readout connections, making RC particularly suitable for physical implementations.

Mathematically, RC is described by the input-driven map $g:U\times X\to X\subset R^K$, where $U$ and $X$ are the input and the reservoir’s state space, respectively. Here, $K$ is considered the dimension of the reservoir’s state. If we feed a discrete-time input sequence $\{\dots ,u_{-1},u_0,u_1,\dots \}$ into the reservoir, the readout reservoir state $x_l$ is represented by the following recurrent relation:

In temporal supervised learning tasks, we are given an input sequence $\{u_1,\dots ,u_L\}$ and the corresponding target sequence $\hat{y}=\{{\hat{y}}_1,\dots ,{\hat{y}}_L\}$, where ${\hat{y}}_k\in R^d$ with $d$ is the output dimension. We consider a parameterized readout map $h_w:X\to R^d$, where the output signal is $y_l=h_w\left(x_l\right)$. The readout map is often taken as a linear combination of the readout reservoir states as $y_l=h_w\left(x_l\right)=w^{\top }{x}_l$. Here, $w$ is the trainable parameter obtained by minimizing the error between $y_l$ and $\widehat{y_l}$ over $l=1,\dots ,L$. For example, we often consider the mean square error

where $\parallel \cdot {\parallel }_2$ denotes the Euclidean norm between two vectors in $R^d$. For training, we add a constant bias term $x_{l,K+1}=1$ to the readout reservoir state $x_l$ and optimize $w$ via the linear regression $\hat{Y}=Xw$. Here, $\hat{Y}=[{\hat{y}}_1\dots {\hat{y}}_L{]}^{\top }$ is the $L\times d$ target matrix and $X$ is the $L\times (K+1)$ matrix that combines the readout reservoir states $x_1,x_2,\dots ,x_L$ of the training data. The optimal value of $w$ is obtained via the Ridge regression in the matrix form ${\hat{w}}^{\top }=(X^{\top }X+\eta I{)}^{-1}{X}^{\top }\hat{Y}$, where $\eta$ is a positive constant for the regularization.

Echo State Network (ESN). ESN is a realization of the input-driven map in RC. It belongs to the concept of artificial recurrent neural network (RNN), as we have a large network with randomly fixed inner and recurrent connections. Consider ESN with $N$ computational nodes and a discrete-time input sequence $\left\{u_l\right\}$, the reservoir state $x_l$ at time step $l$ is described by

where $\text{tanh}(\cdot )$ is the activation function applied on vector $x=(x_1,x_2,\dots ,x_K{)}^{\top }$ as $\text{tanh}\left(x\right)=\left(\tanh \right(x_1),\dots ,\tanh (x_K){)}^{\top }$. Here, $W^{in}$ and $W$ are the input weight matrix and recurrent weight matrix, respectively. In the channel equalization task, the input weight matrix $W^{in}$ is generated from random uniform distribution in $[-1,1]$. We also set the connection probability and the spectral radius of the recurrent weight matrix $W$ to 0.1 and 0.9, respectively.

Learning Quantum Tomography. A quantum device can be described by a function of quantum input $\beta$ and classical control $u$ as $F(u,\beta )$, where we consider the scalar $u$ for ease of notation. Given a sequence of hybrid inputs $(u_1,{\beta }_1)$, $(u_2,{\beta }_2)$, …and a quantum device with a time-dependent behavior, we can describe it using the temporal map $F_l=F\left(\right\{(u_i,{\beta }_i){\}}_{i=1:l})$ of the current and past inputs [33]. We assume that we have full tomography for the corresponding output states of $F$ in the training, where we are given a hybrid input sequence $\left\{\right(u_1,{\beta }_1),\dots ,(u_L,{\beta }_L\left)\right\}$ and the corresponding target sequence $\hat{y}=\{{\hat{y}}_1,\dots ,{\hat{y}}_L\}$. Here, ${\hat{y}}_l$ is the real vector form of $F_l$. If $F_l$ is described by the density matrix, ${\hat{y}}_l$ is formed by stacking the real and imaginary elements of $F\left({\beta }_l\right)$. In the evaluation stage, we are given an input sequence $\left\{\right(u_{L+1},{\beta }_{L+1}),\dots ,(u_{L+T},{\beta }_{L+T}\left)\right\}$ with the corresponding target $\{{\hat{\sigma }}_{L+1},\dots ,{\hat{\sigma }}_{L+T}\}$, where ${\hat{\sigma }}_i=F(u_i,{\beta }_i)$. The output sequence is $\{y_{L+1},\dots ,y_{L+T}\}$, which is rearranged in the matrix form $\{{{\sigma }^{\prime }}_{L+1},\dots ,{{\sigma }^{\prime }}_{L+T}\}$. To obtain the final positive semidefinite matrix ${\sigma }_i$, we project ${\sigma }_i^{\prime }$ onto the spectrahedron such that the trace of ${\sigma }_i$ is equal to 1 and the Frobenius norm of ${\sigma }_i-{\sigma }_i^{\prime }$ is minimized [46, 47].

Tomography learning can be performed with other forms of $F_l$, for example, in the Wigner function representation of continuous variable states. Given a density matrix $\sigma$, the continuous variable states associated with $\sigma$ can be described by the Wigner function

where the integral is evaluated in the whole space and the states $|x_i\pm \frac{y}{2}⟩$ represent continuous position bases. We evaluate Wigner functions on a $61\times 61$ grid of $x_i$ and $p_j$, where we divide the interval $[-3,3]$ into $60$ equal intervals for the range of $x_i$ and $p_j$. The target of continuous variable tomography is to reconstruct these Wigner functions, i.e., the real $61\times 61$ dimensional matrices.

Quantum Switch. In the classical counterpart, a switch is an operation of control that can route a target system through two classical channels $C_A$ and $C_B$ in series following one causal order ($C_A$ then $C_B$) or the other ($C_B$ then $C_A$). The quantum switch is different in that it induces entirely new quantum trajectories where the order of the two operators is indefinite. Technically, a quantum switch includes two quantum channels $N_A$ and $N_B$ to create a new channel $S(N_A,N_B)$, which uses the channels $N_A$ and $N_B$ in an order that is entangled with an independent switch quantum state ${\rho }_s$. The channel $S(N_A,N_B)$ returns the state $(N_A\circ N_B(\rho \left)\right)\otimes |0⟩⟨0|$ if ${\rho }_s=|0⟩⟨0|$ and $(N_B\circ N_A(\rho \left)\right)\otimes |1⟩⟨1|$ if ${\rho }_s=|1⟩⟨1|$. When ${\rho }_s$ is in a superposition of $|0⟩$ and $|1⟩$, the channel returns a correlated state as a result of $N_A$ and $N_B$ acting on $\rho$ in a superposition of two alternative orders.

To describe $S(N_A,N_B)$, we denote the Kraus operators of channels $N_A$ and $N_B$ as $\left\{K_i^{\left(A\right)}\right\}$ and $\left\{K_j^{\left(B\right)}\right\}$, respectively, where $N_A={\sum }_i{K}_i^{\left(A\right)}\rho K_i^{\left(1\right)†}$ and $N_B={\sum }_j{K}_j^{\left(B\right)}\rho K_j^{\left(2\right)†}$. The Kraus operators of $S(N_A,N_B)$ are defined as

The action of the quantum switch is given by

In our study, we consider $N_A$ and $N_B$ as two depolarizing channels with parameters $q_A$ and $q_B$, which are given by

where $D\times D$ is the dimension of $\rho$ and $\{U_i{\}}_{i=1}^{D^2}$, $\{V_j{\}}_{j=1}^{D^2}$ are orthonormal bases of the space of $D\times D$ matrices. Here, we introduce the notation $U_0=\frac{D\sqrt{1-q_A}}{\sqrt{q_A}}I$ and $V_0=\frac{D\sqrt{1-q_B}}{\sqrt{q_B}}I$. We define the extension Kraus operators for $N_A$ and $N_B$ as $K_i^{\left(A\right)}=\frac{\sqrt{q_A}}{D}{U}_i$ for $i=0,1,\dots ,D^2$ and $K_j^{\left(B\right)}=\frac{\sqrt{q_B}}{D}{V}_j$ for $j=0,1,\dots ,D^2$, respectively. We can express the Kraus operators of $S(N_A,N_B)$ as