## Pulse family optimization for parametrized quantum gates using spectral clustering

R.J.P.T. de Keijzer, 1, 2 J.A.C. Snijders, 2, 3 A. R. R. Carvalho, 3 and S.J.J.M.F. Kokkelmans 1, 2

1 Department of Applied Physics, Eindhoven University of Technology,

P. O. Box 513, 5600 MB Eindhoven, The Netherlands

2 Eindhoven Hendrik Casimir Institute, Eindhoven University of Technology,

P. O. Box 513, 5600 MB Eindhoven, The Netherlands ∗

3 Q-CTRL, Berlin, Germany

(Dated: August 2, 2024)

Parametrized gate circuits are used in plentiful applications in the current NISQ era of quantum computing. These parametrized gates are chiefly implemented using analytically found pulse protocols, often yielding suboptimal gate times, and consequently, fidelities. Alternatively, gate optimization algorithms are designed to construct high fidelity pulses for individual, fixed points in continuous parameter space. Gates for intermediate parameters can subsequently be found by some form of interpolation between previously constructed pulses. Nevertheless, it is not guaranteed (as with analytic protocols) that the pulses found by the optimization algorithms belong to the same family of solutions and thus show resemblance. Interpolation between two pulses of differing solution families often leads to high infidelities, as the pulse strays away from the minimum in the parameter/fidelity landscape. In this work, we introduce a spectral clustering method to sort high-fidelity, optimized pulses in families, and interpolating solely between pulses of the same family. Accordingly, interpolations will always approach maximal fidelity. Furthermore, as more than one pulse family is constructed, the parameter space can be partitioned according to which family prevails fidelity-wise. This work provides a meticulous demonstration of our constitutive continuous gate family construction by applying it to a universal gate set for Rydberg and Cat qubits under noise.

## INTRODUCTION

In order for a quantum computer to be able to perform all possible computations, it has to possess the ability to execute a universal gate set, consisting of all possible single qubit rotations and at least one entangling two-qubit operation [1]. The canonical single qubit gate set are the rotations gates over angle θ : R X ( θ , R ) Y ( θ , ) and R Z ( θ ), which form the basis for many variational quantum algorithms (VQAs) [2]. These parametrized gates are executed in physical systems using some control function on the qubits, i.e. laser pulses or electrical currents, specific to each individual parameter. Current era quantum computing systems are in the NISQ-era, where qubits are highly susceptible to noise [3]. This means that it is important to choose the optimal control function to mitigate noise, which is often not a standard analytical protocol, but some faster control pulse as constructed by a pulse optimization algorithm [4-9]. These variational methods often lead to faster, and consequently higher fidelity, pulses.

The problem with these optimization algorithms is that, unlike analytical protocols, they only prescribe pulses for one fixed point in the parameter space at a time, and not for the entire continuous set of parameters, which is especially relevant for VQAs [10-12], where copious different parameter realizations are necessary. Constructing control pulses every time a new unique parameter is required is far too computationally expensive, making it infeasible in the majority of VQAs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-9bc27dbb-a65a-47a9-8cff-deee0c65d2a6.jpg)

Figure 1. Graphical representation of the pulse interpolation problem. A pulse optimization algorithm might find one of several local minimizing pulses for θ i and θ i +1 . If these are corresponding local minima, the pulses belong to the same family (blue or pink) and interpolation yields a pulse for θ i +1 2 / with high fidelity. If the minima do not correspond, interpolation happens between pulses of different families, yielding a pulse for θ i +1 2 / with low fidelity (purple).

The task at hand is to extend a collection of solutions for discrete points in parameter space to the entire continuous set by some form of interpolation. However, this approach has an underlying complication. Assume a set of M subsequent, close points in the gate parameter space { θ i } M i for which one wishes to optimize the pulses. Initially, the pulse optimization algorithm will find a local minimizing pulse for θ 0 , and because the optimization problem is regular, there will exist a corresponding local minimizing pulse for θ 1 as long as the two values for θ are adequately close. In general, there will exist corresponding pulse solutions for each θ i , which we then call a pulse family , as introduced in Ref. [13]. Interpolating within a pulse family will yield pulses for intermediate values of θ , with comparable high fidelities. The problem, however, is that there is no guarantee that for every value of θ i , matching local minima will be found, and thus there is no guarantee that all the optimized pulses found will belong to a single family. Interpolating between two pulse families, will generally yield lower fidelities because the resulting pulses will not be in local minima of the pulse/fidelity landscape, see Fig. 1.

Previous work has been done on the construction of high fidelity continuous gate sets. References [4, 9, 14] construct pulses to mitigate hardware errors based on a parameterized Hamiltonian. These results are extended to the continuous case in [15]. Reference [16] is the first to mention the problem of similarity between pulses, but only related to single parameter optimization. References [13, 17-19] all devise methods to optimize continuous gate sets while maintaining pulse similarity across different parameters for interpolation purposes. In [13, 17], a neural network approach is used, where the networks are trained by randomly sampling the parameter space. The disadvantage of such a black box method is that verification of the resulting pulses is non-trivial. In [18], Trotterization is used to implement continuous parameter exponential matrices, but does not account for time optimization. The work in [19] utilizes Tikhonov regularization and feedforward to enforce pulse similarity, which is partly employed in our method. However, our work is the first to utilize clustering methods from graph theory [20, 21] to find multiple families of gates, which has three main advantages over all previous methods:

- - Foremost, partitioning the parameter space according to where certain pulse families outperform the others has the potential to yield exceedingly higher fidelities than a singular family approach;
- - Finding multiple families allows for the selection of pulses which are best suited to experimental procedures after verification on a setup;
- - The clustering method provides a more well-suited manner of interpolating pulses using Wasserstein-2 distances from optimal transport (see Sec. 3).

The layout of this paper is as follows. Section I describes optimal control methods used to construct optimal pulses for discrete parameter values θ i . Section II prescribes our pulse family clustering algorithm, its relation to other methods, as well as its advantages. Section III details the interpolation methods employed in this work. In Sec. IV, we show initial results for our approach, starting with results on the Wasserstein-2 distance to characterize pulse similarity in Sec. IV A. Section IV B illustrates the increases in fidelities our method has for interpolation. Finally, Sec. V summarizes and presents a future outlook.

## I. PULSE OPTIMIZATION

The application analyzed in this work is creating a universal parametrized gate set of R X ( θ , R ) Y ( θ , R ) Z ( θ ) and R ZZ ( θ ), where θ ∈ [0 , π ]. This will be performed on Rydberg qubits [22, 23], an architecture that recently has become exceedingly mature, and also on Cat qubits [24, 25], which recently are gaining traction within the quantum computing community. This illustrates the versatility of our methods. In this section, we briefly describe the optimization procedure for fixed parameters θ i , which has been well-described in literature before [2, 26, 27]. The evolution of the qubit system density matrix ρ t follows a Lindblad equation [28] of the form

$$\mathbb { m } ^ { \text{is} } \quad & \partial _ { t } \rho _ { t } = - i \left [ H _ { s y s } [ z ( t ) ], \rho _ { t } \right ] + \mathcal { L } ( \rho _ { t } ), \quad \rho ( 0 ) = \rho _ { 0 }, \\ \mathbb { m } ^ { \text{$\mathcal{L}$} } _ { t } \quad & \mathcal { L } ( \rho ) = \sum _ { k } \gamma _ { k } V _ { k } \rho V _ { k } - \frac { 1 } { 2 } \gamma _ { k } \{ V _ { k } ^ { \dagger } V _ { k }, \rho \}, \\ \mathbb { m } ^ { \text{$\mathcal{M}$} } _ { t } \quad & \mathcal { M }$$

where H sys [ z t ( )] is the system Hamiltonian controlled by user-defined, and optimizable pulses z t ( ), and L is the Lindblad operator responsible for the different sources of decoherence in the system. V k and γ k respectively are the jump operators with corresponding strengths, defining the Lindblad operator. The pulses are functions on [0 , T ], where T is the gate end time, bounded by experimental limitations as z min ≤ z t ( ) ≤ z max . The physical interpretation of the pulse type is system specific. In Rydberg systems, we can control the transitions between the | 0 , ⟩ | 1 ⟩ and | r ⟩ states on qubit j using coupling strengths Ω ab,j and detunings ∆ b,j on transition | a ⟩ ↔ | ⟩ b to get z t ( ) ∈ { Ω 01 ,j ( ) t , ∆ 1 ,j ( ) t , Ω 1 r,j ( ) t , ∆ r,j } [2]. Meanwhile, Cat qubits offer control on the single photon drive E j ( ), t the detuning ∆ ( ), interaction j t strength g t ( ), and the two photon drive G j ( ) [24]. t This gives z t ( ) ∈ { E j ( ) t , ∆ ( ) j t , G j ( ) t , g ( ) t } . For the Rydberg and the Cat qubit systems, we introduce characteristic timescales τ Ryd = 1 Ω / max and τ Cat = 1 /K respectively, where Ω max is a maximal coupling strength in the Rydberg system and K is the Kerr non-linearity in the Cat qubit system. Since our methods are agnostic to the qubit architecture and the physical implementation is not integral to understanding our results, we do not further discuss details on the specific systems here. For detailed Hamiltonians, Lindblad operators and a discussion on the pulse end times T , see App. A.

A solution to the Lindblad equation (1) for a specific set of pulses z t ( ) is described by a trace preserving, completely positive operator called a quantum channel E z , s.t. ρ T = E z ( ρ 0 ) [29]. For each specific gate U ∈ { R X ( θ i ) , R Y ( θ i ) , R Z ( θ i ), R ZZ ( θ i ) } to be created, the goal is to optimize the fidelity given by

$$\min _ { z } F ( \mathcal { E } _ { z }, \hat { U } ) \coloneqq \min _ { z } \frac { \sum _ { j } \text{Tr} \left ( U P _ { j } ^ { \dagger } U ^ { \dagger } \mathcal { E } _ { z } \left ( P _ { j } \right ) \right ) + d ^ { 2 } } { 4 ^ { N } d + d ^ { 2 } }, \ ( 2 )$$

as in [30]. Here N is the number of qubits and d is the dimension of the qubit system. For Rydberg systems d = 3 N and for Cat systems d = 2 N . P j are the Pauli matrices (for Rydberg extended to a 3-dimensional system) of which there are 4 N . Using the software package Boulder Opal, we optimize the pulses z t ( ) for the cost function in (2), using the package's standard convergence criteria [6, 31]. For all problems in this work, we optimize M = 20 equidistant predefined angles { θ i } M i in the parameter space [0 , π ].

## II. PULSE CLUSTERING

After having optimized the pulses for discrete fixed parameters { θ i } M i , they need to be partitioned in corresponding solution branches, or families. Here, we detail our method of clustering pulses in a fixed number of families using spectral clustering [21]. This method sorts the nodes V of a weighted graph G = ( V, E, w ) into a fixed number of clusters. E are the edges of the graph, and w are the weights on these edges encoding the similarity of two nodes. To apply this method to the discrete pulse sorting, we define a fully-connected graph, where each node i corresponds to one of the optimized pulses found for a specific parameter θ i , as shown in Fig. 2. The weights are set as w z , z ( 1 2 ) = (dist( z , z 1 2 ) + ϵ ) -1 , where dist is some distance on the space of pulses and ϵ = 10 -4 is a small regularization parameter. As a basic approach, the L 2 distance [32] can be used as

$$L _ { 2 } ( z _ { 1 }, z _ { 2 } ) = \int _ { 0 } ^ { T } | z _ { 1 } ( t ) - z _ { 2 } ( t ) | ^ { 2 } d t.$$

This distance function is often implemented as it is efficient to compute. However, it does not fully serve our needs, since it does not capture shifts on the time axis well. The L 2 distance only compares pulse values pointwise in time (see Fig. 3, where the highest peak shifts, causing a big L 2 distance between the pulses even though they are quite similar in shape). As long as two pulses do not overlap, a large time shift causes the same maximal L 2 distance as a small shift. To remedy this issue, the Wasserstein-2 distance [33] can be used. This distance has its origins in the mathematical framework of optimal transport and defines the cost of moving one distribution onto another, which is more closely related to our intuitive likeness of pulses, see Fig. 3. The important point

Figure 2. To distinguish between different pulse families, pulses are represented as nodes in a graph with the edge weighs between them equal to some distance function d . The spectral clustering algorithm takes as input the distance matrix (right), and creates clusters minimizing the distances between pulses of the same cluster.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-eceb579d-67bb-48fc-8230-9a119c4c5998.jpg)

Figure 3. Example of Wasserstein-2 vs. L 2 distance. a) If two pulses are compared in L 2 distance, small or large time shifts both can result in maximal dissimilarity. The Wasserstein2 distance is able to capture these shifts because it tries to transport one pulse onto the other. b) Example of a coupling γ ∈ Π between pulses z 1 and z 2 , the Wasserstein-2 distance aims to construct a minimal coupling in the sense of (B1).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-d996ca1e-1975-4556-8168-191c4a23299f.jpg)

to understand is that the Wasserstein-2 distance treats the pulses as distributions on the plane R 2 and calculates the cost of moving one onto the other and defines this as the distance [34]. The exact definition of the distance is rather mathematically involved, and therefore deferred to App. B [35].

Using one of the previously defined distances, a similarity matrix consisting of the edge weights w can be constructed, see Fig. 2. This serves as the input for the spectral clustering algorithm. As the spectral clustering algorithm requires a fixed number of clusters, we determine the best possible number of clusters by means of the Elbow method [36]. This heuristic determines the point at which an additional cluster does not capture the pulse differences better, and over-fitting starts. For this work this value is always found to be 3 clusters, but in more complex problems with high dimensional parameters, could definitely be expected to increase.

Figure 4. A family of pulses is found, defined at a subset of the predefined angles { θ i } M i , and is to be extended to all predefined angles. First, ansatzes are created by interpolation if the unassigned angle is within the range of assigned angles, or else set as the highest or lowest found pulse. These ansatzes are optimized using Boulder Opal to extend the family to all predefined angles θ i . Pulses can then be interpolated within a family in order to find high fidelity pulses for the entire parameter space.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-c0430af3-b315-4d9e-a0b1-b1707e417f2c.jpg)

## III. PULSE INTERPOLATION

After the originally optimized pulses have been partitioned into families by the clustering algorithm, the following step is to extend the families to the entire parameter space by means of interpolation. In this work, the parameter space will always be θ ∈ [0 , π ] so that we can refer to angles. For multiparameter gates, an analogous approach can be pursued.

For each family, we want to first construct pulses on all of { θ i } M i , instead of only the subset assigned to it in the clustering. Consider one of the families, there will be a lowest and highest angle assigned to this family in the clustering. Note that a family does not necessarily include only subsequent angles, but might include gaps between disjoint regions. If an angle θ i is not assigned to this family and is in such as gap, we take a linear interpolation as an ansatz and optimize the pulses using Boulder Opal. For fixed angles θ i outside this range, we simply use the lowest or highest angle pulse as an ansatz and optimize. Employing these ansatzes, hopefully a pulse is found for that parameter belonging to the same family (reminiscent of Tikhonov regularization [19]). This leaves us with several extended families of pulses Z j on all predefined parameters { θ i } M i (see Fig. 4).

Using linear interpolation within a single family, new pulses for all parameters in the parameter space can be constructed, finalizing our construction of a continuous gate set for each family. Linear interpolation is computationally very efficient, especially when compared to constructing pulses from the ground up. By interpolating within a family, we expect higher fidelities than if we were to interpolate on the original set of pulses.

.

## IV. RESULTS

## A. Distance function comparison

This section contrasts the effectiveness of clustering families for the Wasserstein-2 and L 2 distances. For arbitrarily optimized pulses, there is no established family structure to validate against. In order to straightforwardly compare the distance measures, we construct J mock families of pulses Z j = { z θ j,i } M i . The pulses within one family are constructed in such a way that pulses z θ i and z θ j look alike for i close to j , as they would in actual applications. The exact details of this are left for App. C. For each angle θ i , a corresponding pulse z true ,i from a family Z true ,i ∈ { Z , ..., Z 1 J } is randomly picked for the final pulses. This mimics the pulses found by the Boulder Opal optimization procedure and accrues a set of pulses as in Fig. 5b, of which we have knowledge on the family structure.

On these mock pulses, the distance metrics are compared by constructing a distance matrix as in Fig. 5, and having the spectral clustering algorithm assign to each pulse i a family Z dist,i ∈ { Z , ...Z 1 J } . The conditional probability P dist of finding a correct match is approximated as

$$\text{$ref$ will} \\ \text{$\text{ily in}$} \\ \text{sansify} \\ \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text                                                                                                                                                                                                         }                                                                                                                                                                                                       
                                                                                                                                                                                                        $                                                                                                                                                                                                       .                                                                                                                                                                                                        -                                                                                                                                                                                                        +                                                                                                                                                                                                        )                                                                                                                                                                                                       ,                                                                                                                                                                                                        "                                                                                                                                                                                                       '                                                                                                                                                                                                       (                                                                                                                                                                                                        ;                                                                                                                                                                                                        :$$

✶ The reason for choosing this figure of merit is that when two pulses are found to be matching, they only lead to high fidelity as long as they belong to the same family. Missing a match because two pulses from the same family are clustered in different families will not lead to faulty interpolations, and is thus less important from a fidelity perspective (nevertheless, the Wasserstein-2 distance also outperforms the L 2 distance in terms of missed matches). Figure 5 shows the comparison between the two distances. Wasserstein-2 can be seen to largely outperform the standard L 2 distance, as it is able to capture the pulse family characteristics better due to its transport based nature. In the rest of this work, the Wasserstein-2 distance is employed for clustering.

Figure 5. Comparison of the Wasserstein-2 and L 2 distance for pulse clustering. a) Box plots of the probability of finding a correct match P (cor. fnd.) | for both distances, data from 2000 trials with M = 20 angles in the interval [0 , π ] randomly selected from three pulse families. b) Example of the authentic clustering of optimized pulse solutions. c) Clusters found by the Wasserstein-2 distance for the pulses from b). d) Clusters found by the L 2 distance for the pulses from b).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-d6d0da8a-e60d-40c1-a478-3c568508ca49.jpg)

## B. Interpolation of gates

pulses from the families in the highlighted areas correspond well and thus result in high fidelities.

To test our clustering method for the interpolation of pulses, we first optimize for M = 20 equidistant angles θ i ∈ [0 , π ] for R X ( θ , R ) Y ( θ , R ) Z ( θ ) and R ZZ ( θ ) for both Rydberg and cat qubits. This allows us to gather results on a universal set of gates on two vastly different systems, showing the versatility and robustness of our method.

The pulses are clustered using our spectral clustering method from Sec. II, where the Elbow method indicates that for all considered applications 3 clusters is the optimal number. The families found are extended to the entire M = 20 predefined angles { θ i } M i using the methods from Sec. III. Lastly, for each family, we interpolate the pulses halfway the predefined angles, and determine their fidelity using (2).

Figure 6 shows the results for the R Y ( θ ) gate on cat qubits. Generally, the R Y gate is the most interesting of the three qubit gates because it requires both amplitude and phase control. In Fig. 6a, we see that when optimizing using Boulder Opal, both the original pulses and the clustered families reach high fidelities, but the interpolation on the original pulses is very bad for almost all angles. For the clustered families, on the other hand, those with good interpolations are found for the entire parameter space. This is highlighted in Fig. 6b, where the parameter space is partitioned according to where each family performs best, leading to great improvements in the interpolation fidelities compared to the original pulses. Figures. 6c and 6d further illustrate this by showing that the original pulses are dissimilar for subsequent angles, resulting in low fidelities, whereas the

Similar results can be seen in Fig. 7 for the 2-qubit R ZZ ( θ ) on a Rydberg system. Here we see good interpolations of the original pulses around θ = π/ 2, which is also retrieved for one of the families. However, for low and high angles the interpolation fidelities become quite low for the original pulses (likely due to the presence of many local minima), but high fidelity pulses ones are constructed for the clustered families.

Figure 8 shows all fidelities of interpolated pulses for the problems considered. Across all parametrized gates for both qubit types, there is an increase of half up to a full order of magnitude in fidelity. Partitioning the parameter space with specific families (Clustered, pink results) results in a big advantage in both mean and variance of the fidelities compared to individual clusters (purple). Here individual clusters would be comparable to a Tikhonov regularization method as in [19], since ansatzes for pulse optimization are used. Along the same lines, especially note the much smaller variances of the fidelities indicating that our partition clustered method yields good pulses over the entire parameter space instead of only in a specific region. The individual family method definitely results in good interpolations for certain regions of parameter space, but multiple families seem to be necessary to patch together a faithful continuous gate set over the entire parameter space. This hails especially true for the cat qubits, as seen from the fidelities in Fig. 8, which is a more complex optimization problem where likely many local minima in parameter/fidelity space exist.

Figure 6. Optimization and interpolation results for R Y ( θ ) gate on a cat qubit system. a) Black, original pulses and interpolations. Colored pulses show optimized and extended clusters plus interpolations, showing matching or improved interpolation infidelities. Vertical lines are min and max angle of each cluster in the original set of pulses. b) Narrowed down results from a) where the best clusters for each regime of angles are shown, highlighting improvements in fidelity. c) Single photon drive E t ( ) for the original set of pulses, showing a jagged interpolation landscape leading to low fidelities. Horizontal lines indicate chosen regime for each cluster in d). d) Clustered pulses for the chosen regimes, showing more regular landscape, necessary for faithful interpolations. Individual clusters are somewhat displaced on z -axis for visibility.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-bd01708d-1cfd-4272-9979-28f2fb1e892e.jpg)

Figure 7. Optimization and interpolation results for a R ZZ ( θ ) gate on a Rydberg neutral atom qubit system. a) Black, original pulses and interpolations, for which interpolations have low fidelities for small and large angles. Colored pulses show optimized and extended clusters plus interpolations, showing matching or improved interpolation infidelities. Vertical lines are min and max angle of each cluster in the original set of pulses. b) Narrowed down results from a) where the best clusters for each regime of angles is chosen, highlighting improvements in fidelity. c) Coupling strength of | 0 ⟩ ↔ | 1 ⟩ transition on qubit 0 Ω 01 0 , ( ) for t the original set of pulses, showing a jagged interpolation landscape leading to low fidelities. Horizontal lines indicate chosen regime for each cluster in d). d) Clustered pulses for the chosen regimes, showing more regular landscape, necessary for faithful interpolations. Individual clusters are somewhat displaced on z -axis for visibility.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-37743e16-9d68-4561-ae85-13ddb6eb7f02.jpg)

.

.

Figure 8. Average and standard variations of infidelities of interpolations halfway between optimized angles, as in Figs. 6 and 7. Non-clustered pulses are the original pulses (cyan), three clustered are shown (purple) as well as the fully clustered method where the parameter space is split up in regimes where the best cluster is assigned (pink). Note the higher mean fidelities as well as much lower standard variations, indicating better interpolations over the entire parameter space. a) Rydberg gates, b) Cat qubit gates.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-31d36391-d105-4072-af58-21a7f802c0ee.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00119v1-figure-7026cec5-0008-43dd-b3f6-65e4581ba826.jpg)

## V. CONCLUSION

This work discusses and analyzes a new method for constructing continuous parameter gate sets based on the clustering of pulse families. This is an important problem as for virtually all NISQ era applications parametrized gates are required, and optimizing for each individual parameter necessary in a problem is extremely computationally expensive and slow. In our method, spectral clustering using a similarity measure based on Wasserstein-2 distances results in multiple families of pulses on which inter-family interpolation can be performed. The advantage of this method is that multiple families of pulses are found from which the best ones can be selected for a specific regime of parameter space. This method will prove to be useful for experimental realization, where one family is possibly easier to implement than another. On all gates and on both qubit systems considered, our method leads to significantly better fidelities and much more consistent interpolations than would be the case for original pulse optimization or Tikhonov based feedforward methods (individual family). In this, we highlight the great performance in cat qubits system, which shows the importance of our method in increasingly more important complex qubit schemes.

In future work, we hope to extend these methods to multiple parameter gates, which might be implemented in more complex VQA or QAOA type problems. Furthermore, we are curious to see the influence of directly optimizing the fidelity via the quantum channel, thus optimizing more directly for the losses, instead of optimizing the unitary evolution and post-processing the influences of the Lindbladian terms. Different parameters have different optimal gate end times, we wish to extend our methods to accommodate for varying gate times instead of one fixed optimal gate time for the largest parameter value ( θ = π in this work). Lastly, we hypothesize that a more sophisticated way of interpolation, potentially based on Wasserstein-2 transport or on the Hamiltonian dynamics, will yield better fidelities for the interpolated pulses.

## ACKNOWLEDGEMENTS

We thank Klaas Wijnsma, Jasper Postema, Julian Teske, Kisa Barkemeyer, Kara Maller, and Oliver Tse for fruitful discussions. This research is financially supported by the Dutch Ministry of Economic Affairs and Climate Policy (EZK), as part of the Quantum Delta NL program, the Horizon Europe programme HORIZON-CL4-2021DIGITAL-EMERGING-01-30 via the project 101070144 (EuRyQa), and by the Netherlands Organisation for Scientific Research (NWO) under Grant No. 680.92.18.05.

## COMPETING INTERESTS

The authors declare no competing interests.

## DATA AVAILABILITY

The data supporting the findings are available from the corresponding author upon reasonable request.

## CODE AVAILABILITY

The code supporting the findings are available from the corresponding author upon reasonable request.

- [1] A. B. D. E. Deutsch and A. Ekert, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences 449 , 669-677 (1995).
- [2] R. J. P. T. de Keijzer, V. E. Colussi, B. Skori´, ˇ c and S. J. J. M. F. Kokkelmans, 'Optimization of the variational quantum eigensolver for quantum chemistry applications,' (2021), arXiv:2102.01781 [quant-ph].
- [3] J. Preskill, Quantum 2 , 79 (2018).
- [4] M. Mohan, R. de Keijzer, and S. Kokkelmans, Phys. Rev. Res. 5 , 033052 (2023).
- [5] S. Jandura and G. Pupillo, Quantum 6 , 712 (2022).
- [6] H. Ball, M. J. Biercuk, A. R. R. Carvalho, J. Chen, M. Hush, L. A. D. Castro, L. Li, P. J. Liebermann, H. J. Slatyer, C. Edmunds, V. Frey, C. Hempel, and A. Milne, Quantum Science and Technology 6 , 044011 (2021).
- [7] E. Knill, D. Leibfried, R. Reichle, J. Britton, R. B. Blakestad, J. D. Jost, C. Langer, R. Ozeri, S. Seidelin, and D. J. Wineland, Phys. Rev. A 77 , 012307 (2008).
- [8] T. Caneva, T. Calarco, and S. Montangero, Phys. Rev. A 84 , 022326 (2011).
- [9] R. de Keijzer, O. Tse, and S. Kokkelmans, Quantum 7 , 908 (2023).
- [10] N. Lacroix, C. Hellings, C. K. Andersen, A. Di Paolo, A. Remm, S. Lazar, S. Krinner, G. J. Norris, M. Gabureac, J. Heinsoo, A. Blais, C. Eichler, and A. Wallraff, PRX Quantum 1 (2020), 10.1103/prxquantum.1.020304.
- [11] C. Grange, M. Poss, and E. Bourreau, 4OR 21 , 363 (2023).
- [12] M. Cerezo, A. Arrasmith, R. Babbush, S. C. Benjamin, S. Endo, K. Fujii, J. R. McClean, K. Mitarai, X. Yuan, L. Cincio, and P. J. Coles, Nature Reviews Physics 3 , 625 (2021).
- [13] F. Sauvage and F. Mintert, Phys. Rev. Lett. 129 , 050507 (2022).
- [14] P. Luchi, F. Turro, S. Quaglioni, X. Wu, V. Amitrano, K. Wendt, J. L. DuBois, and F. Pederiva, The European Physical Journal A 59 , 196 (2023).
- [15] J.-S. Li, J. Ruths, T.-Y. Yu, H. Arthanari, and G. Wagner, Proceedings of the National Academy of Sciences 108 , 1879 (2011), https://www.pnas.org/doi/pdf/10.1073/pnas.1009797108.
- [16] Y. Shi, A. R. Castelli, X. Wu, I. Joseph, V. Geyko, F. R. Graziani, S. B. Libby, J. B. Parker, Y. J. Rosen, L. A. Martinez, and J. L. DuBois, Phys. Rev. A 103 , 062608 (2021).
- [17] F. Preti, T. Calarco, and F. Motzoi, PRX Quantum 3 , 040311 (2022).
- [18] M. Schilling, F. Preti, M. M. Muller, T. Calarco, and F. Motzoi (2024).
- [19] J. D. Chadwick and F. T. Chong, in 2023 IEEE International Conference on Quantum Computing and Engineering (QCE) (IEEE Computer Society, Los Alamitos, CA, USA, 2023) pp. 1286-1294.
- [20] L. Rokach and O. Maimon, 'Clustering methods,' in Data Mining and Knowledge Discovery Handbook , edited by O. Maimon and L. Rokach (Springer US, Boston, MA, 2005) pp. 321-352.
- [21] A. Ng, M. Jordan, and Y. Weiss, Advances in neural information processing systems 14 (2001).
- [22] M. Saffman, Journal of Physics B: Atomic, Molecular and Optical Physics 49 , 202001 (2016).
- [23] M. Morgado and S. Whitlock, AVS Quantum Science 3 , 023501 (2021), https://doi.org/10.1116/5.0036562.
- [24] P. Vikst˚l, a L. Garc´ ıa- ´ lvarez, A S. Puri, and G. Ferrini, 'Quantum approximate optimization algorithm with cat qubits,' (2023), arXiv:2305.05556 [quant-ph].
- [25] A. Grimm, N. E. Frattini, S. Puri, S. O. Mundhada, S. Touzard, M. Mirrahimi, S. M. Girvin, S. Shankar, and M. H. Devoret, Nature 584 , 205 (2020).
- [26] E. S. Matekole, Y.-L. L. Fang, and M. Lin, in Quantum 2.0 Conference and Exhibition (Optica Publishing Group, 2022) p. QM3A.2.
- [27] A. R. R. Carvalho, H. Ball, M. J. Biercuk, M. R. Hush, and F. Thomsen, Phys. Rev. Appl. 15 , 064054 (2021).
- [28] C. A. Brasil, F. F. Fanchini, and R. d. J. Napolitano, Revista Brasileira de Ensino de F´ ısica 35 , 01 (2013).
- [29] D. Kretschmann, in Encyclopedia of Mathematical Physics , edited by J.-P. Fran¸ coise, G. L. Naber, and T. S. Tsun (Academic Press, Oxford, 2006) pp. 424-431.
- [30] M. A. Nielsen, Physics Letters A 303 , 249 (2002).
- [31] Q-CTRL, 'Boulder Opal,' https://q-ctrl.com/boulderopal (2023), [Online].
- [32] K. Yosida, Functional analysis , Vol. 123 (Springer Science &amp; Business Media, 2012).
- [33] F. Santambrogio, Birk¨user, NY a 55 , 94 (2015).
- [34] C. Villani, Topics in optimal transportation , Vol. 58 (American Mathematical Soc., 2021).
- [35] For the Wasserstein-2 distance, pulses are discretized as point clouds with points equidistant along the curve, and Sinkhorn-Knopp [37] is used to identify the minimizing coupling.
- [36] M. Cui et al. , Accounting, Auditing and Finance 1 , 5 (2020).
- [37] P. A. Knight, SIAM Journal on Matrix Analysis and Applications 30 , 261 (2008).
- [38] S. Deffner and S. Campbell, Journal of Physics A: Mathematical and Theoretical 50 , 453001 (2017).
- [39] C. S. Adams, J. D. Pritchard, and J. P. Shaffer, Journal of Physics B: Atomic, Molecular and Optical Physics 53 , 012002 (2019).
- [40] J. Zha, Z. Qin, J. Yan, N. Cao, Q. Wei, and P. Wang, Results in Physics 22 , 103594 (2021).
- [41] H.-g. Lee, Y. Song, and J. Ahn, Phys. Rev. A 96 , 012326 (2017).
- [42] X.-F. Shi and Y. Lu, Phys. Rev. A 104 , 012615 (2021).
- [43] A. Pagano, S. Weber, D. Jaschke, T. Pfau, F. Meinert, S. Montangero, and H. P. B¨chler, Phys. Rev. Res. u 4 , 033019 (2022).
- [44] M. Cosacchi, T. Seidelmann, J. Wiercinski, M. Cygorek, A. Vagov, D. E. Reiter, and V. M. Axt, Phys. Rev. Res. 3 , 023088 (2021).
- [45] H. Goto, Phys. Rev. A 93 , 050301 (2016).

In this appendix we describe the qubit systems resulting in the system Hamiltonians H sys , jump operators V k and control pulses z t ( ). Furthermore, we describe the pulse optimization procedure in more detail.

In order to restrict the parameter space and consequently enhance the likelihood of yielding pulses in the same family, we set the pulses for qubit 1 equal to those of qubit 0 in the R ZZ two qubit case. All pulses are discretized as 30 segment piecewise constant functions, which are smoothened via a Gaussian kernel and are equated to 0 at their beginning and end values, lastly being put through a low-pass filter. For each problem analyzed, the pulse end time T is determined as the end time for which the problem at θ = π has the highest fidelity. The existence of such an optimal time follows from the reasoning that for too large gate times, decoherence effects dominate and decrease fidelity. On the other hand, if the gate time is too small, the quantum speed limit is not reached [38] and there is no way of constructing a minimizing pulse, leading to low fidelities. The quantum speed limit likely has the highest value for θ = π because it least resembles the identity. To ensure we reach the quantum speed limit for all θ , we choose T for all problems as the optimal time for θ = π .

## 1. Rydberg neutral atoms

This section introduces basic Rydberg physics to identify what control pulses can look like for this system, as discussed in Sec. I. In general, this will yield the Hamiltonian and possible Lindblad operators for the systems, as in (1). A Rydberg system is a system of individual neutral atoms trapped in laser optical tweezers, where the electronic states encode for the qubit manifold [22]. For this work, we assume gg qubits [23] such that for each qubit we have the states | 0 ⟩ , | 1 ⟩ and | r ⟩ , leading to find d = 3 N as in (2).

To perform single qubit transitions on qubit j between states | a ⟩ and | b ⟩ , a laser interacts with the atom to realize the Hamiltonian [22, 23, 39]

$$H _ { j } ^ { a b } = \frac { \Omega _ { a b, j } ( t ) } { 2 } \left ( e ^ { i \varphi _ { a b, j } ( t ) } | a \rangle _ { j } \langle b | _ { j } + e ^ { - i \varphi _ { a b, j } ( t ) } | b \rangle _ { j } \langle a | _ { j } \right ) - \Delta _ { b, j } ( t ) | b \rangle _ { j } \langle b | _ { j }.$$

Here, Ω ab,j ( ) denotes the coupling strength, t φ ab,j ( ) the phase of the laser coupled to atom t j , and ∆ b,j ( ) = t ω ab,j ( ) t -˜ ω ab the detuning of the laser frequency ω ab,j ( ) from the energy level difference ˜ t ω ab . In current Rydberg systems, one has exquisite control over Ω and ∆, and less over φ [40], thus we set φ = 0. For our systems, we assume transitions | 0 ⟩ ↔ | 1 ⟩ and | 1 ⟩ ↔ | r . ⟩ This renders control pulses z t ( ) ∈ { Ω 01 ,j ( ) t , ∆ 1 ,j ( ) t , Ω 1 r,j ( ) t , ∆ r,j } . Notice that having both coupling and detuning allows for full control on the Bloch sphere of each individual qubit, allowing for rotational control [41]. For this reason (and to restrict parameter space), we only optimize for Ω 01 and ∆ 1 for the single qubit gates and set the other pulses equal to z t ( ) = 0.

The Rydberg states | r ⟩ are high-energy excited states that have a passive 'always-on' interaction, which is described by a non-controllable Hamiltonian H V [23] as a Van der Waals interaction (VdW) [42]

$$H _ { V } = \sum _ { i = 1 } ^ { m } \sum _ { j = 1 } ^ { m } \frac { C _ { 6 } } { R _ { i j } ^ { 6 } } | r r \rangle _ { i j } \langle r r | _ { i j },$$

where R ij is the interatomic distance and C 6 is an interaction coefficient. For the total Hamiltonian H sys

$$H _ { s y s } = \sum _ { j } H _ { j } ^ { 0 1 } + H _ { j } ^ { 1 r } + H _ { V }.$$

Furthermore, the Rydberg states have a finite lifetime before decaying to | 1 . ⟩ For completion, we also introduce decay from | 1 ⟩ to | 0 , which results in the jump operators ⟩

$$V _ { 1, j } = | 0 \rangle _ { j } \langle 1 | _ { j }, \ V _ { r, j } = | 1 \rangle _ { j } \langle r | _ { j }.$$

In the optimization procedure, Ω max = 10MHz defines the timescale, and the pulses are constrained as | z t ( ) | ≤ Ω max . Furthermore, we set γ , γ 1 r = 0 01Ω . max which will incur fidelity losses. Lastly, C /R 6 6 = 10 3 · Ω max is taken (we only consider up to two qubits in this work). The choices for physical parameters are inspired by [43] and summarized in Table I.

Table I. Optimizable pulses parameters used in the optimizations. The minima and maxima values of each drive are inspired by [43] The columns for each gate denote if the pulse values are real ( R ), complex ( C ) or kept constant (const.). When kept constant, the values are equal to 0. All units in [Ω max ].

| Pulse z ( t )   | min[ z ( t )]   | max[ z ( t )]   | z (0) = z ( T )   | R Z ( ϕ )   | R X ( θ )   | R Y ( φ )   | R ZZ (Θ)   |
|-----------------|-----------------|-----------------|-------------------|-------------|-------------|-------------|------------|
| Ω 01 ( t )      | -1.0            | 1.0             | 0.0               | R           | R           | R           | R          |
| ∆ 1 ( t )       | -1.0            | 1.0             | 0.0               | R           | R           | R           | R          |
| Ω 1 r ( t )     | -1.0            | 1.0             | 0.0               | const.      | const.      | const.      | R          |
| ∆ r ( t )       | -1.0            | 1.0             | 0.0               | const.      | const.      | const.      | R          |
|                 |                 |                 | T [1 / Ω max ]    | 1.36        | 0.73        | 1.36        | 2.83       |

## 2. Cat qubit Hamiltonians

Our discussion of cat qubits largely follows the explanations from [24]. The cat qubit can be realized in a Kerr parametric oscillator with a two-photon pump [25, 44, 45]. In the rotating frame of the two-photon pump frequency ω p , the main system is described by

$$\hat { H } _ { 1 } = - K \hat { a } ^ { \dagger 2 } \hat { a } ^ { 2 } + G _ { j } ( t ) \left ( \hat { a } ^ { \dagger 2 } + \hat { a } ^ { 2 } \right )$$

with K the Kerr non-linearity which is constant and chosen as our timescale, G j the two-photon drive amplitude, ˆ a the annihilation operator and ˆ a † the creation operator. Rewriting the system Hamiltonian gives

$$H _ { 1, j } = - K \left ( \hat { a } ^ { \dagger 2 } - \frac { G _ { j } ( t ) } { K } \right ) \left ( \hat { a } ^ { 2 } - \frac { G _ { j } ( t ) } { K } \right ) + \frac { G _ { j } ( t ) ^ { 2 } } { K }.$$

With the knowledge that ˆ a |± α ⟩ = ± |± α α ⟩ one can see that the degenerate eigenstates of the Hamiltonian are the coherent states | ± α ⟩ with eigenenergy G /K 2 . Since these eigenstates are degenerate, a linear combination

$$| C _ { \alpha } ^ { + } \rangle = N _ { \pm } ( | \alpha \rangle \pm | - \alpha \rangle )$$

of these states is also an eigenstate. The states from (A1) are the so-called cat states, with N ± a normalization constant. The orthogonality of the cat states from (A1) allows for the following encoding

$$| 0 \rangle = \frac { | C _ { \alpha } ^ { + } \rangle + | C _ { \alpha } ^ { - } \rangle } { \sqrt { 2 } }, \ \ | 1 \rangle = \frac { | C _ { \alpha } ^ { + } \rangle - | C _ { \alpha } ^ { - } \rangle } { \sqrt { 2 } },$$

as our computational basis, with | 0 ⟩ and | 1 ⟩ denoting the computational basis states. The corresponding Pauli matrices are defined through the computational basis states, e.g. Y = i | 0 ⟩⟨ 1 | -i | 1 ⟩⟨ 0 . | The computational basis is chosen this way, foremost because single photon-loss (the dominant error mechanism in the system) is exponentially suppressed in α . Because we can not simulate the entire Fock space, we truncate its dimension at 20 states for single qubit gates and at 16 for R ZZ .

The control of these systems is described through the following Hamiltonians

$$H _ { 2, j } = - \Delta _ { j } ( t ) \hat { a } _ { j } ^ { \dagger } \hat { a } _ { j } + E _ { j } ( t ) ( \hat { a } _ { j } e ^ { - i \theta _ { j } } + \hat { a } _ { j } ^ { \dagger } e ^ { i \theta _ { j } } ), \quad H _ { 3 } = g ( t ) ( \hat { a } _ { 1 } \hat { a } _ { 2 } ^ { \dagger } + \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 2 } ),$$

where a detuning ∆ ( ) is created between the two-photon drive and the resonator ∆ = j t ω r -2 ω p , with ω r and ω p the resonator and two-photon pump frequency respectively. θ j the phase of the drive, which is chosen at θ j = 0 (except for R Y ). E j ( ) is t the amplitude of the single-photon drive. g t ( ) is a two-photon exchange between the two resonators. When | g t ( ) | , | E j ( ) t | , | ∆ ( ) j t | ≪ G the computational states are approximately kept in the computational basis. Together with single photon loss, dynamics to outside the qubit manifold is the main source of error. This leaves us with the pulses z t ( ) ∈ { E j ( ) t , ∆ ( ) j t , G j ( ) t } . The total system Hamiltonian takes the form

$$H _ { s y s } = \sum _ { j } H _ { 1, j } + H _ { 2, j } + H _ { 3 }.$$

Single photon losses on qubit j are obviously modelled by an annihilation operator a j , resulting in a jump operator of the form V j = a j .

Table II. Optimizable pulses parameters used in the optimizations. The minima and maxima values of each drive are based on the minima and maxima found in the paper [24]. The columns for each gate denote if the pulse values are real ( R ) or complex ( C ). All units in [ K ].

| Pulse z ( t )   | min[ z ( t )]   | max[ z ( t )]   | z (0) = z ( T )   | R Z ( ϕ )   | R X ( θ )   | R Y ( φ )   | R ZZ (Θ)   |
|-----------------|-----------------|-----------------|-------------------|-------------|-------------|-------------|------------|
| E ( t )         | -0.308          | 0.308           | 0.000             | R           | R           | C           | R          |
| ∆( t )          | 0.000           | 4.000           | 0.000             | R           | R           | R           | R          |
| G ( t )         | -4.000          | 4.000           | 4.000             | R           | R           | R           | R          |
| g ( t )         | 0.000           | 0.154           | 0.000             | n.a.        | n.a.        | n.a.        | R          |
|                 |                 |                 | T [1 /K ]         | 2.06        | 2.06        | 4.03        | 3.00       |

## Appendix B: Wasserstein-2 norm

The Wasserstein-2 distance W z ,z 2 ( 1 2 ) calculates the cost of moving the pulse z 1 onto z 2 , and is more closely related to our intuitive understanding of similarity between pulses than the L 2 distance [33]. To understand this distance, we first define a distribution on a pulse µ z [ ] as

$$\mu [ z ] ( A ) = \frac { 1 } { c [ z ] } \int _ { ( t, y ) \in A } \delta ( y = z [ t ] ) d ( t, y ), \ \ A \subset \mathbb { R } ^ { 2 }.$$

Thus, the measure is over R 2 and only puts weight on the points of the pulse ( t, z ( )). t Here c z [ ] is a normalization constant s.t. µ z [ ]( R 2 ) = 1 and is simply equal to the length of the curve, e.g.

$$c [ z ] = \int _ { 0 } ^ { T } \sqrt { 1 + ( \partial _ { t } z ) ^ { 2 } } d t.$$

We now want to transport the distribution µ z [ 1 ] to the distribution µ z [ 2 ], and for that purpose define a coupling γ ∈ Π( µ z [ 1 ] , µ z [ 2 ]), which is a joint probability measure on R 2 × R 2 whose marginals are µ z [ 1 ] and µ z [ 2 ] [34] (see Fig. 3), e.g.

$$\int _ { \mathbb { R } ^ { 2 } } \gamma ( x, y ) d y = \mu [ z _ { 1 } ] ( x ), \ \int _ { \mathbb { R } ^ { 2 } } \gamma ( x, y ) d x = \mu [ z _ { 2 } ] ( y ).$$

Here, Π( µ z [ 1 ] , ν [ z 1 ]) is the set of all transport couplings between distributions µ z [ 1 ] and µ z [ 2 ]. The Wasserstein-2 distance W z ,z 2 ( 1 2 ) is then defined as the minimal coupling cost

$$W _ { 2 } ( z _ { 1 }, z _ { 2 } ) ^ { 2 } \coloneqq \min _ { \gamma \in \Pi ( \mu [ z _ { 1 } ], \mu [ z _ { 2 } ] ) } \int _ { \mathbb { R } ^ { 2 } \times \mathbb { R } ^ { 2 } } | x - y | ^ { 2 } \, d \gamma ( x, y ).$$

## Appendix C: Pulse Construction

In order to construct a mock family of pulses, we define for θ 0 = 0 a sum of G Gaussians where each Gaussian g has a random amplitude a 0 ,g, 1 , center a 0 ,g, 2 and variance a 0 ,g, 3 chosen according to some uniform random distributions. For subsequent θ i +1 , we then have these coefficients evolve according to

$$a _ { i + 1, g, k } = a _ { i, g, k } + \mu _ { g, k } \Delta \theta + \sigma _ { g, k } \mathcal { N } \sqrt { \Delta \theta },$$

where µ g,k and σ g,k are predefined drifts and variances, respectively. N is a normally distributed random variable and ∆ θ is the step size. These kinds of pulses mimic the pulse optimization solutions we get when constructing pulses for a continuously parametrized gate U θ ( ). When the step size, drift and variances are high, subsequent pulses look less alike and the clustering problem becomes harder to solve, leading to lesser results for both Wasserstein-2 and L 2 distances. Nevertheless, in all tried out configurations we saw Wasserstein-2 outperform L 2 , especially when subsequent pulses are more alike, which can always be achieved by lowering ∆ θ .
Robert de Keijzer Jurgen Snijders André Carvalho Servaas Kokkelmans
Pulse family optimization for parametrized quantum gates using spectral clustering
Parametrized gate circuits are used in plentiful applications in the current
NISQ era of quantum computing. These parametrized gates are chiefly implemented
using analytically found pulse protocols, often yielding suboptimal gate times,
and consequently, fidelities. Alternatively, gate optimization algorithms are
designed to construct high fidelity pulses for individual, fixed points in
continuous parameter space. Gates for intermediate parameters can subsequently
be found by some form of interpolation between previously constructed pulses.
Nevertheless, it is not guaranteed (as with analytic protocols) that the pulses
found by the optimization algorithms belong to the same \textit{family} of
solutions and thus show resemblance. Interpolation between two pulses of
differing solution families often leads to high infidelities, as the pulse
strays away from the minimum in the parameter/fidelity landscape. In this work,
we introduce a \textit{spectral clustering} method to sort high-fidelity,
optimized pulses in families, and interpolating solely between pulses of the
same family. Accordingly, interpolations will always approach maximal fidelity.
Furthermore, as more than one pulse family is constructed, the parameter space
can be partitioned according to which family prevails fidelity-wise. This work
provides a meticulous demonstration of our constitutive continuous gate family
construction by applying it to a universal gate set for Rydberg and Cat qubits
under noise.