We present a path integral formalism for expressing matrix elements of the density matrix of a quantum many-body system between any two coherent states in terms of standard Matsubara action with periodic (antiperiodic) boundary conditions on bosonic (fermionic) fields. We show that this enables us to express several entanglement measures for bosonic and fermionic many-body systems described by a Gaussian action in terms of the Matsubara Green function. We apply this formalism to compute various entanglement measures for the two-dimensional Bose-Hubbard model in the strong-coupling regime, both in the presence and absence of Abelian and non-Abelian synthetic gauge fields, within a strong-coupling mean-field theory. In addition, our method provides an alternative formalism for addressing time evolution of quantum many-body systems, with Gaussian actions, driven out of equilibrium without the use of Keldysh technique. We demonstrate this by deriving analytical expressions of the return probability and the counting statistics of several operators for a class of integrable models represented by free Dirac fermions subjected to a periodic drive in terms of the elements of their Floquet Hamiltonians. We provide a detailed comparison of our method with the earlier, related, techniques used for similar computations, discuss the significance of our results, and chart out other systems where our formalism can be used.
We present a path integral formalism for expressing matrix elements of the density matrix of a quantum many-body system between any two coherent states in terms of standard Matsubara action with periodic (antiperiodic) boundary conditions on bosonic (fermionic) fields. We show that this enables us to express several entanglement measures for bosonic and fermionic many-body systems described by a Gaussian action in terms of the Matsubara Green function. We apply this formalism to compute various entanglement measures for the two-dimensional Bose-Hubbard model in the strong-coupling regime, both in the presence and absence of Abelian and non-Abelian synthetic gauge fields, within a strong-coupling mean-field theory. In addition, our method provides an alternative formalism for addressing time evolution of quantum many-body systems, with Gaussian actions, driven out of equilibrium without the use of Keldysh technique. We demonstrate this by deriving analytical expressions of the return probability and the counting statistics of several operators for a class of integrable models represented by free Dirac fermions subjected to a periodic drive in terms of the elements of their Floquet Hamiltonians. We provide a detailed comparison of our method with the earlier, related, techniques used for similar computations, discuss the significance of our results, and chart out other systems where our formalism can be used.
Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime.
This chapter examines the Schwinger-Keldysh technique, applying path integrals, and the functional method of 2PI effective actions—with the latter serving to construct powerful symmetry-preserving nonperturbative approximations for equilibrium as well as nonequilibrium descriptions of strongly interacting quantum many-body systems. It provides a rigorous mathematical basis for exploring the quantum mechanical basis of thermalization and decoherence. Mesoscopic systems with quantum processes, in particular, play a key role. The domain of nonequilibrium mesoscopic systems is flooded with open questions some of which are solved but still others have provisional or incomplete solutions. A correlation function encodes the correlation between the field operators at different times and is an inference to all thermodynamic quantities. The beauty of the Keldysh time path is the definition of the Green's function, which plays an analogous role in the nonequilibrium formalism similar to the role the causal Green's function plays in the equilibrium theory.
Abstract The functional integral representation of the density matrix at thermal equilibrium in non-relativistic quantum mechanics (QM) with many degrees of freedom, in the grand canonical formulation is introduced. In QM, Hamiltonians H(p,q) can be also expressed in terms of creation and annihilation operators, a method adapted to the study of perturbed harmonic oscillators. In the holomorphic formalism, quantum operators act by multiplication and differentiation on a vector space of analytic functions. Alternatively, they can also be represented by kernels, functions of complex variables that correspond in the classical limit to a complex parametrization of phase space. The formalism is adapted to the description of many-body boson systems. To this formalism corresponds a path integral representation of the density matrix at thermal equilibrium, where paths belong to complex spaces, instead of the more usual position–momentum phase space. A parallel formalism can be set up to describe systems with many fermion degrees of freedom, with Grassmann variables replacing complex variables. Both formalisms can be generalized to quantum gases of Bose and Fermi particles in the grand canonical formulation. Field integral representations of the corresponding quantum partition functions are derived.
Schwinger Keldysh field theory is a widely used paradigm to study nonequilibrium dynamics of quantum many-body systems starting from a thermal state. We extend this formalism to describe nonequilibrium dynamics of quantum systems starting from arbitrary initial many-body density matrices. We show how this can be done for both bosons and fermions, and for both closed and open quantum systems, using additional sources coupled to bilinears of the fields at the initial time, calculating Green's functions in a theory with these sources, and then taking appropriate set of derivatives of these Green's functions with respect to initial sources to obtain physical observables. The set of derivatives depends on the initial density matrix. The physical correlators in a dynamics with arbitrary initial conditions do not satisfy Wick's theorem, even for noninteracting systems. However, our formalism constructs intermediate ``$n\ensuremath{-}\mathrm{particle}$ Green's functions'' that obey Wick's theorem and provide a prescription to obtain physical correlation functions from them. This allows us to obtain analytic answers for all physical many-body correlation functions of a noninteracting system even when it is initialized to an arbitrary density matrix. We use these exact expressions to obtain an estimate of the violation of Wick's theorem and relate it to presence of connected multiparticle initial correlations in the system. We illustrate this new formalism by calculating density and current profiles in many-body fermionic and bosonic open quantum systems initialized to nontrivial density matrices. We have also shown how this formalism can be extended to interacting many-body systems.
One of the most important predictions of quantum mechanics is the existence of quantum entanglement. Since the inception of quantum theory, the study of entanglement has become ubiquitous in several branches of physics, from high energy to condensed matter physics and quantum information. To this avail, different measures of entanglement in quantum systems have been introduced. The most successful ones are the von Neumann and the Rényi entanglement entropies, which have been heavily studied both in lattice models and in quantum field theories (QFTs). A different, more complete characterization of entanglement is provided by the entanglement (or modular) Hamiltonian, defined as the logarithm of the reduced density matrix. One of the key property of this operator is that, under generic conditions, in 1+1-dimensional conformal field theories (CFTs) presents a local and few-body structure. Beyond just being of theoretical interest, this property has been recently used to efficiently construct the ground state of spin chains in quantum simulators. A recently introduced generalisation of the entanglement Hamiltonian is the negativity Hamiltonian, a non-Hermitian operator which is more useful for mixed states and tripartite entanglement. In this thesis we study these quantifications in many-body quantum systems and in QFTs. In the first part of the thesis, we focus on fermionic systems at equilibrium, and we obtain analytical expressions for the entanglement and the negativity Hamiltonians in the massless Dirac fermion QFT, both in presence of boundaries and at finite temperature. We compare these predictions in the continuum QFT with numerical calculations on the lattice, showing that a proper comparison requires a non-trivial continuum limit of the lattice expression, which we detail. We also investigate a lattice realisation of a non-unitary CFT, the bc-ghost theory, showing that the entanglement Hamiltonian presents an additional term with respect to the unitary case. In the second part, we study the entanglement entropies in an integrable QFT, the massless renormalisation group flow from the tricritical to the critical Ising CFT, demonstrating that the entropy is monotonically decreasing under the renormalisation group flow, as expected from the entropic c-theorem. In the third part we consider out-of-equilibrium problems. We study the melting of a domain wall configuration, computing analytically the entanglement Hamiltonian and showing that the insertion of a defect dramatically modifies the entanglement properties, changing the growth of the entropies from logarithmic to linear in time. Finally, we study generic quantum quenches in free fermionic theories, using the celebrated quasiparticle picture to obtain a prediction for the evolution of the entanglement Hamiltonian during the quench.
The core of this thesis is the path-integral formulation of quantum field theory and its ability to describe strongly-coupled quantum many-body systems of finite size. Collective behaviors can be efficiently described in such systems through the implementation of spontaneous symmetry breaking (SSB) in mean-field approaches. However, as the thermodynamic limit does not make sense in finite-size systems, the latter can not exhibit any SSB and the symmetries which are broken down at the mean-field level must therefore be restored. The efficiency of theoretical approaches in the treatment of finite-size quantum systems can therefore be studied via their ability to restore spontaneously broken symmetries. In this thesis, a zero-dimensional O(N) model is taken as a theoretical laboratory to perform such an investigation with many state-of-the-art path-integral techniques: perturbation theory combined with various resummation methods (Padé-Borel, Borel-hypergeometric, conformal mapping), enhanced versions of perturbation theory (transseries derived via Lefschetz thimbles, optimized perturbation theory), self-consistent perturbation theory based on effective actions (auxiliary field loop expansion (LOAF), Cornwall-Jackiw-Tomboulis (CJT) formalism, 4PPI effective action, ...), functional renormalization group (FRG) techniques (FRG based on the Wetterich equation, DFT-FRG, 2PI-FRG). Connections between these different techniques are also emphasized. In addition, the path-integral formalism provides us with the possibility to introduce collective degrees of freedom in an exact fashion via Hubbard-Stratonovich transformations: the effect of such transformations on the aforementioned methods is also examined in detail.