When the two largest eigenvalues of an entanglement or negativity spectrum meet, the gap of the corresponding entanglement or negativity Hamiltonian closes while entropies stay smooth, an event we call a temporal criticality. We ask when such criticalities occur in open fermionic systems, in which order, and whether they precede the geometric decoherence time, at which the negativity starts to fall. We compare loss and gain, which respect fermion parity only on average (weak symmetry) and keep the state Gaussian, with density dephasing, which respects it in every jump (strong symmetry) and makes the state non-Gaussian. We find that Gaussianity, rather than the symmetry, fixes the criterion for a criticality. In a Gaussian state the entanglement gap of a subsystem closes exactly when its two parity sectors carry equal weight, and the negativity gap closes exactly when this happens for the subsystem or its complement, so the negativity criticality never comes after the entanglement one. Uniform balanced loss and gain only shrink the two parities and leave both criticalities in place at every rate, for any quadratic Hamiltonian. In an open Kitaev chain, weak imbalanced loss and gain shift the criticalities by amounts set by the occupations at the cut, and at strong rates they allow a closing near the relaxation time only when an odd number of neighboring sites pass half filling. For dephasing we solve the non-Gaussian dynamics exactly at the point where the bond terms of the chain commute. Weak dephasing delays both criticalities and, away from the chain ends, moves them ahead of the parity balance, where the parity rules of Gaussian states would place them. As the dephasing rate approaches critical damping, both criticalities recede to infinite time, and at strong dephasing the negativity starts to fall long before either gap begins to shrink.
When the two largest eigenvalues of an entanglement or negativity spectrum meet, the gap of the corresponding entanglement or negativity Hamiltonian closes while entropies stay smooth, an event we call a temporal criticality. We ask when such criticalities occur in open fermionic systems, in which order, and whether they precede the geometric decoherence time, at which the negativity starts to fall. We compare loss and gain, which respect fermion parity only on average (weak symmetry) and keep the state Gaussian, with density dephasing, which respects it in every jump (strong symmetry) and makes the state non-Gaussian. We find that Gaussianity, rather than the symmetry, fixes the criterion for a criticality. In a Gaussian state the entanglement gap of a subsystem closes exactly when its two parity sectors carry equal weight, and the negativity gap closes exactly when this happens for the subsystem or its complement, so the negativity criticality never comes after the entanglement one. Uniform balanced loss and gain only shrink the two parities and leave both criticalities in place at every rate, for any quadratic Hamiltonian. In an open Kitaev chain, weak imbalanced loss and gain shift the criticalities by amounts set by the occupations at the cut, and at strong rates they allow a closing near the relaxation time only when an odd number of neighboring sites pass half filling. For dephasing we solve the non-Gaussian dynamics exactly at the point where the bond terms of the chain commute. Weak dephasing delays both criticalities and, away from the chain ends, moves them ahead of the parity balance, where the parity rules of Gaussian states would place them. As the dephasing rate approaches critical damping, both criticalities recede to infinite time, and at strong dephasing the negativity starts to fall long before either gap begins to shrink.
Abstract We investigate the dynamics of the fermionic logarithmic negativity in a free-fermion chain with a localized loss, which acts as a dissipative impurity. The chain is initially prepared in a generic Fermi sea. In the standard hydrodynamic limit of large subsystems and long times, with their ratio fixed, the negativity between two subsystems is described by a simple formula, which depends only on the effective absorption coefficient of the impurity. The negativity grows linearly at short times, then saturating to a volume-law scaling. Physically, this reflects the continuous production with time of entangling pairs of excitations at the impurity site. Interestingly, the negativity is not the same as the Rényi mutual information with Rényi index 1 / 2 , in contrast with the case of unitary dynamics. This reflects the interplay between dissipative and unitary processes. The negativity content of the entangling pairs is obtained in terms of an effective two-state mixed density matrix for the subsystems. Criticality in the initial Fermi sea is reflected in the presence of logarithmic corrections. The prefactor of the logarithmic scaling depends on the loss rate, suggesting a nontrivial interplay between dissipation and criticality.
The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometricdecoherencetime, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-21 entropy -- exactly equal across any bipartition for pure states -- breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail -- residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local Z-dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.
The interplay of unitary evolution and local measurements in many-body systems gives rise to a stochastic state evolution and to measurement-induced phase transitions in the pure state entanglement. In realistic settings, however, this dynamics may be spoiled by decoherence, e.g., dephasing, due to coupling to an environment or measurement imperfections. We investigate the impact of dephasing and the inevitable evolution into a non-Gaussian, mixed state, on the dynamics of monitored fermions. We approach it from three complementary perspectives: (i) the exact solution of the conditional master equation for small systems, (ii) quantum trajectory simulations of Gaussian states for large systems, and (iii) a renormalization group analysis of a bosonic replica field theory. For weak dephasing, constant monitoring preserves a weakly mixed state, which displays a robust measurement-induced phase transition between a critical and a pinned phase, as in the decoherence-free case. At strong dephasing, we observe the emergence of a new scale describing an effective temperature, which is accompanied with an increased mixedness of the fermion density matrix. Remarkably, observables such as density-density correlation functions or the subsystem parity still display scale invariant behavior even in this strongly mixed phase. We interpret this as a signature of gapless, classical diffusion, which is stabilized by the balanced interplay of Hamiltonian dynamics, measurements, and decoherence.
Entanglement entropy scaling links geometry to quantum correlations and diagnoses quantum criticality. Fractal lattices extend this link to noninteger dimensions, where boundary and bulk dimensions are independent. It remains open which dimension governs entropy scaling and how strongly gap closing can enhance it. Nonlocal occupied-state projectors and the lack of translation symmetry obstruct standard methods. We overcome these obstacles and classify free-fermion entropy scaling with bounds set by the hopping across the subsystem boundary: an upper bound from its rank and the spectral gap, and a lower bound from how much it lowers the ground-state energy. With only negative-energy states filled and this energy lowering proportional to the boundary, the entanglement entropy obeys the boundary law SX=Θ(LXdb) under a uniform gap and satisfies LXdb≲SX≲LXdblnLX when the gap closes no faster than algebraically. Here LX is the subsystem linear size and db the boundary dimension. Half occupation of the chiral zero modes gives a mixed-state entropy SXmix with max{SX,zXln2}⩽SXmix⩽SX+zXln2, where zX is the local zero-mode weight; an extensive zX thus makes SXmix follow the bulk dimension.
A closed finite quantum system can return arbitrarily close to its initial state, but how far can a region's entropy depart from its initial value during the evolution? For finite quadratic fermionic systems, we reduce the maximum over all initial states and all times to a static spectral calculation: universal two-sided bounds identify it with B = ‖P − P̄‖1, where P selects the region and P̄ is its spectral dephasing. The criterion includes mixed, non-Gaussian and initially entangled states. In an explicit random Kitaev chain of 2N sites, the disorder-averaged maximum stays bounded uniformly in length away from statistical self-duality. At self-duality it equals (ln 2/6) ln N + o(ln N), including after optimally freezing any fixed number of spectral modes. The same leading term holds throughout each fixed finite-size window |χ|√N ≤ M. A single preparation, vacuum on one half and maximally mixed on the other, already produces logarithmic growth in its long-time averaged local entropy. Cross-boundary spectral channels thus determine both the existence of an entropy budget and its exact critical growth, even when the many-body preparation is unknown.
In this thesis we study entanglement in one-dimensional critical quantum many-body systems and in particular we will focus on disconnected regions. Given a system in a pure state, to quantify the amount of entanglement between a multicomponent subsystem and the rest of the full system we can use as an entanglement measure the renowned entanglement entropy. However, if we are interested in the entanglement shared among the disconnected regions, the entanglement entropy fails to be a good quantifier. The reason is that the state of the subsystem is in general mixed once the rest is traced out, and the entanglement entropy is a good measure of entanglement only for pure states. A good measure of entanglement in mixed states is the logarithmic negativity, which is the quantitative version of Peres' criterion of separability. The main advantage of the negativity with respect to other entanglement measures is the simplicity of its definition in terms of the density matrix describing the quantum state. Since the definition does not require any variational calculus, it is much more easily computable than any other entanglement measure, and therefore we can obtain some results also in complicated settings such as many-body systems. We will be mostly interested in the configuration where the full system is divided into two parts, A and B. If we do not have access to the degrees of freedom in B, we can describe subsystem A through its reduced density matrix, where the degrees of freedom in B have been traced out. The residual subsystem in A will in general be left in a mixed state. Subsystem A is then divided again into two parts, A_1 and A_2, and we will be interested in the entanglement shared by these two components. Knowing the full density matrix, the logarithmic negativity can be easily computed from the eigenvalues of its partial transpose with respect to the degrees of freedom living in one of the two subsystems A_1 or A_2. However, the full density matrix of a many-body state is in general unaccessible, even numerically, since the size of the matrix grows exponentially with the size of the system. If we concentrate on critical systems whose low-energy physics can be described by a quantum field theory, we can resort to its powerful tools to compute the entanglement properties. In particular the main tool is the replica trick. The entanglement entropy is obtained from the moments of the reduced density matrix, while the negativity can be computed from the moments of the partial transpose. However, even the computation of the integer-order moments is not an easy task, and exact analytical results are known only for the simplest quantum field theories. Hence, throughout the thesis we will usually consider conformal field theories, which have an enhanced set of symmetries and therefore allow for some exact computations. In the Introduction, we will try to frame this work by briefly reviewing some of the main topics in the literature of quantum many-body systems and quantum field theories where entanglement plays a crucial role. We will also describe what are the main features that are requested to a `good' measure of entanglement, and we will review some of the main measures that have been considered so far in the context of quantum information. We will define the entanglement negativity and stress its main advantages, as well as its drawbacks. After reviewing some basic facts of conformal field theories, we will give some technical details on the computation of entanglement entropy and logarithmic negativity in general quantum field theories, which will be needed in the rest of the thesis. In Chap. 2 we will start the study of the entanglement of several disjoint disconnected regions by considering the entanglement that the union of these regions share with the rest of the system. We will compute the integer-order Rényi entropies of the subsystem from which the entanglement entropy can be obtained through the replica limit. Unfortunately, the exact analytic continuation to real order is still out of reach and therefore the entanglement entropy cannot be computed exactly. However numerical extrapolations allow to obtain some accurate estimates. In the remaining chapters we will focus on the entanglement shared between two non complementary disjoint regions, through the computations of the entanglement negativity in different settings. In Chap. 3, we will consider a global quantum quench starting from a conformal boundary state and evolving through a conformal evolution. A general formula for the mutual information and the logarithmic negativity for two adjacent and disjoint intervals is given in the spacetime limit. In Chap. 4 we focus on the XY chain and we recover a formula for the partial transpose as a sum of four fermionic auxiliary Gaussian density matrices, by generalizing some previous results for the reduced density matrix of spin systems and for the partial transpose of pure fermionc Gaussian states. Even if the computation of the negativity is still out of reach, we can obtain formulas for the integer-order moments of the partial transpose in terms of the correlation matrices relative to the two components and to the region in between them. In Chap. 5 we will study the entanglement negativity for a free Dirac fermion field. Starting from the lattice results, we will obtain a path integral representation of the partial transpose which can be easily generalized to the continuum. With this representation of the partial transpose, we can construct all its integer-order moments and compute them for the simple conformal field theory of the free fermion. Again, the computation of the negativity could not be accomplished due to technical difficulties in the analytic continuation to real order. The same computation has been extended also to the modular invariant Dirac fermion and the Ising model, and the obtained formulas coincide with the ones already present in the literature. This analysis draws an interesting connection between some terms appearing in the formulas for the moments of the partial transpose (as well as of the R\'enyi entropies) found for the lattice models and the ones found for the corresponding theories describing their scaling limit. Whenever possible, all the analytical results will be checked against numerical calculations performed on simple free chain models. In Chaps. 2 and 3 we will consider bosonic Gaussian states, specifically the harmonic chain in its ground state and out of equilibrium, while in Chaps. 4 and 5 we will consider fermionic systems, specifically the XX and Ising spin chains, and the tight-binding model. Finally, we will draw some conclusions and discuss some open problems.