ℤ2 × ℤ2 Equivariant Quantum Neural Networks: Benchmarking against Classical Neural Networks

Abstract

This paper presents a comparative analysis of the performance of Equivariant Quantum Neural Networks (EQNNs) and Quantum Neural Networks (QNNs), juxtaposed against their classical counterparts: Equivariant Neural Networks (ENNs) and Deep Neural Networks (DNNs). We evaluate the performance of each network with three two-dimensional toy examples for a binary classification task, focusing on model complexity (measured by the number of parameters) and the size of the training dataset. Our results show that the ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ EQNN and the QNN provide superior performance for smaller parameter sets and modest training data samples.

1. Introduction

The rapidly evolving convergence of machine learning (ML) and high-energy physics (HEP) offers a range of opportunities and challenges for the HEP community. Beyond simply applying traditional ML methods to HEP issues, a fresh cohort of experts skilled in both areas is pioneering innovative and potentially groundbreaking approaches. ML methods based on symmetries play a crucial role in improving data analysis as well as expediting the discovery of new physics [1,2]. In particular, classical Equivariant Neural Networks (ENNs) exploit the underlying symmetry structure of the data, ensuring that the input and output transform consistently under the symmetry [3]. ENNs have been widely used in various applications including deep convolutional neural networks for computer vision [4], AlphaFold for protein structure prediction [5], Lorentz equivariant neural networks for particle physics [6], and many other HEP applications [7,8,9,10,11].

Meanwhile, the rise of readily available noisy intermediate-scale quantum computers [12] has sparked considerable interest in using quantum algorithms to tackle high-energy physics problems. Modern quantum computers boast impressive quantum volume and are capable of executing highly complex computations, driving a collaborative effort within the community [13,14] to explore their applications in quantum physics, particularly in addressing theoretical challenges in particle physics. Recent research on quantum algorithms for particle physics at the Large Hadron Collider (LHC) covers a range of tasks, including the evaluation of Feynman loop integrals [15], simulation of parton showers [16] and structure [17], development of quantum algorithms for helicity amplitude assessments [18], and simulation of quantum field theories [19,20,21,22,23,24].

An intriguing prospect in this realm is the emerging field of quantum machine learning (QML), which harnesses the computational capabilities of quantum devices for machine learning tasks. With classical machine learning algorithms already proving effective for various applications at the LHC, it is very natural to explore whether QML can enhance these classical approaches [25,26,27,28,29,30,31,32,33,34]. In recent years, significant development has been made in their quantum counterparts, Equivariant Quantum Neural Networks (EQNNs) [35,36,37,38,39].

In this paper we benchmark the performance of EQNNs against various classical and/or non-equivariant alternatives for three two-dimensional toy datasets, which exhibit a ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ symmetry structure. Such patterns often appear in high-energy physics data, e.g., as kinematic boundaries in the high-dimensional phase space describing the final state [40,41]. By a clever choice of the kinematic variables for the analysis, these boundaries can be preserved in projections onto a lower-dimensional feature space [42,43,44,45,46]. For example, one can form various combinations of possible invariant mass for the generic decay chain considered in Ref. [44], $D\to jC\to j{\ell }_n^{\pm }B\to j{\ell }_n^{\pm }{\ell }_f^{\mp }A$, where Particles A, B, C, and D are hypothetical particles in new physics beyond the standard model of masses $\{m_A,m_B,m_C,m_D\}$, while the corresponding standard model decay products consist of a jet j, a “near” lepton ${\ell }_n^{\pm }$, and a “far” lepton ${\ell }_f^{\pm }$. The two-dimensional (bivariate) distribution $\frac{d^2\Gamma }{d{R}_{ij}d{R}_{kl}}$ shows distributions similar to those in Figure 1, Figure 2 and Figure 3, where $R_{ij}=\frac{m_i^2}{m_j^2}$ is the mass square ratio. Symmetric, anti-symmetric, or non-symmetric structures provide information of particle masses involved in the cascade decays.

In this study, we consider simplified two-dimensional datasets that mimic the data arising in such projections. This setup allows us to focus on the comparison between different methods, avoiding unnecessary issues that may arise when dealing with actual particle physics simulation data such as sampling statistics, parton distribution functions, unknown particle mass spectrum, unknown width, detector effects, etc. We explore EQNNs and benchmark them against classical neural network models. We find that the variational quantum circuits learn the data better with the smaller number of parameters and the smaller training dataset compared to their classical counterparts.

2. Dataset Description

In all three examples, we consider two-dimensional data $(x_1,x_2)$ on the unit square ($-1\le x_i\le 1$). The data points belong to two classes: $y=+1$ (blue points) and $y=-1$ (red points).

• (i). Symmetric case:

  In the first example (Figure 1), the labels are generated by the function

  (1)$\begin{array}{ccc}\hfill y(x_1,x_2)& =& 2H\left(R-\sqrt{{(x_1+1)}^2+{(x_2-1)}^2}\right)\hfill \\ & +& 2H\left(R-\sqrt{{(x_1-1)}^2+{(x_2+1)}^2}\right)-1,\hfill \end{array}$

  where $H\left(x\right)$ is the Heaviside step function and for definiteness we choose $R=1.1$. The function (1) respects a ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ symmetry, where the first ${\mathbb{Z}}_2$ is given by a reflection about the $x_1=x_2$ diagonal

  (2)$x_1\to x_2,x_2\to x_1,y\to y,$

  while the second ${\mathbb{Z}}_2$ corresponds to a reflection about the $x_1=-x_2$ diagonal

  (3)$x_1\to -x_2,x_2\to -x_1,y\to y.$

  This ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ example was studied in Ref. [37] and we shall refer to it as the symmetric case since the y label is invariant.

• (ii). Anti-symmetric case:

  The second example is illustrated in Figure 2. The labels are generated by the function

  (4)$\begin{array}{cc}\hfill y(x_1,x_2)& =H\left(-x_1\right)H\left(-x_2\right)+H\left(-x_1\right)H\left(x_2\right)H(x_1+x_2)\hfill \\ \hfill & -H\left(x_1\right)H\left(x_2\right)+H\left(x_1\right)H\left(-x_2\right)H(x_1+x_2).\hfill \end{array}$

  The first ${\mathbb{Z}}_2$ is still realized as in (2). However, this time, the labels are flipped under a reflection along the $x_1=-x_2$ diagonal:

  (5)$x_1\to -x_2,x_2\to -x_1,y\to -y,$

  which is why we shall refer to this case as anti-symmetric.

• (iii). Fully anti-symmetric case:

  The last example is depicted in Figure 3. The labels are generated by the function

  (6)$\begin{array}{cc}\hfill y(x_1,x_2)& =H\left(x_1\right)H\left(x_2\right)(2H\left(x_1-x_2\right)-1)+H\left(-x_1\right)H\left(-x_2\right)(2H\left(x_2-x_1\right)-1)\hfill \\ \hfill & +H\left(-x_1\right)H\left(x_2\right)H\left(x_1+x_2\right)(2H\left(R-\sqrt{{(x_1+1)}^2+{(x_2-1)}^2}\right)-1)\hfill \\ \hfill & +H\left(-x_1\right)H\left(x_2\right)H\left(-x_1-x_2\right)(1-2H\left(R-\sqrt{{(x_1+1)}^2+{(x_2-1)}^2}\right))\hfill \\ \hfill & +H\left(x_1\right)H\left(-x_2\right)H\left(x_1+x_2\right)(1-2H\left(R-\sqrt{{(x_1-1)}^2+{(x_2+1)}^2}\right))\hfill \\ \hfill & +H\left(x_1\right)H\left(-x_2\right)H\left(-x_1-x_2\right)(2H\left(R-\sqrt{{(x_1-1)}^2+{(x_2+1)}^2}\right)-1),\hfill \end{array}$

  where $H\left(x\right)$ is the Heaviside step function and for definiteness we choose $R=1$. In this case, the labels are flipped under both reflections along the $x_1=-x_2$ diagonal as well as the $x_1=x_2$ diagonal, which is why we shall refer to this case as fully anti-symmetric. As we will see later, it is straightforward to incorporate both symmetric and anti-symmetric properties in variational quantum circuits, while it is not obvious how to consider the anti-symmetric case in the classical neural networks.

3. Network Architectures

To assess the importance of embedding the symmetry in the network, and to compare the classical and quantum versions of the networks, we study the performance of the following four different architectures: (i) Deep Neural Network (DNN), (ii) Equivariant Neural Network (ENN), (iii) Quantum Neural Network (QNN), and (iv) Equivariant Quantum Neural Network (EQNN). In each case, we adjust the hyperparameters to ensure that the number of network parameters is roughly the same.

• (i). Deep Neural Networks:

  In our DNN, for the symmetric (anti-symmetric) case, we use one (two) hidden layer(s) with four neurons. For both types of classical networks, we use the softmax activation function, Adam optimizer, and a learning rate of $0.1$. We use the binary cross-entropy for both the DNN and ENN.

• (ii). Equivariant Neural Networks:

  A given map $f:x\in X\to f\left(x\right)\in Y$ between an input space X and an output space Y is said to be equivariant under a group G if it satisfies the following relation:

  (7)$f\left(g_{\mathrm{in}}\left(x\right)\right)=g_{\mathrm{out}}\left(f\left(x\right)\right),$

  where $g_{\mathrm{in}}$ ($g_{\mathrm{out}}$) is a representation of a group element $g\in G$ acting on the input (output) space. In the special case when $g_{\mathrm{out}}$ is the trivial representation, the map is called invariant under the group G, i.e., a symmetry transformation acting on the input data x does not change the output of the map. The goal of ENNs, or equivariant learning models in general, is to design a trainable map f which would always satisfy Equation (7). In tasks where the symmetry is known, such equivariant models are believed to have an advantage in terms of the number of parameters and training complexity. Several studies in high-energy physics have attempted to use classical equivariant neural networks [6,47,48,49,50]. Our ENN model utilizes four ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ symmetric copies for each data point, which are fed into the input layer, followed by one equivariant layer with three (two) neurons and one dense layer with four (four) neurons in the symmetric (anti-symmetric) case.

• (iii). Quantum Neural Networks:

  For the QNN, we utilize the one-qubit data-reuploading model [51], as shown in Figure 4, with depth four (eight) for the symmetric (anti-symmetric and fully anti-symmetric) case, using the angle embedding and three parameters at each depth. This choice leads to a similar number of parameters as in the classical networks. We use the Adam optimizer and the loss

  (8)$L_{QNN}=y(1-|⟨\psi |O_1{|\psi ⟩\left|\right)}^2+(1-y)(1-|⟨\psi |O_2|\psi ⟩{\left|\right)}^2$

  for any choice of two orthogonal operators $O_1$ and $O_2$ (see Ref. [52] for more details.). In this paper, we use

  (9)$O_1=\left(\begin{array}{cc}\hfill 1& 0\\ \hfill 0& 0\end{array}\right),O_2=\left(\begin{array}{cc}\hfill 0& 0\\ \hfill 0& 1\end{array}\right).$

  for all three datasets considered in this paper.

• (iv). Equivariant Quantum Neural Networks.

  In EQNN models, symmetry transformations acting on the embedding space of input features are realized as finite-dimensional unitary transformations $U_g$, $g\in G$. Consider the simplest case where one trainable operator $U(\theta ,x)$ acts on a state $|\psi ⟩$: $U(\theta ,x)|\psi ⟩$. If for a symmetry transformation $U_g$, the condition

  (10)$U(\theta ,x)U_g|\psi ⟩=U_{g}U(\theta ,x)|\psi ⟩,$

  is satisfied, then the operator U is equivariant, i.e., the equivariant gate should commute with the symmetry. In general, the $U_g$ operators on the two sides of Equation (10) do not necessarily have to be in the same representation but are often assumed so for simplicity. The output of a QNN is the measurement of the expectation value of the state with respect to some observable O. If the gates are equivariant and we apply some symmetry transformation $U_g$, then this is equivalent to measuring the observable $U_g^{†}O{U}_g$. Hence, if O commutes with the symmetry $U_g$, the model as a whole would be invariant under $U_g$, which is the case in our symmetric example. Otherwise the model is equivariant, as in our anti-symmetric example.

Our EQNN uses the two-qubit quantum circuit depicted in Figure 5 for depth 1. This circuit is repeated five (ten) times with different parameters for the symmetric (anti-symmetric and fully anti-symmetric). The two $R_Z$ gates embed $x_1$ and $x_2$, respectively. The $R_X$ gates share the same parameter (${\theta }_1$) and the $R_{ZZ}$ gate uses another parameter (${\theta }_2$). The invariant model (for the symmetric case) uses the same observable O for both classes in the data. In the anti-symmetric case, we use two different observables $O_1$ and $O_2$ that correspond to each label. They transform into one another under reflection $g_r$, i.e., $U_{g_r}^{†}{O}_1{U}_{g_r}=O_2$.

In the symmetric case, we use binary cross-entropy loss, assuming the true label y is either 0 or 1,

(11)$L_{EQNN}^{\mathrm{symm}}=ylog\left(\left|⟨\psi |O|\psi ⟩\right|\right)+(1-y)log\left(1-\left|⟨\psi |O|\psi ⟩\right|\right).$

The observables O and the reflection $U_{g_r}$ along $x_1=-x_2$ are defined as follows:

(12)$O={\displaystyle \frac{1}{4}}\left(\begin{array}{cccc}1& 1& 1& 1\\ 1& 1& 1& 1\\ 1& 1& 1& 1\\ 1& 1& 1& 1\end{array}\right),U_{g_r}=\left(\begin{array}{cccc}0& 0& 0& 1\\ 0& 0& 1& 0\\ 0& 1& 0& 0\\ 1& 0& 0& 0\end{array}\right).$

In the anti-symmetric and fully anti-symmetric cases, we used the same loss as in QNN

(13) $L_{EQNN}^{\mathrm{anti}-\mathrm{symm}}=y(1-|⟨\psi |O_1{|\psi ⟩\left|\right)}^2+(1-y)(1-|⟨\psi |O_2|\psi ⟩{\left|\right)}^2.$

For the anti-symmetric case $O_1$ ($O_2$) is the observable corresponding to $y=1$ ($y=0$)

(14)$O_1={\displaystyle \frac{1}{4}}\left(\begin{array}{cccc}\hfill 1& 1& 1\hfill & -1\\ \hfill 1& 1& 1\hfill & -1\\ \hfill 1& 1& 1\hfill & -1\\ \hfill -1& -1& -1\hfill & 1\end{array}\right),O_2={\displaystyle \frac{1}{4}}\left(\begin{array}{cccc}\hfill 1& -1& -1\hfill & -1\\ \hfill -1& 1& 1\hfill & 1\\ \hfill -1& 1& 1\hfill & 1\\ \hfill -1& 1& 1\hfill & 1\end{array}\right).$

For the fully anti-symmetric case, we use another set of observables, so one will transform into the other with reflection along any of the two diagonals. They are given as follows:

(15)$O_1={\displaystyle \frac{1}{4}}\left(\begin{array}{cccc}\hfill 1& -1& 1\hfill & 1\\ \hfill -1& 1& -1\hfill & 1\\ \hfill 1& -1& 1\hfill & -1\\ \hfill 1& 1& -1\hfill & 1\end{array}\right),O_2={\displaystyle \frac{1}{4}}\left(\begin{array}{cccc}\hfill 1& 1& -1\hfill & 1\\ \hfill 1& 1& -1\hfill & -1\\ \hfill -1& -1& 1\hfill & 1\\ \hfill 1& -1& 1\hfill & 1\end{array}\right).$

Since it is anti-symmetric with respect to each of the diagonals, the result is invariant if both reflections are applied. It is difficult to build a classical equivariant neural network using these anti-symmetries since classical equivariant models are built based on the assumption that the target is invariant under certain transformations. When discussing the theory of classical equivariant machine learning models, the models that transform non-trivially under the symmetry group are often discussed mathematically but rarely implemented in code. For our classical model on partially anti-symmetric data, we only implemented the invariant part of the symmetry (${\mathbb{Z}}_2$) and ignored the anti-symmetric portion of the data. While it may not be impossible to consider such asymmetric cases in classical neural networks, implementation can be quite involved.

On the other hand, it is straightforward to build quantum equivariant models. For this purpose, we would only need to exploit the transformation properties of the observables. If one observable transforms to the other under the transformation of interest (reflection along the diagonal in this case), then measurement made on one observable is equivalent to the measurement of the other observable given the transformed input.

We can consider equivariant quantum models with anti-symmetric transformation from the point of view of representation theory. The fully invariant (symmetric) case can be considered as the model transform under the trivial representation of the group, where all the transformations defined by the group do not change the output of the model. The asymmetric (either anti-symmetric or fully anti-symmetric) cases that we considered here can be interpreted as transforms under some other (one-dimensional) representation of the group, where some transformations change the output of the model to its opposite value, while other transformations do not change the output.

4. Results

The left panels in Figure 6 show the receiver operating characteristic (ROC) curves for each network with $N_{\mathrm{train}}=200$ and $N_{\mathrm{test}}=2000$ samples for the symmetric (top), anti-symmetric (middle), and fully anti-symmetric (bottom) dataset. The results for the DNN, ENN, QNN, and EQNN are shown in (green, dotted), (yellow, dotdashed), (red, dashed), and (blue, solid), respectively. As expected, networks with an equivariance structure (EQNN and ENN) improve the performance of the corresponding networks (QNN and DNN) without the symmetry. We also observe that quantum networks perform better than the classical analogs. In the legends, numerical values followed by network acronyms represent the number of parameters used for each network. For the symmetric example, the EQNN uses only 10 parameters; thus, for fair comparison, we constructed the other networks with $\mathcal{O}$(10) parameters as well. For the anti-symmetric example, we use 20 parameters for the EQNN.

The evolution of the accuracy during training and testing is shown in the right panels of Figure 6. The accuracy converges faster (after only 5 epochs) for the QNN and EQNN in comparison to their classical counterparts (10–20 epochs). The same color-scheme is used, but this time, solid curves represent training accuracy, while dashed curves show test accuracy.

To further quantify the performance of our quantum networks, in Figure 7, we show the AUC (Area under the ROC Curve) as a function of the number of parameters (left panels) with a fixed size of the training data ($N_{\mathrm{train}}=200$), and as a function of the number of training samples (right panels) with a fixed number of parameters ($N_{\mathrm{params}}=20$). The top, middle, and bottom panels show results for the symmetric, anti-symmetric, and fully anti-symmetric dataset. As the number of parameters increases, the performance of all networks improves. All AUC values become similar when $N_{\mathrm{params}}\approx 20$ ($N_{\mathrm{params}}\approx 40$) for the symmetric (anti-symmetric) case. As shown in the bottom panels, the performances of all networks become comparable to each other for both examples once the size of the training data reaches ∼400, except for the fully anti-symmetric case. We observe that from the top panel to the bottom panel, the relative improvement from the QNN to the EQNN grows, indicating the importance of symmetry implementation on the network. Similar relative improvement exists from the DNN to the QNN, emphasizing the importance of quantum algorithms. Note that the ENN curves are missing in the bottom panel of both Figure 6 and Figure 7. This is due to the non-trivial implementation of the anti-symmetric property in classical ENNs.

Finally Table 1 shows the accuracy of the DNN for the fully anti-symmetric dataset. The different rows and columns represent different choices of the number of parameters and the number of training samples, respectively. These numbers are compared against those in right-bottom panel of Figure 7. The EQNN achieves 0.95 accuracy with 20 parameters and 200 training samples, while the DNN requires more parameters and/or more training samples.

5. Conclusions

In this paper, we examined the performance of Equivariant Quantum Neural Networks and Quantum Neural Networks, compared against their classical counterparts, Equivariant Neural Networks and Deep Neural Networks, considering two toy examples for a binary classification task. Our study demonstrates that EQNNs and QNNs outperform their classical counterparts, particularly in scenarios with fewer parameters and smaller training datasets. This highlights the potential of quantum-inspired architectures in resource-constrained settings. This point has been emphasized in a similar study recently in Ref. [35], which showed that an EQNN outperforms the non-equivariant one in terms of generalization power, especially with a small training set size. We note a more significant enhancement in the performance of an EQNN and QNN compared to an ENN and DNN, particularly evident in the anti-symmetric example rather than the symmetric one. This underscores the robustness of quantum algorithms. The code used for this study is publicly available at https://github.com/ZhongtianD/EQNN/tree/main (accessed on 7 March 2024).

While our current study has primarily focused on an EQNN with discrete symmetries, it is crucial to acknowledge the significant role that continuous symmetries, such as Lorentz symmetry or gauge symmetries, play in particle physics. In our future research, we aim to compare an EQNN with continuous symmetries against classical neural networks. Exploring more complex datasets with high-dimensional features is another direction we plan to pursue. However, handling such examples would necessitate an increase in the number of network parameters, prompting an investigation into related issues like overparameterization, barren plateaus, and others.

Author Contributions

Conceptualization, Z.D.; methodology, M.C.C., G.R.D., Z.D., R.T.F., S.G., D.J., K.K., T.M., K.T.M., K.M. and E.B.U.; software, Z.D.; validation, M.C.C., G.R.D., Z.D., R.T.F., T.M. and E.B.U.; formal analysis, Z.D.; investigation, M.C.C., G.R.D., Z.D., R.T.F., T.M. and E.B.U.; resources, Z.D., K.T.M. and K.M.; data curation, G.R.D., S.G. and T.M.; writing—original draft preparation, Z.D.; writing—review and editing, S.G., D.J., K.K., K.T.M. and K.M.; visualization, Z.D.; supervision, S.G., D.J., K.K., K.T.M. and K.M.; project administration, S.G., D.J., K.K., K.T.M. and K.M.; funding acquisition, S.G. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The dataset used in this analysis was sampled from Equations (1), (4) and (6).

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:

API	Application Processing Interface
AUC	Area Under the Curve
DNN	Deep Neural Network
ENN	Equivariant Neural Network
EQNN	Equivariant Quantum Neural Network
HEP	High-Energy Physics
LHC	Large Hadron Collider
MDPI	Multidisciplinary Digital Publishing Institute
ML	Machine Learning
NN	Neural Network
QML	Quantum Machine Learning
QNN	Quantum Neural Network
ROC	Receiver Operating Characteristic

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Table

Figure 1. Pictorial illustration of the first dataset used in this study—the symmetric case (1).

Figure 2. Pictorial illustration of the second dataset used in this study—the anti-symmetric case (4).

Figure 3. Pictorial illustration of the third dataset used in this study—the fully anti-symmetric case (6).

Figure 4. Illustration of the quantum circuit used for QNN at depth 2. This circuit is repeated up to depth four (five) times with different parameters for the symmetric (anti-symmetric and fully anti-symmetric) case. The data points [Forumla omitted. See PDF.] are loaded via angle embedding with rotation gates, followed by another rotation [Forumla omitted. See PDF.] with arbitrary angle parameters.

Figure 5. Illustration of the quantum circuit used for EQNN at depth 1. This circuit is repeated five (ten) times with different parameters for the symmetric (anti-symmetric and fully anti-symmetric) case. The data points [Forumla omitted. See PDF.] are loaded via angle embedding with two [Forumla omitted. See PDF.] gates, [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The remaining circuits are parameterized by [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.].

Figure 6. ROC (left) and accuracy (right) curves for the symmetric (top), anti-symmetric (middle), and fully anti-symmetric (bottom) example.

Figure 7. AUC as a function of the number of parameters (left) for fixed [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.], and as a function of [Forumla omitted. See PDF.] (right) with a fixed number of parameters as shown in the legend, for the symmetric (top), anti-symmetric (middle), and fully anti-symmetric example (bottom).

Table 1

Accuracy of DNN for the fully anti-symmetric dataset. The different rows (columns) represent different choices of the number of parameters (the number of training samples).

${\mathit{N}}_{\mathbf{params}}$\${\mathit{N}}_{\mathbf{train}}$	100	200	300	400	500	600	700	800	900
105	0.764	0.855	0.879	0.963	0.973	0.981	0.981	0.982	0.988
85	0.669	0.743	0.804	0.953	0.951	0.978	0.986	0.946	0.981
67	0.587	0.722	0.695	0.946	0.886	0.9632	0.975	0.944	0.980
51	0.624	0.655	0.856	0.926	0.908	0.876	0.846	0.974	0.986
37	0.596	0.696	0.639	0.782	0.747	0.816	0.849	0.922	0.952

References

1. Shanahan, P.; Terao, K.; Whiteson, D. Snowmass 2021 Computational Frontier CompF03 Topical Group Report: Machine Learning. arXiv; 2022; arXiv: 2209.07559

2. Feickert, M.; Nachman, B. A Living Review of Machine Learning for Particle Physics. arXiv; 2021; arXiv: 2102.02770

3. Cohen, T.; Welling, M. Group Equivariant Convolutional Networks. Proceedings of the 33rd International Conference on Machine Learning; New York, NY, USA, 20–22 June 2016; Balcan, M.F.; Weinberger, K.Q. Volume 48, pp. 2990-2999.

4. Krizhevsky, A.; Sutskever, I.; Hinton, G.E. ImageNet Classification with Deep Convolutional Neural Networks. Proceedings of the Advances in Neural Information Processing Systems; Lake Tahoe, NV, USA, 3–6 December 2012; Pereira, F.; Burges, C.; Bottou, L.; Weinberger, K. Curran Associates, Inc.: Sydney, Australia, 2012; Volume 25.

5. Jumper, J.M.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Zídek, A.; Potapenko, A. et al. Highly accurate protein structure prediction with AlphaFold. Nature; 2021; 596, pp. 583-589. [DOI: https://dx.doi.org/10.1038/s41586-021-03819-2]

6. Bogatskiy, A.; Anderson, B.; Offermann, J.; Roussi, M.; Miller, D.; Kondor, R. Lorentz Group Equivariant Neural Network for Particle Physics. Proceedings of the 37th International Conference on Machine Learning; Virtual, 13–18 July 2020; Volume 119, pp. 992-1002.

7. Cohen, T.S.; Weiler, M.; Kicanaoglu, B.; Welling, M. Gauge Equivariant Convolutional Networks and the Icosahedral CNN. arXiv; 2019; 2, arXiv: 1902.04615

8. Boyda, D.; Kanwar, G.; Racanière, S.; Rezende, D.J.; Albergo, M.S.; Cranmer, K.; Hackett, D.C.; Shanahan, P.E. Sampling using SU(N) gauge equivariant flows. Phys. Rev. D; 2021; 103, 074504. [DOI: https://dx.doi.org/10.1103/PhysRevD.103.074504]

9. Favoni, M.; Ipp, A.; Müller, D.I.; Schuh, D. Lattice Gauge Equivariant Convolutional Neural Networks. Phys. Rev. Lett.; 2022; 128, 032003. [DOI: https://dx.doi.org/10.1103/PhysRevLett.128.032003]

10. Dolan, M.J.; Ore, A. Equivariant Energy Flow Networks for Jet Tagging. Phys. Rev. D; 2021; 103, 074022. [DOI: https://dx.doi.org/10.1103/PhysRevD.103.074022]

11. Bulusu, S.; Favoni, M.; Ipp, A.; Müller, D.I.; Schuh, D. Generalization capabilities of translationally equivariant neural networks. Phys. Rev. D; 2021; 104, 074504. [DOI: https://dx.doi.org/10.1103/PhysRevD.104.074504]

12. Preskill, J. Quantum Computing in the NISQ era and beyond. Quantum; 2018; 2, 79. [DOI: https://dx.doi.org/10.22331/q-2018-08-06-79]

13. Feynman, R.P. Simulating physics with computers. Int. J. Theor. Phys.; 1982; 21, pp. 467-488. [DOI: https://dx.doi.org/10.1007/BF02650179]

14. Georgescu, I.M.; Ashhab, S.; Nori, F. Quantum Simulation. Rev. Mod. Phys.; 2014; 86, 153. [DOI: https://dx.doi.org/10.1103/RevModPhys.86.153]

15. Ramírez-Uribe, S.; Rentería-Olivo, A.E.; Rodrigo, G.; Sborlini, G.F.R.; Vale Silva, L. Quantum algorithm for Feynman loop integrals. JHEP; 2022; 5, 100. [DOI: https://dx.doi.org/10.1007/JHEP05(2022)100]

16. Bepari, K.; Malik, S.; Spannowsky, M.; Williams, S. Quantum walk approach to simulating parton showers. Phys. Rev. D; 2022; 106, 056002. [DOI: https://dx.doi.org/10.1103/PhysRevD.106.056002]

17. Li, T.; Guo, X.; Lai, W.K.; Liu, X.; Wang, E.; Xing, H.; Zhang, D.B.; Zhu, S.L. Partonic collinear structure by quantum computing. Phys. Rev. D; 2022; 105, L111502. [DOI: https://dx.doi.org/10.1103/PhysRevD.105.L111502]

18. Bepari, K.; Malik, S.; Spannowsky, M.; Williams, S. Towards a quantum computing algorithm for helicity amplitudes and parton showers. Phys. Rev. D; 2021; 103, 076020. [DOI: https://dx.doi.org/10.1103/PhysRevD.103.076020]

19. Jordan, S.P.; Lee, K.S.M.; Preskill, J. Quantum Algorithms for Fermionic Quantum Field Theories. arXiv; 2014; arXiv: 1404.7115

20. Preskill, J. Simulating quantum field theory with a quantum computer. PoS; 2018; LATTICE2018, 024. [DOI: https://dx.doi.org/10.22323/1.334.0024]

21. Bauer, C.W.; de Jong, W.A.; Nachman, B.; Provasoli, D. Quantum Algorithm for High Energy Physics Simulations. Phys. Rev. Lett.; 2021; 126, 062001. [DOI: https://dx.doi.org/10.1103/PhysRevLett.126.062001]

22. Abel, S.; Chancellor, N.; Spannowsky, M. Quantum computing for quantum tunneling. Phys. Rev. D; 2021; 103, 016008. [DOI: https://dx.doi.org/10.1103/PhysRevD.103.016008]

23. Abel, S.; Spannowsky, M. Quantum-Field-Theoretic Simulation Platform for Observing the Fate of the False Vacuum. PRX Quantum; 2021; 2, 010349. [DOI: https://dx.doi.org/10.1103/PRXQuantum.2.010349]

24. Davoudi, Z.; Linke, N.M.; Pagano, G. Toward simulating quantum field theories with controlled phonon-ion dynamics: A hybrid analog-digital approach. Phys. Rev. Res.; 2021; 3, 043072. [DOI: https://dx.doi.org/10.1103/PhysRevResearch.3.043072]

25. Mott, A.; Job, J.; Vlimant, J.R.; Lidar, D.; Spiropulu, M. Solving a Higgs optimization problem with quantum annealing for machine learning. Nature; 2017; 550, pp. 375-379. [DOI: https://dx.doi.org/10.1038/nature24047]

26. Blance, A.; Spannowsky, M. Unsupervised event classification with graphs on classical and photonic quantum computers. JHEP; 2020; 21, 170. [DOI: https://dx.doi.org/10.1007/JHEP08(2021)170]

27. Wu, S.L.; Chan, J.; Guan, W.; Sun, S.; Wang, A.; Zhou, C.; Livny, M.; Carminati, F.; Di Meglio, A.; Li, A.C. et al. Application of quantum machine learning using the quantum variational classifier method to high energy physics analysis at the LHC on IBM quantum computer simulator and hardware with 10 qubits. J. Phys. G; 2021; 48, 125003. [DOI: https://dx.doi.org/10.1088/1361-6471/ac1391]

28. Blance, A.; Spannowsky, M. Quantum Machine Learning for Particle Physics using a Variational Quantum Classifier. JHEP; 2021; 2, 212. [DOI: https://dx.doi.org/10.1007/JHEP02(2021)212]

29. Abel, S.; Blance, A.; Spannowsky, M. Quantum optimization of complex systems with a quantum annealer. Phys. Rev. A; 2022; 106, 042607. [DOI: https://dx.doi.org/10.1103/PhysRevA.106.042607]

30. Wu, S.L.; Sun, S.; Guan, W.; Zhou, C.; Chan, J.; Cheng, C.L.; Pham, T.; Qian, Y.; Wang, A.Z.; Zhang, R. et al. Application of quantum machine learning using the quantum kernel algorithm on high energy physics analysis at the LHC. Phys. Rev. Res.; 2021; 3, 033221. [DOI: https://dx.doi.org/10.1103/PhysRevResearch.3.033221]

31. Chen, S.Y.C.; Wei, T.C.; Zhang, C.; Yu, H.; Yoo, S. Hybrid Quantum-Classical Graph Convolutional Network. arXiv; 2021; arXiv: 2101.06189

32. Terashi, K.; Kaneda, M.; Kishimoto, T.; Saito, M.; Sawada, R.; Tanaka, J. Event Classification with Quantum Machine Learning in High-Energy Physics. Comput. Softw. Big Sci.; 2021; 5, 2. [DOI: https://dx.doi.org/10.1007/s41781-020-00047-7]

33. Araz, J.Y.; Spannowsky, M. Classical versus quantum: Comparing tensor-network-based quantum circuits on Large Hadron Collider data. Phys. Rev. A; 2022; 106, 062423. [DOI: https://dx.doi.org/10.1103/PhysRevA.106.062423]

34. Ngairangbam, V.S.; Spannowsky, M.; Takeuchi, M. Anomaly detection in high-energy physics using a quantum autoencoder. Phys. Rev. D; 2022; 105, 095004. [DOI: https://dx.doi.org/10.1103/PhysRevD.105.095004]

35. Chang, S.Y.; Grossi, M.; Saux, B.L.; Vallecorsa, S. Approximately Equivariant Quantum Neural Network for p4m Group Symmetries in Images. Proceedings of the 2023 International Conference on Quantum Computing and Engineering; Bellevue, WA, USA, 17–22 September 2023.

36. Nguyen, Q.T.; Schatzki, L.; Braccia, P.; Ragone, M.; Coles, P.J.; Sauvage, F.; Larocca, M.; Cerezo, M. Theory for Equivariant Quantum Neural Networks. arXiv; 2022; arXiv: 2210.08566

37. Meyer, J.J.; Mularski, M.; Gil-Fuster, E.; Mele, A.A.; Arzani, F.; Wilms, A.; Eisert, J. Exploiting Symmetry in Variational Quantum Machine Learning. PRX Quantum; 2023; 4, 010328. [DOI: https://dx.doi.org/10.1103/PRXQuantum.4.010328]

38. West, M.T.; Sevior, M.; Usman, M. Reflection equivariant quantum neural networks for enhanced image classification. Mach. Learn. Sci. Technol.; 2023; 4, 035027. [DOI: https://dx.doi.org/10.1088/2632-2153/acf096]

39. Skolik, A.; Cattelan, M.; Yarkoni, S.; Bäck, T.; Dunjko, V. Equivariant quantum circuits for learning on weighted graphs. npj Quantum Inf.; 2023; 9, 47. [DOI: https://dx.doi.org/10.1038/s41534-023-00710-y]

40. Kim, I.W. Algebraic Singularity Method for Mass Measurement with Missing Energy. Phys. Rev. Lett.; 2010; 104, 081601. [DOI: https://dx.doi.org/10.1103/PhysRevLett.104.081601] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20366926]

41. Franceschini, R.; Kim, D.; Kong, K.; Matchev, K.T.; Park, M.; Shyamsundar, P. Kinematic Variables and Feature Engineering for Particle Phenomenology. Rev. Mod. Phys.; 2022; 95, 045004. [DOI: https://dx.doi.org/10.1103/RevModPhys.95.045004]

42. Kersting, N. On Measuring Split-SUSY Gaugino Masses at the LHC. Eur. Phys. J. C; 2009; 63, pp. 23-32. [DOI: https://dx.doi.org/10.1140/epjc/s10052-009-1063-6]

43. Bisset, M.; Lu, R.; Kersting, N. Improving SUSY Spectrum Determinations at the LHC with Wedgebox Technique. JHEP; 2011; 5, 095. [DOI: https://dx.doi.org/10.1007/JHEP05(2011)095]

44. Burns, M.; Matchev, K.T.; Park, M. Using kinematic boundary lines for particle mass measurements and disambiguation in SUSY-like events with missing energy. JHEP; 2009; 5, 094. [DOI: https://dx.doi.org/10.1088/1126-6708/2009/05/094]

45. Debnath, D.; Gainer, J.S.; Kim, D.; Matchev, K.T. Edge Detecting New Physics the Voronoi Way. EPL; 2016; 114, 41001. [DOI: https://dx.doi.org/10.1209/0295-5075/114/41001]

46. Debnath, D.; Gainer, J.S.; Kilic, C.; Kim, D.; Matchev, K.T.; Yang, Y.P. Detecting kinematic boundary surfaces in phase space: Particle mass measurements in SUSY-like events. JHEP; 2017; 6, 092. [DOI: https://dx.doi.org/10.1007/JHEP06(2017)092]

47. Bogatskiy, A.; Hoffman, T.; Miller, D.W.; Offermann, J.T. PELICAN: Permutation Equivariant and Lorentz Invariant or Covariant Aggregator Network for Particle Physics. arXiv; 2022; arXiv: 2211.00454

48. Hao, Z.; Kansal, R.; Duarte, J.; Chernyavskaya, N. Lorentz group equivariant autoencoders. Eur. Phys. J. C; 2023; 83, 485. [DOI: https://dx.doi.org/10.1140/epjc/s10052-023-11633-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37303461]

49. Buhmann, E.; Kasieczka, G.; Thaler, J. EPiC-GAN: Equivariant point cloud generation for particle jets. SciPost Phys.; 2023; 15, 130. [DOI: https://dx.doi.org/10.21468/SciPostPhys.15.4.130]

50. Batatia, I.; Geiger, M.; Munoz, J.; Smidt, T.; Silberman, L.; Ortner, C. A General Framework for Equivariant Neural Networks on Reductive Lie Groups. arXiv; 2023; arXiv: 2306.00091

51. Pérez-Salinas, A.; Cervera-Lierta, A.; Gil-Fuster, E.; Latorre, J.I. Data re-uploading for a universal quantum classifier. Quantum; 2020; 4, 226. [DOI: https://dx.doi.org/10.22331/q-2020-02-06-226]

52. Ahmed, S. Data-Reuploading Classifier. Available online: https://pennylane.ai/qml/demos/tutorial_data_reuploading_classifier (accessed on 7 March 2024).