All Publications

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.

We study the relaxation of the local ferromagnetic order in the transverse field quantum Ising chain with power-law decaying interactions ∼ 1/rα. We prepare the system in the GHZ state and study the time evolution of the probability distribution function (PDF) of the order parameter within a block of when quenching the transverse field. The model is known to support long range order at finite temperature for α ≤ 2.0. In this regime, quasi-localized topological magnetic defects are expected to strongly affect the equilibration of the full probability distribution. We highlight different dynamical regimes where gaussification mechanism may be slowed down by confinement and eventually breaks. We further study the PDF dynamics induced by changing the effective dimensionality of the system; we mimic this by quenching the range of the interactions. As a matter of fact, the behavior of the system crucially depends on the value of α governing the unitary evolution.

Nonequilibrium 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 an abrupt change in the scaling law of the bipartite entanglement entropy, thus suggesting the existence of different nonequilibrium regimes. However, our understanding of how these regimes appear and whether they survive in the thermodynamic limit is much less established. Here we investigate these questions on a one-dimensional quadratic fermionic model: this allows us to reach system sizes relevant in the thermodynamic sense. We show that local projective measurements induce a qualitative modification of the time growth of the entanglement entropy which changes from linear to logarithmic. However, in the stationary regime, the logarithmic behavior of the entanglement entropy does not survive in the thermodynamic limit and, for any finite value of the measurement rate, we numerically show the existence of a single area-law phase for the entanglement entropy. Finally, exploiting the quasiparticle picture, we further support our results by analyzing the fluctuations of the stationary entanglement entropy and its scaling behavior.

In many cases, neural networks can be mapped into tensor networks with an exponentially large bond dimension. Here, we compare different sub-classes of neural network states, with their mapped tensor network counterpart for studying the ground state of short-range Hamiltonians. We show that when mapping a neural network, the resulting tensor network is highly constrained and thus the neural network states do in general not deliver the naive expected drastic improvement against the state-of-the-art tensor network methods. We explicitly show this result in two paradigmatic examples, the 1D ferromagnetic Ising model and the 2D antiferromagnetic Heisenberg model, addressing the lack of a detailed comparison of the expressiveness of these increasingly popular, variational ansätze.

We present an unconstrained tree-tensor-network approach to the study of lattice gauge theories in two spatial dimensions, showing how to perform numerical simulations of theories in the presence of fermionic matter and four-body magnetic terms, at zero and finite density, with periodic and open boundary conditions. We exploit the quantum-link representation of the gauge fields and demonstrate that a fermionic rishon representation of the quantum links allows us to efficiently handle the fermionic matter while finite densities are naturally enclosed in the tensor network description. We explicitly perform calculations for quantum electrodynamics in the spin-one quantum-link representation on lattice sizes of up to 16×16 sites, detecting and characterizing different quantum regimes. In particular, at finite density, we detect signatures of a phase separation as a function of the bare mass values at different filling densities. The presented approach can be extended straightforwardly to three spatial dimensions.

In spin chains with local unitary evolution preserving the magnetization Sz, the domain-wall state typically "melts."At large times, a nontrivial magnetization profile develops in an expanding region around the initial position of the domain wall. For nonintegrable dynamics, the melting is diffusive, with entropy production within a melted region of size t. In contrast, when the evolution is integrable, ballistic transport dominates and results in a melted region growing linearly in time, with no extensive entropy production: The spin chain remains locally in states of zero entropy at any time. Here we show that, for the integrable spin-1/2 XXZ chain, low-energy quantum fluctuations in the melted region give rise to an emergent Luttinger liquid which, remarkably, differs from the equilibrium one. The striking feature of this emergent Luttinger liquid is its quasiparticle charge (or Luttinger parameter K), which acquires a fractal dependence on the XXZ chain anisotropy parameter Δ.

We study the relaxation of the local ferromagnetic order in the quantum Ising chain in a slant field with both longitudinal and transverse components. After preparing the system in a fully polarised state, we analyse the time evolution of the entire probability distribution function (PDF) of the magnetisation within a block of l spins. We first analyse the effect of confinement on the Gaussification of the PDF for large l, showing that the melting of initial order is suppressed when the longitudinal field is aligned to initial magnetisation while it is speed up when it is in the opposite direction. Then we study the thermalisation dynamics. In the paramagnetic region, the PDF quickly shows thermal features. Conversely, in the ferromagnetic phase, when confinement takes place, the relaxation suffers a typical slowing down which depends on the interplay between the strength of the longitudinal field, the density of excitations, and the direction of the initial polarisation. Even when the initial magnetisation is aligned oppositely to the longitudinal field, confinement prevents thermalisation in the accessible timescale, as it is neatly bared by the PDF.

Confinement of excitations induces quasilocalized dynamics in disorder-free isolated quantum many-body systems in one spatial dimension. This occurrence is signaled by severe suppression of quantum correlation spreading and of entanglement growth, long-time persistence of spatial inhomogeneities, and long-lived coherent oscillations of local observables. In this work, we present a unified understanding of these dramatic effects. The slow dynamical behavior is shown to be related to the Schwinger effect in quantum electrodynamics. We demonstrate that it is quantitatively captured for long-time scales by effective Hamiltonians exhibiting Stark localization of excitations and weak growth of the entanglement entropy for arbitrary coupling strength. This analysis explains the phenomenology of real-time string dynamics investigated in a number of lattice gauge theories, as well as the anomalous dynamics observed in quantum Ising chains after quenches. Our findings establish confinement as a robust mechanism for hindering the approach to equilibrium in translationally invariant quantum statistical systems with local interactions.