## Modeling stochastic eye tracking data: A comparison of quantum generative adversarial networks and Markov models

Shailendra Bhandari shailendra.bhandari@oslomet.no Department of Computer Science, OsloMet - Oslo Metropolitan University, N-0130 Oslo, Norway

Pedro Lencastre pedrog@oslomet.no Department of Computer Science, OsloMet - Oslo Metropolitan University, N-0130 Oslo, Norway

## ABSTRACT

We explore the use of quantum generative adversarial networks QGANsformodelingeyemovementvelocity data. We assess whether the advanced computational capabilities of QGANs can enhance the modeling of complex stochastic distribution beyond the traditional mathematical models, particularly the Markov model. The findings indicate that while QGANs demonstrate potential in approximating complex distributions, the Markov model consistently outperforms in accurately replicating the real data distribution. This comparison underlines the challenges and avenues for refinement in time series data generation using quantum computing techniques. It emphasizes the need for further optimization of quantum models to better align with real-world data characteristics.

## Pedro Lind

pedrolin@oslomet.no Department of Computer Science, OsloMet - Oslo Metropolitan University, N-0130 Oslo Norway; Simula Research Laboratory, Numerical Analysis and Scientific Computing, 0164 Oslo, Norway generator when it produces data that is not convincing. Throughout the training process, the generator continually strives to enhance its ability to create increasingly convincing forgeries. Simultaneously, the discriminator aims to improve its ability to discern between real and fake data accurately. Therefore, a GAN is a battle between two adversaries, the generator and the discriminator. The point of balance in this dynamic battle is achieved when the generator is capable of producing forgeries that are indistinguishable from the original training data, leaving the discriminator with a mere 50% certainty in distinguishing real from fake outputs.

## KEYWORDS

Quantum Generative Adversarial Networks (QGANs), Markov models, Eye-tracking data

## ACMReference Format:

Shailendra Bhandari, Pedro Lencastre, and Pedro Lind. 2024. Modeling stochastic eye tracking data: A comparison of quantum generative adversarial networks and Markov models. In Genetic and Evolutionary Computation Conference (GECCO '24 Companion), July 14-18, 2024, Melbourne, VIC, Australia. ACM, New York, NY, USA, 8 pages. https://doi.org/10.1145/3638530.3664134

## 1 INTRODUCTION

The generative adversarial networks (GANs) [11, 12] are the type of artificial intelligence algorithms that consist of two neural networks: the generator and the discriminator . The generator is trained to create realistic data, which are then used as negative instances for the discriminator to learn from. This process is also known as generative modeling, a branch of machine learning that involves training a model to produce new data similar to a given dataset. Meanwhile, the discriminator improves its ability to differentiate between authentic and generated data, providing feedback to the

Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for third-party components of this work must be honored. For all other uses, contact the owner/author(s).

GECCO '24, July 14-18, 2024, Melbourne, VIC, Australia

© 2024 Copyright held by the owner/author(s).

ACM ISBN 979-8-4007-0495-6/24/07

https://doi.org/10.1145/3638530.3664134

GANs become tremendously interesting and challenging topics in Machine Learning [20]. It is beyond the boundary line of computational creativity, showcasing increasingly amazing examples each year [23]. While GANs have shown remarkable power and interest, it is also limited by various challenges: difficulties in achieving stable training of the GAN [9, 27], the vanishing gradient, and the mode collapse [1]. Efforts from both the academic community and the commercial sector have been made to address these problems [8, 16, 25, 34, 36]. The exploration of Quantum Deep Learning models has been gaining popularity recently, followed by the concept of quantum supremacy, which Google has showcased [2], along with the potential for a quantum edge in machine learning. Recently, discussions around Quantum Generative Adversarial Networks (QGANs) mainly focused on applications like efficient data handling or seeking benefits over traditional classical models [38]. An example of a successful QGAN involves using a manageable number of quantum bits to train the GAN to recognize and replicate random patterns. The motivation behind quantumenabled research is the potential quantum advantage that could partially alleviate the computational complexity issue and search for the best solution to the mode collapse and vanishing gradient problem.

In the context of integrating quantum principles into neural networks, it is important to recognize foundational work that has set the stage for contemporary advancements. One of the important works is by Behera et al. [3], who developed a recurrent quantum neural network (RQNN) model to describe eye movements in tracking moving targets. Their work utilized a nonlinear Schrödinger wave equation to mediate the neural responses, providing a foundation for understanding how quantum mechanics can be integrated with neural processing. This early exploration into quantum neural networks is important. Our approach to QGANs aims to enhance

Figure 1: Generative adversarial network models workflow: the generator generates data samples ( 𝑔 𝑡 ) to imitate the realworld data and tries to fool the discriminator. The discriminator differentiates the generated and the training data samples by training both the generator and discriminator alternatively until the loss converges towards the Nash equilibrium.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00673v1-figure-87921be5-92de-4288-ad26-cac135ffa1d4.jpg)

the generative capabilities of neural networks through the incorporation of quantum bits, quantum gates, and the fundamental properties of quantum mechanics: superposition and entanglement.

The main goal of this paper is to use the QGAN and test its performance in reproducing stochastic trajectories and compare them with some benchmarks from stochastic modeling. Previously, Lencastre et al. [21] tested the classical GANs'performance, comparing with a mathematical model, namely a Markov chain. They found that GANs struggle to capture rare events and cross-feature relations and are unable to create faithful synthetic data successfully. Here, we test the QGANs performance compared with the previous classical GANs and mathematical models. One of the problems of the classical GAN is the vanishing gradient problem when generating discrete data [21] to overcome this issue we explore the ability of quantum computing, which naturally has advantages to equip GANs with the ability to deal with discrete data.

## 2 THE QUANTUM GAN MODEL

In our study, the utilized GAN framework comprises two neural networks: the generator and the discriminator, as illustrated in Figure 1. These networks undergo alternating training phases. Consider a classical training data set 𝑋 = { 𝑥 0 , ..., 𝑥 𝑠 -1 } drawn from unknown time series distribution. The generator ( 𝐺 𝜃 ) receives a random noise vector 𝑧 and then produces the generated sample 𝐺 𝑧 ( ) . The discriminator 𝐷 𝜙 is trained to distinguish between the training data 𝑥 and the generated data 𝐺 𝑧 ( ) . The parameters of 𝐷 𝜙 are updated in order to maximize [33]:

$$\mathbb { E } _ { x \sim P _ { d } ( x ) } [ \log D _ { \phi } ( x ) ] + \mathbb { E } _ { z \sim P _ { z } ( z ) } [ \log ( 1 - D _ { \phi } ( G _ { \theta } ( z ) ) ) ] \, \ \ ( 1 )$$

where 𝑃 𝑑 ( 𝑥 ) is the real-time series distribution and the 𝑃 𝑧 ( 𝑧 ) is the distribution of the input noise. The output 𝐷 𝜙 ( 𝑥 ) can be explained as the probability that 𝐷 𝜙 thinks the sample 𝑥 is real. The goal of the discriminator here is to maximize the probability of correctly classifying to make 𝐷 𝜙 a better adversary so that 𝐺 𝜃 has to try harder to fool the 𝐷 𝜙 . Similarly, the parameters of the 𝐺 𝜃 are updated to maximize:

$$\mathbb { E } _ { z \sim P _ { z } ( z ) } [ \log ( D _ { \phi } ( G _ { \theta } ( z ) ) ) ] \,,$$

to convince 𝐷 𝜙 that the generates samples 𝐺 𝑧 ( ) are real. This process is illustrated in Figure 1.

The goal of optimizing classical GANs can be approached from several perspectives. In this study, we adopt the non-saturating loss function [7], which is also implemented in the original GAN publication's code [10]. The generator loss function is given by:

$$L _ { G } ( \theta ) = - \mathbb { E } _ { z \sim P _ { z } ( z ) } [ \log ( D _ { \phi } ( G _ { \theta } ( z ) ) ) ] \,, \quad \ \ ( 3 )$$

which aims at maximizing the likelihood that the generator creates samples that are labeled as real data samples. Also, the discriminator's loss function is given by the expression:

$$\text{real} _ { \ } l _ { D } ( \phi ) = \mathbb { E } _ { x \sim P _ { d } ( x ) } \left [ \log D _ { \phi } ( x ) \right ] + \mathbb { E } _ { z \sim P _ { z } ( z ) } \left [ \log ( 1 - D _ { \phi } ( G _ { \theta } ( z ) ) ) \right ] \,, \\ \text{in} \, \infty \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dists \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dcl \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dCL \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dcd \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \d.. \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \d. \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ddots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad d. \end{$$

which aims at maximizing the likelihood that the discriminator labels training data samples as training data samples and generated data samples as generated data samples. In practice, the expected values are approximated by batches of size 𝑚

$$L _ { G } ( \phi, \theta ) = - \frac { 1 } { m } \sum _ { i = 1 } ^ { m } [ \log D _ { \phi } ( G _ { \theta } ( z ^ { ( i ) } ) ) ],$$

and

$$\underset { \mathbb { I } y, } { \mathrm r e } \underset { \mathbb { I } y, } { \mathrm l } L _ { D } ( \phi, \theta ) = \frac { 1 } { m } \sum _ { i = 1 } ^ { m } [ \log D _ { \phi } ( x ^ { ( i ) } ) + \log ( 1 - D _ { \phi } ( G _ { \theta } ( z ^ { ( i ) } ) ) ) ], \quad ( 6 )$$

for 𝑥 ′ in 𝑋 and 𝑧 ′ sampled from 𝑃 𝑑 ( 𝑥 ) .

## 2.1 Quantum generator

The model we designed is a hybrid architecture (quantum-classical) where the input is in classical form with a classical discriminator and generates an output by measurement of a state from the parametrized quantum circuit. Specifically, it is a quantum GAN that generates a classical discrete time series, which is composed of a parametrized quantum circuit as a generator and a classical feed-forward neural network as a discriminator.

The parameterized quantum circuit is then expressed as 𝐺 ( 𝜽 ) , where 𝜽 = { 𝜃 1 , . . . , 𝜃 𝑘 } represents the set of parameters which can be tuned as necessary. This circuit is operated by the single qubit 𝑅 𝑌 and 𝑅 𝑍 rotation gates, along with the two-qubit Controlled-NOT (CX) gates, all operating on an initial state [4]. When applying the rotation operators on each qubit, it is expressed as:

$$\prod _ { i = 1 } ^ { N } R _ { Z } ^ { i } ( \theta _ { l, 2 } ^ { i } ) R _ { Y } ^ { i } ( \theta _ { l, 1 } ^ { i } ),$$

Figure 2: The quantum generator circuit in variational form with 𝐿 layers acting on 𝑛 qubits. Each layer in the circuit is composed of single qubit rotation gates ( 𝑅 𝑍 , 𝑅 𝑌 ) and two-qubit controlled-phase gates.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00673v1-figure-3d692cba-d0fa-4ea3-9eb3-268bcd7d844d.jpg)

where 𝑙 denotes the 𝑙 𝑡ℎ layer and 𝑖 denotes the 𝑖 𝑡ℎ qubit. 𝑅 𝑧 and 𝑅 𝑌 are the rotation gates given by following expressions:

$$R _ { Z } ( \theta ) = e ^ { - i \theta \sigma _ { z } / 2 } = \cos \frac { \theta } { 2 } I - i \sin \frac { \theta } { 2 } \sigma _ { z } = \begin{bmatrix} e ^ { - i \theta / 2 } & 0 \\ 0 & e ^ { i \theta / 2 } \end{bmatrix} \, \\ R _ { Y } ( \theta ) = e ^ { - i \theta \sigma _ { y } / 2 } = \cos \frac { \theta } { 2 } I - i \sin \frac { \theta } { 2 } \sigma _ { y } = \begin{bmatrix} \cos \frac { \theta } { 2 } & - i \sin \frac { \theta } { 2 } \\ i \sin \frac { \theta } { 2 } & \cos \frac { \theta } { 2 } \end{bmatrix}. \quad \text{or eq} \\ \sum \quad \sim \. \ \, \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots$$

The Controlled-NOT gate is used to create entanglement between the qubits, and this process is expressed as:

$$\prod _ { i = 1 } ^ { N } C U _ { ( i \, m o d \, N ) + 1 } ^ { i } \,, \quad \ \ ( 9 ) \quad \ g$$

where 𝑖 denote the controlled qubit and ( 𝑖 𝑚𝑜𝑑 𝑁 ) + 1 represent the target qubit. The number of single-qubit gates represents the number of parameters in the generative quantum circuit in a circuit which is 3 𝑁 per layer.

The parametrized quantum circuit for generating 𝑁 -𝑏𝑖𝑡 samples consists of 𝑁 qubit and 𝑈 𝜃 ( 𝐿 ) layers as shown in Figure 2. The input quantum state | 𝜓 𝑖𝑛 ⟩ is initialized to | 0 ⟩ ⊗ 𝑁 and passed through 𝐿 layers of unitary operator. All the qubits outcomes of the circuit are measured on the computational basis at the end of the circuit. The measurement outcome is then collected to form 𝑁 -bit sample 𝑋 . Overall, the quantum generator is trained to transform 𝑁 -qubit input to 𝑁 -qubit output state [38]:

$$\ m p u a \, \infty \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a\, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ R r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ n o r \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \,
 \text{with} a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q ua \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q ua \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q ua \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ n o r \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r g u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r g u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r g u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a \, \ r r q u a$$

where samples are drawn by measuring the output state | 𝑔 𝜃 ⟩ in computational basis with measurement outcomes | 𝑗 ⟩ , 𝑗 ∈ { 0 1 , ...., 2 𝑁 -1 . } 𝑁 is the number of qubits. In addition, the term 𝑝 𝑗 is the sampling probability associated with the basis state | 𝑥 𝑗 ⟩ . For 𝑚 data samples ( 𝑔 𝑙 ) from the quantum generator and 𝑚 randomly chosen training data samples ( 𝑥 𝑙 ), where 𝑙 = 1 , ...𝑚 , the generator and discriminator are optimized with a respective loss function.

The generator's goal is to generate data that can fool the discriminator. The generator loss function for a data batch of size 𝑚 is expressed as:

$$L _ { G } ( \phi, \theta ) = - \frac { 1 } { m } \sum _ { i = 1 } ^ { m } [ \log D _ { \phi } ( g ^ { l } ) ], \quad \quad ( 1 1 )$$

or equivalently,

$$L _ { G } ( \phi, \theta ) = \sum _ { j = 0 } ^ { 2 ^ { N } - 1 } p _ { \theta } ^ { j } \log ( D _ { \phi } ( \varrho ^ { j } ) ), \quad \quad ( 1 2 )$$

where 𝑝 𝑗 𝜃 = | ⟨ 𝑗 | 𝑔 𝜃 ⟩| 2 . The loss function of the discriminator is given by

$$L _ { D } ( \phi, \theta ) = \frac { 1 } { m } \sum _ { l = 1 } ^ { m } [ l o g D _ { \phi } ( x ^ { l } ) + l o g ( 1 - D _ { \phi } ( g ^ { l } ) ) ]. \quad ( 1 3 )$$

Gradient-based optimization techniques can enhance convergence speed, particularly near local optima in a convex region, when compared to methods that do not utilize gradient information [14]. This method for calculating analytical gradients [6, 24, 37, 38] for the variational circuit is as discussed below. The parameters 𝜽 can be updated with gradient-based methods that require the evaluation of

$$\frac { \partial \mathcal { L } _ { G } ( \phi, \theta ) } { \partial \theta ^ { i, l } } = - \sum _ { j = 1 } ^ { m } \frac { \partial p _ { \theta } ^ { j } } { \partial \theta _ { i, l } } \log ( D _ { \phi } ( g ^ { j } ) ). \quad \quad ( 1 4 )$$

Eq. (14) can be further evaluated based on Ref. [6]

$$\frac { \partial p _ { j } ^ { \theta } } { \partial \theta ^ { i, l } } = \frac { 1 } { 2 } \left ( p _ { \theta _ { + } ^ { i, j } } ^ { j } - p _ { \theta _ { - } ^ { i, j } } ^ { j } \right ),$$

with 𝜃 𝑖,𝑙 ± = 𝜃 𝑖,𝑙 ± 𝜋 2 𝑒 𝑖,𝑙 and 𝑒 𝑖,𝑙 denoting the ( 𝑖, 𝑙 ) -unit vector of the respective parameter space.

The selection of parameters 𝜽 is important, especially when 𝑙 &gt; 1. A key factor to consider is the circuit depth, as increasing it enhances the complexity of the quantum circuit. Utilizing a circuit depth greater than 1 is advantageous for training on complex datasets, as deeper circuits can capture more intricate structural patterns. Thus, it is essential to design a circuit that is both deep and rich in parameters. This choice is driven by the need to effectively capture and represent the complex distribution characteristics of the data. By increasing both the number of parameters and the

circuit's depth, we can significantly enhance the generator's ability to accurately model the data. Such a strategy is pivotal for quantum generators tasked with processing and generating data distributions of considerable complexity.

## 2.2 The discriminator

The discriminator model, a classical feed-forward neural network, is constructed and implemented using PyTorch's [28] neural network module. It comprises a three-layer LSTM [15] (Long Short-Term Memory) network with 128 hidden units and bidirectional processing. This LSTM feeds into a series of linear layers that further process the data. A dropout of 0.3 is applied after the first linear transformation to prevent overfitting. The output of the discriminator is a single scalar value representing the probability that the input data is from the real distribution, obtained through a sigmoid activation function.

The discriminator is trained using AMSGRAD [30] with a learning rate of 0.002 and a momentum of 0.999. Using the first and the second momentum terms, AMASGRAD is a robust optimization technique for non-stationary objective functions and noisy gradients [17]. The loss function employed is binary cross-entropy (BCE) and is suitable for binary classification tasks such as distinguishing between real and fake data sequences. During training, the discriminator evaluates both real and generated data sequences, updating its weights to minimize the loss function. The training stability is maintained using a gradient penalty on a discriminator's loss function [9, 18]. The analytic computation of the quantum generator loss function gradients is implemented based on Equations. (14) and (15). The training procedure is iteratively conducted for a specified number of epochs, with both generator and discriminator losses monitored to ensure the model's convergence.

The adversarial training of the quantum GAN described above is illustrated in Algorithm 1. The training of the GAN iterates for the specified number of epochs, until the convergence of the loss function. At each epoch, the loss function optimizes alternatively to parameters 𝜃 and 𝜙 .

## 3 MARKOV-CHAIN MODELS OF TIME-SERIES

We generate a two-dimensional time series based on a Markov transition matrix 𝑇 , where each element 𝑇 𝑖 𝑗 denotes the transition probability from state 𝑖 to state 𝑗 . Mathematically, a time series 𝑋 𝑡 is said to follow a Markov process if it fulfills a Markov property:

$$P r ( X _ { t = j } = \hat { x } _ { j } | X _ { t = j - 1 } = \hat { x } _ { j - 1 }, \dots, X _ { t = 0 } = \hat { x } _ { 0 } ) = & & \quad \text{and} \, t \\ P r ( X _ { t = j } = \hat { x } _ { j } | X _ { t = j - 1 } = \hat { x } _ { j - 1 } ), & & \quad \text{and} \, t$$

for all positive integers j, where capital letters mean stochastic variables at different time steps and lowercase letters are the respective values of those variables.

In a Markov model, time-series generation is based on computing the conditional probability 𝑃𝑟 ( 𝑋 𝑡 = 𝑛 + 1 = 𝑥 𝑛 + 1 | 𝑋 𝑡 = 𝑛 = 𝑥 𝑛 ) . This is estimated empirically with the Gaussian estimation kernel 𝐾 as follows:

$$P r ( X _ { t = n + 1 } = x _ { n + 1 } \, | \, X _ { t = n } = x _ { n } ) = \frac { P r ( X _ { t = n + 1 } = x _ { n + 1 }, X _ { t = n } = x _ { n } ) } { P r ( X _ { t = n } = x _ { n } ) }, \quad \text{ of dal} \\ ( 1 7 ) \quad \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        }$$

Algorithm 1: QGAN training for eye-tracking data

- 1. Input: Eye-tracking velocity distribution,

hyperparameters;

- 2. Initialization:
- 2.1. Load and preprocess eye-tracking data;
- 2.2. Set seed for reproducibility across quantum and classical computations;
- 2.3. Initialize quantum circuit for the generator with specified qubits and an EfficientSU2 ansatz;
- 2.4. Initialize a PyTorch-based classical neural network that represents the classical discriminator;
- 2.5. Define adversarial loss function for training;
- 2.6. Define optimizer for both generator and discriminator according to Equations (5) and (6) with specified learning rates and hyperparameters;
- 3. Training QGAN (generator and discriminator);

## for each epoch do

for each batch in training data do

- a. Generate fake data using the quantum generator; b. Compute discriminator loss on both real and generated data, backpropagate and update
- discriminator;
- c. Generate fresh fake data and compute generator loss against the discriminator's feedback, backpropagate, and update generator;
- d. Record losses for analysis and model tuning;

## end

- e. Save generator and discriminator losses for the current epoch;

## end

- 4. Output: Generator and discriminator losses; visualize training progression with generator and discriminator losses;

with

$$\text{IS} \\ \text{kov} \\ \text{tion} \\ s \, X _ { t } \\ r t y \colon \\ \text{s} \, X _ { t } \\ \text{try} \colon \\ \text{end} & \text{$\omega$ had $K$ had $\omega$} \, \infty$$

and the kernel 𝐾 defined as

$$K \left ( \frac { x _ { n } - \hat { x } _ { i } } { h } \right ) = \frac { 1 } { \sqrt { 2 \pi } } \exp \left ( - \frac { 1 } { 2 } \left ( \frac { x _ { n } - \hat { x } _ { i } } { h } \right ) ^ { 2 } \right ). \quad \ \ ( 1 9 )$$

The bandwidth ℎ of the Gaussian estimation kernel is computed following Silverman's rule [32]:

$$h = \left ( \frac { 4 \hat { \sigma } ^ { 5 } } { 3 \hat { N } } \right ) ^ { \frac { 1 } { 3 } } \approx 1. 0 6 \hat { \sigma } ( \hat { N } - 1 ) ^ { - 1 / 5 }, \quad \quad ( 2 0 )$$

where ˆ is the standard deviation of the sample and 𝜎 𝑁 is the number of data points in our sample.

For empirical data analysis, the conditional probability can be represented by a transition matrix 𝑇 of dimension 𝑁 𝑠 × 𝑁 𝑠 with

Modeling stochastic eye tracking data: A comparison of quantum generative adversarial networks and Markov models GECCO '24, July 14-18, 2024, Melbourne, VIC, Australia entries

## 5 RESULTS

$$T _ { i j } = P r ( X _ { t = n + 1 } \in [ k _ { i }, k _ { i + 1 } ) | X _ { n } \in [ k _ { j }, k _ { j + 1 } ) ), \quad ( 2 1 ) \quad \ h e _ { m }.$$

with 𝑖, 𝑗 ∈ 𝑁,𝑖, 𝑗 ∈ [ 0 , 𝑁 𝑠 ) and 𝑘 𝑚 &gt; 𝑘 𝑛 for 𝑚 &gt; 𝑛 .

The efficacy of the proposed method is contingent upon three critical factors: (i) the robustness of the Markov assumption [31]; (ii) the number of states 𝑁 𝑠 , which directly influences computational efficiency and scales with 𝑁 4 𝑠 for bi-dimensional processes; and (iii) the sample size 𝑁 , which affects the bandwidth ℎ and thereby determines the spatial resolution of the model. We implemented this process using the open source code accessible via GitHub [22].

## 4 DATA AND THE STATISTICAL MEASURE

We used the Eye-tracking data 1 which was gathered at Oslo Metropolitan University utilizing the advanced eye-link Duo device, capable of reaching up to 2000 Hz but was adjusted to 200 Hz for this study. The measurements, recorded in screen pixels, were taken as participants searched for specific targets in images from the book 'Where's Waldo?' [13]. Each of the eight selected images was viewed for two minutes, a duration not expected to suffice for finding all targets but intended to keep participants focused. The data from eye-tracking measurements were preprocessed and utilized to train a QGAN. Initially, the velocity data for both left and right eyes are calculated by finding the Euclidean distance between consecutive position points and then dividing by the time interval, set at 1/200 seconds, to convert this distance into velocity. This velocity data is then aggregated into a structured format suitable for feeding into the discriminator. The data undergoes resampling to a fixed interval of 10 seconds, aggregating measurements by their mean within these intervals to reduce noise and temporal variability. The resultant dataset is then normalized using a MinMaxScaler to ensure the data fits within the required operational range of 0 and 1 [5]. The normalized data is transformed into sequences of a specified length (100 data points in this case), which are then used to create training batches. These sequences were fed into the discriminator, where it assessed them and learned to discern features that distinguish real data from the outputs generated by the quantum generator. This process aided the discriminator in enhancing its accuracy in identifying authentic samples.

The performance of the QGAN model is measured by the JensenShannon (JS) divergence [26]. This metric, a symmetrized variant of the Kullback-Leibler (KL) divergence [19], offers a symmetric distance measure between probability distributions. The JS divergence [35] is defined by

$$D _ { J S } ( P | Q ) = \frac { 1 } { 2 } D _ { K L } ( P | M ) + \frac { 1 } { 2 } D _ { K L } ( Q | M ), \quad \ \ ( 2 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where 𝑃 and 𝑄 are distributions, and 𝑀 = 1 2 ( 𝑃 + 𝑄 ) . The KL divergence is standard for assessing distributional similarity which enhances maximum likelihood estimates. The JS divergence, while preserving these properties, is more intuitive as it assesses the approximation of synthetic distributions to empirical ones. Therefore JS divergence in this case can be relevant for discriminators within QGANs to distinguish synthetic data from the generator.

1 All data collected was anonymized and follows the ethical requirements from the Norwegian Agency for Shared Services in Education and Research (SIKT), under the application with Ref. 129768.

The implementation of QGAN was done using PyTorch [28] and IBM's Qiskit [29] simulator 0.45. The discriminator's LSTM network is built with linear layers and a dropout mechanism distinguished between real and generated data sequences. The quantum generator is modeled by a parametrized quantum circuit and was optimized to replicate the complex distribution of the eye movement velocities. For the QGAN applied to eye-tracking data, the discretization and representation in quantum states was done using 3 and 4 qubits, each offering a different resolution of data representation. With 3 qubits, the velocity data is mapped onto 2 3 = 8 discrete levels, resulting in the quantum state:

$$G _ { \theta } \left | \psi _ { i n } \right \rangle = \sum _ { i = 0 } ^ { 7 } \sqrt { p _ { \theta } ^ { i } } | i \rangle,$$

where | 𝑖 ⟩ represents the 𝑖 -th discretized state of the velocity data. 𝑝 𝑖 𝜃 are the associated probabilities, learned through the training of the quantum generator as elaborated in Eq. (10). The training of our QGAN model was carried out over 300 epochs, with data batches of 500 samples each.

Figures 3 illustrate the histograms of log-transformed real versus generated data for three and four-qubit configurations across five layers. The use of log transformation ensures a focus on the distribution's dynamics rather than its absolute scale to facilitate a clearer understanding of the model's performance across varying complexities.

The evaluation of configurations of 3 and 4 qubits presents a nuanced understanding of the models' capabilities in reproducing real eye movement velocity data. Tab. 1 shows the statistical analysis for the generated data across different circuit layers including the Markovs model and the real data. The data generated by the QGANs exhibits a trend of increased skewness and kurtosis with increasing layers in quantum circuits, indicating a divergence from the real data's distribution especially in distribution tails. The generated distributions to the actual data are further quantified through the JSD values, as shown in Tab. 2. For 3 qubits QGAN performance is relatively stable across layers. For 4 qubits the second and third layers show optimal JSD, indicating minimal divergence points. Notably the Markov model exhibits superior performance with the lowest JSD, highlighting its efficiency in achieving minimal divergence compared to the QGAN-generated data.

The comparison to a mathematical Markov model, with its significantly lower JSD value, highlights the distinct challenges faced in quantum data generation. Figure 4 illustrates the comparative histograms of log-transformed eye movement velocities between real data and the generated data by a Markov model. The visual comparison and the lower JSD value clearly emphasize the effectiveness of Markov models in simulating complex data distributions closely resembling real data.

## 6 DISCUSSION AND CONCLUSIONS

The comprehensive analysis of QGANs alongside Markov models, particularly in modeling stochastic eye movement velocity, underscores the promising aspects of quantum computing in handling complex datasets. However, when we put side by side with the mathematical Markov model, our analysis paints a comprehensive

Figure 3: Comparative histograms of log-transformed real and generated eye movement velocity data across different circuit layers (1 - 5) for three and four-qubit QGAN. Each subplot illustrates the distribution of real data (in blue) with the corresponding generated data (in magenta) at respective circuit layers. The use of log transformation ensures a focus on the distribution's dynamics rather than its absolute scale to facilitate a clearer understanding of the model's performance across varying complexities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00673v1-figure-5c29360d-9fcf-43b5-a0d5-257dc151d94e.jpg)

Table 1: Comparison of generated data for 3 and 4 qubits with real data

Table 2: Jensen-Shannon divergence for 3 and 4 qubits across depths compared with a Markov model.

| Layers [3 Qubits]   | Mean       | Std Dev    | Skewness   | Kurtosis   |
|---------------------|------------|------------|------------|------------|
| 1 - Generated       | 7.1196e-05 | 0.00131    | 26.0450    | 933.951    |
| 2 - Generated       | 7.1196e-05 | 0.00147    | 30.7205    | 1253.51    |
| 3 - Generated       | 7.1196e-05 | 0.00120    | 41.7435    | 2740.014   |
| 4 - Generated       | 7.1196e-05 | 0.00180    | 54.2916    | 4773.553   |
| 5 - Generated       | 7.1196e-05 | 0.00183    | 60.87      | 6174.75    |
| Layers [4 Qubits]   | Mean       | Std Dev    | Skewness   | Kurtosis   |
| 1 - Generated       | 7.1196e-05 | 0.0012     | 22.964     | 674.157    |
| 2 - Generated       | 7.1196e-05 | 0.0013     | 28.720     | 1117.93    |
| 3 - Generated       | 7.1196e-05 | 0.0014     | 32.326     | 1743.243   |
| 4 - Generated       | 7.1196e-05 | 0.0015     | 37.681     | 2377.453   |
| 5 - Generated       | 7.1196e-05 | 0.0015     | 38.528     | 2641.427   |
| Real data           | 0.001735   | 0.04855    | 9.19       | 128.26     |
| Markov model        | 0.000371   | 0.00044752 | 3.493265   | 19.6143    |

picture of the present state of synthetic data generation via quantum methods.

| Layers (L)   | 3 Qubits   |    4 Qubits |
|--------------|------------|-------------|
| 1            | 0.08940    | 0.080354    |
| 2            | 0.08895    | 0.080198    |
| 3            | 0.08885    | 0.080153    |
| 4            | 0.08872    | 0.080196    |
| 5            | 0.08859    | 0.080102    |
| Markov model |            | 0.000340215 |

We have shown that by increasing the number of qubits as well as the number of layers, QGANs improve their performance in reproducing eye-gaze trajectories. However, despite QGANs' innovative nature and their capability to discern intricate eye gaze data patterns, Markov models outshine them in closely mirroring the actual data distributions.

This observation not only reinforces the previously identified challenges faced by AI in creating synthetic data-a topic already explored within classical GANs [21] research, but it also broadens

Figure 4: Comparative histograms of log-transformed eye movement velocities: real versus Markov model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00673v1-figure-4cf24b2a-1769-494d-be30-46e61fecc8b8.jpg)

this understanding to include the emerging domain of quantum computing. Such findings underscore the critical need for ongoing advancements in quantum algorithm development. The goal is clear: to overcome these obstacles and fully harness quantum computing's potential to produce synthetic datasets that closely resemble their real-world analogs.

## 7 ACKNOWLEDGMENTS

Weextend our sincere gratitude to Mr. Ramesh Uprety for his invaluable contributions and insightful discussions on the architecture of GANs, hyperparameter tuning strategies, and careful data-cleaning processes. We are deeply appreciative of his collaborative spirit and expert advice, which have significantly enriched this work. This work was funded by the Research Council of Norway under grant number 335940 for the project 'Virtual-Eye'.

## REFERENCES

- [1] Martin Arjovsky, Soumith Chintala, and Léon Bottou. 2017. Wasserstein GAN. arXiv:1701.07875 [stat.ML]
- [2] Frank Arute, Kunal Arya, Ryan Babbush, Dave Bacon, Joseph C. Bardin, Rami Barends, Rupak Biswas, Sergio Boixo, Fernando G. S. L. Brandao, David A. Buell, Brian Burkett, Yu Chen, Zijun Chen, Ben Chiaro, Roberto Collins, William Courtney, Andrew Dunsworth, Edward Farhi, Brooks Foxen, Austin Fowler, Craig Gidney, Marissa Giustina, Rob Graff, Keith Guerin, Steve Habegger, Matthew P. Harrigan, Michael J. Hartmann, Alan Ho, Markus Hoffmann, Trent Huang, Travis S. Humble, Sergei V. Isakov, Evan Jeffrey, Zhang Jiang, Dvir Kafri, Kostyantyn Kechedzhi, Julian Kelly, Paul V. Klimov, Sergey Knysh, Alexander Korotkov, Fedor Kostritsa, David Landhuis, Mike Lindmark, Erik Lucero, Dmitry Lyakh, Salvatore Mandrà, Jarrod R. McClean, Matthew McEwen, Anthony Megrant, Xiao Mi, Kristel Michielsen, Masoud Mohseni, Josh Mutus, Ofer Naaman, Matthew Neeley, Charles Neill, Murphy Yuezhen Niu, Eric Ostby, Andre Petukhov, John C. Platt, Chris Quintana, Eleanor G. Rieffel, Pedram Roushan, Nicholas C. Rubin, Daniel Sank, Kevin J. Satzinger, Vadim Smelyanskiy, Kevin J. Sung, Matthew D. Trevithick, Amit Vainsencher, Benjamin Villalonga, Theodore White, Z. Jamie Yao, Ping Yeh, Adam Zalcman, Hartmut Neven, and John M. Martinis. 2019. Quantum supremacy using a programmable superconducting processor. Nature 574, 7779 (01 10 2019), 505-510. https://doi.org/10.1038/s41586-019-1666-5
- [3] Laxmidhar Behera, Indrani Kar, and Avshalom C. Elitzur. 2005. A Recurrent Quantum Neural Network Model to Describe Eye Tracking of Moving Targets. Foundations of Physics Letters 18 (2005), 357-370. https://doi.org/10.1007/s10702005-7125-6
- [4] M. P. Cuéllar, M. C. Pegalajar, L. G. B. Ruiz, and C. Cano. 2023. Time Series Forecasting with Quantum Neural Networks. In Advances in Computational Intelligence , Ignacio Rojas, Gonzalo Joya, and Andreu Catala (Eds.). Springer Nature Switzerland, Cham, 666-677.
- [5] Lucas B.V. de Amorim, George D.C. Cavalcanti, and Rafael M.O. Cruz. 2023. The choice of scaling technique matters for classification performance. Applied Soft Computing 133 (2023), 109924. https://doi.org/10.1016/j.asoc.2022.109924
- [6] Edward Farhi and Hartmut Neven. 2018. Classification with Quantum Neural Networks on Near Term Processors. arXiv:1802.06002 [quant-ph]
- [7] William Fedus, Mihaela Rosca, Balaji Lakshminarayanan, Andrew M. Dai, Shakir Mohamed, and Ian Goodfellow. 2018. Many Paths to Equilibrium: GANs Do Not Need to Decrease a Divergence At Every Step. arXiv:1710.08446 [stat.ML]
- [8] Siddhant Garg and Goutham Ramakrishnan. 2020. Advances in Quantum Deep Learning: An Overview. arXiv:2005.04316 [quant-ph]
- [9] W. Gerych, K. Hickey, T. Hartvigsen, L. Buquicchio, A. Alajaji, K. Chandrasekaran, H. Mansoor, E. Agu, and E. Rundensteiner. 2023. Stabilizing Adversarial Training for Generative Networks. In 2023 IEEE International Conference on Big Data (BigData) . IEEE Computer Society, Los Alamitos, CA, USA, 5223-5232. https: //doi.org/10.1109/BigData59044.2023.10386654
- [10] Ian Goodfellow, Jean Pouget-Abadie, Mehdi Mirza, Bing Xu, David Warde-Farley, Sherjil Ozair, Aaron Courville, and Yoshua Bengio. 2014. Generative Adversarial Nets. In Advances in Neural Information Processing Systems , Z. Ghahramani, M. Welling, C. Cortes, N. Lawrence, and K.Q. Weinberger (Eds.), Vol. 27. Curran Associates, Inc., USA. https://proceedings.neurips.cc/paper\_files/paper/2014/ file/5ca3e9b122f61f8f06494c97b1afccf3-Paper.pdf
- [11] Ian Goodfellow, Jean Pouget-Abadie, Mehdi Mirza, Bing Xu, David Warde-Farley, Sherjil Ozair, Aaron Courville, and Yoshua Bengio. 2020. Generative adversarial networks. Commun. ACM 63, 11 (10 2020), 139-144. https://doi.org/10.1145/ 3422622
- [12] Ian J. Goodfellow, Jean Pouget-Abadie, Mehdi Mirza, Bing Xu, David WardeFarley, Sherjil Ozair, Aaron Courville, and Yoshua Bengio. 2014. Generative Adversarial Networks. arXiv:1406.2661 [stat.ML]
- [13] Martil Handford. 2007. Where's Waldo . Candlewick Press, Somerville, MA, USA.
- [14] Aram W. Harrow and John C. Napp. 2021. Low-Depth Gradient Measurements Can Improve Convergence in Variational Hybrid Quantum-Classical Algorithms. Phys. Rev. Lett. 126 (Apr 2021), 140502. Issue 14. https://doi.org/10. 1103/PhysRevLett.126.140502
- [15] Sepp Hochreiter and Jürgen Schmidhuber. 1997. Long Short-Term Memory. Neural Computation 9, 8 (11 1997), 1735-1780. https: //doi.org/10.1162/neco.1997.9.8.1735 arXiv:https://direct.mit.edu/neco/articlepdf/9/8/1735/813796/neco.1997.9.8.1735.pdf
- [16] Ling Hu, Shu-Hao Wu, Weizhou Cai, Yuwei Ma, Xianghao Mu, Yuan Xu, Haiyan Wang, Yipu Song, Dong-Ling Deng, Chang-Ling Zou, and Luyan Sun. 2019. Quantum generative adversarial learning in a superconducting quantum circuit. Science Advances 5, 1 (2019), eaav2761. https://doi.org/10.1126/sciadv.aav2761 arXiv:https://www.science.org/doi/pdf/10.1126/sciadv.aav2761
- [17] Diederik P. Kingma and Jimmy Ba. 2017. Adam: A Method for Stochastic Optimization. arXiv:1412.6980 [cs.LG]
- [18] Naveen Kodali, James Hays, Jacob D. Abernethy, and Zsolt Kira. 2018. On Convergence and Stability of GANs. https://api.semanticscholar.org/CorpusID: 37428828
- [19] Solomon Kullback and R. A. Leibler. 1951. On Information and Sufficiency. Annals of Mathematical Statistics 22, 1 (1951), 79-86. https://doi.org/10.1214/ aoms/1177729694
- [20] M. R. Pavan Kumar and Prabhu Jayagopal. 2020. Generative adversarial networks: a survey on applications and challenges. International Journal of Multimedia Information Retrieval 10 (2020), 1 - 24. https://api.semanticscholar.org/CorpusID: 226335714
- [21] Pedro Lencastre, Marit Gjersdal, Leonardo Rydin Gorjão, Anis Yazidi, and Pedro G. Lind. 2023. Modern AI versus century-old mathematical models: How far can we go with generative adversarial networks to reproduce stochastic processes? Physica D: Nonlinear Phenomena 453 (2023), 133831. https://doi.org/10.1016/j. physd.2023.133831
- [22] Pedro Lencastre, Marit Gjersdal, Leonardo Rydin Gorjão, Anis Yazidi, and Pedro G. Lind. 2023. Modern AI versus century-old mathematical models: How far can we go with generative adversarial networks to reproduce stochastic processes? Physica D: Nonlinear Phenomena 453 (2023), 133831. https://github.com/134f/ Physica-D-Markov-model
- [23] Yiyi Liao, Katja Schwarz, Lars Mescheder, and Andreas Geiger. 2020. Towards Unsupervised Learning of Generative Models for 3D Controllable Image Synthesis. arXiv:1912.05237 [cs.CV]
- [24] Jin-Guo Liu and Lei Wang. 2018. Differentiable learning of quantum circuit Born machines. Phys. Rev. A 98 (12 2018), 062324. Issue 6. https://doi.org/10.1103/ PhysRevA.98.062324
- [25] Sirui Lu, Lu-Ming Duan, and Dong-Ling Deng. 2020. Quantum adversarial machine learning. Phys. Rev. Res. 2 (8 2020), 033212. Issue 3. https://doi.org/10. 1103/PhysRevResearch.2.033212

- [26] Frank Nielsen. 2019. On the Jensen-Shannon Symmetrization of Distances Relying on Abstract Means. Entropy 21, 5 (2019), 485. https://doi.org/10.3390/ e21050485
- [27] Pekka Orponen. 1994. Computational complexity of neural networks: a survey.
- [28] Adam Paszke, Sam Gross, Francisco Massa, Adam Lerer, James Bradbury, Gregory Chanan, Trevor Killeen, Zeming Lin, Natalia Gimelshein, Luca Antiga, Alban Desmaison, Andreas Köpf, Edward Yang, Zach DeVito, Martin Raison, Alykhan Tejani, Sasank Chilamkurthy, Benoit Steiner, Lu Fang, Junjie Bai, and Soumith Chintala. 2019. PyTorch: An Imperative Style, High-Performance Deep Learning Library. arXiv:1912.01703 [cs.LG]
- [29] Qiskit contributors. 2023. Qiskit: An Open-source Framework for Quantum Computing. https://doi.org/10.5281/zenodo.2573505
- [30] Sashank J. Reddi, Satyen Kale, and Sanjiv Kumar. 2019. On the Convergence of Adam and Beyond. arXiv:1904.09237 [cs.LG]
- [31] Ana Russo, Frank Raischel, and Pedro G. Lind. 2013. Air quality prediction using optimal neural networks with stochastic variables. Atmospheric Environment 79 (2013), 822-830. https://doi.org/10.1016/j.atmosenv.2013.07.072
- [32] B.W. Silverman. 1998. Density Estimation for Statistics and Data Analysis (1st ed.). Routledge, New York, NY, USA. https://doi.org/10.1201/9781315140919
- [33] Haozhen Situ, Zhimin He, Yuyi Wang, Lvzhou Li, and Shenggen Zheng. 2020. Quantum generative adversarial network for generating discrete distribution. Information Sciences 538 (2020), 193-208. https://doi.org/10.1016/j.ins.2020.05.127
- [34] Samuel A. Stein, Betis Baheri, Daniel Chen, Ying Mao, Qiang Guan, Ang Li, Shuai Xu, and Caiwen Ding. 2022. QuClassi: A Hybrid Deep Neural Network Architecture based on Quantum State Fidelity. arXiv:2103.11307 [quant-ph]
- [35] Lilian Weng. 2019. From GAN to WGAN. arXiv:1904.08994 [cs.LG]
- [36] Nathan Wiebe, Ashish Kapoor, and Krysta M. Svore. 2015. Quantum Deep Learning. arXiv:1412.3489 [quant-ph]
- [37] Jinfeng Zeng, Yufeng Wu, Jin-Guo Liu, Lei Wang, and Jiangping Hu. 2019. Learning and inference on generative adversarial quantum circuits. Phys. Rev. A 99 (05 2019), 052306. Issue 5. https://doi.org/10.1103/PhysRevA.99.052306
- [38] C Zoufal, A Lucchi, and S Woerner. 2019. Quantum Generative Adversarial Networks for learning and loading random distributions. npj Quantum Information 5 (2019), 103. https://doi.org/10.1038/s41534-019-0223-2
Shailendra Bhandari Pedro Lincastre Pedro Lind
Modeling stochastic eye tracking data: A comparison of quantum generative adversarial networks and Markov models
We explore the use of quantum generative adversarial networks QGANs for
modeling eye movement velocity data. We assess whether the advanced
computational capabilities of QGANs can enhance the modeling of complex
stochastic distribution beyond the traditional mathematical models,
particularly the Markov model. The findings indicate that while QGANs
demonstrate potential in approximating complex distributions, the Markov model
consistently outperforms in accurately replicating the real data distribution.
This comparison underlines the challenges and avenues for refinement in time
series data generation using quantum computing techniques. It emphasizes the
need for further optimization of quantum models to better align with real-world
data characteristics.