We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $⟨A(t)B{⟩}_{\beta }\to ⟨A{⟩}_{\beta }⟨B{⟩}_{\beta }$, thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

We consider the optimization of a finite-time Carnot engine characterized by small dissipations. We bound the power with a simple inequality and show that the optimal strategy is to perform small cycles around a given working point, which can be, thus, chosen optimally. Remarkably, this optimal point is independent of the figure of merit combining power and efficiency that is being maximized. Furthermore, for a general class of dissipative dynamics the maximal power output becomes proportional to the heat capacity of the working substance. Since the heat capacity can scale supraextensively with the number of constituents of the engine, this enables us to design optimal many-body Carnot engines reaching maximum efficiency at finite power per constituent in the thermodynamic limit.

The occurrence of an inverse current, where the sign of the induced current is opposite to the applied force, is a highly counterintuitive phenomenon. We show that inverse currents in coupled transport (ICC) of energy and particle can occur in a one-dimensional interacting Hamiltonian system when its equilibrium state is perturbed by coupled thermodynamic forces. This seemingly paradoxical result is possible due to the self-organization occurring in the system in response to the applied forces.

We study a method to simulate quantum many-body dynamics of spin ensembles using measurement-based feedback. By performing a weak collective measurement on a large ensemble of two-level quantum systems and applying global rotations conditioned on the measurement outcome, one can simulate the dynamics of a mean-field quantum kicked top, a standard paradigm of quantum chaos. We analytically show that there exists a regime in which individual quantum trajectories adequately recover the classical limit, and show the transition between noisy quantum dynamics to full deterministic chaos described by classical Lyapunov exponents. We also analyze the effects of decoherence, and show that the proposed scheme represents a robust method to explore the emergence of chaos from complex quantum dynamics in a realistic experimental platform based on an atom-light interface.

By using highly entangled states, quantum metrology guarantees a precision impossible with classical measurements. Unfortunately such states can be very susceptible to noise, and it is a great challenge of the field to maintain quantum advantage in realistic conditions. In this Letter we investigate the practicality of graph states for quantum metrology. Graph states are a natural resource for much of quantum information, and here we characterize their quantum Fisher information for an arbitrary graph state. We then construct families of graph states which approximately achieves the Heisenberg limit, we call these states bundled graph states. We demonstrate that bundled graph states maintain a quantum advantage after being subjected to independent and identically distributed dephasing or finite erasures. This shows that these graph states are good resources for robust quantum metrology. We also quantify the number of $n$ qubit stabilizer states that are useful as a resource for quantum metrology.

We analyze the dynamics of an initially trapped cloud of interacting quantum particles on a lattice under a linear (Stark) potential. We reveal a dichotomy: initially trapped interacting systems possess features typical of both many-body-localized and thermalizing systems. We consider both fermions ($t\text{−}V$ model) and bosons (Bose-Hubbard model). For the zero and infinite interaction limits, both systems are integrable: we provide analytic solutions in terms of the moments of the initial cloud shape and clarify how the recurrent dynamics (many-body Bloch oscillations) depends on the initial state. Away from the integrable points, we identify and explain the timescale at which Bloch oscillations decohere.

We use a near quantum limited detector to experimentally track individual quantum state trajectories of a driven qubit formed by the hybridization of a waveguide cavity and a transmon circuit. For each measured quantum coherent trajectory, we separately identify energy changes of the qubit as heat and work, and verify the first law of thermodynamics for an open quantum system. We further establish the consistency of these results by comparison with the master equation approach and the two-projective-measurement scheme, both for open and closed dynamics, with the help of a quantum feedback loop that compensates for the exchanged heat and effectively isolates the qubit.

A generalization of quantum discord to multipartite systems is proposed. A key feature of our formulation is its consistency with the conventional definition of discord in bipartite systems. It is by construction zero only for systems with classically correlated subsystems and is a non-negative quantity, giving a measure of the total nonclassical correlations in the multipartite system with respect to a fixed measurement ordering. For the tripartite case, we show that the discord can be decomposed into contributions resulting from changes induced by nonclassical correlation breaking measurements in the conditional mutual information and tripartite mutual information. The former gives a measure of the bipartite nonclassical correlations and is a non-negative quantity, while the latter is related to the monogamy of the nonclassical correlations.

We propose a time-of-arrival operator in quantum mechanics by conditioning on a quantum clock. This allows us to bypass some of the problems of previous proposals, and to obtain a Hermitian time of arrival operator whose probability distribution arises from the Born rule and which has a clear physical interpretation. The same procedure can be employed to measure the “time at which some event happens” for arbitrary events (and not just specifically for the arrival time of a particle).

A remarkable feature of quantum many-body systems is the orthogonality catastrophe that describes their extensively growing sensitivity to local perturbations and plays an important role in condensed matter physics. Here we show that the dynamics of the orthogonality catastrophe can be fully characterized by the quantum speed limit and, more specifically, that any quenched quantum many-body system, whose variance in ground state energy scales with the system size, exhibits the orthogonality catastrophe. Our rigorous findings are demonstrated by two paradigmatic classes of many-body systems—the trapped Fermi gas and the long-range interacting Lipkin-Meshkov-Glick spin model.

The Rényi entanglement entropy in quantum many-body systems can be viewed as the difference in free energy between partition functions with different trace topologies. We introduce an external field $\lambda$ that controls the partition function topology, allowing us to define a notion of nonequilibrium work as $\lambda$ is varied smoothly. Nonequilibrium fluctuation theorems of the work provide us with statistically exact estimates of the Rényi entanglement entropy. This framework also naturally leads to the idea of using quench functions with spatially smooth profiles, providing us a way to average over lattice scale features of the entanglement entropy while preserving long distance universal information. We use these ideas to extract universal information from quantum Monte Carlo simulations of SU($N$) spin models in one and two dimensions. The vast gain in efficiency of this method allows us to access unprecedented system sizes up to $192\times 96$ spins for the square lattice Heisenberg antiferromagnet.

We prove that matrix-product unitaries with on-site unitary symmetries are completely classified by the (chiral) index and the cohomology class of the symmetry group $G$, provided that we can add trivial and symmetric ancillas with arbitrary on-site representations of $G$. If the representations in both system and ancillas are fixed to be the same, we can define symmetry-protected indices (SPIs) which quantify the imbalance in the transport associated to each group element and greatly refines the classification. These SPIs are stable against disorder and measurable in interferometric experiments. Our results lead to a systematic construction of two-dimensional Floquet symmetry-protected topological phases beyond the standard classification, and thus shed new light on understanding nonequilibrium phases of quantum matter.

As a promising candidate for exhibiting quantum computational supremacy, Gaussian boson sampling (GBS) is designed to exploit the ease of experimental preparation of Gaussian states. However, sufficiently large and inevitable experimental noise might render GBS classically simulable. In this work, we formalize this intuition by establishing a sufficient condition for approximate polynomial-time classical simulation of noisy GBS—in the form of an inequality between the input squeezing parameter, the overall transmission rate, and the quality of photon detectors. Our result serves as a nonclassicality test that must be passed by any quantum computational supremacy demonstration based on GBS. We show that, for most linear-optical architectures, where photon loss increases exponentially with the circuit depth, noisy GBS loses its quantum advantage in the asymptotic limit. Our results thus delineate intermediate-sized regimes where GBS devices might considerably outperform classical computers for modest noise levels. Finally, we find that increasing the amount of input squeezing is helpful to evade our classical simulation algorithm, which suggests a potential route to mitigate photon loss.

A powerful perspective in understanding nonequilibrium quantum dynamics is through the time evolution of its entanglement content. Yet apart from a few guiding principles for the entanglement entropy, to date, much less is known about the refined characteristics of entanglement propagation. Here, we unveil signatures of the entanglement evolving and information propagating out of equilibrium, from the view of the entanglement Hamiltonian. We investigate quantum quench dynamics of prototypical Bose-Hubbard model using state-of-the-art numerical technique combined with conformal field theory. Before reaching equilibrium, it is found that a current operator emerges in the entanglement Hamiltonian, implying that entanglement spreading is carried by particle flow. In the long-time limit the subsystem enters a steady phase, evidenced by the dynamic convergence of the entanglement Hamiltonian to the expectation of a thermal ensemble. Importantly, the entanglement temperature in steady state is spatially independent, which provides an intuitive trait of equilibrium. These findings not only provide crucial information on how equilibrium statistical mechanics emerges in many-body dynamics, but also add a tool to exploring quantum dynamics from the perspective of the entanglement Hamiltonian.

It is now experimentally possible to entangle thousands of qubits, and efficiently measure each qubit in parallel in a distinct basis. To fully characterize an unknown entangled state of $n$ qubits, one requires an exponential number of measurements in $n$, which is experimentally unfeasible even for modest system sizes. By leveraging (i) that single-qubit measurements can be made in parallel, and (ii) the theory of perfect hash families, we show that all $k$-qubit reduced density matrices of an $n$ qubit state can be determined with at most $e^{\mathcal{O}(k)}{\log }^2(n)$ rounds of parallel measurements. We provide concrete measurement protocols which realize this bound. As an example, we argue that with near-term experiments, every two-point correlator in a system of 1024 qubits could be measured and completely characterized in a few days. This corresponds to determining nearly 4.5 million correlators.

The chain rule for the classical relative entropy ensures that the relative entropy between probability distributions on multipartite systems can be decomposed into a sum of relative entropies of suitably chosen conditional distributions on the individual systems. Here, we prove a chain rule inequality for the quantum relative entropy. The new chain rule allows us to solve an open problem in the context of asymptotic quantum channel discrimination: surprisingly, adaptive protocols cannot improve the error rate for asymmetric channel discrimination compared to nonadaptive strategies.

We consider current-current deformations that generalize $T\overline{T}$ ones, and show that they may be also introduced for integrable spin chains. In analogy with the integrable QFT setup, we define the deformation as a modification of the $S$ matrix in the Bethe equations. Using results by Bargheer, Beisert and Loebbert we show that the deforming operator is composite and constructed out of two currents on the lattice; its expectation value factorizes like for $T\overline{T}$. Such a deformation may be considered for any combination of charges that preserve the model’s integrable structure.