All Publications

We introduce a novel hybrid approach combining tensor network methods with the stabilizer formalism to address the challenges of simulating many-body quantum systems. By integrating these techniques, we enhance our ability to accurately model unitary dynamics while mitigating the exponential growth of entanglement encountered in classical simulations. We demonstrate the effectiveness of our method through applications to random Clifford T-doped circuits and random Clifford Floquet dynamics. This approach offers promising prospects for advancing our understanding of complex quantum phenomena and accelerating progress in quantum simulation.

We investigate the robustness of the many-body localized (MBL) phase to the quantum-avalanche instability by studying the dynamics of a localized spin chain coupled to a T=∞ thermal bath through its leftmost site. By analyzing local magnetizations we estimate the size of the thermalized sector of the chain and find that it increases logarithmically slowly in time. This logarithmically slow propagation of the thermalization front allows us to lower-bound the slowest thermalization time, and find a broad parameter range where it scales fast enough with the system size that MBL is robust against thermalization induced by avalanches. The further finding that the imbalance - a global quantity measuring localization - thermalizes over a timescale that is exponential both in disorder strength and system size is in agreement with these results.

In this paper, we investigate the relationship between entanglement and nonstabilizerness (also known as magic) in matrix product states (MPSs). We study the relation between magic and the bond dimension used to approximate the ground state of a many-body system in two different contexts: full state of magic and mutual magic (the nonstabilizer analog of mutual information, thus free of boundary effects) of spin-1 anisotropic Heisenberg chains. Our results indicate that obtaining converged results for nonstabilizerness is typically considerably easier than entanglement. For full state magic at critical points and at sufficiently large volumes, we observe convergence with 1/χ2, with χ being the MPS bond dimension. At small volumes, magic saturation is so quick that, within error bars, we cannot appreciate any finite-χ correction. Mutual magic also shows a fast convergence with bond dimension, whose specific functional form is however hindered by sampling errors. As a byproduct of our study, we show how Pauli-Markov chains (originally formulated to evaluate magic) resets the state of the art in terms of computing mutual information for MPS. We illustrate this last fact by verifying the logarithmic increase of mutual information between connected partitions at critical points. By comparing mutual information and mutual magic, we observe that, for connected partitions, the latter is typically scaling much slower - if at all - with the partition size, while for disconnected partitions, both are constant in size.

We present a novel classical algorithm designed to learn the stabilizer group - namely, the group of Pauli strings for which a state is a ±1 eigenvector - of a given matrix product state (MPS). The algorithm is based on a clever and theoretically grounded biased sampling in the Pauli (or Bell) basis. Its output is a set of independent stabilizer generators whose total number is directly associated with the stabilizer nullity, notably a well-established nonstabilizer monotone. We benchmark our method on T-doped states randomly scrambled via Clifford unitary dynamics, demonstrating very accurate estimates up to highly entangled MPS with bond dimension χ∼103. Our method, thanks to a very favorable scaling O(χ3), represents the first effective approach to obtain a genuine magic monotone for MPS, enabling systematic investigations of quantum many-body physics out of equilibrium.

Conserved-charge densities are very special observables in quantum many-body systems as, by construction, they encode information about the dynamics. Therefore, their evolution is expected to be of much simpler interpretation than that of generic observables and to return universal information on the state of the system at any given time. Here, we study the dynamics of the fluctuations of conserved U(1) charges in systems that are prepared in charge-asymmetric initial states. We characterize the charge fluctuations in a given subsystem using the full-counting statistics of the truncated charge and the quantum entanglement between the subsystem and the rest resolved to the symmetry sectors of the charge. We show that, even though the initial states considered are homogeneous in space, the charge fluctuations generate an effective inhomogeneity due to the charge-asymmetric nature of the initial states. We use this observation to map the problem into that of charge fluctuations on inhomogeneous, charge-symmetric states and treat it using a recently developed space-time duality approach. Specializing the treatment to interacting integrable systems we combine the space-time duality approach with generalized hydrodynamics to find explicit predictions.

Repeated measurements can induce entanglement phase transitions in the dynamics of quantum systems. Interacting models, both chaotic and integrable, generically show a stable volume-law entangled phase at low measurement rates that disappears for free, Gaussian fermions. Interactions break the Gaussianity of a dynamical map in its unitary part, but non-Gaussianity can be introduced through measurements as well. By comparing the entanglement and non-Gaussianity structure of different protocols, we propose a single particle indicator of the measurement-induced phase transition, and we use it to argue in favor of the stability of the transition when non-Gaussianity is purely provided by measurements.

Nonstabilizerness, also colloquially referred to as magic, is a resource for advantage in quantum computing and lies in the access to non-Clifford operations. Developing a comprehensive understanding of how nonstabilizerness can be quantified and how it relates to other quantum resources is crucial for studying and characterizing the origin of quantum complexity. In this work, we establish a direct connection between nonstabilizerness and entanglement spectrum flatness for a pure quantum state. We show that this connection can be exploited to efficiently probe nonstabilizerness even in the presence of noise. Our results reveal a direct connection between nonstabilizerness and entanglement response, and define a clear experimental protocol to probe nonstabilizerness in cold atom and solid-state platforms.

We study the thermalization of the transverse field Ising chain with a power law decaying interaction ∼ 1/rα following a global quantum quench of the transverse field in two different dynamical regimes. The thermalization behavior is quantified by comparing the full probability distribution function (PDF) of the evolving states with the corresponding thermal state given by the canonical Gibbs ensemble (CGE). To this end, we used the matrix product state (MPS)-based Time Dependent Variational Principle (TDVP) algorithm to simulate both real time evolution following a global quantum quench and the finite temperature density operator. We observe that thermalization is strongly suppressed in the region with strong confinement for all interaction strengths α, whereas thermalization occurs in the region with weak confinement.

The dynamics of a quantum system, undergoing unitary evolution and continuous monitoring, can be described in term of quantum trajectories. Although the averaged state fully characterizes expectation values, the entire ensemble of stochastic trajectories goes beyond simple linear observables, keeping a more attentive description of the entire dynamics. Here we go beyond the Lindblad dynamics and study the probability distribution of the expectation value of a given observable over the possible quantum trajectories. The measurements are applied to the entire system, having the effect of projecting the system into a product state. We develop an analytical tool to evaluate this probability distribution at any time t. We illustrate our approach by analyzing two paradigmatic examples: a single qubit subjected to magnetization measurements, and a free hopping particle subjected to position measurements.

We introduce a novel approach to evaluate the nonstabilizerness of an N-qubits matrix product state (MPS) with bond dimension χ. In particular, we consider the recently introduced stabilizer Rényi entropies (SREs). We show that the exponentially hard evaluation of the SREs can be achieved by means of a simple perfect sampling of the many-body wave function over the Pauli string configurations. The sampling is achieved with a novel MPS technique, which enables us to compute each sample in an efficient way with a computational cost O(Nχ3). We benchmark our method over randomly generated magic states, as well as in the ground-state of the quantum Ising chain. Exploiting the extremely favorable scaling, we easily have access to the nonequilibrium dynamics of the SREs after a quantum quench.

Owing to its probabilistic nature, a measurement process in quantum mechanics produces a distribution of possible outcomes. This distribution - or its Fourier transform known as full counting statistics (FCS) - contains much more information than say the mean value of the measured observable, and accessing it is sometimes the only way to obtain relevant information about the system. In fact, the FCS is the limit of an even more general family of observables - the charged moments - that characterize how quantum entanglement is split in different symmetry sectors in the presence of a global symmetry. Here we consider the evolution of the FCS and of the charged moments of a U(1) charge truncated to a finite region after a global quantum quench. For large scales these quantities take a simple large-deviation form, showing two different regimes as functions of time: while for times much larger than the size of the region they approach a stationary value set by the local equilibrium state, for times shorter than region size they show a nontrivial dependence on time. We show that, whenever the initial state is also U(1) symmetric, the leading order in time of FCS and charged moments in the out-of-equilibrium regime can be determined by means of a space-time duality. Namely, it coincides with the stationary value in the system where the roles of time and space are exchanged. We use this observation to find some general properties of FCS and charged moments out of equilibrium, and to derive an exact expression for these quantities in interacting integrable models. We test this expression against exact results in the Rule 54 quantum cellular automaton and exact numerics in the XXZ spin-1/2 chain.

In quadratic fermionic models, we determine a quantum correction to the work statistics after both a sudden quench and a time-dependent driving. Such a correction lies in the noncommutativity of the initial quantum state and the time-dependent Hamiltonian, and is revealed via the Kirkwood-Dirac quasiprobability (KDQ) approach to two-times correlators. Thanks to the latter, one can assess the onset of nonclassical signatures in the KDQ distribution of work, in the form of negative and complex values that no classical theory can reveal. By applying these concepts on the one-dimensional transverse-field Ising model, we relate nonclassical behaviors of the KDQ statistics of work in correspondence of the critical points of the model. Finally, we also prove the enhancement of the extracted work in nonclassical regimes where the noncommutativity takes a role.

Non-equilibrium dynamics of many-body quantum systems under the effect of measurement protocols is attracting an increasing amount of attention. It has been recently revealed that measurements may induce different non-equilibrium regimes and an abrupt change in the scaling-law of the bipartite entanglement entropy. However, our understanding of how these regimes appear, how they affect the statistics of local quantities and, finally whether they survive in the thermodynamic limit, is much less established. Here we investigate measurement-induced phase transitions in the Quantum Ising chain coupled to a monitoring environment. In particular we show that local projective measurements induce a quantitative modification of the out-of-equilibrium probability distribution function of the local magnetization. Starting from a GHZ state, the relaxation of the paramagnetic and the ferromagnetic order is analysed. In particular we describe how the probability distributions associated to them show different behaviour depending on the measurement rate.

Here, we focus on the problem of minimizing complex classical cost functions associated with prototypical discrete neural networks, specifically the paradigmatic Hopfield model and binary perceptron. We show that the adiabatic time evolution of QA can be efficiently represented as a suitable Tensor Network. This representation allows for simple classical simulations, well-beyond small sizes amenable to exact diagonalization techniques. We show that the optimized state, expressed as a Matrix Product State (MPS), can be recast into a Quantum Circuit, whose depth scales only linearly with the system size and quadratically with the MPS bond dimension. This may represent a valuable starting point allowing for further circuit optimization on near-term quantum devices.

A large ongoing research effort focuses on variational quantum algorithms (VQAs), representing leading candidates to achieve computational speed-ups on current quantum devices. The scalability of VQAs to a large number of qubits, beyond the simulation capabilities of classical computers, is still debated. Two major hurdles are the proliferation of low-quality variational local minima, and the exponential vanishing of gradients in the cost-function landscape, a phenomenon referred to as barren plateaus. In this work, we show that by employing iterative search schemes, one can effectively prepare the ground state of paradigmatic quantum many-body models, also circumventing the barren plateau phenomenon. This is accomplished by leveraging the transferability to larger system sizes of a class of iterative solutions, displaying an intrinsic smoothness of the variational parameters, a result that does not extend to other solutions found via random-start local optimization. Our scheme could be directly tested on near-term quantum devices, running a refinement optimization in a favorable local landscape with nonvanishing gradients.

We study the dynamics of the quantum Dyson hierarchical model in its paramagnetic phase. An initial state made by a local excitation of the paramagnetic ground state is considered. We provide analytical predictions for its time evolution, solving the single-particle dynamics on a hierarchical network. A localization mechanism is found, and the excitation remains close to its initial position at arbitrary times. Furthermore, a universal scaling among space and time is found that is related to the algebraic decay of the interactions as r-1-σ. We compare our predictions to numerics, employing tensor network techniques, for large magnetic fields, discussing the robustness of the mechanism in the full many-body dynamics.

We consider the two-dimensional quantum Ising model, in absence of disorder, subject to periodic imperfect global spin flips. We show by a combination of exact diagonalization and tensor-network methods that the system can sustain a spontaneously broken discrete time-translation symmetry. Employing careful scaling analysis, we show the feasibility of a two-dimensional discrete time-crystal (DTC) prethermal phase. Despite an unbounded energy pumped into the system, in the high-frequency limit, a well-defined effective Hamiltonian controls a finite-temperature intermediate regime, wherein local time averages are described by thermal averages. As a consequence, the long-lived stability of the DTC relies on the existence of a long-range ordered phase at finite temperature. Interestingly, even for large deviations from the perfect spin flip, we observe a nonperturbative change in the decay rate of the order parameter, which is related to the long-lived stability of the magnetic domains in 2D.

By taking inspiration from the backflow transformation for correlated systems, we introduce a tensor network Ansatz which extends the well-established matrix product state representation of a quantum many-body wave function. This structure provides enough resources to ensure that states in dimensions larger than or equal to one obey an area law for entanglement. It can be efficiently manipulated to address the ground-state search problem by means of an optimization scheme which mixes tensor-network and variational Monte Carlo algorithms. We benchmark the Ansatz against spin models both in one and two dimensions, demonstrating high accuracy and precision. We finally employ our approach to study the challenging S=1/2 two-dimensional (2D) J1-J2 model, demonstrating that it is competitive with the state-of-the-art methods in 2D.

Discrete time crystals represent a paradigmatic nonequilibrium phase of periodically driven matter. Protecting its emergent spatiotemporal order necessitates a mechanism that hinders the spreading of defects, such as localization of domain walls in disordered quantum spin chains. In this work, we establish the effectiveness of a different mechanism arising in clean spin chains: the confinement of domain walls into "mesonic"bound states. We consider translationally invariant quantum Ising chains periodically kicked at arbitrary frequency, and we discuss two possible routes to domain-wall confinement: longitudinal fields and interactions beyond nearest neighbors. We study the impact of confinement on the order-parameter evolution by constructing domain-wall-conserving effective Hamiltonians and analyzing the resulting dynamics of domain walls. On the one hand, we show that for arbitrary driving frequency, the symmetry-breaking-induced confining potential gets effectively averaged out by the drive, leading to deconfined dynamics. On the other hand, we rigorously prove that increasing the range R of spin-spin interactions Ji,j beyond nearest neighbors enhances the order-parameter lifetime exponentially in R. Our theory predictions are corroborated by a combination of exact and matrix-product-state simulations for finite and infinite chains, respectively. The long-lived stability of spatiotemporal order identified in this work does not rely on Floquet prethermalization nor on eigenstate order, but rather on the nonperturbative origin of vacuum-decay processes. We point out the experimental relevance of this new mechanism for stabilizing a long-lived time-crystalline response in Rydberg-dressed spin chains.

We develop a general approach to compute the symmetry-resolved Rényi and von Neumann entanglement entropies (SREE) of thermodynamic macrostates in interacting integrable systems. Our method is based on a combination of the thermodynamic Bethe ansatz and the Gärtner-Ellis theorem from large deviation theory. We derive an explicit simple formula for the von Neumann SREE, which we show to coincide with the thermodynamic Yang-Yang entropy of an effective macrostate determined by the charge sector. Focusing on the XXZ Heisenberg spin chain, we test our result against iTEBD calculations for thermal states, finding good agreement. As an application, we provide analytic predictions for the asymptotic value of the SREE following a quantum quench.