
2
surements [89,90]. This enables numerical investigation
of the transition between QEC and QZ phases in higher
dimensions. In particular, Ref. [91–93] studied the (2+1)-
dimensional systems but found conflicting results, likely
due to the limited system sizes in their computations or
the choice of observable used to locate the transition.
Therefore, it is desirable to employ large-scale numeri-
cal simulations and clarify the phase diagram in higher
dimensions for stabilizer circuits.
This paper revisits the MIPT in (d+ 1)-dimensional
stabilizer circuits. We utilize large-scale numerics up to
N= 32768 qubits to unveil a precise characterization of
the critical exponents of the MIPT for d= 1,2, as well
as consider the previously unexplored case d= 3. We
provide a coherent picture of measurement-induced crit-
icality by probing the three aspects of dynamics of the
circuits that change between QEC and QZ phases: the
entanglement of the pure state, the purification of ini-
tially maximally mixed state as well as the structure of
the wave function along quantum trajectories. We test
the robustness of our results by considering different ob-
servables and various spatiotemporal architectures of the
circuits. Furthermore, we investigate the bulk-boundary
critical properties by considering ancilla qubits coupled
to our circuits. In the spirit of Ref. [94], we suggest that
the resulting statistical field theory at the transition be-
tween QEC and QZ phases is that of a perturbed perco-
lation model, and hence the critical point is conformal.
This ansatz is compatible with our numerical findings.
We show that the critical exponents for (1+1)D circuits
are close but distinct from the percolation transition in
2dimensions, whereas the results for (d+ 1)-dimensional
circuits for d= 2,3are compatible with the percolation
transition in (d+ 1) spatial dimensions, see Table. Ifor
the summary of the results.
The rest of this work is structured as follows. First,
in Sec. II we give a brief overview of our approach and
summarize the obtained critical exponents for MIPT in
(d+ 1)-dimensional circuits. Sec. III describes the ar-
chitectures of stabilizer circuits that are employed in our
investigations of MIPTs. In Sec. IV we detail the entan-
glement entropy measures and extract the universal prop-
erties of the MIPT through finite-size scaling. Sec. Vcon-
tains results about the purification aspect of MIPT while
Sec. VI describes the wave-function structure change at
the transition between QEC and QZ phases. Sec. VII
is devoted to the investigation of the bulk-boundary
features of the MIPT through coupling with ancillary
qubits. Finally, Sec. IV discusses our results and illus-
trates the conformal character of the critical point, while
Sec. IX concludes the manuscript with some outlooks. In
the Appendices, we detail the most technical aspects of
our work.
II. OVERVIEW OF THE RESULTS
Let us first present the rationale of our analysis
and summarize our main findings. After motivating
the choice of observables for the quantum trajectories,
we briefly discuss how they describe the measurement-
induced phases and transitions, highlighting the key re-
sults. In summary, we establish that (d+ 1) spacetime
dimensional stabilizer circuits have a MIPT with univer-
sal behavior close to a percolation field theory in (d+ 1)
spatial dimensions. In particular, for d= 1 the stabilizer
circuit critical exponents are distinct and within five er-
ror bars from the 2D percolation exact results, while for
d= 2 (d= 3)all our estimates are compatible within one
error bar with those of 3D (4D) percolation I. In addition,
we report evidence suggesting an emergence of a confor-
mal field theory right at the MIPT in (d+1)-dimensional
stabilizer circuits for d= 2 and d= 3.
a. Observables of interest and quantum trajectories
This manuscript considers quantum circuits build from:
(i) random unitary gates (drawn uniformly from the Clif-
ford group), (ii) projective measurements onto the local
qubit basis. In particular, a circuit realization Km(with
implicit depth t, also referred to as time) stems from
the choice of Clifford gates, the spacetime location of
measurements, and the measurement outcomes. Start-
ing from the initial state ρ0, the final state is given by
ρm=Kmρ0K†
m
tr(K†
mKmρ0).(1)
We denote the collective average by Em. The trajectory
nature of Eq. (1) translates to statistical properties of an
observable Ξon the state ρm. If Ξis linear in the density
matrix, it is easy to see that
Em[Ξ[ρm]] = Ξ[Em(ρm)].(2)
Notably, the above equation does not hold when Ξis a
non-linear function of the density matrix. Consider for
instance the purity Ξ[ρ] = tr(ρ2), and a pure state quan-
tum trajectory. Since the average state ρ=Em(ρm)is
mixed, we have 1 = Em[Ξ[ρm]] 6= Ξ[ρ]≤1. In particular,
the average state evolves through a quantum channel and
is insensible of the trajectory registry m. In other words,
the statistical structure of quantum trajectories contains
more information (including measurement-induced criti-
cality) than the average state. We consider entropic ob-
servables: being non-linear in the density matrix, they
allow us to investigate the various aspects of the MIPT.
Crucially, as discussed in the remaining, they are ef-
ficiently computable for stabilizer states, allowing ex-
tensive numerical characterization of MIPTs in (d+ 1)-
dimensional stabilizer circuits. In Sec. III we detail the
properties of stabilizer circuits and specify the circuital
architectures treated in this manuscript.
b. Entanglement transition. For a given density ma-
trix ρ, the von Neumann entropy is defined as
S(ρ)≡ −tr(ρlog2ρ)(3)