
Entanglement transitions with free fermions
Joseph Merritt and Lukasz Fidkowski
Department of Physics, University of Washington, Seattle, WA 98195-1560, USA
(Dated: October 2022)
We use Majorana operators to study entanglement dynamics under random free fermion unitary
evolution and projective measurements in one dimension. For certain choices of unitary evolution,
namely those which swap neighboring Majorana operators, and measurements of neighboring Ma-
jorana bilinears, one can map the evolution to the statistical model of completely packed loops
with crossings (CPLC) and study the corresponding phase diagram. We generalize this model us-
ing the language of fermionic Gaussian states to a general free fermion unitary evolution acting on
neighboring Majorana operators, and numerically compute its phase diagram. We find that both
the Goldstone and area law phases persist in this new phase diagram, but with a shifted phase
boundary. One important qualitative aspect of the new phase boundary is that even for the case of
commuting measurements, the Goldstone phase persists up to a finite non-zero measurement rate.
This is in contrast with the CPLC, in which non-commuting measurements are necessary for real-
izing the Goldstone phase. We also numerically compute the correlation length critical exponent at
the transition, which we find to be near to that of the CPLC, and give a tentative symmetry based
explanation for some differences in the phase transition line between the CPLC and generalized
models.
I. INTRODUCTION
Recently the study of ‘hybrid’ quantum circuits, in-
volving both unitary dynamics and projective measure-
ments, has received a great deal of attention [1–3]. By
focusing on the ensemble of quantum trajectories of pure
states defined by the various measurement outcomes one
can study new types of non-equilibrium phase transitions,
with the canonical example being the ‘entanglement tran-
sition’. In the entanglement transition, the ensemble-
averaged entanglement entropy changes from scaling as
an area law to scaling as a volume law as the measure-
ment rate is decreased. A closely related concept is that
of a ‘purification’ transition, where instead of pure quan-
tum state trajectories one studies the purifying behavior
of an initial maximally mixed state.
In the special case of free fermion dynamics [4] - i.e.
unitaries which are exponentials of bilinears of the cre-
ation and annihilation operators, and measurements only
of Fock space mode occupation numbers - the mixed (or
‘volume law’) phase is known to be unstable to any non-
zero measurement rate [5]. However, such free fermion
dynamics can still accommodate an interesting phase
transition from a purifying phase to a so-called ‘Gold-
stone’ phase [6]. The latter still exhibits purifying be-
havior, but on time scales parametrically longer in system
size. More precisely, in the Goldstone phase the entropy
for a system of length Lafter a time of order Lscales like
log L(whereas in the purifying phase this entropy would
be close to 0, i.e. the state would have approximately pu-
rified long before this time scale was reached). This phase
transition can be realized in a specific free fermion model,
one that can be solved [7] by exact mapping to a known
statistical mechanical model, the completely packed loop
model with crossings (CPLC) [6].
A natural question one may ask is, how generic is the
CPLC phase diagram in the context of free fermion hy-
brid dynamics? In other words, do the area law and
Goldstone phases persist when the dynamics is deformed
slightly away from the specific point dual to the CPLC
model? Is the phase transition still continuous, and is it
in the same universality class as that of the CPLC model?
In this work, we investigate these questions by extend-
ing the CPLC-dual free fermion model to a more general
family of free fermion models. Specifically, following [7],
the CPLC-dual model has a convenient description in
terms of a one dimensional chain of Majorana fermions,
as follows: the unitary gates swap a neighboring pair of
Majoranas (γ1→γ2, γ2→ −γ1), and the measurements
measure the occupation number of a free fermion mode
defined by a neighboring pair of Majoranas. These gates
are implemented with respect to one pairing of Majo-
ranas and its complementary pairing in an alternating
fashion, as described in detail below. The phase diagram
is a function of two parameters, pand q, which control,
respectively, the rate of measurements and the asymme-
try between the two complementary pairings of Majo-
ranas. Our generalized model replaces the swap gate,
which may be thought of as a π
2rotation in the SO(2)
that rotates γ1into γ2, by a rotation by a random an-
gle inside this SO(2). The measurement gates are as in
the CPLC model, and the phase diagram is once again a
function of the two parameters pand q.
Our generalized model no longer admits an easily solv-
able statistical mechanical dual, although in principle
some statistical mechanical dual should exist, as dis-
cussed below. To study it, we therefore instead leverage
the free fermion nature of the dynamics to perform effi-
cient Monte Carlo simulations using the Gaussian state
formalism. The essential feature of both the CPLC and
our generalized models which makes this possible is the
fact that, for a particular quantum trajectory, the many-
body quantum state of 2NMajoranas remains Gaussian,
meaning that it can be efficiently encoded in a correla-
arXiv:2210.05681v3 [cond-mat.str-el] 15 Nov 2022