for every Gaussian state |Gi, we can choose a canonical basis ξ
a
≡ (q
i
, p
i
), such that
G ≡ 1. This means the bilinear form G does not contain information that is invariant
under changing the canonical basis or put simply: “All Gaussian states look the same if
we can choose the right basis for each individual state.” This changes of course, if we have
two Gaussian states |Gi and |
˜
Gi in the same system
12
and force ourselves to represent
the two-point functions G and
˜
G with respect to the same canonical basis. Again, we
can choose a basis, such that G ≡ 1, but we will not be able to accomplish the same
for
˜
G. The remaining freedom of choosing a canonical basis is described by the group
U(N) = Sp(2N, R) ∩ SO(2N ) consisting of canonical transformation (i.e., MΩM
|
= Ω)
that simultaneously orthogonal with respect to G (i.e., MGM
|
= G). The invariant
information about the relation between the original state |ψi and the transformed state
|
˜
ψi is completely captured by the eigenvalues of the relative covariance matrix
13
∆
a
b
=
˜
G
ac
g
cb
with g = G
−1
, (3.12)
i.e., G
ac
g
cb
= δ
a
b
. In particular, any quantities that depend on the two states in a
Sp(2N, R)-invariant way, e.g., their inner product,
14
can be computed purely from ∆.
This will apply to the complexity provided that we choose a geometry that is Sp(2N, R)-
invariant, e.g., we do not introduce penalty factors which conflict with the group structure.
For our Bogoliubov transformation (3.5), we have spec(∆) = (e
2r
, e
−2r
). We say that |
˜
ψi
arises from a one-mode squeezing of |ψi with squeezing parameter r. For bosonic Gaussian
states, understanding one-mode squeezing is the key to relate any two states. That is, for
any two bosonic Gaussian states |ψi and |
˜
ψi with N degrees of freedom, there exists a
normal mode basis (q
1
, ··· , q
N
, p
1
, ··· , p
N
), such that |
˜
ψi is the result of N independent
one-mode squeezing operations in the N different normal modes [39, 51, 52]. This is related
to the Iwasawa (or KAN) decomposition of Sp(2N, R), e.g., see [56, 57].
3.2 Two fermions
We now turn to the case of fermionic Gaussian states. In this case, the space of Gaussian
states for N fermionic degrees of freedom is given by the quotient M
f,N
= O(2N)/U(N),
which has dimension N(N − 1), e.g., [52]. Of course, this space is a small submanifold
within the full 2
N
-dimensional Hilbert space H of the fermionic system. Further, it is not
preserved by general unitary transformations U(2
N
) acting on H, but only the subgroup
O(2N) corresponding to Bogoliubov transformations. That is, the most straightforward
way to think of characterizing the fermionic Gaussian states is in terms of the annihilation
and creation operators. With N fermionic pairs (a
i
, a
†
i
) satisfying {a
i
, a
†
j
} = δ
ij
, the corre-
sponding Gaussian state is again defined by a
i
|ψi = 0 and the Bogoliubov transformations
mixing these fermionic operators map Gaussian states to Gaussian states.
12
Of course, this is the situation where we are examining circuit complexity of states since we have both
the target state and the reference state.
13
Note that one could have just as easily defined
b
∆ = G ˜g with ˜g =
˜
G
−1
. However, one then has
b
∆ = ∆
−1
and due to the fact that ∆ is symplectic, the two have the same spectrum. This is discussed in more detail
in section 5.
14
For bosonic states, we find the simple formula |hG|
˜
Gi|
2
= det
√
2∆
1/4
√
1+∆
derived in [52].
– 11 –