
2
This work uses the Hilbert-Schmidt metric ∥A∥=
ptr(AA†), as it seems both most appropriate and ac-
cessible to analytical considerations. As a strength mea-
sure, one can choose the metric to satisfy desirable
properties based on the application [16]. It should lead
to a function on unitaries that satisfy: (i) f(uA⊗uB) =
0, (ii) f(UV)≤f(U) + f(V)and (iii) f(U⊗I) = f(U).
The first two properties are satisfied by all strength
measures, and we show that if we pick the Hilbert-
Schmidt metric, the third property, known as stability, is
satisfied by KD(U), see Appendix A. A particular appli-
cation of interest is in circuit complexity; one can show
that the strength measure on the operator norm of a
bipartite gate can be used to understand its capacity
to create entanglement when acting upon pure states,
based on recent arguments in Ref. [23].
Thus far there are partial results concerning KD(U)
for d=2 and no known results for d>2. Apart from
deriving several new results concerning general bounds
and exact evaluations of this measure, we show that
the set of dual-unitaries are maximally and equally away
from product unitaries for the case of qubits, d=2, and
present strong numerical evidence that this is the case
even for d>2. This motivated the caricature in Fig. (1)
which indicates the special place of dual-unitaries in the
space of bipartite operators.
These considerations are intimately related to an in-
teresting property that is easy to see for qubits, but also
appears to hold in general: for any bipartite unitary U
in Hd⊗ Hd, there exists a pair of maximally entangled
states |Φ1⟩and |Φ2⟩such that U|Φ1⟩=|Φ2⟩. A dy-
namical map is devised that converges to such a pair
of maximally entangled states for any U, which can be
used to find the closest local unitaries to dual-unitary
gates. A closely related procedure works surprisingly
well also for general (non-dual) unitaries, which allows
for detailed numerical explorations of KD(U). We find
KD(U)analytically for several special families of uni-
taries which also verifies the numerical procedure.
I. GENERAL CONSIDERATIONS AND BOUNDS FOR
KD(U)
Let the operator Schmidt decomposition of an arbi-
trary two-qudit unitary gate Ube:
U=
d2
∑
i=1pλimA
i⊗mB
i. (2)
Here the set {mA
i,i=1, ··· ,d2}forms an orthonormal
operator basis in subspace Athat is
trmA
imA†
j=dδij,
and mBis a similar set in subspace B. The operator
Schmidt coefficients λkare chosen to be in decreasing
order, 1 ≥λ1≥ ··· ≥ λd2≥0, and from the unitarity
of Uit follows that ∑d2
k=1λk=1.
The Schmidt decomposition provides measures of
operator entanglement, for example one such is the lin-
ear entropy
E(U) = 1−
d2
∑
k=1
λ2
k. (3)
This range is 0 ≤E(U)≤E(S) = (d2−1)/d2, and is
0 iff Uis a product unitary (in which case KD(U) = 0).
The maximum value is attained when for all k,λk=
1/d2, and are the operator equivalents of Bell states.
The swap operator S(S|kl⟩=|lk⟩) achieves this value
and is an important example of a maximally entangled
unitary operator. However, Sis by far not the only uni-
tary operator to achieve this value. If an unitary Uis
such that E(U) = E(S), this may also be taken as the
definition of dual-unitary operators. Restriction of an
unitary Uto this special set will be generically denoted
as Udual, that is E(Udual) = 1. One of our goal is to
calculate KD(Udual), but we first turn to general state-
ments and bounds.
For the distance measure KD(U)in Eq. (1), for the rest
of the paper we use the Hilbert-Schmidt metric. Define
K∗
D(U):=∥U−mA
1⊗mB
1∥=q2d2−2d2pλ1. (4)
As the Schmidt decomposition already provides the
nearest product operator (for a proof in the case of
states, see [17]), this quantity is the distance to the near-
est product operator, removing the unitarity constraint
from the local operators in Eq. (1). As KD(U)is the
distance to a more constrained set, the following lower-
bound follows:
KD(U)≥K∗
D(U). (5)
Alternatively, the definition in Eq. (1) implies
K2
D(U) = min
uA,uB∥U−uA⊗uB∥2
=min
uA,uB2d2−2Re htr(U†(uA⊗uB))i
=2d2−max
uA,uB
2tr(U†(uA⊗uB)).
(6)
The second equality follows from the fact that the
phase can be absorbed by the local unitaries. Ex-
pand local unitaries uAand uBin the orthonormal
bases from the Schmidt decomposition of Uin Eq. (2),
uA=∑d2
i=1αimA
i,uB=∑d2
i=1βimB
i, where ∑i|αi|2=
∑i|βi|2=1. This leads to tr(U†(uA⊗uB))=
d2d2
∑
i=1pλiαiβi≤d2pλ1
d2
∑
i=1|αi||βi| ≤ d2pλ1. (7)