Entanglement fluctuations associated with the Schrödinger evolution of wave functions offer a unique perspective on various fundamental issues ranging from quantum thermalization to state preparation in quantum devices. Very recently, a subset of the present authors showed that, in a class of free-fermion lattice models and interacting spin chains, entanglement dynamics enters into a new regime at long time, with entanglement probes displaying persistent temporal fluctuations, whose statistics falls into the seemingly disparate paradigm of mesoscopic fluctuations in condensed-matter physics. This motivates us to revisit here entanglement dynamics of a canonical bosonic model in many-body physics, i.e., a coupled harmonic oscillator chain. We find that when the system is driven out of equilibrium, the long-time entanglement dynamics exhibits strictly the same statistical behaviors as those of free-fermion models. Specifically, irrespective of entanglement probes and microscopic parameters, the statistical distribution of entanglement fluctuations is flanked by asymmetric tails: sub-Gaussian for upward fluctuations and sub-Gamma for downward; moreover, the variance exhibits a crossover from the scaling ∼ 1 / L to ∼ L A 3 / L 2 as the subsystem size L A increases ( L the total system size). This insensitivity to the particle statistics, dubbed boson-fermion universality, is contrary to the common wisdom that statistical phenomena of a many-body nature depend strongly on particle statistics. Together with our previous work, the present study indicates rich fluctuation phenomena in long-time entanglement dynamics that await in-depth exploration.
Entanglement fluctuations associated with the Schrödinger evolution of wave functions offer a unique perspective on various fundamental issues ranging from quantum thermalization to state preparation in quantum devices. Very recently, a subset of the present authors showed that, in a class of free-fermion lattice models and interacting spin chains, entanglement dynamics enters into a new regime at long time, with entanglement probes displaying persistent temporal fluctuations, whose statistics falls into the seemingly disparate paradigm of mesoscopic fluctuations in condensed-matter physics. This motivates us to revisit here entanglement dynamics of a canonical bosonic model in many-body physics, i.e., a coupled harmonic oscillator chain. We find that when the system is driven out of equilibrium, the long-time entanglement dynamics exhibits strictly the same statistical behaviors as those of free-fermion models. Specifically, irrespective of entanglement probes and microscopic parameters, the statistical distribution of entanglement fluctuations is flanked by asymmetric tails: sub-Gaussian for upward fluctuations and sub-Gamma for downward; moreover, the variance exhibits a crossover from the scaling ∼ 1 / L to ∼ L A 3 / L 2 as the subsystem size L A increases ( L the total system size). This insensitivity to the particle statistics, dubbed boson-fermion universality, is contrary to the common wisdom that statistical phenomena of a many-body nature depend strongly on particle statistics. Together with our previous work, the present study indicates rich fluctuation phenomena in long-time entanglement dynamics that await in-depth exploration.
Understanding fluctuation phenomena plays a dominant role in the development of many-body physics. The time evolution of entanglement is essential to a broad range of subjects in many-body physics, ranging from exotic quantum matter to quantum thermalization. Stemming from various dynamical processes of information, fluctuations in entanglement evolution differ conceptually from out-of-equilibrium fluctuations of traditional physical quantities. Their studies remain elusive. Here we uncover an emergent random structure in the evolution of the many-body wavefunction in two classes of integrable-either interacting or noninteracting-lattice models. It gives rise to out-of-equilibrium entanglement fluctuations which fall into the paradigm of mesoscopic fluctuations of wave interference origin. Specifically, the entanglement entropy variance obeys a universal scaling law in each class, and the full distribution displays a sub-Gaussian upper and a sub-Gamma lower tail. These statistics are independent of both the system's microscopic details and the choice of entanglement probes, and broaden the class of mesoscopic universalities. They have practical implications for controlling entanglement in mesoscopic devices.
On finite system sizes, monitored one-dimensional free-fermion chains display a broad crossover in entanglement scaling as the measurement rate increases, from sub-extensive behavior at weak monitoring toward an area law at strong monitoring. Analytical field theory and recent large-scale simulations indicate that this apparent change is not a finite-rate transition in the thermodynamic limit. I test, using trajectory-resolved correlation-matrix simulations (chains up to L=96 and up to 128 Born-rule trajectories per parameter point), whether two quantities built from the single-particle entanglement spectrum are useful finite-size diagnostics: the bipartite particle-number fluctuation FA=∑kνk(1−νk) and a participation-style effective number of entangling modes, MA=exp[−∑kwklnwk], with wk∝νk(1−νk). Both quantities track the crossover. However, for trajectory- and time-averaged steady-state values, M is nearly a deterministic function of F: a pooled curve explains 99.7\% of its variance across all rates, so M carries little independent information. Its main advantage is statistical: its trajectory-to-trajectory coefficient of variation is up to a factor of ∼2 smaller than that of the entropy or F under strong monitoring. Subsystem-scaling comparisons at L=48--96 also show that a pure logarithmic law is not statistically preferred under weak monitoring; a small residual quasi-extensive component and a drifting logarithmic coefficient are instead consistent with a crossover rather than a critical phase. Thus, F remains the natural experimentally motivated companion to the entropy, whereas M is best viewed as a variance-reduced numerical summary of essentially the same information.
Motivated by recent development of the concept of the disorder operator and its relation with entanglement entropy in bosonic systems, here we show the disorder operator successfully probes many aspects of quantum entanglement in fermionic many-body systems. From both analytical and numerical computations in free and interacting fermion systems in 1D and 2D, we find the disorder operator and the entanglement entropy exhibit similar universal scaling behavior, as a function of the boundary length of the subsystem, but with subtle yet important differences. In 1D they both follow the \log{L} log L scaling behavior with the coefficient determined by the Luttinger parameter for disorder operator, and the conformal central charge for entanglement entropy. In 2D they both show the universal L\log L L log L scaling behavior in free and interacting Fermi liquid states, with the coefficients depending on the geometry of the Fermi surfaces. However at a 2D quantum critical point with non-Fermi-liquid state, extra symmetry information is needed in the design of the disorder operator, so as to reveal the critical fluctuations as does the entanglement entropy. Our results demonstrate the fermion disorder operator can be used to probe quantum many-body entanglement related to global symmetry, and provide new tools to explore the still largely unknown territory of highly entangled fermion quantum matter in 2 or higher dimensions.
Out-of-equilibrium quantum many-body systems stand at the forefront of modern theoretical physics, addressing fundamental questions on thermalization, transport, and universal phenomena. These fields, already well established in condensed matter, statistical physics, quantum optics, and quantum information theory, are progressively gaining even greater relevance with the rapid development of quantum technologies and simulation, which inherently operate in dynamical regimes. In recent years, the traditional paradigm of unitary quantum evolution has been expanded to include measurements, opening new directions in out-of-equilibrium physics. At the core of these advances lie measurement-induced phase transitions (MIPTs), which have emerged as a new class of dynamical critical phenomena characterizing the general behavior of monitored quantum dynamics. When external monitoring intertwines with unitary evolution, many-body quantum correlations change their structure, giving rise to distinct entanglement phases of matter. This discovery has sparked enormous interest in MIPTs, leading to significant advances in open quantum systems, entanglement theory, and more broadly quantum complexity. Despite much progress, a full understanding of monitored many-body dynamics is far from complete, leaving several open questions on the nature of MIPTs, their experimental observability, and the possibilities offered by measurements to enhance control over synthetic quantum matter. These issues persist due to the intrinsic complexity of the problem and the lack of efficient tools to study it, mainly caused by the stochastic character of monitored evolution. This thesis addresses these challenges by expanding the investigation of measurement-induced phenomena in new settings and introducing innovative probes of entanglement and many-body quantum complexity for MIPTs. A core question we investigate is the role of symmetries, non-ergodicity, and especially integrability in measurement-induced criticality, which dramatically affect the non-equilibrium phases. We further explore how these phenomena extend beyond bipartite quantum correlations to multipartite entanglement and quantum non-stabilizerness, highlighting the non-trivial interplay between measurements and complexity notions rooted in quantum information theory. Finally, we focus on the compelling problem of decoherence, modeling how noise spoils entanglement structures. These findings, supported by advanced numerical simulations and theoretical analysis, deepen the current understanding of entanglement, complexity, and integrability in monitored quantum many-body systems, offering new perspectives on their rich behavior. In parallel, we address the experimental problem of dissipation in MIPTs, which is of key relevance for practical implementations. We anticipate the present investigation to foster future research on the nature of monitored dynamical critical phenomena and, more broadly, the applications of measurements in quantum state evolution.
One way to look for complex behaviours in many‐body quantum systems is to let the number N of degrees of freedom become large and focus upon collective observables. Mean‐field quantities scaling as tend to commute, whence complexity at the quantum level can only be inherited from complexity at the classical level. Instead, fluctuations of microscopic observables scale as and exhibit collective Bosonic features, typical of a mesoscopic regime half‐way between the quantum one at the microscopic level and the classical one at the level of macroscopic averages. Here, we consider the mesoscopic behaviour emerging from an infinite quantum spin chain undergoing a microscopic dissipative, irreversible dynamics and from global states without long‐range correlations and invariant under lattice translations and dynamics. We show that, from the fluctuations of one site spin observables whose linear span is mapped into itself by the dynamics, there emerge bosonic operators obeying a mesoscopic dissipative dynamics mapping Gaussian states into Gaussian states. Instead of just depleting quantum correlations because of decoherence effects, these maps can generate entanglement at the collective, mesoscopic level, a phenomenon with no classical analogue that embodies a peculiar complex behaviour at the interface between micro and macro regimes.
With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.