No Arabic abstract
In this work we propose to simulate many-body thermodynamics of infinite-size quantum lattice models in one, two, and three dimensions, in terms of few-body models of only O(10) sites, which we coin as quantum entanglement simulators (QESs). The QES is described by a temperature-independent Hamiltonian, with the boundary interactions optimized by the tensor network methods to mimic the entanglement between the bulk and environment in a finite-size canonical ensemble. The reduced density matrix of the physical bulk then gives that of the infinite-size canonical ensemble under interest. We show that the QES can, for instance, accurately simulate varieties of many-body phenomena, including finite-temperature crossover and algebraic excitations of the one-dimensional spin liquid, the phase transitions and low-temperature physics of the two- and three-dimensional antiferromagnets, and the crossovers of the two-dimensional topological system. Our work provides an efficient way to explore the thermodynamics of intractable quantum many-body systems with easily accessible systems.
Many-body localization is a striking mechanism that prevents interacting quantum systems from thermalizing. The absence of thermalization behaviour manifests itself, for example, in a remanence of local particle number configurations, a quantity that is robust over a parameter range. Local particle numbers are directly accessible in programmable quantum simulators, in systems of cold atoms even in two spatial dimensions. Yet, the classical simulation aimed at building trust in quantum simulations is highly challenging. In this work, we present a comprehensive tensor network simulation of a many-body localized systems in two spatial dimensions using a variant of an iPEPS algorithm. The required translational invariance can be restored by implementing the disorder into an auxiliary spin system, providing an exact disorder average under dynamics. We can quantitatively assess signatures of many-body localization for the infinite system: Our methods are powerful enough to provide crude dynamical estimates for the transition between localized and ergodic phases. Interestingly, in this setting of finitely many disorder values, which we also compare with simulations involving non-interacting fermions and for which we discuss the emergent physics, localization emerges in the interacting regime, for which we provide an intuitive argument, while Anderson localization is absent.
Quantum coherence quantifies the amount of superposition a quantum state can have in a given basis. Since there is a difference in the structure of eigenstates of the ergodic and many-body localized systems, we expect them also to differ in terms of their coherences in a given basis. Here, we numerically calculate different measures of quantum coherence in the excited eigenstates of an interacting disordered Hamiltonian as a function of the disorder. We show that quantum coherence can be used as an order parameter to detect the well-studied ergodic to many-body-localized phase transition. We also perform quantum quench studies to distinguish the behavior of coherence in thermalized and localized phases. We then present a protocol to calculate measurement-based localizable coherence to investigate the thermal and many-body localized phases. The protocol allows one to investigate quantum correlations experimentally in a non-destructive way, in contrast to measures that require tracing out a subsystem, which always destroys coherence and correlation.
We consider quantum nonlinear many-body systems with dissipation described within the Caldeira-Leggett model, i.e., by a nonlocal action in the path integral for the density matrix. Approximate classical-like formulas for thermodynamic quantities are derived for the case of many degrees of freedom, with general kinetic and dissipative quadratic forms. The underlying scheme is the pure-quantum self-consistent harmonic approximation (PQSCHA), equivalent to the variational approach by the Feynman-Jensen inequality with a suitable quadratic nonlocal trial action. A low-coupling approximation permits to get manageable PQSCHA expressions for quantum thermal averages with a classical Boltzmann factor involving an effective potential and an inner Gaussian average that describes the fluctuations originating from the interplay of quanticity and dissipation. The application of the PQSCHA to a quantum phi4-chain with Drude-like dissipation shows nontrivial effects of dissipation, depending upon its strength and bandwidth.
Non-Hermtian (NH) Hamiltonians effectively describing the physics of dissipative systems have become an important tool with applications ranging from classical meta-materials to quantum many-body systems. Exceptional points, the NH counterpart of spectral degeneracies, are among the paramount phenomena unique to the NH realm. While realizations of second-order exceptional points have been reported in a variety of microscopic models, higher-order ones have largely remained elusive in the many-body context, as they in general require fine tuning in high-dimensional parameter spaces. Here, we propose a microscopic model of correlated fermions in three spatial dimensions and demonstrate the occurrence of interaction-induced fourth-order exceptional points that are protected by chiral symmetry. We demonstrate their stability against symmetry breaking perturbations and investigate their characteristic analytical and topological properties.
Quantum many-body systems (QMBs) are some of the most challenging physical systems to simulate numerically. Methods involving approximations for tensor network (TN) contractions have proven to be viable alternatives to algorithms such as quantum Monte Carlo or simulated annealing. However, these methods are cumbersome, difficult to implement, and often have significant limitations in their accuracy and efficiency when considering systems in more than one dimension. In this paper, we explore the exact computation of TN contractions on two-dimensional geometries and present a heuristic improvement of TN contraction that reduces the computing time, the amount of memory, and the communication time. We run our algorithm for the Ising model using memory optimized x1.32x large instances on Amazon Web Services (AWS) Elastic Compute Cloud (EC2). Our results show that cloud computing is a viable alternative to supercomputers for this class of scientific applications.