Entanglement, noise, and the cumulant expansion

Universal Features of Entanglement Entropy in the Honeycomb Hubbard Model

The entanglement entropy is a unique probe to reveal universal features of strongly interacting many-body systems. In two or more dimensions these features are subtle, and detecting them numerically requires extreme precision, a notoriously difficult task. This is especially challenging in models of interacting fermions, where many such universal features have yet to be observed. In this Letter we tackle this challenge by introducing a new method to compute the Rényi entanglement entropy in auxiliary-field quantum Monte Carlo simulations, where we treat the entangling region itself as a stochastic variable. We demonstrate the efficiency of this method by extracting, for the first time, universal subleading logarithmic terms in a two-dimensional model of interacting fermions, focusing on the half-filled honeycomb Hubbard model at T=0. We detect the universal corner contribution due to gapless fermions throughout the Dirac semi-metal phase and at the Gross-Neveu-Yukawa critical point, where the latter shows a pronounced enhancement depending on the type of entangling cut. Finally, we observe the universal Goldstone mode contribution in the antiferromagnetic Mott insulating phase.

Entanglement Hamiltonian of Interacting Fermionic Models

Recent numerical advances in the field of strongly correlated electron systems allow the calculation of the entanglement spectrum and entropies for interacting fermionic systems. An explicit determination of the entanglement (modular) Hamiltonian has proven to be a considerably more difficult problem, and only a few results are available. We introduce a technique to directly determine the entanglement Hamiltonian of interacting fermionic models by means of auxiliary field quantum Monte Carlo simulations. We implement our method on the one-dimensional Hubbard chain partitioned into two segments and on the Hubbard two-leg ladder partitioned into two chains. In both cases, we study the evolution of the entanglement Hamiltonian as a function of the physical temperature.

Out-of-equilibrium protocol for Rényi entropies via the Jarzynski equality

In recent years entanglement measures, such as the von Neumann and the Rényi entropies, provided a unique opportunity to access elusive features of quantum many-body systems. However, extracting entanglement properties analytically, experimentally, or in numerical simulations can be a formidable task. Here, by combining the replica trick and the Jarzynski equality we devise an alternative effective out-of-equilibrium protocol for measuring the equilibrium Rényi entropies. The key idea is to perform a quench in the geometry of the replicas. The Rényi entropies are obtained as the exponential average of the work performed during the quench. We illustrate an application of the method in classical Monte Carlo simulations, although it could be useful in different contexts, such as in quantum Monte Carlo, or experimentally in cold-atom systems. The method is most effective in the quasistatic regime, i.e., for a slow quench. As a benchmark, we compute the Rényi entropies in the Ising universality class in 1+1 dimensions. We find perfect agreement with the well-known conformal field theory predictions.